# Estimating the volume and age of water stored in global lakes using a geo-statistical approach

Mathis Loïc Messager, Bernhard Lehner, Günther Grill, Irena Nedeva, Oliver Schmitt

► To cite this version:

Mathis Loïc Messager, Bernhard Lehner, Günther Grill, Irena Nedeva, Oliver Schmitt. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 2016, 7 (1), pp.13603. hal-04151085

HAL Id: hal-04151085

https://hal.science/hal-04151085

Submitted on 4 Jul 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

Received 16 Feb 2016 | Accepted 19 Oct 2016 | Published 15 Dec 2016

DOI: 10.1038/ncomms13603

OPEN

# Estimating the volume and age of water stored in global lakes using a geo-statistical approach

Mathis Loïc Messager<fnref n="1" />,†, Bernhard Lehner<fnref n="1" />, Günther Grill<fnref n="1" />, Irena Nedeva<fnref n="1" /> & Oliver Schmitt<fnref n="1" />

Lakes are key components of biogeochemical and ecological processes, thus knowledge about their distribution, volume and residence time is crucial in understanding their properties and interactions within the Earth system. However, global information is scarce and inconsistent across spatial scales and regions. Here we develop a geo-statistical model to estimate the volume of global lakes with a surface area of at least 10 ha based on the surrounding terrain information. Our spatially resolved database shows 1.42 million individual polygons of natural lakes with a total surface area of $2.67 \times 10^6$ km² (1.8% of global land area), a total shoreline length of $7.2 \times 10^6$ km (about four times longer than the world's ocean coastline) and a total volume of $181.9 \times 10^3$ km³ (0.8% of total global non-frozen terrestrial water stocks). We also compute mean and median hydraulic residence times for all lakes to be 1,834 days and 456 days, respectively.

<footnote n="1">¹Department of Geography, McGill University, 805 Sherbrooke Street West, Montreal, Quebec, Canada H3A 0B9. †Present address: School of Aquatic and Fishery Sciences, 1122 Boat Street NE, Box 355020, University of Washington, Seattle, Washington 98195-5020, USA. Correspondence and requests for materials should be addressed to B.L. (email: bernhard.lehner@mcgill.ca).</footnote>
The role of lakes in the global hydrological and biogeochemical cycles is intimately tied to their geometric characteristics of surface area, depth, stored water volume and shoreline length. In addition, the rate of water flowing into and out of lakes depends on their location within the river network, which then defines their hydraulic residence (or turnover) time, that is, the average time that water spends in a lake. Spatially explicit knowledge of all these parameters is crucial for understanding and modelling a wide variety of Earth system processes and interactions with the environment, including hydrological budgets¹; carbon or methane exchange rates²; sediment trapping³; heat fluxes and coupled weather and climate effects⁴; dissolved silica retention⁵; the cycling of pollutants and nutrients⁶; as well as associated ecological processes such as lake productivity⁷; species richness⁸; food chain dynamics⁹; and inland fishery yields¹⁰.

In the absence of globally consistent data describing the multitude of lake parameters, statistical models have been developed that use readily available lake characteristics such as surface area or perimeter as proxies for interpolating large-scale processes¹¹, yet these simplifications pose critical limitations. For instance, despite recent efforts to revise global carbon models, most estimates of fluxes in and out of inland waters default to multiplying an average flux by the total surface area of lakes in a region¹² resulting in wide confidence intervals and high uncertainties.

Data on the volume and depth of lakes on a global scale are scarce and inconsistent. Previous estimates of the total global volume of water contained in lakes ranged from 166 to $275 \times 10^3 \text{ km}^3$ (refs 13–16), with more recent figures converging to $176–180 \times 10^3 \text{ km}^3$ (refs 17–19). However, these estimates are not spatially explicit, are based on incomplete data sets of lake distribution or use simple extrapolation methods that rely on a limited set of variables. Moreover, our attention typically focuses on large lakes, resulting in the oversight of small lentic water bodies when analysing continental- to global-scale systems. Various authors showed the importance of small lakes for processes ranging from evaporation to sediment trapping, greenhouse gas emissions, catchment interactions, lake mixing, diagenetic reactions or aquatic habitat conservation²⁰⁻²³. In particular, while large lakes might dominate processes driven by volume or surface area due to their prevalence at a global scale, small lakes contribute more to the total aquatic-terrestrial interface than large lakes²⁴. It remains inconclusive how different size classes might relate in terms of residence time.

The gap in explicit lake volume data is largely due to the fact that our ability to produce bathymetric maps from limnological surveys is severely constrained by a combination of both technical and operational challenges and the sheer amount of water bodies on Earth. Despite increasingly precise and accurate geo-positioning technologies, satellite imagery and computational cartography²⁵, the bathymetric map creation process through acoustic profiling and interpolation remains a time- and cost-intensive method, as it requires extensive field work by qualified technicians followed by lengthy data analyses. Furthermore, while optical remote sensing methods have been widely used for bathymetric mapping of coastal benthic habitats for several decades²⁶ these techniques are still limited to shallow environments and favourable water conditions²⁷.

In this article, we use a literature review in combination with a novel geo-statistical model to produce a consistent estimate of lake volumes at the global scale for virtually all lakes with a surface area of at least 10 ha (0.10 km²). The main assumption of this model is that the land surface topography surrounding a given lake can be used as a predictor of lake bathymetry. The model is trained and validated on 12,150 existing records of lake

depth and then applied to a new global data set—termed HydroLAKES—that was compiled as part of this project and contains the shoreline polygons of 1.43 million individual lakes and reservoirs. Hydraulic residence time, that is, the average age of lake water, is estimated for each lake using discharge estimates derived from a global hydrological model at high spatial resolution. It should be noted that throughout this article, every effort is made to distinguish ‘natural lakes’ from ‘human-made reservoirs’; yet the general term ‘lake’ is used in all instances where this distinction is considered non-critical.

## Results

**Global lake abundance and characteristics.** HydroLAKES distinguishes 1.43 million individual polygons of natural lakes and human-made reservoirs with a surface area of at least 10 ha. Natural lakes cover a total of $2.67 \times 10^6 \text{ km}^2$ or 1.8% of global land area (Table 1), while large human-made reservoirs (that is, all reservoir polygons $\ge$10 ha from the Global Reservoir and Dam database²⁸) add another $0.26 \times 10^6 \text{ km}^2$ or 0.2% of global land area.

The shoreline length of all natural lake polygons combined stretches for a total of $7.2 \times 10^6 \text{ km}$, plus another $0.5 \times 10^6 \text{ km}$ for large reservoirs. Given the observation that shoreline length is dominated by smaller lake size classes (Table 1; also confirmed by Winslow *et al.*²⁴) we assume that these estimates are significantly underestimating global shoreline length as lakes below 10 ha are not included. Also, our results depend on map scale (in the case of HydroLAKES 1:100,000 and coarser; see Methods) that defines how finely the shorelines are resolved, and it has been observed that a doubling of measurement resolution will cause shoreline length to increase by 15% (ref. 24). In comparison, the global ocean coastline length has been calculated to measure $1.6 \times 10^6 \text{ km}$ at a consistent map scale of 1:250,000 (ref. 29). Thus, the shoreline length of lakes and reservoirs $\ge$10 ha is estimated to be roughly four times longer than the global ocean shoreline, even if scaling effects are accounted for.

According to our study, all natural lakes of at least 10 ha in size contain a total of $181.9 \times 10^3 \text{ km}^3$ of water, representing 0.8% of total global non-frozen terrestrial water stocks, which is equivalent to 1.7% of all the predominantly fresh groundwater on Earth¹⁹. When large human-made reservoirs are added, the storage volume increases to $187.9 \times 10^3 \text{ km}^3$. Of this combined amount, 40.2% is contained in the Caspian Sea and another 33.2% in the next four largest lakes, while total reservoir storage represents 3.2% (Fig. 1 and Table 1). If spread over the world’s landmass, the total water volume would form a layer ~1.26 m in depth, representing 1.7 times the amount of precipitation over land each year³⁰. While the mean depth of all lakes, computed as the sum of lake depths divided by the total number of lakes, is only 3.9 m due to the high number of small lakes, the global area-weighted mean depth, calculated as the total volume divided by the total area of all lakes and reservoirs, is estimated to be 64.2 m. Grouping all natural lakes into bins using logarithmic area and volume size classes (Table 1 and Fig. 2a, respectively) shows that the largest lakes contribute the most to global lake volume. In contrast, no clear size pattern emerges for the contribution to total surface area.

As has been proposed in the literature²⁸,³¹, we observe a roughly 10-fold increase in the number of natural lakes from one logarithmic size class (using surface area) to the next smaller one, confirming that the global distribution of natural lakes and their sizes is scale-invariant down to at least 0.1 km² for surface area, and can be approximated with a Pareto distribution model (Fig. 3). Using this model to extrapolate the next smaller size class of lakes, we estimate that there are 21.2 million natural lakes with
<table><caption>Table 1 | Distribution and morphometric characteristics of lakes worldwide.*</caption><thead><tr><th rowspan="2">Spatial unit</th><th rowspan="2">Number of lakes (10<sup>3</sup>)</th><th rowspan="2">Area (10<sup>3</sup> km<sup>2</sup>)</th><th rowspan="2">Limnicity (% lake area)</th><th rowspan="2">Shoreline length (10<sup>3</sup> km)</th><th rowspan="2">Average shoreline development<sup>†</sup></th><th rowspan="2">Volume (10<sup>3</sup> km<sup>3</sup>)</th><th rowspan="2">Average depth (m)</th><th colspan="2">Residence time (years)</th></tr><tr><th>Mean</th><th>Median</th></tr></thead><tbody><tr><td colspan="10"><strong>Size class (km<sup>2</sup>)</strong></td></tr><tr><td>0.1-1</td><td>1,241.2</td><td>348.4</td><td>0.2</td><td>3,637.6</td><td>1.6</td><td>1.3</td><td>3.5</td><td>4.4</td><td>1.2</td></tr><tr><td>1-10</td><td>165.1</td><td>411.0</td><td>0.3</td><td>2,022.3</td><td>2.2</td><td>2.5</td><td>5.4</td><td>8.2</td><td>1.4</td></tr><tr><td>10-100</td><td>13.4</td><td>331.6</td><td>0.2</td><td>862.4</td><td>3.7</td><td>3.8</td><td>10.4</td><td>17.5</td><td>2.0</td></tr><tr><td>100-1,000</td><td>1.22</td><td>313.5</td><td>0.2</td><td>412.3</td><td>6.1</td><td>6.5</td><td>19.6</td><td>66.1</td><td>4.1</td></tr><tr><td>1,000-10,000</td><td>0.115</td><td>313.3</td><td>0.2</td><td>158.5</td><td>7.8</td><td>10.0</td><td>32.3</td><td>81.4</td><td>4.2</td></tr><tr><td>&gt; 10,000<sup>‡</sup></td><td>0.018</td><td>959.1</td><td>0.6</td><td>79.0</td><td>5.8</td><td>157.8</td><td>139.6</td><td>102.4</td><td>12.3</td></tr><tr><td colspan="10"><strong>World</strong></td></tr><tr><td>Natural lakes ≥10 ha</td><td>1,421.0</td><td>2,676.8</td><td>1.8</td><td>7,172.2</td><td>1.7</td><td>181.9</td><td>3.8</td><td>5.0</td><td>1.2</td></tr><tr><td>Incl. large reservoirs<sup>§</sup></td><td>1,427.7</td><td>2,926.7</td><td>2.0</td><td>7,661.9</td><td>1.7</td><td>187.9</td><td>3.9</td><td>5.0</td><td>1.2</td></tr><tr><td>Natural lakes ≥1 ha</td><td>21,152.4<sup>||</sup></td><td>3232.2</td><td>2.2</td><td>NA</td><td>NA</td><td>182.9<sup>*</sup></td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td colspan="10"><strong>Continent</strong></td></tr><tr><td>North America<sup>††</sup></td><td>991.9</td><td>1,229.5</td><td>5.1</td><td>4,990.9</td><td>1.8</td><td>36.6</td><td>3.7</td><td>4.5</td><td>1.4</td></tr><tr><td>Europe<sup>#</sup></td><td>280.7</td><td>781.2</td><td>3.4</td><td>1,264.5</td><td>1.5</td><td>103.8</td><td>4.6</td><td>7.4</td><td>1.6</td></tr><tr><td>Asia<sup>**</sup></td><td>66.2</td><td>274.8</td><td>0.9</td><td>391.7</td><td>1.6</td><td>7.3</td><td>2.9</td><td>5.9</td><td>0.3</td></tr><tr><td>South America</td><td>53.8</td><td>103.7</td><td>0.6</td><td>296.7</td><td>1.6</td><td>3.1</td><td>3.3</td><td>1.6</td><td>0.2</td></tr><tr><td>Africa</td><td>15.2</td><td>232.0</td><td>0.8</td><td>120.1</td><td>1.7</td><td>30.6</td><td>2.5</td><td>3.5</td><td>0.3</td></tr><tr><td>Oceania<sup>‡‡</sup></td><td>13.2</td><td>55.7</td><td>0.7</td><td>108.3</td><td>1.6</td><td>0.4</td><td>2.8</td><td>7.7</td><td>0.8</td></tr><tr><td colspan="10"><strong>Countries with most lakes<sup>§§</sup></strong></td></tr><tr><td>Canada</td><td>879.8</td><td>856.5</td><td>8.6</td><td>4,498.1</td><td>1.8</td><td>12.6</td><td>3.7</td><td>4.6</td><td>1.5</td></tr><tr><td>Russia</td><td>201.2</td><td>667.4</td><td>4.0</td><td>832.9</td><td>1.4</td><td>102.2</td><td>4.5</td><td>9.5</td><td>2.5</td></tr><tr><td>USA</td><td>102.5</td><td>340.3</td><td>3.6</td><td>427.1</td><td>1.6</td><td>23.5</td><td>3.4</td><td>3.3</td><td>0.8</td></tr><tr><td>China</td><td>23.8</td><td>81.0</td><td>0.9</td><td>144.7</td><td>1.6</td><td>1.0</td><td>3.4</td><td>6.5</td><td>0.4</td></tr><tr><td>Sweden</td><td>22.6</td><td>34.3</td><td>7.7</td><td>119.9</td><td>1.6</td><td>0.5</td><td>4.8</td><td>2.7</td><td>1.1</td></tr><tr><td>Brazil</td><td>20.9</td><td>31.4</td><td>0.4</td><td>122.3</td><td>1.8</td><td>0.2</td><td>2.5</td><td>0.6</td><td>0.1</td></tr><tr><td>Norway</td><td>20.0</td><td>13.9</td><td>4.3</td><td>86.7</td><td>1.7</td><td>0.3</td><td>7.0</td><td>2.1</td><td>0.8</td></tr><tr><td>Argentina</td><td>13.6</td><td>27.9</td><td>1.0</td><td>67.7</td><td>1.5</td><td>0.6</td><td>2.2</td><td>3.6</td><td>0.7</td></tr><tr><td>Kazakhstan</td><td>12.4</td><td>61.9</td><td>2.3</td><td>84.6</td><td>1.5</td><td>0.4</td><td>1.9</td><td>16.7</td><td>3.5</td></tr><tr><td>Australia</td><td>11.4</td><td>49.5</td><td>0.6</td><td>92.8</td><td>1.6</td><td>0.1</td><td>2.3</td><td>9.0</td><td>1.2</td></tr><tr><td colspan="10"><strong>Largest lakes by volume</strong></td></tr><tr><td>Caspian Sea</td><td>NA</td><td>377.0</td><td>NA</td><td>15.8</td><td>7.3</td><td>75.6</td><td>201</td><td>295.4</td><td>NA</td></tr><tr><td>Lake Baikal</td><td>NA</td><td>32.0</td><td>NA</td><td>2.7</td><td>4.2</td><td>23.6</td><td>739</td><td>374.6</td><td>NA</td></tr><tr><td>Lake Tanganyika</td><td>NA</td><td>32.8</td><td>NA</td><td>2.1</td><td>3.3</td><td>18.9</td><td>577</td><td>402.6</td><td>NA</td></tr><tr><td>Lake Superior</td><td>NA</td><td>81.8</td><td>NA</td><td>5.2</td><td>5.2</td><td>12.0</td><td>147</td><td>132.5</td><td>NA</td></tr><tr><td>Lake Malawi</td><td>NA</td><td>29.5</td><td>NA</td><td>1.7</td><td>2.8</td><td>7.7</td><td>261</td><td>218.5</td><td>NA</td></tr><tr><td>Lake Michigan</td><td>NA</td><td>57.7</td><td>NA</td><td>2.9</td><td>3.4</td><td>4.9</td><td>84</td><td>82.0</td><td>NA</td></tr><tr><td>Lake Huron</td><td>NA</td><td>59.4</td><td>NA</td><td>8.9</td><td>10.3</td><td>3.6</td><td>60</td><td>12.3</td><td>NA</td></tr><tr><td>Lake Victoria</td><td>NA</td><td>67.2</td><td>NA</td><td>7.4</td><td>8.1</td><td>2.6</td><td>41</td><td>50.4</td><td>NA</td></tr><tr><td>Great Bear Lake</td><td>NA</td><td>30.5</td><td>NA</td><td>5.3</td><td>8.6</td><td>2.2</td><td>72</td><td>130.3</td><td>NA</td></tr><tr><td>Kara-Bogaz-Gol</td><td>NA</td><td>18.7</td><td>NA</td><td>1.0</td><td>2.0</td><td>1.9</td><td>101</td><td>NA</td><td>NA</td></tr></tbody></table>

NA, not applicable.

*Unless noted otherwise, data refer to all natural lakes ≥ 10 ha (0.1 km²) contained in the HydroLAKES database including regulated (natural) lakes and small unreported reservoirs.

†Shoreline development measures the degree of deviation of a lake's surface shape from a circle; a value of 1 implies a perfectly circular shape and higher values indicate increased shoreline sinuosity.

‡Includes all lakes 10,000-100,000 km² and Caspian Sea.

§All reservoir polygons ≥10 ha from GRanD database²⁸.

||Extrapolated using a Pareto distribution model for lakes 0.01-0.35 km².

¶Extrapolated using a Pareto distribution model for lakes 0.00001-0.0005 km³.

#Includes all of Russia.

**Includes Middle East and Turkey.

††Includes Mexico, the Caribbean and Central America.

‡‡Includes Australia, New Zealand, Micronesia, Melanesia and Polynesia.

§§International lakes are assigned to only one country based on the location of the lake outlet; the Caspian Sea is assigned to Russia (and thus Europe).

a surface area of at least 0.01 km² (1 ha), representing a total of 3.23 × 10⁶ km² or 2.2% of global land area. Similarly, we extrapolate that there are 25.4 million natural lakes with a volume of at least 0.00001 km³ (the equivalent of a lake with an area of 0.01 km² and a depth of 1 m), storing a total water volume of 182.9 × 10³ km³.

In terms of regional distribution, areas with high lake densities are nearly all encompassed within the extent of the last glacial maximum in Northern Canada and Scandinavia as well as in

some areas of Alaska and Russia (Fig. 4a), with some smaller hotspots found in the alpine zones of the Andes, Rockies and Himalayas. Divergent plate boundaries with deep tectonic lakes such as observed in the African Rift Valley, as well as large floodplains such as in the Amazon Basin or in coastal China are also places where lakes are prominent features of the landscape. The global distribution of lake volume (not shown) follows the same general trends as lake area, though alpine zones and rift valleys with deep lakes tend to hold more water than areas where
<figure><img src="image_3.png" /><figcaption>**Figure 1** | Global distribution of water volume stored in lakes and reservoirs with a surface area of at least 10 ha. Total volume is $187.9 \times 10^3 \text{ km}^3$. Data for large lakes are empirical, while volumes of medium and small lakes are modelled. Data for large and medium human-made reservoirs are from the Global Reservoir and Dam (GRanD) database²⁸. Distinction between fresh and saline water is only available in the empirical data for large lakes.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2** | Distribution of global lakes sizes and hydraulic residence times. (a) Total volume of natural lakes in km³ and surface area in km² across volume size classes by continent and for the world. (b) The frequency distribution of hydraulic residence times by continent and for the world. Colour legend in a refers to both panels. For definition of continents see Table 1. Note that Europe includes the Caspian Sea and Lake Baikal, the two most voluminous lakes on Earth.</figcaption></figure>

shallower lakes were formed by continental glaciers (for example, on the Canadian Shield) or by fluviatile processes (see lake depth distribution on Fig. 4b).

**Hydraulic residence times.** The mean lake residence time computed for all natural lakes with a surface area of at least 10 ha,

calculated as the sum of individual lake residence times divided by the total number of lakes, is 1,834 days (about 5.0 years), while the median is 456 days. The discrepancy between these statistics stems from the skewed frequency distribution of residence times, which is dominated by smaller lakes with residence times below the average (Fig. 2b). Some lakes show much longer residence times; this generally applies to large lakes located in areas with
<figure><img src="image_5.png" /><figcaption>**Figure 3** | Global size and volume distributions of natural lakes using a Pareto model. Distributions are plotted as the total number of global lakes larger than a given surface area (a) or volume (b) derived from data contained in HydroLAKES. Yellow points represent data that were not included for fitting the log-log regression. Red triangles represent extrapolated values based on the log-log regression. See Methods for further explanations.</figcaption></figure>

little discharge or small catchment areas such as Lake Baikal with an estimated residence time of 375 years (reported as 321 years³²). Others display very short residence times; this may indicate lakes situated within a larger river channel that are flushed rather rapidly by river discharge; but may also include some errors due to incorrect co-registration of lake outfalls to the river network in our model (such as small oxbow lakes being assigned to the main-stem river).

Overall, we find a trend of shorter residence times for smaller lakes and longer residence times for larger lakes, yet there is significant scatter among this result, and long residence times can occur throughout the entire lake size spectrum. Finally, the average age of all lake water, calculated as the volume-weighted average of all residence times, is 255 years. This estimate, however, is strongly dominated by the influence of the most voluminous lakes (Table 1).

## Discussion

This study investigated the relationship between lake geometry, surrounding landscape characteristics and lake bathymetry. Our geo-statistical model confirmed earlier studies in that, beyond lake surface area, the variables with the most predictive power are based on expressions of the variation in elevation around the lakes. We found that the depth of lakes is generally best reflected by the slope of their immediate surroundings (that is, within 100 m of the shoreline), which is in agreement with earlier studies, such as that of Sobek et al.³³ who used a similar approach based on the maximum slope within 50 m of the shoreline. Despite its larger extent, the statistical performance of our global model is comparable to that of regional studies, including a recent model developed for 55,000 lakes in Quebec, Canada³⁴ (see Methods for more model comparisons).

The limited amount of globally consistent data on the geological age of the landscape and the various morphological processes involved in shaping lakes constrained the ability of our model to differentiate between different lake types. This lack of data places an inherent limit on the degree of variance in lake depths that can be explained. There is an exceptional diversity of lakes in terms of formation processes, including fluvial, aeolian, landslide or volcanic origins, which highly influence their

morphometry³⁵,³⁶. Ideally, each of these types should be treated by customized sub-models. Furthermore, different bathymetry may be due to differential rates of sedimentation, which in turn depend on variations in geology, climate and vegetation in the associated upstream catchment.

As with any model, errors in the underpinning data sets represent a substantial source of uncertainty (Supplementary Discussion). Accordingly, identifying whether the over- and underestimations found in the model predictions are due to reference data errors or modelling errors is an arduous task. Given the uncertainties, we caution from interpreting volume estimates for single lakes, as individual errors can be very large. Rather, our results are designed for use in regional to global scale studies, where errors are reduced due to data aggregation. We believe that progress in the resolution and accuracy of elevation data set holds promise for the improvement of future global lake depth estimates.

Previous estimates of the number and total surface area of lakes have varied widely (Fig. 5 and Supplementary Table 1), with total surface area ranging from about 2 to 5 million km². Some of this variation can be attributed to different definitions and inclusion or exclusion of certain lake sizes or types (for example, fresh versus saline or permanent versus intermittent). Our new estimate of 2.67 × 10⁶ km² for 1.42 million natural lakes ≥ 10 ha, and 3.23 × 10⁶ km² for 21.2 million natural lakes ≥ 1 ha, respectively, are generally exceeding older estimates indicating the comprehensiveness of HydroLAKES. At the same time our results stay below some higher estimates of the past decade³¹,³⁷, which we believe is due to our enhanced ability to differentiate natural lakes from reservoirs, as well as variations in the statistical extrapolation of smaller lakes.

More recently, Verpoorter et al.³⁸ proposed much higher global estimates of lakes and lake area (Fig. 5 and Supplementary Table 1) and attributed the soaring number of small lakes to improvements in the resolution of satellite sensors. In this regard, however, a significant caveat lies in the difficulty to distinguish lakes from other open water features including rivers, reservoirs and inundated floodplains or wetlands. Producing 'lake-only' data sets typically requires extensive supervised classification and manual adjustments³⁹. As HydroLAKES underwent thorough manual corrections to detect and remove fluvial features as well as
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4** | Patterns of global lake distribution. (a) Lake area density (limnicity) calculated as percent area covered by lakes within a 25 km radius. (b) Average depth of all lakes within a 25 km radius, weighted by their partial area within that radius. Both maps include reservoirs from the Global Reservoir and Dam (GRanD) database²⁸.</figcaption></figure>

large reservoirs (see Methods), our reversal to lower 'lake-only' numbers may indicate a possible inflation of global lake estimates through remote sensing imagery due to misclassifications or differences in defining lakes versus rivers and wetlands. Also, as opposed to the polygon structure of HydroLAKES, spectral classification of raster-based remote sensing data³⁸,⁴⁰ cannot easily distinguish individual functional water bodies, and is therefore better suited to represent continuous surface water masks rather than discrete objects.

Furthermore, current remote sensing data sets mostly represent a snapshot of surface water on Earth at a given time, and thus may fail to consistently capture intermittent water bodies or historic maximum water extents. In contrast, topographic maps often integrate knowledge over much longer time periods to delimit lakes. However, this apparent limitation of remote sensing imagery also provides an opportunity for future

advancements as it allows for the dynamic monitoring of water extent and storage variations over time using multi-temporal data sets⁴¹,⁴².

Our total volume estimate for natural lakes ≥10 ha of 181.9 × 10³ km³ is similar to those found in the literature of the past decades¹⁷⁻¹⁹ (Supplementary Table 1). Estimates for the largest global lakes seem particularly robust. For example, our data matches the estimate of Tamrazyan¹⁶ of 160,600 km³ for lakes over 6,000 km² nearly exactly. Without the 10 most voluminous lakes, we estimate that the remaining world lakes contain 28,991 km³ of water, while Shiklomanov and Rodda¹⁹, and Ryanzhin et al.¹⁸ estimated 23,265 and 26,465 km³, respectively. It remains important to note, however, that we rely on depth estimates from literature for all lakes larger than 500 km² (see Methods), thus the amount of modelled lake volume accounts for only 12,224 km³, or about 6.5% of total global lake
<figure><img src="image_8.png" /><figcaption>**Figure 5** | Comparison of the total number of global lakes between previous estimates and this study. Lake numbers are shown by logarithmic size class for lakes with a surface area of at least 10 ha. The asterisk (*) indicates that the numbers were extrapolated from canonical data sets with log-log regressions between the number of lakes and their surface area (as shown in Fig. 3). The number of lakes over 10 km² for Lehner and Döll⁷⁰ and Downing et al.³¹ are both obtained from the Global Lakes and Wetlands Database (GLWD)⁷⁰. All data refer to natural lakes except Verpoorter et al.³⁸ which include reservoirs.</figcaption></figure>

volume. Also, our volume estimate does not take seasonal fluctuations into account.

This study provides the first global data set where each lake is individually assigned a residence time, which is a required metric for understanding and managing lentic ecosystems and the effect of anthropogenic disturbance on inland waters in general⁴³. Despite the mentioned limitations and errors, we believe that our estimates provide a reasonable first-order proxy of residence times for regional assessments, yet require further validation and adjustment when investigating individual lakes. Our new finding of 5 years for the mean global hydraulic residence time of all lakes ≥10 ha falls within the range provided by L'vovich¹⁴ who calculated the 'total rate of exchange of surface water on land' as 7 years. On the other hand, Korzoun et al.⁴⁴ approximated global lake residence time to be 17 years in their widely cited work¹⁹,⁴⁵. It should be noted, however, that these previous estimates only relate aggregated global or regional freshwater stocks and fluxes and are thus conceptually different to our new method of calculating lake-specific residence times by taking the actual location of each lake within the drainage network into account. Our lower estimate is evidence of the fact that the surface-water-land interface is dominated by lakes with smaller volumes (see Table 1), resulting in a relatively faster cycling of water through these systems than if they were lumped and related to regional flows.

In conclusion, we believe that the findings of this study support an enhanced understanding of the role of lakes in the Earth

system by defining a robust and reliable global database of lake distribution and characterization. For the field of global limnology to achieve its full potential, data on lake area, volume and depth are 'of rudimentary importance'⁴⁶. The model presented here aims to be a stepping-stone towards providing this information. It relies on estimating lake bathymetry based on literature as well as a novel statistical approach that was derived by combining and analysing numerous existing lake data sets in conjunction with surrounding topography data. Despite these advancements, continued efforts to conduct actual bathymetric measurements remain the most vital component for future improvements in global lake volume estimation.

The model results presented here contribute to the Hydro-LAKES database, a spatially resolved object-oriented baseline data set where each uniquely identified lake polygon comes with a suite of morphometric attributes: surface area, perimeter, mean depth, volume and residence time. As every lake polygon is linked via its pour point to the river network of the HydroSHEDS database⁴⁷, a variety of additional characteristics are readily available for each lake's catchment area, such as average discharge, or can be derived by incorporating auxiliary data layers, such as landscape characteristics or population statistics.

Beyond providing a mapping repository, HydroLAKES also offers a structured framework to which additional information can be added in support of subsequent analyses. We consider it particularly promising to combine the object-oriented features and current attributes of the HydroLAKES database with ecological information such as species and biodiversity distributions and/or with the dynamic capabilities of multi-temporal remote sensing imagery to monitor seasonal variations. These combinations would offer a new perspective on eco-hydrological assessments regarding the state and future of global lakes and could establish a baseline for lake monitoring as required by global inland water ecological research programs⁴⁸. We hope that HydroLAKES will promote this endeavor as a new source of information for scientists, managers and decision-makers to be used in biogeochemical models, river network routing, environmental planning or ecological applications at the global scale.

## Methods

**Topographic data.** All topographic calculations were performed in a Geographic Information System (GIS) environment on the elevation data provided by EarthEnv-DEM90 (ref. 49), a composite global 3 arc-second (~90 m) digital elevation model spanning from 60°S to 83°N. Slope data were computed based on Horn's method⁵⁰ with latitudinal corrections for the distortion in the XY spacing of geographic coordinates by approximating the geodesic distance between cell centres.

**Empirical data on lake depths and volumes.** Bathymetric and/or depth data for 12,150 natural lakes were obtained from 14 existing national and international data sets (Supplementary Table 2), the majority being from the USA, Canada and Europe; and own compilations based on about 40 peer-reviewed articles. Data that were originally provided as bathymetric lines or in raster format were converted into volumetric and depth estimates using 3D GIS procedures. Supplementary Fig. 1 provides a map overview of the spatial distribution of the reference data.

Of the 12,150 lake records, 7,049 were used as training data for developing the model of lake depth, and 5,101 were used as independent validation data. To achieve the best possible model performance, higher quality data was preferred for model training while the validation data were either of limited reliability or from regions that were already well represented in the training data set. In particular, we refrained from training our model with data from the Global Lake Database (GLDB)⁵¹, a commonly applied global lake depth data set for parameterization in numerical weather prediction and climate models, as only few GLDB depth estimates contain a reference; we therefore perceived its quality as not consistent and difficult to judge. Nevertheless, we included GLDB as an independent validation data set because of its ability to cover lakes outside of our training data set, serving as a benchmark to test model performance in untrained regions.

The Svenskt Vattenarkiv (SVAR) from the Swedish Meteorological and Hydrological Institute was split in halves to serve for both training and validation purposes by random sampling from logarithmic lake size classes (0.1–1, 1–10,
10–100 and 100–500 km²). In cases of overlap between SVAR and GLDB, priority was given to SVAR due to its higher data quality, and the corresponding record of GLDB was removed. Finally, the Quebec data set by Heathcote et al.³⁴ was used as independent validation data as it provided similar lake types as the Ontario training data set. The inclusion of extensive training and validation data from British Columbia, Ontario, Quebec and Sweden was particularly important to adequately account for lakes in boreal regions, which represent the largest portion of world lakes.

In addition, large lakes—here defined as those over 500 km²—were given special attention as their geomorphologic characteristics, including depth, are often determined by highly unique and complex interactions of different formation processes and geologic constraints. For these 347 lakes, we refrained from modelling their volume and instead used empirical mean depth estimates from about 170 additional sources, the great majority being peer-reviewed or governmental and institutional documents including a global compilation by Herdendorf⁵². Finally, the storage capacities of 6,797 large reservoirs or regulated lakes were added mostly from the Global Reservoir and Dam (GRanD) database (version 1.1)²⁸.

**Discharge data.** To compute the individual hydraulic residence time of each lake, we used estimates of long-term (1971–2000) average 'naturalized' discharge provided by the state-of-the-art global integrated water model WaterGAP (model version 2.2 as of 2014)⁵³. The data were spatially downscaled from their original 0.5 degree pixel resolution (~50 km at the equator) to the 15 arc-second (~500 m) resolution of the global HydroSHEDS river network⁴⁷ using geo-statistical approaches⁵⁴.

**HydroLAKES polygon layer.** To apply the lake volume model, a global database of all lakes with a surface area of at least 10 ha was generated, comprising the shoreline polygons of 1,427,688 individual water bodies. This novel data set aims to be as comprehensive and consistent as possible at a global scale and contains both freshwater and saline lakes, including the Caspian Sea, as well as human-made reservoirs and regulated lakes. We created the HydroLAKES database by compiling, correcting and unifying several near-global and regional data sets (Supplementary Table 3), foremost the SRTM Water Body Data⁵⁵ for regions from 56°S to 60°N, and CanVec⁵⁶ for most North American lakes.

If the original data were provided in raster format, they were first vectorized using boundary smoothing procedures to create polygons. Other main processing steps in the creation of HydroLAKES included manual identification and removal of river and wetland polygons; removal of duplicates and overlapping polygons; dissolving of segmented polygons into individual lake entities; correction of corrupt or incorrect polygon geometry; removal of small islands (<3 ha) within lakes; smoothing of water body shorelines to reduce inconsistencies between data sets of different initial resolution; and establishing a 10 ha (0.10 km²) cut-off based on lake surface area. A more detailed description of the production steps is provided in a technical documentation that is distributed together with the database.

The resolution of the underpinning source data ranged from 1:24,000 to 1:1 million for the data in vector format, and from 30 to 250 m pixels for the data in raster format. Due to these inconsistencies in scale and the various polygon transformations, including smoothing and generalization steps during the map consolidation process, the resulting resolution of the global HydroLAKES database cannot be strictly defined. However, regional comparisons with maps at a variety of known resolutions, as well as tests using shoreline scaling laws²⁴ suggest the following scales as best approximation: about 1:100,000 for Canada and Alaska (that is, accounting for two thirds of global lakes); about 1:250,000 for Europe and all areas below 60°N outside of Canada (that is, accounting for most of the global landmass); and about 1:1 million for the remaining areas (that is, northern Russia and Greenland). Therefore, the resulting map scale is estimated to be between 1:100,000 and 1:250,000 for most lakes globally, with some coarser ones at 1:1 million.

A spatial co-registration between HydroLAKES and the river network of the HydroSHEDS database⁴⁷ was established by linking each lake to the most downstream river pixel that drains the lake. This pour point is typically near the lake's shoreline but can also fall inside a polygon for terminal lakes in endorheic basins. HydroLAKES was also combined with the GRanD database²⁸ to identify and flag 6,797 polygons as large reservoirs or regulated lakes.

**Topographic lake depth.** Before applying geo-statistical models, we investigated the deterministic estimation of lake volumes using geospatial tools to explicitly simulate the bathymetry inside each lake (Supplementary Fig. 2a). The underlying assumption of this approach is that a lake's bathymetry is typically formed by the same geophysical processes as the surrounding landscape³⁶ and that therefore the terrain slope around a lake can be extrapolated into the lake to represent its shape (Supplementary Fig. 2b). Each lake pixel's depth can then be calculated as the product of the tangent of the slope and the distance from the shore. In the simplest case, this can be achieved through linear expansion of the slope into the middle of the lake. In reality, however, the bottom of many lakes tends to flatten towards the centre to form a bowl-like shape due to the accumulation of sediments. This can be approximated with more advanced extrapolation or splining techniques. We tested

a variety of GIS-based methods including power functions where depth was derived as [(distance from shore)⁸ · tan(slope)] using values of 0.90, 0.95 and 0.97 for the exponent x.

Furthermore, we compared different approaches on how to sample and allocate the available shoreline slope values, following two basic concepts (Supplementary Fig. 2c). In the first concept, each pixel in the lake is assigned the slope value of the closest pixel on shore. In the second concept, all pixels in the lake are assigned the average slope value of all pixels within a given distance around the lake. While the former approach is more spatially explicit by incorporating local terrain variations as well as the geometric shape of the water body, the latter characterizes the landscape in a more general way with a single slope proxy representing the surrounding topography. For the first approach, we controlled for potential outliers and data artefacts that can affect local slope calculations by applying smoothing filters; that is, we calculated the mean slope within a 3 × 3, 5 × 5 and 7 × 7 neighbourhood and ultimately selected the 3 × 3 kernel. For the second approach, we tested a variety of different buffer distances around the lake, ranging from a single pixel width to 1,000 m. Finally, all pixel depths of a lake were averaged to calculate the mean 'topographic depth' ($D_t$) per lake.

To constrain computing efforts, we conducted initial tests to identify the best method of determining $D_t$ using empirical reference depth data for 173 European lakes located in a variety of geomorphologic settings, mostly clustered in the Baltic countries, the Alps and Ireland. Comparisons between different extrapolation methods and settings were performed through visual examination of the data as well as goodness-of-fit tests. Although the various models performed differently for the test regions, the gains in predictive power from applying different extrapolation functions and settings did not follow a clear pattern. Thus we selected the simplest approach of linearly extending the slope of the closest shoreline pixel to the centre of the lake and applied this method to all 7,049 reference lakes of the training data set to estimate their mean topographic depth.

Supplementary Fig. 3 shows the resulting scatterplot between predicted and reference mean depths for different regions. While we observe a good overall trend, there is also a clear tendency in the model to overestimate lake depth. We believe that at least part of this error can be explained by our neglect to represent sedimentation effects through adjusting (flattening) of the simulated lake bottom. However, as global sedimentation rates are not readily available, we refrained from further investigating deterministic or region-specific solutions and instead continued to use 'topographic depth' only as a predictor variable within statistical models.

**Geo-statistical model development.** Most statistical models that aim to predict lake volumes have used lake surface area as a proxy, often relying on a power relationship between the two parameters⁵⁷,⁵⁸. However, in analogy to the concept of 'topographic lake depth' as explained above, it can be generally postulated that deep lakes preferentially occur in mountainous regions while plains are indicative of shallower lakes. Based on this principle, multiple studies drew statistical relationships between the topography surrounding a lake and its depth using a variety of terrain indices for the prediction, such as elevation ranges and slopes³³,⁵⁹⁻⁶². Most recently, Heathcote et al.³⁴ developed an empirical model to estimate the mean and maximum lake depths of over 55,000 lakes in Quebec, Canada, using a multiple linear regression based on surface area and an index of topographic variation within a buffer zone around the lake as the predictor variables.

Following the same conceptual approaches, we tested a large variety of methods and equations, including those suggested in previous studies, to estimate the volume of our 7,049 empirical lake records. In particular, we investigated a multitude of regression equations of different types and complexities to predict lake depth using an array of 41 predictor variables that describe topographic landscape and lake geometry characteristics, including slopes and elevation ranges in varying buffer distances from the shoreline; lake surface area; perimeter; shape factors; and shoreline development (that is, the ratio of the length of the lake shoreline to the perimeter of a circle of equal area to that of the lake³⁶). As part of the model development process, we also tested the novel proxy of 'topographic lake depth' as an alternative predictor variable.

Following the approach of Håkanson and Peters⁵⁹, different lake size classes were distinguished (that is, 0.1–1, 1–10, 10–100 and 100–500 km²) and specific regression coefficients were derived for each class. The separation into groups allows for the use of scale-appropriate, customized predictor variables per size class and limits the extent to which large lakes dominate parameter results by avoiding unequal probability sampling.

After comparing a large range of models and settings (for more details see below), we identified a multi-variable statistical model as the best performing method for predicting lake depth and volume at a global scale. The model produces an estimate of mean depth for each lake that is then multiplied by the lake polygon surface area to obtain lake volume.

**Model comparisons and selection.** A main goal of the model selection process was to develop a parsimonious model for the data at hand, while avoiding the risk of overfitting. Many of the investigated 41 predictor variables were significantly correlated to the mean lake depth as well as to each other. This stems from the fact that a large portion of these candidate predictors represent the same characteristic,
yet applied to different buffer sizes around the lakes (for example, elevation range within buffers of 1 and 2 km width). Thus, to avoid multicollinearity, only a limited number of variables were tested simultaneously in the multiple regression models. Multicollinearity was assessed using the variance inflation factor, together with observations of coefficient parameters (for example, parameters of unlikely sign or magnitude) following Belsley *et al.*<sup>63</sup>

After testing both simpler (for example, using surface area as the only predictor variable) and more complex methods (for example, non-linear models or classification and regression trees), statistical model types ranging from simple linear to multiple log-linear regressions were chosen for further investigation. The following six main model types represent typical stages of model complexity (see Supplementary Table 4 for equations): Model 1: a simple linear regression using surface area as the only predictor; Model 2: a simple linear regression using mean topographic lake depth as the only predictor; Model 3: a simple linear regression using surface area together with a topographic variable as suggested by Heathcote *et al.*<sup>34</sup> (that is, the difference between lake surface area and mean landscape elevation within a buffer width equal to 25% of the diameter of a circle that represents the lake area); Model 4: group-specific multiple regressions using the same predictors as Model 3; Model 5: group-specific multiple regressions using surface area together with terrain and lake shape variables as predictors; Model 6: group-specific multiple regressions using surface area together with topographic lake depth and other terrain and shape variables as predictors.

We compared a large number of different instances of these six model types using visual examinations and interpretation of the resulting scatterplots (Supplementary Figs 4 and 5). While the significance of tested predictor variables was a requirement for their inclusion in a model, their performance and ultimately their selection was determined by a variety of statistical indices (Supplementary Table 5) including the adjusted correlation coefficient (R<sup>2</sup>); the root-mean-square error (RMSE); the mean absolute error (MAE); and the symmetric mean absolute percent error (SMAPE = 100 × 1/N × (observed − predicted value)/(observed value + predicted value)/2)).

Both Models 3 and 4 using surface area and the topographic variable from Heathcote *et al.*<sup>34</sup> as the predictors performed well, particularly for small lakes, yet were not as satisfactory as Models 5 and 6, in part as they yielded more skewed predictions for large lakes. Model types 5 and 6 showed the best performances and delivered overall equivalent results. Model 5 was ultimately selected for application due to its simpler calculation method (that is, it does not require the GIS-based calculation of topographic depth for every lake), and its best-fit equations were determined as presented in Supplementary Table 4. It should be noted, however, that despite its inability to enhance our current model performance we believe that the new proxy of 'topographic lake depth' can potentially improve bathymetry estimates if combined with regional information regarding lake morphological types, as it incorporates local terrain and shape characteristics.

Model 5 was able to reduce the SMAPE for the predicted mean depth of all lakes from 64.1% (simplest model type 1) to 47.4% while increasing R<sup>2</sup> from 0.16 to 0.51. Also, Model 5 yielded the least biased estimate of volume for the entire data set—reduced from 43% to <1%, respectively. Improvements in predictions were particularly strong for larger lakes but remained marginal for smaller ones. It should be noted that Models 3 and 4 performed slightly better than Model 5 when compared for Quebec data only, indicating their adequacy in the region for which their topographic variable was originally developed<sup>34</sup>.

As a general observation, all models performed much better when analysing volume estimates (maximum R<sup>2</sup> of 0.91) as compared to depth estimates (maximum R<sup>2</sup> of 0.51). This can be explained by volume being derived through multiplication of average depth with lake surface area—the latter being an observed value with little uncertainty—which makes volume a highly constrained variable. Nevertheless, while the predictive power of our model may be better described by its ability to estimate average depth, the ultimate goal of the desired model was to provide volume estimates for which high correlations were achieved.

**Validation and uncertainty of statistical model results.** To validate our results and to quantify the uncertainty of the regression models as well as the potential bias inherent in the sample of lakes used for model training, we performed a two-pronged validation analysis. First, a bootstrapping analysis was performed on the multiple regression models and lake size classes using 10,000 replicates following Fox and Weisberg<sup>64</sup> (results are presented in Supplementary Table 5). This analysis confirmed that the uncertainty and potential bias in the regression coefficients, performance indices and total predicted volume are generally within acceptable ranges, yet with less predictive power for smaller lakes in comparison with larger ones. However, we found that variability identified by resampling with replacement generally increased from smaller to larger lakes. This can be attributed to the paucity of data for large lakes and the increasing range of depths that larger lakes can have. Moreover, as lake size increases, the link between the mean depth of a lake and its surrounding landscape becomes more tenuous in terms of geomorphology. For the smallest lakes in our data set (0.1–1 km<sup>2</sup>), we found the 95% confidence intervals of the SMAPE and the linear regression R<sup>2</sup> to be 46.9–49.8% and 0.20–0.27, respectively, whereas for the largest lakes (100–500 km<sup>2</sup>), the 95% confidence interval of the SMAPE was 45.6–70.2% with much higher R<sup>2</sup> values of 0.59–0.83.

The second validation approach consisted of comparing the estimates of our best performing Model 5 against the independent empirical lake volume data of the

5,101 lakes in our validation data set (see Supplementary Table 2 for data sources). Overall the model performance was satisfying for unobserved lakes given the breadth of world lakes analysed (Supplementary Figs 1 and 6). Spatially, the model tends to underestimate the depth and volume of deep lakes in young mountainous regions such as the Andes or the European Alps, whereas lakes in low-relief areas of northwestern Russia, Finland and Sweden tend to be overestimated. No clear bias is visible for other regions of the world.

When computed for the entire validation data set, the SMAPE between predicted and reference volume was 48.8%, that is, nearly identical to the SMAPE of the training data (47.4%). The modelled mean depths explained 46.3% of the variance in reference mean depths and 91.6% of the variance in reference volumes, but underestimated the total volume of lakes in the validation data by 16.5% (predicted = 1,503 km<sup>3</sup>, observed = 1,800 km<sup>3</sup>). This can be attributed primarily to the abundance of deep alpine lakes in the validation data set, which the model tends to underestimate. When considering only those 190 lakes outside of North America and Europe for which validation data were available, and which were under-represented in the training data, the SMAPE increased to 64.4% and the R<sup>2</sup> of the regressions for mean depth and volume were 66.8% and 86.5%, respectively. Here the model underestimated volume by 28.3% (predicted = 549 km<sup>3</sup>, observed = 765 km<sup>3</sup>), again largely due to underestimation of deep lakes in alpine areas of Argentina and the South Island of New Zealand.

**Volume calculations.** We used Model 5 and its regression coefficients as shown in Supplementary Table 4 to predict the log<sub>10</sub> of mean depths for all 1.42 million natural lakes smaller than 500 km<sup>2</sup> in the HydroLAKES database. Lake mean depth was then obtained by back-transforming the log<sub>10</sub> of the predicted mean depths and adding a bias-correction for log–log regressions as suggested by Ferguson<sup>65</sup> and applied by Heathcote *et al.*<sup>34</sup>:

$$D_{\text{bias-corrected}} = 10^{\log_{10}\bar{D}} \times \exp(2.65 \times s^2) \quad (1)$$

where D<sub>bias-corrected</sub> is the bias-corrected predicted mean depth, log<sub>10</sub> D̄ is the uncorrected predicted mean depth calculated by applying Model 5, s<sup>2</sup> is the residual variance from Model 5 and 2.65 is a constant.

**Extrapolation of smaller lakes.** To extrapolate the number, area and volume of lakes smaller than 10 ha in size, we fitted a line to the full empirical cumulative distribution of all natural lake polygons that are contained in the HydroLAKES database using least-square regression. We found that a power law model, specifically the Pareto distribution, provides a satisfactory fit for both surface area and volume (Fig. 3). Least-square regression was shown to result in minimal biases when applied to a full empirical cumulative distribution with large amounts of data and an appropriate threshold in the lower and upper tail<sup>66,67</sup>. We thus truncated our data at the lower end at 0.35 km<sup>2</sup> and 0.0005 km<sup>3</sup>, respectively, to exclude increasingly incomplete data records; and at the upper end at 10,000 km<sup>2</sup> and 1,000 km<sup>3</sup>, respectively, to exclude increasingly random large lakes from the statistical assessment.

As it has been suggested that other models, in particular a log-normal distribution, could equally describe lower-truncated lake size distributions<sup>68</sup>, we also fitted log-normal and exponential distribution functions to the empirical data using maximum-likelihood techniques. However, the power law (Pareto) distribution model provided a superior fit to the data.

Using surface area as the size criterion, the trend found in the Pareto model indicates that the global distribution of natural lakes is scale-invariant, that is, there is a 10-fold increase in the number of lakes from one logarithmic size class to the next. In fact, the coefficient of the linear fit is nearly equal to −1 (Fig. 3a):

$$N = 164,698 A^{-1.05434} \quad (2)$$

where N is the number of lakes larger than area A (km<sup>2</sup>) (R<sup>2</sup> = 0.999). We find a similar relationship when using volume as the size class criterion (Fig. 3b):

$$N = 1,596.9 V^{-0.84038} \quad (3)$$

where N is the number of lakes larger than volume V (km<sup>3</sup>) (R<sup>2</sup> = 0.999). A flattening of the curve is observed for lakes below 0.0005 km<sup>3</sup> due to the truncation of our lake data set at 10 ha.

Based on these scaling relationships, we derived estimates for the number of smaller lakes as well as their total surface area and volume using Pareto class means as suggested by Downing *et al.*<sup>31</sup> To curtain uncertainty, and given that McDonald *et al.*<sup>37</sup> found a departure from the power law distribution for lake sizes below ~0.001 km<sup>2</sup>, we capped our extrapolation at a lake size of 0.01 km<sup>2</sup> to avoid a deterioration of the results. We combined all existing records contained in HydroLAKES for lakes ≥0.35 km<sup>2</sup> and ≥0.0005 km<sup>3</sup>, respectively, with the extrapolated results for lakes below these thresholds.

**Hydraulic residence time calculation and validation.** We calculated lake residence time as the ratio between lake volume and the estimated discharge rate at the lake's pour point, assuming that all lakes are well-mixed and that evaporation and seepage are negligible. For lakes with high evaporation rates or losses to groundwater, the flux term would increase and residence times would become shorter, while lakes that do not destratify seasonally, such as Lake Baikal or
Tanganyika, may contain much older water than the residence time suggests. Nevertheless, the residence time as calculated here—also called 'theoretical water residence time' by some limnologists due to the described shortcomings—is a widely accepted proxy⁴⁵. Despite inconsistencies in these definitions, we achieved a good correlation ($R^2$ = 0.57) between our modelled residence times and documented values of 374 natural lakes provided by the European Waterbase⁶⁹ and Kalff⁴⁵.

**Data availability.** The HydroLAKES database is offered for free for scientific and educational applications at http://www.hydrosheds.org. All figures were produced with data from the HydroLAKES database. All data used in the development of the HydroLAKES polygon database are publicly available and cited in Supplementary Table 3. All data used in the training and validation of the model of lake depth are listed in Supplementary Table 2. Data from public sources used in the training and validation of the model are available from the corresponding author upon request. Some data that support the findings of this study are subject to availability restrictions (for example, Ontario Ministry of Natural Resources; Quebec lake depth data³⁴) and were used under license for the current study, and so are not publicly available.

## References

1. Muller Schmied, H. *et al.* Sensitivity of simulated global-scale freshwater fluxes and storages to input data, hydrological model structure, human water use and calibration. *Hydrol. Earth. Syst. Sci.* **18**, 3511–3538 (2014).

2. Bastviken, D., Tranvik, L. J., Downing, J. A., Crill, P. M. & Enrich-Prast, A. Freshwater methane emissions offset the continental carbon sink. *Science* **331**, 50 (2011).

3. Downing, J. A. *et al.* Sediment organic carbon burial in agriculturally eutrophic impoundments over the last century. *Global Biogeochem. Cycles* **22**, GB1018 (2008).

4. Balsamo, G. *et al.* On the contribution of lakes in predicting near-surface temperature in a global weather forecasting model. *Tellus A* **64**, 15829 (2012).

5. Harrison, J. A., Frings, P. J., Beusen, A. H. W., Conley, D. J. & McCrackin, M. L. Global importance, patterns, and controls of dissolved silica retention in lakes and reservoirs. *Global Biogeochem. Cycles* **26**, GB2037 (2012).

6. Nöges, T. Relationships between morphometry, geographic location and water quality parameters of European lakes. *Hydrobiologia* **633**, 33–43 (2009).

7. Staehr, P. A., Baastrup-Spohr, L., Sand-Jensen, K. & Stedmon, C. Lake metabolism scales with lake morphometry and catchment conditions. *Aquat. Sci.* **74**, 155–169 (2012).

8. Dodson, S. I., Arnott, S. E. & Cottingham, K. L. The relationship in lake communities between primary productivity and species richness. *Ecology* **81**, 2662–2679 (2000).

9. Post, D. M., Pace, M. L. & Hairston, N. G. Ecosystem size determines food-chain length in lakes. *Nature* **405**, 1047–1049 (2000).

10. De Graaf, G., Bartley, D., Jorgensen, J. & Marmulla, G. The scale of inland fisheries, can we do better? Alternative approaches for assessment. *Fisheries Manag. Ecol.* **1**, 64–70 (2012).

11. Winslow, L. A., Read, J. S., Hanson, P. C. & Stanley, E. H. Does lake size matter? Combining morphology and process modeling to examine the contribution of lake classes to population-scale processes. *Inland Waters* **5**, 7–14 (2015).

12. Seekell, D. A., Carr, J. A., Gudas, C. & Karlsson, J. Upscaling carbon dioxide emissions from lakes. *Geophys. Res. Lett.* **41**, 7555–7559 (2014).

13. Lvovich, M. I. in *Proceedings of the Reading Symposium (World Water Balance) Vol. 2*, Publ. 93, 401–415 (IASH-UNESCO-WMO, Gentbrugge-Paris-Geneva, 1970).

14. L'vovich, M. I. *World Water Resources and Their Future* (American Geophysical Union, 1979).

15. Nace, R. L. *Are We Running Out of Water? Report No. 536* (US Geological Survey, 1967).

16. Tamrazyan, G. Total lake water resources of the planet. *Bull. Geol. Soc. Finland* **46**, 23–27 (1974).

17. Meybeck, M. in *Physics and Chemistry of Lakes* (eds Lerman, A., Imboden, D. M. & Gat, J. R.) 1–35 (Springer, 1995).

18. Ryanzhin, S. V., Subetto, D. A., Kochkov, N. V., Akhmetova, N. S. & Veinmeister, N. V. Polar lakes of the World: current data and status of investigations. *Water Resour.* **37**, 427–436 (2010).

19. Shiklomanov, I. A. & Rodda, J. C. *World Water Resources at the Beginning of the Twenty-First Century* Vol. 13 (Cambridge University Press, 2003).

20. Downing, J. A. Emerging global role of small lakes and ponds: little things mean a lot. *Limnetica* **1**, 9–24 (2010).

21. Holgerson, M. A. & Raymond, P. A. Large contribution to inland water CO₂ and CH₄ emissions from very small ponds. *Nat. Geosci.* **9**, 222–226 (2016).

22. Smith, S., Renwick, W., Bartley, J. & Buddemeier, R. Distribution and significance of small, artificial water bodies across the United States landscape. *Sci. Total Environ.* **299**, 21–36 (2002).

23. Weyhenmeyer, G. A. in *European Large Lakes Ecosystem Changes and Their Ecological and Socioeconomic Impacts* (eds Nöges, T. *et al.*) 105–110 (Springer, 2008).

24. Winslow, L. A., Read, J. S., Hanson, P. C. & Stanley, E. H. Lake shoreline in the contiguous United States: quantity, distribution and sensitivity to observation resolution. *Freshw. Biol.* **59**, 213–223 (2014).

25. Dost, R. J. J. & Mannaerts, C. M. M. in *Proceedings of the ESRI International User Conference: GIS, Geography in action*, August 4-8 (Redmond, 2008).

26. Goodman, J. A., Purkis, S. J. & Phinn, S. R. *Coral Reef Remote Sensing: A Guide for Mapping, Monitoring and Management* (Springer Science & Business Media, 2013).

27. Legleiter, C. J., Roberts, D. A. & Lawrence, R. L. Spectrally based remote sensing of river bathymetry. *Earth. Surf. Proc. Land.* **34**, 1039–1059 (2009).

28. Lehner, B. *et al.* High-resolution mapping of the world's reservoirs and dams for sustainable river-flow management. *Front. Ecol. Environ.* **9**, 494–502 (2011).

29. Burke, L. *et al.* *Pilot Analysis of Global Ecosystems: Coastal Ecosystems* (World Resources Institute, 2000).

30. Oki, T. & Kanae, S. Global hydrological cycles and world water resources. *Science* **313**, 1068–1072 (2006).

31. Downing, J. A. *et al.* The global abundance and size distribution of lakes, ponds, and impoundments. *Limnol. Oceanogr.* **51**, 2388–2397 (2006).

32. O'Sullivan, P. & Reynolds, C. S. *The Lakes Handbook: Limnology and Limnetic Ecology* Vol. 1 (John Wiley & Sons, 2008).

33. Sobek, S., Nisell, J. & Fölster, J. Predicting the volume and depth of lakes from map-derived parameters. *Inland Waters* **1**, 177–184 (2011).

34. Heathcote, A. J., del Giorgio, P. A., Prairie, Y. T. & Brickman, D. Predicting bathymetric features of lakes from the topography of their surrounding landscape. *Can. J. Fish. Aquat. Sci.* **72**, 1–8 (2015).

35. Branstrator, D. K. in *Lake Ecosystem Ecology: A Global Perspective* (ed. Likens, G. E.) 191–202 (Academic Press, 2009).

36. Hutchinson, G. E. *A Treatise on Limnology* Vol. 1 (Wiley, 1957).

37. McDonald, C. P., Rover, J. A., Stets, E. G. & Striegl, R. G. The regional abundance and size distribution of lakes and reservoirs in the United States and implications for estimates of global lake extent. *Limnol. Oceanogr.* **57**, 597–606 (2012).

38. Verpoorter, C., Kutser, T., Seekell, D. A. & Tranvik, L. J. A global inventory of lakes based on high-resolution satellite imagery. *Geophys. Res. Lett.* **41**, 6396–6402 (2014).

39. Emmerton, C. A., Lesack, L. F. W. & Marsh, P. Lake abundance, potential water storage, and habitat distribution in the Mackenzie River Delta, western Canadian Arctic. *Water Resour. Res* **43**, W05419 (2007).

40. Carroll, M. L., Townshend, J. R., DiMiceli, C. M., Noojipady, P. & Sohlberg, R. A. A new global raster water mask at 250 m resolution. *Int. J. Digital Earth* **2**, 291–308 (2009).

41. Gao, H. Satellite remote sensing of large lakes and reservoirs: from elevation and area to storage. *Wiley Interdisciplinary Reviews Water* **2**, 147–157 (2015).

42. Papa, F. *et al.* Interannual variability of surface water extent at the global scale, 1993–2004. *J. Geophys. Res.* **115**, D12111 (2010).

43. Jorgensen, S. E. The application of models to find the relevance of residence time in lake and reservoir management. *J. Limnol.* **62**, 16–20 (2003).

44. Korzoun, V. I. *et al.* in *Studies and Reports in Hydrology* (UNESCO, 1978).

45. Kalff, J. *Limnology: Inland Water Ecosystems* Vol. 592 (Prentice Hall, 2002).

46. Downing, J. A. Plenary lecture global limnology: Up-scaling aquatic services and processes to planet Earth. *Verh. Internat. Verein. Limnol.* **30**, 1149–1166 (2009).

47. Lehner, B., Verdin, K. & Jarvis, A. New global hydrography derived from spaceborne elevation data. *EOS Trans. Am. Geophys. Union* **89**, 93–94 (2008).

48. Goodman, K. J., Parker, S. M., Edmonds, J. W. & Zeglin, L. H. Expanding the scale of aquatic sciences: the role of the National Ecological Observatory Network (NEON). *Freshw. Sci.* **34**, 377–385 (2015).

49. Robinson, N., Regetz, J. & Guralnick, R. P. *EarthEnv-DEM90: A nearly-global, void-free, multi-scale smoothed, 90m digital elevation model from fused ASTER and SRTM data*. *ISPRS J. Photogramm. Remote Sens.* **87**, 57–67 (2014).

50. Horn, B. K. P. Hill shading and the reflectance map. *Proc. IEEE* **69**, 14–47 (1981).

51. Choulga, M., Kourzeneva, E., Zakharova, E. & Doganovsky, A. Estimation of the mean depth of boreal lakes for use in numerical weather prediction and climate modelling. *Tellus A* **66**, 21295 (2014).

52. Herdendorf, C. E. Large lakes of the world. *J. Great Lakes Res.* **8**, 379–412 (1982).

53. Döll, P., Kaspar, F. & Lehner, B. A global hydrological model for deriving water availability indicators: model tuning and validation. *J. Hydrol.* **270**, 105–134 (2003).

54. Lehner, B. & Grill, G. Global river hydrography and network routing: baseline data and new approaches to study the world's large river systems. *Hydrol. Process.* **27**, 2171–2186 (2013).

55. Slater, J. A. *et al.* The SRTM data 'finishing' process and products. *Photogramm. Eng. Remote Sens.* **72**, 237–247 (2006).
56. Natural Resources Canada. CanVec Hydrography: Waterbody Features. Version 12.0. Data. Available at ftp://ftp2.cits.rncan.gc.ca/pub/canvec/. Accessed May 2013.

57. Hayashi, M. & Van der Kamp, G. Simple equations to represent the volume-area-depth relations of shallow wetlands in small topographic depressions. *J. Hydrol.* **237**, 74–85 (2000).

58. Ryanzhin, S. V. New estimates for global surface area and volume of natural world lakes. *Doklady Earth Sciences c/c of Doklady-Akademija Nauk* **401**, 253–257 (2005).

59. Håkanson, L. & Peters, R. H. *Predictive Limnology: Methods for Predictive Modelling* (SPB Academic Pub.; Distributors for the USA and Canada, Kugler Publications, 1995).

60. Hollister, J. W. & Milstead, W. B. Using GIS to estimate lake volume from limited data. *Lake Reserv. Manage* **26**, 194–199 (2010).

61. Hollister, J. W., Milstead, W. B. & Urrutia, M. A. Predicting maximum lake depth from surrounding topography. *PLoS ONE* **6**, e25764 (2011).

62. Pistocchi, A. & Pennington, D. European hydraulic geometries for continental SCALE environmental modelling. *J. Hydrol.* **329**, 553–567 (2006).

63. Belsley, D. A., Kuh, E. & Welsch, R. E. *Regression Diagnostics: Identifying Influential Data and Sources of Collinearity* (John Wiley & Sons, 2005).

64. Fox, J. & Weisberg, S. *An R Companion to Applied Regression* (Sage, 2010).

65. Ferguson, R. I. River loads underestimated by rating curves. *Water Resour. Res.* **22**, 74–76 (1986).

66. Clauset, A., Shalizi, C. R. & Newman, M. E. J. Power-law distributions in empirical data. *SIAM Rev.* **51**, 661–703 (2009).

67. White, E. P., Enquist, B. J. & Green, J. L. On estimating the exponent of power-law frequency distributions. *Ecology* **89**, 905–912 (2008).

68. Seekell, D. A. & Pace, M. L. Does the Pareto distribution adequately describe the size-distribution of lakes? *Limnol. Oceanogr.* **56**, 350–356 (2011).

69. European Environment Agency. ECRINS (European Catchments and Rivers Network System): Lakes. Version 1 (2012).

70. Lehner, B. & Döll, P. Development and validation of a global database of lakes, reservoirs and wetlands. *J. Hydrol.* **296**, 1–22 (2004).

## Acknowledgements

This work was supported by funding provided by the Natural Sciences and Engineering Research Council of Canada (NSERC) as part of the Canadian Network for Aquatic Ecosystem Services CNAES (NETGP 417353-11) and through two NSERC

Undergraduate Student Research Awards; and by McGill University, Montreal, Canada. We thank the data-contributing organizations and scholars that provided their individual databases for use in the model development and/or incorporation in HydroLAKES. In addition, we thank the many colleagues and individuals who have dedicated time and energy to the study of lakes and contributed in various ways to this research.

## Author contributions

M.L.M. and B.L. led the work, designed the study, methodology and modelling approach, and collected data; M.L.M. performed the analysis and interpretation, and wrote the initial draft of the manuscript; M.L.M. and B.L. revised the manuscript; G.G. made substantial contributions to the study design and methodology; I.N. compiled, corrected and unified data sets for the creation of the HydroLAKES polygon database; O.S. performed a pilot study designed by B.L.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Messager, M. L. et al. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. *Nat. Commun.* **7**, 13603 (doi: 10.1038/ncomms13603 (2016)).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2016