# Lake-TopoCat: a global lake drainage topology and catchment database

Md Safat Sikder, Jida Wang, George H. Allen, Yongwei Sheng, Dai Yamazaki, Chunqiao Song, Meng Ding, Jean-François Cretaux, Tamlin M. Pavelsky

► To cite this version:

Md Safat Sikder, Jida Wang, George H. Allen, Yongwei Sheng, Dai Yamazaki, et al. Lake-TopoCat: a global lake drainage topology and catchment database. Earth System Science Data, 2023, 10.5194/essd-15-3483-2023. hal-04720449

HAL Id: hal-04720449

https://hal.science/hal-04720449v1

Submitted on 4 Oct 2024

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>

# Lake-TopoCat: a global lake drainage topology and catchment database

Md Safat Sikder<fnref n="1" />, Jida Wang<fnref n="2" />,<fnref n="1" />, George H. Allen<fnref n="3" />, Yongwei Sheng<fnref n="4" />, Dai Yamazaki<fnref n="5" />, Chunqiao Song<fnref n="6" />,<fnref n="7" />, Meng Ding<fnref n="1" />, Jean-François Crétaux<fnref n="8" />, and Tamlin M. Pavelsky<fnref n="9" />

**Correspondence:** Md Safat Sikder (msikder@ksu.edu) and Jida Wang (jidawang@ksu.edu)

Received: 11 December 2022 – Discussion started: 18 January 2023

Revised: 10 May 2023 – Accepted: 17 June 2023 – Published: 8 August 2023

**Abstract.** Lakes and reservoirs are ubiquitous across global landscapes, functioning as the largest repository of liquid surface freshwater, hotspots of carbon cycling, and sentinels of climate change. Although typically considered lentic (hydrologically stationary) environments, lakes are an integral part of global drainage networks. Through perennial and intermittent hydrological connections, lakes often interact with each other, and these connections actively affect water mass, quality, and energy balances in both lacustrine and fluvial systems. Deciphering how global lakes are hydrologically interconnected (or the so-called “lake drainage topology”) is not only important for lake change attribution but also increasingly critical for discharge, sediment, and carbon modeling. Despite the proliferation of river hydrography data, lakes remain poorly represented in routing models, partially because there has been no global-scale hydrography dataset tailored to lake drainage basins and networks. Here, we introduce the global Lake drainage Topology and Catchment database (Lake-TopoCat), which reveals detailed lake hydrography information with careful consideration of possible multifurcation. Lake-TopoCat contains the outlet(s) and catchment(s) of each lake; the interconnecting reaches among lakes; and a wide suite of attributes depicting lake drainage topology such as upstream and downstream relationship, drainage distance between lakes, and a priori drainage type and connectivity with river networks. Using the HydroLAKES v1.0 (Messager et al., 2016) global lake mask, Lake-TopoCat identifies ~1.46 million outlets for ~1.43 million lakes larger than 10 ha and delineates 77.5 × 10⁶ km<fnref n="2" /> of lake catchments covering 57% of the Earth’s landmass except Antarctica. The global lakes are interconnected by ~3 million reaches, derived from MERIT Hydro v1.0.1 (Yamazaki et al., 2019), stretching a total distance of ~10 × 10⁶ km, of which ~80% are shorter than 10 km. With such unprecedented lake hydrography details, Lake-TopoCat contributes towards a globally coupled lake–river routing model. It may also facilitate a variety of limnological applications such as attributing water quality from lake scale to basin scale, tracing inter-lake fish migration due to changing climate, monitoring fluvial–lacustrine connectivity, and improving estimates of terrestrial carbon fluxes. Lake-TopoCat is freely accessible at https://doi.org/10.5281/zenodo.7916729 (Sikder et al., 2023).

<footnote n="1">¹Department of Geography and Geospatial Sciences, Kansas State University, Manhattan, KS, USA</footnote>

<footnote n="2">²Department of Geography and Geographic Information Science, University of Illinois Urbana-Champaign, Urbana, IL, USA</footnote>

<footnote n="3">³Department of Geosciences, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA</footnote>

<footnote n="4">⁴Department of Geography, University of California, Los Angeles, CA, USA</footnote>

<footnote n="5">⁵Institute of Industrial Science, The University of Tokyo, Tokyo, Japan</footnote>

<footnote n="6">⁶Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing, China</footnote>

<footnote n="7">⁷University of Chinese Academy of Sciences, Nanjing (UCASNJ), Nanjing, China</footnote>

<footnote n="8">⁸Laboratoire d'Études en Géophysique et Océanographie Spatiales (LEGOS), Centre National d'Études Spatiales (CNES), Toulouse, France</footnote>

<footnote n="9">⁹Department of Earth, Marine and Environmental Sciences, University of North Carolina, Chapel Hill, NC, USA</footnote>
# 1 Introduction

Natural lakes, ponds, and human-made reservoirs, hereafter "lakes", store the largest amount of liquid freshwater on Earth's surface (Oki and Kanae, 2006; Abbott et al., 2019). Although widely perceived as water stores and lentic systems, lakes are often connected to each other through river networks and are inherent components of the global drainage system (Fergus et al., 2017; Gardner et al., 2019). The water balance in a lake and its water quality (e.g., turbidity, dissolved oxygen, acidity, and temperature) reflect climatic, hydrological, geologic, and land use characteristics of the local drainage catchment (Yang et al., 2022). The limnological properties of one lake also affect others via the transfers of water mass (e.g., Huziy and Sushama, 2017), sediments, nutrients (e.g., Stieglitz et al., 2003), and energy along connecting rivers. Since lakes sequester a large amount of organic and inorganic carbon (Tranvik et al., 2009; Mendonça et al., 2017), deciphering how global lakes are connected through drainage networks will complement routing models, which currently emphasize rivers (Leibowitz et al., 2018; Liu et al., 2022), to better constrain the terrestrial carbon cycling (Cardille et al., 2007). The drainage paths between lakes, including their distances, gradients, connectivity, and climatic zonation, also determine the migration pattern for lacustrine species, which is crucial for monitoring ecosystem health and services under climate change (Woolway et al., 2020). For these reasons, a hydrography dataset tailored to global lakes, which provides fine spatial details of the boundary of each lake catchment and how lakes are topologically interconnected, is overdue and has the potential for high hydrological and ecological significance.

In recent years, significant improvements have been achieved in global lake mapping and inventories. A few prominent examples are HydroLAKES (v1.0), which inventories 1.4 million lakes larger than 10 ha (Messager et al., 2016); GLAKES comprising 3.4 million lakes with maximum surface area larger than 3 ha (Pi et al., 2022); and the UCLA circa-2000 and circa-2015 Global Lake Inventories with more than 9 million lake polygons larger than 0.4 ha (Sheng et al., 2016). These water masks, unfortunately, emphasize individual lake entities and provide little to no metadata depicting drainage hydrography (such as catchment boundaries and hydrological connectivity) among the lakes. The lack of drainage hydrography has restricted application of these refined lake inventories to watershed-scale processes and has left lake–river integration a challenging task.

Only a few attempts have been made to integrate global lakes into drainage networks. In the customized format of HydroBASINS (Lehner and Grill, 2013), for example, lake polygons were directly clipped into the river sub-basin polygons, in order to register the lakes to the topological (upstream and downstream) coding structure of HydroBASINS.

Nevertheless, since the original sub-basins in HydroBASINS were generated without consideration of lake presence, there was no guarantee that the sub-basin pour points are aligned with the lake inlets or outlets. This inconsistency leads to cases where a sub-basin polygon contains multiple lakes or a lake polygon contains the pour point of the most downstream intersecting sub-basin. In such cases, the drainage topology among lakes is incompletely depicted, and the intersecting sub-basins may underrepresent or overshoot the spatial domains of the lake catchments. Similar issues were also acknowledged in the Global Lake area, Climate, and Population dataset (GLCP), where HydroBASINS was also used to match lakes to the surrounding sub-basins (Meyer et al., 2020). Han et al. (2020, 2023) developed integrated lake–river routing products for Canada and North America using their “BasinMaker” GIS (geographic information system) tool, in which they redefine the flow direction to derive consistent lake–river topological connections.

In HydroLAKES (Messager et al., 2016) and the most recent LakeATLAS dataset (Lehner et al., 2022), the outlet and the upstream drainage area associated with each lake were derived from the 15 arcsec HydroSHEDS hydrography dataset (Lehner et al., 2008). Their methods of lake outlet delineation, however, excluded the possibility of lake bifurcation or multifurcation (i.e., a lake flowing out in two or more directions), and the 15 arcsec resolution (about 450 m at the Equator) of HydroSHEDS might be too coarse to derive reliable catchment boundaries for many small lakes in HydroLAKES (with a minimum size of 10 ha). In addition, lake catchments in HydroLAKES or LakeATLAS were provided with only area values not geometric boundaries, meaning that spatially explicit applications of such catchment information remain inconvenient. In fact, despite the proliferation of river basin datasets (e.g., Lehner and Grill, 2013; Lin et al., 2021), there has been no global-scale catchment data tailored specifically for lakes. To our knowledge, only a few regional lake watershed datasets are available so far, and examples are the National Hydrography Dataset Plus Version 2 (NHDPlusV2) (McKay et al., 2012), the Lake-Catchment (Lake-Cat) dataset (Hill et al., 2018), and the LAGOS-US LOCUS dataset (Cheruvelil et al., 2021) for the USA; the Comprehensive Data set for China's Lake Basins (CODCLAB) for 767 large lakes (> 10 km²) in China (T. Chen et al., 2022); the lake–river routing product of Canada (Han et al., 2020) and North America (Han et al., 2023); and several lake catchment datasets for the endorheic Tibetan Plateau (e.g., Liu et al., 2020, 2021; Yan et al., 2019).

More detailed configurations of lake drainage connectivity are also limited to regional scales. For the USA, Schmadel et al. (2018) and Gardner et al. (2019) integrated lakes into river networks as depicted in the NHDPlusV2. In these studies, lakes were considered part of the drainage system only if the lake polygon was directly intersected by river channels.
While this is reasonable in many cases, hydrological connec-
tivity among lakes can be established through intermittent
channels, ephemeral flow (such as via “fill and spill”), and
subsurface flow (Leibowitz et al., 2016; McDonnell et al.,
2021), which are not always visible or channelized in the
NHDPlusV2 vector networks. Alternatively, multitemporal
spectral images have been applied to better capture the vari-
ability in intermittent connectivity. For instance, Vanderhoof
et al. (2016, 2017) used Landsat images to examine the pat-
terns of wetland connections in relation to wetland arrange-
ments and surface water expansion across the Prairie Pot-
hole region of the USA. Tan et al. (2019) used Landsat and
MODIS imagery to quantify the surface water connectivity to
understand the complex surface water dynamics of Poyang
Lake in China. More recently, Dolan et al. (2021) applied
Landsat-observed sediments in water to understand the vari-
ation in lake-to-channel connectivity and the impact on ice
phenology within the Colville Delta, Alaska. The complete-
ness and consistency of the observed lake–river connectivity
are, however, affected by the image quality and resolution.
This is particularly problematic when channel visibility is
obstructed by riparian canopies. In addition, although higher-
resolution images may improve the mapping of lake–river
connections (e.g., Wu and Lane, 2017), optical images alone
are not yet sufficient to decipher complex drainage topology
among global lakes.

We argue that, while lake–river connectivity may vary
through time, a prerequisite for monitoring such important
dynamics is to sort out the drainage topology (i.e., upstream
and downstream relationships) among global lakes based on
their potential connecting paths. This lake topology, albeit
temporarily static, ensures that all lakes on the continental
surface are registered to the drainage system, and this sys-
tem will serve as a useful a priori network for examining the
temporal variability in surface water connection. This topol-
ogy will also facilitate a more thorough inclusion of lakes
into global hydrological models. For example, Bowling and
Lettenmaier (2010) showed that lakes in the Arctic region
can store up to 80% of the snowmelt water each year, ul-
timately reducing the spring peak of river flows in that re-
gion. Considering the unique roles of lakes, numerous studies
have incorporated lakes into hydrological and/or hydrody-
namic models to improve the accuracy of simulated stream-
flow (e.g., Han et al., 2020; Lin et al., 2015), water temper-
ature (Tokuda et al., 2021), and regional climate (Huziy and
Sushama, 2017). Lakes can be registered manually into the
routing networks in the case of small areas, where the num-
ber of lakes is manageable (e.g., Lin et al., 2015). On the
other hand, usually only large lakes have been resolved in
continental- or global-scale modeling networks (e.g., Tokuda
et al., 2021), due to the constraints of model and hydrography
data resolutions, computing power, and additional workload.

To address these gaps, we introduce the global Lake
drainage Topology and Catchment database or Lake-
TopoCat. Lake-TopoCat offers an unprecedented detail

of global lake hydrography, including spatially resolved
perimeters of the lake catchments; drainage outlets from
the lakes; perennial, intermittent, and potential drainage
reaches between the lakes; and, importantly, the topolog-
ical relationships among the lake, catchment, and reach
structures. To tackle the above-described challenges, we ap-
plied the MERIT Hydro high-resolution raster hydrography
maps (Yamazaki et al., 2019). Our algorithm also included
a mechanism to capture possible lake bifurcation. This pa-
per describes the initial version of Lake-TopoCat, which
was constructed for nearly all 1.4 million lakes in Hydro-
LAKES v1.0; however, our algorithm is generic and ap-
plicable to other lake masks (see the “Database availabil-
ity and versioning” section). We name the initial version
“HydroLAKES-TopoCat v1.1”, where the first part of the
name (HydroLAKES) indicates the input lake mask, and the
second part of the name (TopoCat v1.1) reflects the ver-
sion of the algorithm and other specifics used to construct
TopoCat, and this naming convention will be followed for
future versioning and updates. For conciseness, we use the
generic name Lake-TopoCat or TopoCat to refer to this ini-
tial version in the remaining text. For user convenience, the
following sections first provide an overview of the compo-
nents, structure, and functionality of Lake-TopoCat (Sect. 2),
before elaborating detailed input data sources and methodol-
ogy (Sect. 3). They are then followed by quality assessment
and technical validation (Sect. 4) and database availability
(Sect. 5). We conclude with a discussion on Lake-TopoCat’s
potential applications and limitations (Sect. 6).

## 2 Data description and structure

Lake-TopoCat consists of five feature components, each with
multiple attributes depicting lake drainage relationships. The
five features are (1) lake boundaries (polygons), (2) lake
outlets (points), (3) unit catchment boundaries (polygons)
defining the drainage areas between cascading (i.e., imme-
diately upstream and downstream; see definition in Sect. 2.2
and 2.3) lake outlets, (4) inter-lake reaches (lines) defining
the drainage networks that connect the lake outlets to the in-
land sinks or the ocean, and (5) lake-network basins (poly-
gons) that define the entirety of the drainage area containing
each inter-lake network (i.e., a complete basin from the head-
water to an inland sink or the ocean for all basins contain-
ing lakes). An example of the feature components is given
for a focal region (about 6300 km²) in the Canadian Shield
(Fig. 1). The attribute tables for each of the feature compo-
nents are explained in Table 1. Each of the components and
their associated attributes are described in further detail be-
low.

2.1 Lake boundaries

Lake boundaries in TopoCat are geometrically the same as
the HydroLAKES lake polygons, except that the former in-
<table><caption>**Table 1.** Attribute definitions in each of the five feature components of TopoCat. Here, dd, m a.s.l., and n/a are abbreviations of decimal degrees, meters above sea level, and not applicable, respectively.</caption><thead><tr><th>Feature</th><th>Attribute</th><th>Description</th><th>Unit</th></tr></thead><tbody><tr><td rowspan="15">Lakes:<br/>Individual lake polygons<br/>as in HydroLAKES. File<br/>name: Lakes_pfaf_xx, where<br/>"pfaf_xx" indicates the IDs<br/>of the 68 Pfafstetter level-2<br/>basins for data organization<br/>(see Sect. 3.2).</td><td>Hylak_id</td><td>HydroLAKES ID for this lake</td><td>n/a</td></tr><tr><td>Lake_area</td><td>Lake water surface area</td><td>km<sup>2</sup></td></tr><tr><td>Outlet_n</td><td>Count of outlets of this lake (&gt; 1 indicating multifurcation)</td><td>n/a</td></tr><tr><td>D_hylak_id</td><td>List of IDs of the next (i.e., directly connected) downstream lakes</td><td>n/a</td></tr><tr><td>D_lake_n</td><td>Count of the next downstream lakes</td><td>n/a</td></tr><tr><td>D_lak_ntot</td><td>Count of all downstream lakes to the sink or ocean, excluding this lake</td><td>n/a</td></tr><tr><td>U_lake_n</td><td>Count of the next upstream lakes</td><td>n/a</td></tr><tr><td>U_lak_ntot</td><td>Count of all upstream lakes from the headwater, excluding this lake</td><td>n/a</td></tr><tr><td>Cat_a_lake</td><td>Area of the total upstream catchment(s) (accumulative from the headwater) for this lake; in the case of multifurcation, this value aggregates the areas of the upstream catchments for all outlets of this lake</td><td>km<sup>2</sup></td></tr><tr><td>Lake_type</td><td>Drainage type of this lake in relation to the inter-lake reach network: headwater, flow-through, endorheic (terminal), or coastal</td><td>n/a</td></tr><tr><td>Lake_order</td><td>Strahler order of the lake in the inter-lake reach network</td><td>n/a</td></tr><tr><td>Laktyp_mhv</td><td>Drainage type of this lake in relation to rivers in MERIT Hydro-Vector: isolated, headwater, inflow-headwater, flow-through, endorheic (terminal), or coastal</td><td>n/a</td></tr><tr><td>Lperm_glcp</td><td>Average percentage of the lake area being permanent water, derived from the GLCP dataset; -99.99 indicates null value</td><td>n/a</td></tr><tr><td>Basin_id</td><td>ID of the lake-network basin (see below) this lake belongs to</td><td>n/a</td></tr><tr><td rowspan="15">Outlets:<br/>Outlet or pour points from<br/>each of the lakes. There are<br/>multiple outlets from a multi-<br/>furcation lake. File name:<br/>Outlets_pfaf_xx.</td><td>Outlet_id</td><td>ID of this lake outlet, which is different from Hylak_ID as there can be more than one outlet for a multifurcation lake</td><td>n/a</td></tr><tr><td>Hylak_id</td><td>ID of the associated HydroLAKES lake</td><td>n/a</td></tr><tr><td>Outlet_lat</td><td>Latitude at this outlet</td><td>dd</td></tr><tr><td>Outlet_lon</td><td>Longitude at this outlet</td><td>dd</td></tr><tr><td>Outlet_elv</td><td>Elevation at this outlet (based on hydrologically adjusted MERIT DEM, as for other elevation values)</td><td>m a.s.l.</td></tr><tr><td>D_out_id</td><td>ID of the next downstream outlet</td><td>n/a</td></tr><tr><td>D_hylak_id</td><td>ID of the next downstream lake</td><td>n/a</td></tr><tr><td>D_lak_ntot</td><td>Count of all downstream lakes to the sink or ocean, excluding this lake</td><td>n/a</td></tr><tr><td>D_reach_id</td><td>ID of the connecting downstream reach</td><td>n/a</td></tr><tr><td>D_slope</td><td>Hydraulic gradient (slope) from this outlet to the next downstream outlet or sink</td><td>n/a</td></tr><tr><td>D_distance</td><td>Drainage distance from this outlet to the next downstream outlet or sink</td><td>m</td></tr><tr><td>D_dst_sink</td><td>Total drainage distance to the sink or ocean</td><td>m</td></tr><tr><td>U_lake_n</td><td>Count of the next upstream lakes</td><td>n/a</td></tr><tr><td>U_dist_avg</td><td>Mean distance from the next upstream outlets</td><td>m</td></tr><tr><td>U_dist_min</td><td>Minimum distance from the next upstream outlets</td><td>m</td></tr><tr><td>U_dist_max</td><td>Maximum distance from the next upstream outlets</td><td>m</td></tr><tr><td>U_lak_ntot</td><td>Count of all upstream lakes from the headwater, excluding this lake</td><td>n/a</td></tr><tr><td>Cat_area</td><td>Area of the associated unit catchment</td><td>km<sup>2</sup></td></tr><tr><td>Cat_a_tot</td><td>Area of the entire upstream drainage basin, i.e., from the headwater to this outlet</td><td>km<sup>2</sup></td></tr><tr><td>Out_type</td><td>Drainage type of this lake outlet in relation to the inter-lake reach network: headwater, flow-through, terminal, or coastal</td><td>n/a</td></tr><tr><td>Out_order</td><td>Strahler order of the outlet in the inter-lake reach network</td><td>n/a</td></tr><tr><td>Basin_id</td><td>ID of the lake-network basin this outlet belongs to</td><td>n/a</td></tr></tbody></table>
<table><caption>**Table 1. Continued.**</caption><thead><tr><th>Feature</th><th>Attribute</th><th>Description</th><th>Unit</th></tr></thead><tbody><tr><td rowspan="7">Unit catchments:<br>Boundaries of the drainage area<br>between each lake outlet and<br>the next upstream outlet(s) or<br>the drainage area from the<br>headwater to the outlet if there<br>are no lakes upstream. File<br>name: Catchments_pfaf_xx.</td><td>Outlet_id</td><td>ID of the associated outlet, also used as the ID of the unit catchment</td><td>n/a</td></tr><tr><td>Hylak_id</td><td>ID of the associated lake</td><td>n/a</td></tr><tr><td>D_out_id</td><td>ID of the next downstream outlet</td><td>n/a</td></tr><tr><td>D_hylak_id</td><td>ID of the next downstream lake</td><td>n/a</td></tr><tr><td>Cat_area</td><td>Area of this unit catchment</td><td>km<sup>2</sup></td></tr><tr><td>Cat_type</td><td>Drainage type of this unit catchment in relation to the inter-lake reach network (same as in Outlet_type)</td><td>n/a</td></tr><tr><td>Basin_id</td><td>ID of the lake-network basin this unit catchment belongs to</td><td>n/a</td></tr><tr><td rowspan="11">Inter-lake reaches:<br>Drainage networks inter-<br>connecting the lake outlets to<br>the next downstream lake out-<br>lets, inland sinks, or the ocean.<br>File name: Reaches_pfaf_xx.</td><td>Reach_id</td><td>ID of this inter-lake reach</td><td>n/a</td></tr><tr><td>D_reach_id</td><td>ID of the next downstream reach</td><td>n/a</td></tr><tr><td>D_out_id</td><td>ID of the connected downstream outlet</td><td>n/a</td></tr><tr><td>D_hylak_id</td><td>ID of the connected downstream lake</td><td>n/a</td></tr><tr><td>U_out_id</td><td>ID of the connected upstream outlet</td><td>n/a</td></tr><tr><td>U_hylak_id</td><td>ID of the connected upstream lake</td><td>n/a</td></tr><tr><td>U_lak_ntot</td><td>Count of the total upstream lakes from the headwater</td><td>n/a</td></tr><tr><td>Start_elv</td><td>Elevation at the starting node of this reach</td><td>m a.s.l.</td></tr><tr><td>End_elv</td><td>Elevation at the ending node of this reach</td><td>m a.s.l.</td></tr><tr><td>Rch_length</td><td>Reach geodesic length</td><td>m</td></tr><tr><td>Rch_slope</td><td>Hydraulic gradient (slope) of this reach</td><td>n/a</td></tr><tr><td>Rchint_mhv</td><td>Proportion of this inter-lake reach (in % length) overlapped by river channels in MERIT Hydro-Vector, implying<br>the likelihood of this reach being more perennial or ephemeral</td><td>n/a</td></tr><tr><td rowspan="7">Lake-network basins:<br>Entire drainage basins from the<br>headwater to the ocean or in-<br>land sinks, defined by each<br>of the inter-lake reach net-<br>works. Lakes in each lake basin<br>are topologically related. File<br>name: Basins_pfaf_xx.</td><td>Rch_order</td><td>Strahler order of this reach</td><td>n/a</td></tr><tr><td>Basin_id</td><td>ID of the lake-network basin this reach belongs to</td><td>n/a</td></tr><tr><td>Basin_id</td><td>ID of the lake-network basin</td><td>n/a</td></tr><tr><td>Basin_type</td><td>Type of this lake-network basin: endorheic (draining to a terminal lake or an inland sink) or exorheic (draining<br>to the ocean)</td><td>n/a</td></tr><tr><td>Basin_area</td><td>Area of this lake-network basin</td><td>km<sup>2</sup></td></tr></tbody></table>

cludes new attributes informing on lake drainage relations. About 99.95% of the ~1.4 million HydroLAKES polygons were located within or on the MERIT Hydro dataset boundary. These lakes were used to construct TopoCat (Fig. 2a and b). For consistency, we used the original HydroLAKES lake IDs (“Hylak_id”) to index lakes and their associated topological attributes (Table 1). This allows users to link TopoCat features or attributes to HydroLAKES, as well as to other databases derived from HydroLAKES, such as GLCP (Meyer et al., 2020), LakeATLAS (Lehner et al., 2022), the Global Lake Evaporation Volume (GLEV) dataset (Zhao et al., 2022), and the GLObal lakes Bathymetry dataset (GLOBathy) (K hazaei et al., 2022).

The attributes associated with each lake feature include the counts of the next (i.e., first or directly connected) upstream and downstream lakes (“U_lake_n” and “D_lake_n”), the ID(s) of the next downstream lake(s) (“D_hylak_id”), the count of all upstream lakes which are connected in a cas-

cade network from the headwater to this lake (“U_lak_ntot”), and the count of all downstream lakes from this lake to the drainage sink or ocean (“D_lak_ntot”). To take into account lake bifurcation or multifurcation, we also report the count of outlets for each lake (“Outlet_n”; see Sect. 2.2), with values larger than 1 indicating possible multifurcation (see the example in Fig. 3). The attributes also include the area of the entire lake catchment (“Cat_a_lake”), which is defined as the upstream drainage area from the headwater to the lake (note this is different from the unit catchment areas; see Sect. 2.2). For a multifurcation lake, the lake catchment aggregates the upstream catchments associated with all lake outlets (see Sect. 2.2 and 2.3). More explanation for each of the attributes is given in Table 1.

TopoCat also provides the drainage type of each lake (hereafter lake type) based on the position of the lake in relation to the associated drainage network (Fig. 2a and b). Following the terminology in lake surface hydrologic posi-
<figure><img src="image_3.png" /><figcaption>**Figure 1.** An example of TopoCat in the Canadian Shield. This example covers an area of ~6300 km² with ~850 lakes represented. Displayed features are lake outlets, inter-lake reaches (with width representing the abundance or count of the upstream lakes), lake unit catchments, and lake boundaries (lake boundaries are the same as HydroLAKES; Messager et al., 2016). All lakes in this region belong to the same lake-network basin shown in the inset.</figcaption></figure>

tion (Martin and Soranno, 2006; illustrated in Fig. 2a and b), our considered lake types include (1) headwater if the lake spatially concurs with the river origin; (2) inflow-headwater if the lake is fed by a river without any upstream lakes; (3) flow-through if the lake is located on the pathway of a river drained from an upstream lake; (4) terminal or endorheic, meaning that the lake is located at an inland sink or terminal; (5) coastal, indicating the lake is immediately connected to the ocean (such as a lagoon) without any evident outflow river; and (6) isolated (or seepage) if the lake is disconnected from the surficial river network.

We acknowledge that the reliability of the configured lake type depends on the accuracy and completeness of the global river networks, and some of the “isolated” lakes may actually be connected through non-channelized surface runoff and/or subsurface flow (associated with surficial terrain). In addition, since lake–river connectivity can be intermittent, an isolated lake may also be connected back to the net-

work during high-flow conditions. For these reasons, we offer two a priori lake type attributes to allow for more flexibility (Table 1). The first is “Lake_type”, where lake drainage types were configured based on the drainage paths connecting the entirety of global lakes (hereafter “inter-lake connecting reaches”), which we delineated from the 3 arcsec MERIT Hydro hydrography dataset (Sect. 2.4). In other words, the inter-lake connecting reaches are composed of a set of global drainage networks, where drainage density is determined by the observed spatial variability of lake density. Lake_type ensures that all global lakes are topologically connected and eliminates any isolated and inflow-headwater lakes from lake drainage configuration. The numbers of identified headwater, flow-through, endorheic, and coastal lakes as Lake_type are 57.3%, 42.5%, 0.1%, and 0.1% of the total lakes, respectively. The other attribute is “Laktyp_mhv”. Lake types in the Laktyp_mhv attribute were configured based on MERIT Hydro-Vector (Lin et al., 2021), a high-resolution global river
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Global lakes and their drainage positions in TopoCat. (a) Lake types ("Lake_type" attribute) configured based on the hydrologic positions in the inter-lake reaches (see Sect. 2.4) of TopoCat. (b) Lake types ("Laktyp_mhv" attribute) configured based on the hydrologic position in global perennial and intermittent rivers provided in MERIT Hydro-Vector (Lin et al., 2021). Pie charts in both figures show the global lake area composition by drainage type. Here, lake boundaries are the same as HydroLAKES (Messager et al., 2016).</figcaption></figure>

network dataset which accounts for variable drainage density across the land surface (see Sect. 3.1.2 for more details). This attribute represents how global lakes are related to the perennial and intermittent rivers. About 4.5% of the global lakes are categorized as isolated, meaning that they are disconnected from any rivers in MERIT Hydro-Vector (Fig. 2b). The drainage types for the other lakes are the same as those in attribute Lake_type. The numbers of identified isolated,

headwater, inflow-headwater, flow-through, endorheic, and coastal lakes as Laktyp_mhv are 31.5%, 15.0%, 20.8%, 32.5%, 0.1%, and 0.1% of the total lakes, respectively.

The pie chart in Fig. 2a shows the global lake area composition according to the lake drainage types based on the Lake_type attribute. Here, flow-through lakes cover the largest proportion of the global lake surface area (~67.1%), although only ~42.5% of the global lakes were identified as
<figure><img src="image_5.png" /><figcaption>**Figure 3.** An example of a multifurcation lake (Horsetooth Reservoir in Colorado, USA) with four verified drainage outlets (lake boundaries are the same as HydroLAKES; Messager et al., 2016, and map background from Esri, Maxar, Earthstar Geographics, and the GIS User Community).</figcaption></figure>

flow-through. Usually, flow-through lakes are larger in size than any of the upstream types (i.e., isolated, headwater, and inflow-headwater). Therefore, covering the largest lake surface area with relatively low frequency was expected. Similarly, the pie chart in Fig. 2b shows the global lake area composition according to the lake drainage types based on Laktyp_mhv. Considering that isolated lakes have the potential to be connected to the other lakes, 15.5% of the isolated, 16.3% of the headwater, and 12.7% of the inflow-headwater lakes

from Laktyp_mhv become flow-through lakes in Lake_type, leading to a 30.7% increase in the flow-through lake frequency. The dominant proportion for flow-through lakes further accentuates the intrinsic roles of lakes in the global drainage system and the necessity of building this global-scale lake topology dataset. Endorheic or terminal lakes also cover a significant proportion of the global lake area (~21%), mainly due to large terminal saline lakes such as the Caspian Sea, the Great Salt Lake, and Lake Balkhash,
as well as the cluster of terminal lakes across the endorheic
Tibetan Plateau (Fig. 2a and b).

As many lakes are highly dynamic in inundation extent and level, so are their drainage topology and connectivity. While our drainage relations were based on the static Hydro-LAKES lake mask and should be treated as a priori topology, we provide another attribute, “Lperm_glcp”, to indicate the seasonal or intra-annual stability of water area in each lake, thus implying how the drainage relationship is subject to change. Lperm_glcp quantifies the average percentage of the lake area being permanently inundated, calculated using 6-year (2010–2015) statistics from the GLCP database (Meyer et al., 2020). A higher value of Lperm_glcp indicates that this lake is more stable or permanent, whereas a lower value indicates that the lake is more seasonal, unstable, or intermittent.

## 2.2 Lake outlets

The outlet features, in principle, represent the locations of the maximum flow accumulation pouring from each lake to each of its next downstream lakes or drainage sink. If a lake is terminal, the outlet refers to the drainage sink. Although the majority of lakes have one outlet only, lake bifurcation or even multifurcation does exist in both natural lakes and artificial reservoirs (see the example in Fig. 3). To provide a more realistic representation of global lake topology, we leveraged the detailed lake boundaries in HydroLAKES and the high-accuracy hydrography information in MERIT Hydro (see Sect. 3.1.2) to allow for possible multifurcation. The method is elaborated on in Sect. 3.2. However, it is important to note that only water body maps were used to develop the MERIT Hydro, without considering lake masks as an input. Therefore, in many lakes, flow directions within the lake might be influenced by adjacent topography and causing multiple outlets. Since multifurcation was considered, the count of outlets in TopoCat is larger than the count of lakes. We identified 1459201 outlets for 1426967 lakes, where 29190 lakes (~2% of the global lakes) drain to multiple destinations (Fig. 4).

The outlet attributes contain the topology information detailing how the outlets, lakes, and catchments are hydrologically connected to each other. They include the unique ID of the outlet ("Outlet_id"); the associated HydroLAKES lake ID ("Hylak_id"); and the drainage relationship with other outlets and lakes, such as the IDs of the next downstream outlets ("D_out_id"), the IDs of the next downstream lakes ("D_hylak_id"), and the count of the next upstream lakes ("U_lake_n").

The attributes also contain the area of the unit catchment associated with each outlet ("Cat_area"). Here a unit catchment (also see Sect. 2.3) is defined as the drainage area between any two cascading outlets (i.e., from an outlet to its next upstream outlet) or the drainage area from the headwater to this outlet if there is no lake upstream. In addition,

we also report the area of the entire upstream catchment (“Cat_a_tot”), defined as the drainage area from the headwater to this outlet, the count of lakes within the entire upstream catchment (“U_lak_ntot”), and the count of lakes within the entire downstream catchment to the inland sink or ocean (“D_lake_ntot”).

Similar to lake types, the type of each lake outlet ("Outlet_type") was assigned based on the drainage position in our delineated inter-lake reach network (Fig. 4). For lakes with only one outlet, the outlet type equals the corresponding lake type (Lake_type). For multifurcation lakes with disagreeing outlet types, the outlet type associated with the most downstream position took precedence in representing the lake type. For instance, a multi-outlet lake might have a headwater outlet and a flow-through outlet at the same time, which are flowing in two different directions. Here, the flow-through outlet has at least one upstream lake, and the headwater outlet has no upstream lake. This mainly occurs due to the limitation of the D8 flow-direction, which cannot depict flow-bifurcation. In this case, the lake was treated to be on the drainage pathway and thus assigned to be a flow-through lake.

We also provided useful proximity information that measures how global lakes are close to one another in terms of drainage distance, i.e., the geodesic distance along our delineated inter-lake reaches (Sect. 2.4). These attributes include "D_distance", defining the drainage distance from each lake outlet to the next downstream outlet; "D_dist_sink", measuring the cumulative drainage distance from the lake outlet to the inland sink or the ocean; and the maximum, minimum, and mean drainage distances to the next upstream outlets ("U_dist_max", "U_dist_min", "U_dist_avg").

## 2.3 Unit catchments

The unit catchment features depict the spatial boundaries of unit catchments (see definition in Sect. 2.2) associated with each of the lake outlets. Note that unit catchments are not based on a one lake to one catchment relationship but are based on a one outlet to one catchment relationship. This was to ensure that the catchment boundaries were given at a spatial detail consistent with the lake outlets or drainage destinations. This way, the count of unit catchments is equal to the count of lake outlets, and bifurcation or multifurcation lakes have multiple unit catchments.

Attributes of each unit catchment include the ID of the associated outlet ("Outlet_id"); the ID of the associated lake ("Hylak_id"); the IDs of the next downstream outlet and lake, respectively ("D_out_id" and "D_hylak_id"); and the area of this unit catchment ("Cat_area"). We also reported the drainage type of each unit catchment in relation to the inter-lake reach network ("Cat_type"), which is identical to the drainage type of the associated outlet (Outlet_type). The pie chart in Fig. 5 summarizes the composition of global lake unit catchment areas in terms of their drainage types,
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Global map of lake outlets and drainage types in TopoCat. Outlet types are based on drainage positions in the inter-lake reach network. The schematic diagram shows the outlet type definitions. The pie chart shows the frequencies of different types of outlets.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Global map of lake unit catchments in TopoCat. The catchment boundaries define the drainage areas between cascading lake outlets (i.e., one catchment for each lake outlet), and when there is no lake further upstream, the catchment defines the drainage area from the headwater to the lake outlet. The pie chart shows the composition of global lake catchment area by outlet drainage type.</figcaption></figure>

where flow-through catchments cover the largest proportion (~79%) in area, followed by headwater catchments and terminal catchments. It is important to note that a terminal unit catchment is not always equivalent to an endorheic basin. The former refers to the unit catchment of a terminal lake, which is the drainage area between the sink and the next up-

stream lake(s), whereas the latter refers to the entire land-locked drainage basin from the headwater to the sink (i.e., the cases where "Basin_type" = "endorheic" in "Lake-network basins" features).

It is worth noting that although the unit catchments were offered at the outlet level, when needed, users can easily dis-
solve (i.e., aggregate) them based on their lake IDs ("Hy-
lak_ID") to form the local catchment boundaries for each
lake (i.e., the drainage area from the upstream lake(s), if
any, to this lake). In total, the delineated unit catchments
in TopoCat cover about 77.5 × 10⁶ km², which is about
26.4 times that of the total lake area and about 57% of the
Earth's land mass excluding Antarctica (Fig. 5).

## 2.4 Inter-lake reaches

The inter-lake reach features contain a set of detailed reach
networks that tie any lake in the world to its upstream and
downstream lakes (if any) through the outlets. Together with
the topology attributes (Table 1), each inter-lake reach net-
work depicts how the lakes in the headwater drain to other
lakes downstream in a cascade pattern, which eventually
reaches a drainage terminus, i.e., an inland sink or the ocean.
Because these networks are determined by lake distribution,
their drainage density varies with lake density across the
continental surface. Such inter-lake networks can be differ-
ent from other river datasets such as MERIT Hydro-Vector,
which were not tailored for lake-determined drainage den-
sity. More specifically, each reach in an inter-lake network
starts from the outlet of one lake and ends with the outlet
of the next downstream lake, but along the way, the reach
is often segmented further by the confluences with tributary
reaches from other lakes. If a lake is not coastal or endorheic
and does not have a downstream lake, we kept the outflow-
ing reach from this lake all the way to the ocean or inland
sink. In other words, a reach network always stretches from
the most upstream lakes to the drainage terminus even when
there is no terminal lake. About 3 million inter-lake reaches
were generated among the 1.46 million outlets (Fig. 6). The
total length of these reaches is about 10 × 10⁶ km.

Main attributes of each reach include the unique reach ID ("Reach_id"); the IDs of the connected upstream outlet and lake, respectively ("U_out_id" and "U_hylak_id"); the count of all upstream lakes from the headwater ("U_lak_ntot"); the ID of the next downstream reach ("D_reach_id"); the IDs of the connected downstream outlet and lake, respectively ("D_out_id" and "D_hylak_id"); the geodesic distance and slope of the reach ("Rch_length" and "Rchslope"); and the Strahler stream order (Strahler, 1957) of the reach based on the inter-lake network ("Rch_order"). Figure 6 shows the global distribution of our delineated inter-lake reach networks, where different colors illustrate the total abundance (count) of lakes upstream of each reach (as in the attribute "U_lak_ntot"). For example, the abundance of upstream lakes per reach in some of the most lake-rich regions, such as the Canadian Shield and Scandinavia, can be several orders of magnitude higher than the rest of the world.

As previously noted, the inter-lake reach networks aim
to define the potential drainage paths connecting global
lakes and reservoirs. These reaches are not always peren-
nial or intermittent rivers but can be highly temporally dy-

namic, non-perennial, and non-channelized. Similar to the
lake type assignment, we here give a rough inference of
the type of each reach based on its relationship with the
river channels in MERIT Hydro-Vector. We reported the
proportion of each inter-lake reach ("Rchint_mhv") over-
lapped by river channels in MERIT Hydro-Vector. The over-
lapping proportion indicates "probably river" if an inter-
lake reach is completely contained by the river network in
MERIT Hydro-Vector (i.e., "Rchint_mhv" = 100%), "par-
tially river" if the inter-lake reach partially overlaps MERIT
Hydro-Vector (i.e., "Rchint_mhv" < 100% and > 0%), and
"other drainage path" if the inter-lake reach is beyond the
spatial extent of the river network in MERIT Hydro-Vector
(i.e., "Rchint_mhv" = 0%).

## 2.5 Lake-network basins

For user convenience, we provided another feature file:
"Lake-network basins". Each lake-network basin, in princi-
ple, defines the entire drainage area of an inter-lake reach net-
work, which stretches from the headwater to an inland sink
or the ocean. Nevertheless, our inter-lake reaches were con-
structed based on lake outlets (refer back to Sect. 2.4). If a
multifurcation lake infringes more than one reach network
that drains to different termini, these networks were consid-
ered related to each other through this shared lake, and their
drainage areas were merged into one lake-network basin.

This logic resulted in 47 340 lake-network basins in the
world as shown in Fig. 7. Among them, endorheic basins (as
indicated by the "Basin_type" attribute) account for 5.1 %
by count and 18 % by area of all lake-network basins. These
endorheic basins cover ~ 15.4 % of global surface excluding
Antarctica. This proportion is smaller than the conventional
size of endorheic basins, which is about a fifth of the global
landmass excluding Antarctica (Wang et al., 2018; Wang,
2020). This is because a lake-network basin delineated here
must contain at least one lake, and some of the conventional
endorheic basins, such as with extreme aridity in the Sahara,
the Arabian Desert, and the Gobi Desert, are arheic and have
no lake presence (according to HydroLAKES), which led
to their exclusion from our delineated lake-network basins.
But, overall, the spatial pattern of the lake-network basins is
consistent with the river basins, signifying again that lakes
are an integrated component of the global river system. Im-
portantly, these lake-network basins define how the millions
of global lakes are partitioned by their drainage dependence,
which can potentially help users gauge or optimize parallel
processing (e.g., for global lake modeling and coupled lake-
river modeling). For this reason, we also assigned the ID of
the affiliated lake-network basin ("Basin_id") to each of the
other TopoCat features.
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Global map of inter-lake reaches in TopoCat. Reach colors illustrate the accumulative lake abundance (count) upstream from each inter-lake reach. Land background is displayed in black to enhance the visual transition from lower to higher upstream lake abundance (akin to increasing stream order). The histogram on the lower-left corner shows the frequency of upstream lake abundance.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Global map of lake-network basins in TopoCat. A lake-network basin is defined as the entire drainage area from the headwater to the ocean or the inland sink, which contains a hydrologically independent inter-lake network. Gray indicates basins having no lake presence. The pie chart illustrates the composition of global lake basin areas by their drainage type (endorheic and exorheic).</figcaption></figure>
# 3 Methods

## 3.1 Input data sources

The presented initial version of Lake-TopoCat was developed based on HydroLAKES v1.0 (Messager et al., 2016) primarily using the 3 arcsec MERIT Hydro (v1.0.1) hydrography dataset (Yamazaki et al., 2019). Other auxiliary data sources include the annual permanent and total lake areas from GLCP (Meyer et al., 2020) to calculate the average permanent water proportion in each lake and the river flowlines in MERIT Hydro-Vector (Lin et al., 2021) to assist in lake and reach type configurations. More details of each of the input datasets are given below.

### 3.1.1 Lake mask and water extent

Constructing global lake topology requires two essential types of data: a global lake inventory with detailed lake boundaries and a DEM (or hydrography) dataset depicting accurate drainage directions. For the former, we selected the widely used HydroLAKES, which contains about 1.4 million lake polygons larger than 10 ha (0.1 km²) compiled from various digital and remote sensing products (Messager et al., 2016). HydroLAKES is, to our knowledge, one of the most comprehensive global-scale lake and reservoir inventories publicly available. Associated with each lake boundary, HydroLAKES also provides the outlet point and a number of attributes such as lake mean depth, lake volume, and catchment area. The lake outlet points and catchment areas in HydroLAKES were derived from the 15 arcsec (~500 m at the Equator) HydroSHEDS database (Lehner et al., 2008) without consideration of bifurcation, and they were later used for comparison with the outlets and lake catchment areas in TopoCat (see Sect. 4.2 for details).

In addition to the lake mask, we used the annual water area statistics reported in the GLCP dataset (Meyer et al., 2020) to calculate the average percentage of the total lake area being permanently inundated for each HydroLAKES polygon. GLCP harmonized HydroLAKES, HydroBASINS, and several climate and socioeconomic datasets to provide a number of attributes such as annual average precipitation, temperature, and water extent for each water body in HydroLAKES. For example, Meyer et al. (2020) obtained the annual water extents of each lake from 1995 to 2015 from the Joint Research Centre (JRC) Global Surface Water (GSW) dataset (Pekel et al., 2016). They also reported the annual permanent and seasonal water areas per lake with the assistance of the JRC-GSW dataset. However, GLCP was unable to distinguish between ice and snow due to the limitation of the water identification algorithm of Pekel et al. (2016). Therefore, the permanently and seasonally frozen lake areas were assigned as “no water” and “entirely seasonal water” areas, respectively.

### 3.1.2 Hydrography

We considered two global hydrography datasets for building lake drainage topology: HydroSHEDS (Lehner et al., 2008) and MERIT Hydro (Yamazaki et al., 2019). The HydroSHEDS hydrography database is available at three different resolutions. Limited by the spatial coverage of SRTM DEM (Farr et al., 2007), the highest-resolution HydroSHEDS v1.1 at 3 arcsec (~90 m at the Equator) is missing the landmass north to 60° N, where 50.5% of the global lakes (by count) are distributed. Although the 15 arcsec version covers the pan-Arctic region above 60° N, the accuracy of the hydrography is significantly inferior due to a poorer quality of the underlying HYDRO1k DEM (EROS Center, 2018). On the other hand, MERIT Hydro v1.0.1 was developed consistently based on the 3 arcsec, high-accuracy MERIT DEM (Yamazaki et al., 2017) with a complete coverage of the global landmass (excluding Antarctica). Compared with SRTM DEM, MERIT DEM eliminated major error components including absolute bias, stripe noise, speckle noise, and tree canopy height, which often distort surface drainage directions and fragment drainage paths. The production of MERIT Hydro also employed a “stream-burning” technique on MERIT DEM. This technique applied multi-source waterbody datasets to derive flow directions that were more realistically aligned with the observed drainage networks. Small dummy depressions were then removed by an inland basin detection algorithm to further smooth the drainage continuity. For these reasons, we selected MERIT Hydro to configure lake outlets, catchments, and interconnecting drainage paths. As one of the supplementary layers, MERIT Hydro also produced a hydrologically adjusted DEM, where elevations were ensured to comply with the condition that the downstream elevation should not be higher than the upstream elevation. This hydrologically adjusted DEM was used to calculate the drainage gradient (slope) of each inter-lake reach.

Furthermore, we used the vector-based global river network dataset, MERIT Hydro-Vector (Lin et al., 2021), to assign the type of each lake based on its drainage position in the network (attribute “Laktyp_mhv”). In this dataset, the spatial variability of drainage density was considered when the networks were delineated from the original MERIT Hydro. A machine learning technique was employed to optimize the drainage density per watershed (with a flexible channelization threshold down to 1 km²) in relationship with various hydroclimate factors. About 58 million river reaches were delineated, with a total length of ~75 × 10⁶ km. It is important to note that the point of MERIT Hydro-Vector was not to simply vectorize MERIT Hydro to the finest spatial detail (as something done by Hydrography90m; Amatulli et al., 2022) but to simulate spatially variable river networks so that the vectorized channels resemble the actual densities of perennial and/or intermittent rivers in the world. Although the temporal dynamics of channel heads due to climate variability
were not accounted for, and the resolvable headwater sizes
are constrained by the spatial resolution of MERIT DEM,
we considered the flowlines in MERIT Hydro-Vector as one
of the best available global representations of perennial and
intermittent river networks.

## 3.2 Dataset development and algorithm

The TopoCat constructing algorithm is illustrated by the
schematic diagram in Fig. 8. In brief, the algorithm started
with an iterative process that used MERIT Hydro's flow ac-
cumulation layer to identify all possible outlets for each Hy-
droLAKES lake (i.e., pour points to all possible outbound
directions from this lake) and each of their associated unit
catchments. The connecting reaches among the lake outlets
were then delineated using the flow direction layer. The ma-
jority of attributes for lake outlets, unit catchments, and inter-
lake reaches (Table 1) were then generated using the hydrog-
raphy, hydrologically adjusted DEM, and the spatial rela-
tionships among the delineated features. This resulted in the
"preliminary TopoCat" database, which consists of the fea-
ture geometries and topological attributes in the finest spatial
detail possible, based on all identified lake outlets. For user
convenience, we also release the preliminary TopoCat (with
the affix "_prelim" in the file names) together with the final
cleaned-up version (as explained in Sect. 2). The definitions
of attributes are consistent in both versions (Table 1).

Next, the preliminary TopoCat was simplified based on the
principle of “one outlet per drainage destination (i.e., down-
stream lake) per lake”, which formed the final cleaned-up
version of TopoCat. Conceptually, if a lake has multiple pre-
liminary outlets draining to the same lake, a representative
outlet was selected to be the one with the maximum flow
accumulation. Accordingly, all preliminary catchments of a
lake that drain to the same destination were dissolved to
one final unit catchment, and the subset of the preliminary
reaches connecting the representative lake outlets were ex-
tracted as the final reaches. The existing attributes were then
recalculated according to the updated network. The lake-
network basins were delineated using the terminal points of
all the available inter-lake reach networks. In the last step, a
few additional lake outlet attributes were populated and as-
signed to the associated lake polygons, and MERIT Hydro-
Vector was used to determine lake types (“Laktyp_mhv”) ac-
cording to their hydrological position in the river networks.

To enable parallel processing, we partitioned the MERIT
Hydro dataset into multiple Pfafstetter level-2 basins. A total
number of 61 Pfafstetter level-2 basins were obtained from
HydroBASINS (Lehner and Grill, 2013), which were orig-
inally generated from the HydroSHEDS hydrography data
(Lehner et al., 2008). Lin et al. (2021) redefined these level-
2 basin boundaries to make them compatible with MERIT
Hydro. This preprocessing, however, only considered water-
sheds larger than 25 km², so we further updated the 61 level-2
basins by annexing any smaller watersheds to the nearest ma-

jor watershed. In addition, a number of islands and archipel-
gos scattered across the Pacific, Atlantic, and Indian oceans
were covered by MERIT Hydro but excluded from the level-
2 basins. We carefully grouped these islands in such a way
that the area of the rectangular box that covers each group
is not larger than the largest extent of the original 61 level-2
basins. This approach creates another seven clusters, form-
ing a total of 68 Pfafstetter level-2 basins or subdivisions
in the world (Fig. 9). These 68 regions cover all pixels in
MERIT Hydro and were used to define the subset domains
of our TopoCat product. For user convenience, we provided
the polygon boundaries of these 68 Pfafstetter regions in the
file named "Pfaf2_regions".

The abovementioned algorithm was implemented in six
specific steps. All steps were automated by Python scripts
with the ArcPy geoprocessing package, except step 4 which
required manual quality control for the identified terminal
outlets (see Sect. 3.2.4). This manual quality control took
approximately 25–30 h by an expert, and much of it was a
one-time process that will require no replication even if the
input lake mask is updated. We describe each of the steps in
the sections below.

### 3.2.1 Delineating preliminary lake outlets and unit catchments

Although HydroLAKES provides lake outlet points derived
from the 15 arcsec HydroSHEDS data, we regenerated the
lake outlets using the higher-resolution (3 arcsec) MERIT
Hydro dataset. Different from the outlets available in Hy-
droLAKES, we were not constrained by the “one outlet for
one lake” assumption in order to capture possible lake mul-
tifurcation. However, we had to assume that multifurcation
can only occur at the lake level not the outlet level, be-
cause pixel-based bifurcation was not assumed in MERIT
Hydro. Typically, the outlet location of a lake can be identi-
fied by intersecting the lake shoreline with the flow direction
layer and then locating the intersecting cell(s) flowing out-
ward from the lake (Tokuda et al., 2021). However, there is
no guarantee that the lake outlets always fall exactly at the
shoreline cells. An example is an endorheic lake, where the
drainage outlet (the sink in this case) is located inside the
lake boundary. Other examples are some of the lagoon poly-
gons in HydroLAKES, which can overshoot the terrestrial
boundary of MERIT Hydro due to resolution or precision in-
consistency between the lake and hydrography datasets. To
address these issues, we adopted a flow-accumulation-based
approach, which identified the lake outlet as the location of
maximal flow accumulation within the lake boundary.

The flow-accumulation-based approach was employed it-
eratively until all possible outlets for each lake were iden-
tified. As illustrated in Fig. 10a, the cell where flow accu-
mulation is the maximum within the lake boundary was first
identified as the primary outlet, and the unit catchment as-
sociated with the primary outlet (hereafter "primary catch-
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Schematic diagram of the algorithm used to construct TopoCat. Rounded corner rectangles, rectangles, and oval shapes represent input data, intermediate products, and final TopoCat products, respectively.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** A total of 68 partitioning regions (Pfafstetter level-2 basins) in the world for TopoCat production and subset organization.</figcaption></figure>

ment”) was delineated using flow direction. Logically speaking, the drainage basin of a lake should completely contain the lake boundary. If the primary catchment covers only part of the lake polygon, we assumed that possible multifurcation may occur for this lake. In this case, the second iteration was triggered to identify the next maximum flow outlet (i.e., the secondary outlet) from the lake cells outside of the primary catchment, and the unit catchment associated with the secondary outlet was delineated (Fig. 10b). The iteration was repeated until the entire lake area was covered by its unit catchment(s). The iteration sequence, which corresponds to a descending order of the outlet flow accumulation, was as-

signed to each outlet to keep track of the multiple outlets of the same lake. Locations of all possible lake outlets and their associated unit catchments were delineated through this iterative approach. Each outlet and the associated unit catchment were labeled with the same unique outlet ID and associated lake ID. We termed them as “preliminary outlets and unit catchments” to make a distinction from the representative outlets and their associated unit catchments in our final lake topology (see Sect. 3.2.3).

Despite topological redundancy, there are multiple benefits to retrieving all preliminary outlets using the above-described iterative process. First, this process was a neces-
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Schematics of the iterative process for outlet and catchment delineations (a–c) and the improvement in lake connectivity and catchment completeness due to this iterative process (d, e). Lakes are shown as striped polygons, lake unit catchments as color-filled polygons, lake outlets as red dots, and inter-lake reaches as blue lines. (a) Delineated primary lake outlets (i.e., no iteration) with their associated unit catchments and connecting reaches are shown. Note that lake 1 and lake 2 are not contained by their unit catchments. (b) Delineated multiple preliminary outlets and their associated unit catchments through the iterative process are shown. (c) The final selected outlets with their dissolved unit catchments and connecting reaches are shown (described in Sect. 3.2.3). Lake 1 and lake 2 are now contained by their own unit catchments. (d) Lake topology and unit catchments generated by the non-iterative process are shown. Lake 3 drains to lake 6 without intersecting the outlet of lake 5, and the unit catchment of lake 5 appears incomplete. (e) Lake topology and unit catchments delineated by the iterative process, with improved inter-lake drainage connectivity and unit catchment completeness. Here, both outlets of lake 5 have been maintained in the final TopoCat dataset to conserve the topological connectivity between lakes 3, 5, and 6. Here, lake boundaries are the same as HydroLAKES (Messager et al., 2016).</figcaption></figure>

sary preparation for capturing lake multifurcation (see the example in Fig. 3). Second, identifying all lake outlets ensured a proper portrayal of the complete lake catchment (Fig. 10b and c). Third, it also ensured that the lake drainage topology was consistently configured. In the example of Fig. 10e, both lake 3 and lake 4 first drain to lake 5 which then drains to lake 6. Without the iteration, the secondary outlet of lake 5 connecting lake 3 and lake 6 would have been missing, meaning that in the configured topology the directly connected downstream lake from lake 3 would have been misidentified as lake 6.

### 3.2.2 Building preliminary lake topology and inter-lake reaches

Lake topology defines how different lakes are hydrologically connected to each other. In this step, such topological infor-

mation, together with the inter-lake reach networks, was established using the preliminary lake outlets and MERIT Hydro. Given each preliminary outlet, the next downstream outlet was traced by routing from this outlet pixel based on the MERIT Hydro flow directions. The identified downstream outlet ID and its corresponding lake ID ("D_out_id" and "D_hylak_id") were saved in the attributes of both preliminary outlet and catchment features.

During this routing process, the drainage path connecting the preliminary outlets was delineated pixel by pixel until the path reached an inland sink or the ocean. The traced paths were converted into vector format to create an inter-lake reach network, where individual reaches were segmented by lake outlets, drainage termini (sinks or the ocean), and the confluences of the reaches (Fig. 10b). For each reach, the next downstream reach ID ("D_reach_id") was also configured to depict the drainage topology of this inter-lake
<figure><img src="image_13.png" /><figcaption>**Figure 11.** An example of QA/QC for endorheic lakes. (a) Derived TopoCat features surrounding an endorheic lake with partial water extent in HydroLAKES. The incomplete polygon for this endorheic lake does not contain its drainage sink (black dot) in MERIT Hydro. As a result, the derived lake outlet before QA/QC (red dot labeled by “endorheic check”) is on the incomplete lake boundary, and the associated terminal catchment was delineated incompletely. (b) Improvement in the terminal unit catchment and drainage topology after QA/QC, where the “endorheic check” outlet was automatically snapped to the sink, and the updated terminal catchment became complete. Here, lake boundaries are the same as HydroLAKES (Messager et al., 2016), and the map background is from Esri, Maxar, Earthstar Geographics, and the GIS User Community in both figures.</figcaption></figure>

boundary, leading to a misassignment of the outlet type as “headwater” or “flow-through” (rather than “endorheic”). Moreover, the catchment delineated based on the misplaced outlet did not enclose the endorheic basin. To tackle this issue, an intermediate outlet type “endorheic-check” was assigned to any outlet that was directly connected to a MERIT sink for visual quality assurance and quality control (QA/QC) in the next step.

### 3.2.4 Quality control for endorheic lakes

This was the only step in the TopoCat construction that involved manual intervention. In the QA/QC process, each outlet temporarily labeled as “endorheic-check” in the drainage type attribute was visually inspected against high-resolution Google Earth and Esri satellite images. With expert knowledge, this visual inspection aimed at judging whether this labeled outlet is a true outlet of an open lake or an “artifact” pour point of an incomplete endorheic lake (as in the example of Fig. 11a). The outlets verified as the latter scenario were flagged for automated correction. To minimize manual

workload and avoid further complications, we performed this QA/QC on the final (rather than preliminary) outlets.

In the following automated correction, each flagged outlet was snapped to its connected sink (the snapped outlet in Fig. 11b). If there are multiple outlets from the same incomplete endorheic lake polygon (which all drain to a single sink), the outlet associated with the largest total catchment area (“Cat_a_tot”) was selected for this snapping process. If there are fragmented HydroLAKES polygons representing the same endorheic lake (i.e., with different HydroLAKES IDs), the outlet associated with the largest polygon was selected for snapping. This automated snapping was performed on the preliminary outlets, and the unit catchments were re-delineated based on the updated outlets (drainage sinks in this case). Steps 2 and 3 (Sect. 3.2.2 and 3.2.3) were then repeated to accommodate these corrections in all preliminary and final outlet, catchment, and reach attributes. It is worth noting that this QA/QC did not aim at dissolving the fragmented endorheic lake polygons in HydroLAKES. This was to maintain the consistency that every HydroLAKES polygon has at least one outlet and one unit catchment. However, since there was only one outlet snapped to the sink per en-
endorheic lake, only the snapped outlet and its associated unit catchment were corrected to an endorheic drainage type. In other words, for an endorheic lake that has multiple Hydro-LAKES polygons, only the polygon with a snapped outlet was labeled "endorheic".

A total of 1261 outlets were identified as "endorheic-check". Among them, 291 were confirmed to be pour points of incomplete endorheic lake polygons, and 14 of these 291 outlets shared drainage sinks with other outlets. The locations of the other 277 outlets were snapped to their connected drainage sinks, and the drainage catchments were re-delineated. As exemplified in Fig. 11, the result of this QA/QC was a complete endorheic catchment and corrected topological attributes despite an incomplete lake extent in the original HydroLAKES dataset.

### 3.2.5 Configuring reach types, stream orders, outlet orders, and lake-network basins

The inter-lake reaches in TopoCat depict feasible drainage connections between lakes, whereas the river networks in MERIT Hydro-Vector (Lin et al., 2021) represent, in principle, the actual perennial or intermittent channels (refer back to Sect. 3.1.2). Given this logic, we applied MERIT Hydro-Vector to obtain a priori information regarding whether an inter-lake reach in TopoCat is likely a river channel or a non-channelized drainage connection between the lakes. The proportion of each inter-lake reach that intersects MERIT Hydro-Vector was calculated in percentage of length and saved in the attribute "Rchint_mhv" of the inter-lake reaches.

The Strahler stream order (attribute "Rch_order") was also calculated for each inter-lake reach. The order of each lake outlet ("Out_order") was then assigned as the stream order of the connected upstream reach (whose end node coincurs with this outlet). Therefore, the headwater outlets always have a zero outlet order. If multiple reaches converge to the same outlet, the maximum stream order among the confluence reaches was assigned as the outlet order. Finally, the lake-network basins were delineated using the terminal points (inland sinks or coastal outlets) of all inter-lake reach networks. On rare occasions, a multifurcation lake may infringe more than one reach network; when this occurred, the multiple basins encompassing these reach networks were dissolved into one basin. The "Basin_id" was then assigned to all associated lake outlets, unit catchments, and inter-lake reaches in both the preliminary and final TopoCat.

### 3.2.6 Aggregating attributes at the lake level

The TopoCat attributes for each lake were configured by aggregating the values of the associated outlet attributes. They include the counts of directly connected upstream and downstream lakes ("U_lake_n" and "D_lake_n"), the counts of all upstream lakes from the headwater and all downstream lakes to the terminus ("U_lak_ntot" and "D_lak_ntot"), the

count of outlets for each lake ("Outlet_n", > 1 indicating multifurcation), and the area of the entire lake catchment ("Cat_a_lake", a dissolve of all unit catchments of and upstream from this lake). Lake orders ("Lake_order") were directly obtained from the outlet orders, and in the case of a multifurcation lake, its lake order was based on the maximum of its outlet orders. The ID of the lake-network basin that each lake belongs to was also assigned according to the basin ID of the associated lake outlet(s).

The drainage type of each lake was configured in this step. As explained in Sect. 2.1, we provided two lake drainage types depending on the applied drainage networks: "Lake_type" based on our own inter-lake reaches and "Laktyp_mhv" based on MERIT Hydro-Vector. For most lakes with a single outlet, the value of Lake_type was equivalent to the outlet type, but for each multifurcation lake, the value was sorted out by considering the drainage hierarchy among its outlets. As described in Sect. 2.2, precedence was given to the outlet type corresponding to the most downstream drainage position. In such cases, "endorheic" and "coastal" outlets have the most weight in determining their lake types, followed by "flow-through" and "headwater" outlets.

The values of "Laktyp_mhv" were modified from "Lake_type" values with reference to MERIT Hydro-Vector (Fig. 2b). MERIT Hydro-Vector contains river networks with variable drainage densities (hereafter "VDD networks"), as well as river networks with a 25 km² constant drainage density ("CDD networks"). While the VDD networks are the state-of-the-art global river networks available today, the river topology in VDD networks is occasionally disconnected, and these networks do not cover Greenland or parts of other circum-Arctic regions. To address the problems, we independently used both VDD and CDD networks to identify lake drainage positions and then merged the results (i.e., VDD and CDD lake types) to jointly determine Laktyp_mhv.

Specifically, for any lake where the "Lake_type" value was already the most downstream possible, i.e., terminal or coastal, the same type was assigned to its VDD lake type as well. Lakes that are geometrically off the VDD networks were updated to be isolated. For the rest of the lakes, the counts of inlets and outlets were obtained through intersecting the lakes with the VDD networks. Note that the count of directly connected upstream lakes in the VDD networks could be different from that based on our inter-lake reach networks. If a lake has at least one directly connected upstream lake, at least one inlet, and also at least one outlet, the lake was identified as a "flow-through" type. The lake was identified as "inflow-headwater" if the lake has at least one inlet and at least one outlet but no directly connected upstream lake. If the lake has only outlets but no directly connected upstream lake or inlet, the lake was identified as "headwater" type. The same process was repeated for the CDD networks to identify the CDD lake types. The final value of "Laktyp_mhv" was then assigned by comparing the VDD and CDD lake types. Similar to the Lake_type attribute.
tribute, we maintained the preference hierarchy towards more
downstream types, where “endorheic” and “coastal” have
the most preference, followed by “flow-through”, “inflow-
headwater”, “headwater”, and “isolated” types.

Finally, the average proportion of permanent water area ("Lperm_glcp") for each lake was calculated using the annual permanent water area and the annual total lake water area (integrating both permanent and seasonal water) obtained from GLCP. GLCP provides the seasonal, permanent, and total water areas for each HydroLAKES water body per year from 1995 to 2015. Since the purpose of "Lperm_glcp" is to represent expected seasonal stability of each lake, we balanced two factors in selecting the optimal time window. On the one hand, some lakes exhibited significant expanding or shrinking trends. If the time window is too long, the proportion of permanent water area for these lakes would be dominated by the signal of long-term trends rather than seasonality. On the other hand, if the time window is too short, the proportion of permanent water might be biased by short-term climate variability such as the El Niño–Southern Oscillation (ENSO). To avoid these issues, we selected a 6-year time window from 2010 to 2015, which covers at least one major La Niña event (2010–2011) and one major El Niño event (2015–2016) (Wang et al., 2018), thus balancing out both wet and dry anomalies to represent a more objective inundation pattern. We then calculated "Lperm_glcp" (in %) for each lake by averaging the ratios between annual permanent water area and annual total water area during the selected 6-year period.

4 Quality assessment and validation

The quality of TopoCat was first assessed by comparing its
lake catchment areas and lake outlet locations with those
in HydroLAKES or LakeATLAS with reference to high-
accuracy regional hydrography datasets: LakeCat (Hill et al.,
2018) and NHDPlusV2 (McKay et al., 2012) for the conti-
nental USA (CONUS). Since the lake polygons, lake out-
lets, and catchment or watershed area in LakeATLAS are
the same as in HydroLAKES, we only refer to LakeATLAS
hereafter, but results are also true for HydroLAKES. This
comparison is followed by a more detailed validation against
the lake topology acquired from NHDPlusV2 and LakeCat
in a lake-dense area of Minnesota. Finally, the performance
of our delineated endorheic basins in TopoCat was assessed
against endorheic basins delineated by Liu et al. (2021) in
the Tibetan Plateau for the years 2000 and 2018. Details of
the reference datasets, quality assessment, and validation are
given in the following sections.

4.1 Validation data

Detailed and accurate information on lake topology and
basins is lacking on a global scale, so we benchmarked
TopoCat against two high-quality regional hydrographical

networks. NHDPlusV2.1 (McKay et al., 2012), a value-
added version of the 1:100 000 NHD dataset (USGS, 2001),
was used to validate our derived lake catchment areas in the
continental USA and lake topology in part of the Prairie
Pothole region in the USA, one of the world's most lake-
dense areas. NHDPlusV2.1 contains catchment boundaries
for NHD on-network waterbodies and their drainage topol-
ogy based on flow directions. Here, on-network waterbodies
are those which intersect the NHD flowline network. For the
waterbodies off the NHD network, we obtained catchments
and drainage topology from the LakeCat dataset (Hill et al.,
2018). Together, these datasets represent the most up-to-date
and comprehensive lake catchment and drainage topology in-
formation for the USA.

In addition, we used the watershed boundaries on the in-
ner Tibetan Plateau produced by Liu et al. (2021) to vali-
date the endorheic lake catchments for the same region in
TopoCat. The inner Tibetan Plateau hosts the world's largest
cluster of terminal lakes. Delineating lake catchments in arid
and endorheic basins can be challenging for at least two rea-
sons. First, small pseudo depressions due to errors or noise
in the elevation data often coexist with actual inland basins,
and eliminating them from the hydrography data is difficult
(Yamazaki et al., 2019). Second, many terminal lakes on the
Tibetan Plateau are rapidly expanding due to a wetting and
warming climate (W. Chen et al., 2022; Wang et al., 2018;
Yao et al., 2018). Such lake dynamics have led to ongoing
drainage reorganizations such as basin merging and annexa-
tion, meaning that the catchments derived from a static lake
mask (such as HydroLAKES) are prone to boundary migra-
tion through time. Despite the availability of other inland lake
catchment data for this region (e.g., Yan et al., 2019), the
dataset from Liu et al. (2021) is, to our knowledge, the only
one that considered the ongoing drainage changes on the Ti-
betan Plateau due to lake inundation dynamics. They used a
lake-oriented approach (Liu et al., 2020), together with the
calibration of multi-temporal lake mappings, to refine the
basin boundaries initially derived from MERIT DEM. The
dataset offers the endorheic lake catchments for the years
2000 and 2018 (hereafter "Liu21 endorheic basins 2000 and
2018"), with total catchment numbers of 434 and 421, re-
spectively.

4.2 Comparison with LakeATLAS

We hypothesize that the overall accuracy of lake outlets and
catchment areas in TopoCat, which were derived from the
3 arcsec MERIT Hydro, outperforms those in LakeATLAS
derived from the 15 arcsec HydroSHEDS. In addition, the
lake catchments in TopoCat were delineated using an iter-
ative process to consider possible multifurcation, so our de-
rived lake catchment boundaries are likely more accurate and
complete. To test this hypothesis, we benchmarked the to-
tal lake catchment areas (i.e., the entire drainage areas from
the headwater to each lake outlet) derived from LakeATLAS
<figure><img src="image_14.png" /></figure>

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

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

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

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

**Figure 12.** Improvements in TopoCat catchment areas, outlet locations, and catchment-to-lake area ratio in comparison with the information in LakeATLAS. (a) Comparison of total lake catchment areas in between TopoCat and LakeCat (benchmark). Only lake polygons from both datasets overlapping at least 90% of the area were considered for this comparison (n = 3991). (b) Same as Fig. 12a except the comparison was made between LakeATLAS and LakeCat. (c) Locations of lake outlets in TopoCat and LakeATLAS, with respect to the outlets derived from NHD (benchmark). Here, NHD outlets were located using the intersecting points between outbound NHD flowlines and the NHD waterbodies showing at least 90% area overlap with lake polygons in LakeATLAS. Outlets of 2596 NHD waterbodies were identified from 3991 overlapped lakes (lake boundaries are the same as HydroLAKES; Messager et al., 2016; and the map background is from Esri, Maxar, Earthstar Geographics, and the GIS User Community). (d) Offsets (in geodesic distance) of LakeATLAS vs. TopoCat outlets from NHD outlets. Points above the diagonal line indicate that TopoCat outlets are closer to NHD outlets than LakeATLAS outlets are. (e) Comparison between the cumulative distributions of the catchment-area to lake-area ratio calculated from TopoCat and LakeATLAS. About 5% of the lakes in LakeATLAS have larger lake areas than their catchment areas, whereas in TopoCat such unreasonable cases are less than 0.04%.
and TopoCat against LakeCat. Following the method performed by Lehner et al. (2022) for validating LakeATLAS, we identified 3991 pairs of lakes in between TopoCat and LakeCat (lake polygons in the latter acquired from NHD-PlusV2) which have at least a 90% overlapping area per pair. Values of the lake catchment area attribute in LakeCat were then spatially joined to TopoCat. Since TopoCat was developed using the same lake polygons as in LakeATLAS, these lake pairs are exactly the same as those used for LakeATLAS validation (Lehner et al., 2022), and the only difference is the lake catchment areas used for comparison with those in LakeCat. As shown in Fig. 12a, the agreement of lake catchment areas between TopoCat and LakeCat (Fig. 12a; log-log least-square linear regression $R^2 = 0.72$, symmetric mean absolute percentage error (SMAPE) = 36%, and *n* = 3991) is overall higher than that between LakeATLAS and LakeCat (Fig. 12b; log-log least-square linear regression $R^2 = 0.63$, SMAPE = 51%, and *n* = 3991). Figure 12c compares the cumulative distributions (in relative frequency) of the catchment-to-lake area ratio calculated from TopoCat and LakeATLAS. In TopoCat, the catchment-to-lake area ratios for 80% of the global lakes range from 2 to 70. In LakeATLAS, however, about 5% of the lakes have surface areas larger than their own catchments. This logical error was reduced to less than 0.04% in TopoCat, probably due to the use of higher-accuracy hydrography data and a careful consideration of lake multifurcation, leading to more complete catchment boundaries.

We also compared the geodesic offsets of LakeATLAS and TopoCat outlets from the benchmark lake outlets. For simplicity, only the primary outlets in TopoCat were used for this comparison in the case of multifurcation lakes. For benchmarking, we intersected NHD flowlines with the 3991 NHD-PlusV2 lake polygons identified in the previous step. The intersecting point between each lake polygon and the outbound NHD flowline was considered the benchmark outlet (hereafter “the NHD outlet”). To be consistent with the other datasets, the outlet with the largest upstream catchment area was selected in the case of multiple NHD outlets for a single lake. A total of 2569 NHD outlets were generated across CONUS, and the average offset to LakeATLAS and TopoCat outlets is 577 (± 1402) m and 393 (± 1395) m, respectively, indicating an overall better accuracy of lake outlet locations in TopoCat. Figure 12c shows a few examples in northern Wisconsin, where the TopoCat outlets are in much closer proximity to the NHD outlets than the LakeATLAS outlets. In ~76% of the 2569 cases, the offsets between TopoCat and NHD outlets are ~361 (± 735) m smaller than the offsets between LakeATLAS and NHD outlets (Fig. 12d), confirming that the overall quality of the outlet locations in TopoCat is superior to that provided in LakeATLAS.

## 4.3 Validation of lake topology

Validating lake topology is challenging as there are only a few regional hydrography datasets that consider lake components, and additionally, lake topology is not readily available in these datasets, and retrieving it requires a substantial amount of preprocessing. In this study, we benchmarked our TopoCat lake topology against the reference data in a lake-dense area of Minnesota (Fig. 13), located at the intersection of the Canadian Shield (the Laurentian Plateau) and the Prairie Pothole region. We selected this validation site for two reasons. The first reason was the availability of high-quality hydrography data (in this case, NHDPlusV2 and LakeCat). The second reason is that the accuracy of TopoCat is expected to be relatively low in this region given the difficulty of using gridded hydrography to portray the geological and hydrological complexity in such lake-rich areas. Thus, our rationale is that if the performance of TopoCat is acceptable in this lake-rich area, it may work equally well or better for the rest of the world.

The TopoCat and benchmark datasets were preprocessed to enable a comparison at equivalent levels. First, the waterbodies in NHDPlusV2 were intersected with the lake polygons in TopoCat within the validation site to retrieve the common lakes in both datasets. A linking ID between the common lakes was established. In occasional cases where one TopoCat lake covers multiple NHDPlusV2 lakes, the link was established for the most downstream NHDPlusV2 lake among the ones intersecting TopoCat. This way, the topological orders among the linked subset of NHDPlusV2 lakes were still maintained. A total of 427 linked lake pairs were established between NHDPlusV2 and TopoCat, where 112 intersect the NHDPlusV2 flowline network (hereafter referred to as “on-network lakes”). Topological orders of the on-network lakes were obtained from the NHDPlusV2 dataset using the topological orders of their associated catchments. Several on-network NHDPlusV2 lakes intersect no lake in TopoCat. In such cases, we assumed that the lake is absent in TopoCat (i.e., the HydroLAKES polygons), so such lakes were removed from the assessment, and the topological orders of the remaining on-network lakes were updated accordingly.

The topological orders of the off-network lakes in NHD-PlusV2 (315 in total) were obtained from the LakeCat database (i.e., from associated lake catchment orders). Similar to the preprocessing for on-network lakes, any off-network lake intersecting no lake in TopoCat was removed from the assessment, followed by an update of the topological order. The updated off-network lake topology was then appended to the updated on-network ones to form a completely preprocessed validation or benchmark lake topology (hereafter “NHDPlusV2 topology”). It is noteworthy to mention that only on-network waterbodies were burned into the DEM in development of the NHDPlusV2 hydrography dataset (i.e., flow direction, flow accumulation). On
<figure><img src="image_19.png" /><figcaption>**Figure 13.** Skill of TopoCat lake topology in the lake-rich Prairie Pothole region. (a) Comparison with the benchmark lake topology derived from NHDPlusV2 and LakeCat (lake boundaries are the same as HydroLAKES; Messager et al., 2016). Lines illustrate the accuracy of detecting the next downstream lake, whereas lake colors indicate the accuracy of detecting the total count of upstream lakes. (b) The skill of detecting the next downstream lake is shown. The total number of lakes that correctly (and incorrectly) identified the downstream lakes are shown on the bars. (c) Distribution of lake area according to the accuracy of detecting the total upstream lake count.</figcaption></figure>

the other hand, LakeCat directly used the flow direction of NHDPlusV2 to derive the topology of off-network lakes, instead of developing a hydrologically adjusted flow direction through burning off-network lakes in the DEM. Therefore, the reference topology of the on-network lakes is more reliable than the off-network lakes. Akin to the benchmark data processing, TopoCat lakes intersecting no lakes in NHDPlusV2 were excluded from the assessment, and all attributes in TopoCat were updated to accommodate these changes.

The updated TopoCat topology was compared with NHDPlusV2 to assess the skill in identifying the downstream lake, as well as the skill in identifying the upstream lake abundance (i.e., the total count of upstream lakes from the headwater to each lake). The spatial distributions of these skills are shown in Fig. 13a. Here, thick lines represent correct connections made by TopoCat, and thin lines represent incorrect connections. About 80% of the on-network lakes in TopoCat are correctly connected to the next downstream lake, and for off-network lakes the skill is ~55% (Fig. 13b). Different colors of lakes show the skill in capturing the upstream lake abundance (Fig. 13a). We considered TopoCat to be skillful

when the difference in the total upstream lake count between TopoCat and the benchmark data is within ±1. Given this tolerance, about 76% of the lakes in TopoCat capture the total count of upstream lakes, which cover ~88% of the lake surface area within this region (Fig. 13c).

## 4.4 Validation of endorheic basins

Since the Tibetan Plateau is one of the most endorheic-lake-rich regions in the world, we selected this region to validate the endorheic basins derived in TopoCat (as in the lake-network basin layer), using the benchmark data “Liu21 endorheic basins 2000 and 2018” (Liu et al., 2021; refer back to Sect. 4.1). As a preparation of catchment delineation, the raw DEM is usually processed to fill any trivial depressions or artifact pits before generating hydrologically realistic flow direction and flow accumulation. However, during this pit-filling process, true inland depressions may also be removed from the DEM (Grimaldi et al., 2007), which makes it challenging to delineate endorheic catchments properly. Liu et al. (2021) used a lake-oriented approach to delineate endorheic basins on the Tibetan Plateau. In their method, individual
terminal lakes were first mapped from spectral images, and their associated endorheic catchments were then delineated accordingly. A total of 253 endorheic basins were selected from TopoCat on the studied Tibetan Plateau. Unlike Liu et al. (2021), the topography modification in MERIT Hydro was not tailored for each observed endorheic lake. Instead, Yamazaki et al. (2019) assumed that a significant volume change in topography was required to merge a true inland depression to its adjacent basin. After trial-and-error tests, MERIT Hydro adopted a threshold modification volume (2500000 m³) to automatically separate true depressions from the dummy ones. Although this automatic approach has proven to be effective at a global scale, it occasionally failed to detect true depressions on the Tibetan Plateau, ultimately connecting some of the actual endorheic basins to their neighboring catchments. This type of error particularly occurs for small endorheic basins, where an effective pit-fill volume is smaller than the threshold modification volume.

The boundaries of endorheic basins in TopoCat are visually compared with those of Liu21 endorheic basins 2018 (the validation benchmark) in Fig. 14a. The numbers of basins in TopoCat and the benchmark data on the endorheic Tibetan Plateau are 259 and 245, respectively, and in 13 cases, the latter covers multiple basins in TopoCat. In such cases, we merged the TopoCat basins enclosed by the corresponding benchmark basin to make the comparison feasible. The filling colors in Fig. 14a illustrate the critical success index (CSI) value for each TopoCat endorheic basin, which is here calculated as

$$ CSI = \frac{H}{H + M + F}, \qquad (1) $$

where H is the area overlapped by TopoCat and the benchmark basins, M is the benchmark basin area not overlapped by the corresponding TopoCat basin, and F is the TopoCat basin area not overlapped by the corresponding benchmark basin. Although the input hydrography and lake mask data are different in the two datasets, their basin boundaries show a fairly good agreement. The overall CSI value of all the 245 endorheic basins is ~0.85, with the individual basin CSI values ranging from 0.03 to 0.99. The largest difference occurs in the endorheic basin of Pangong Lake, which was misidentified as part of the exorheic Indus Basin in TopoCat. Our visual examination attributed this omission error to the automated process of dummy depression removal in MERIT Hydro, where the adopted threshold modification volume was too large for this basin.

Figure 14b further shows the scatter-plot comparison between the TopoCat and benchmark endorheic basins. Here TopoCat basins were validated against Liu21 endorheic basins for both year 2000 (log-log least-square regression R² = 0.86, SMAPE = 17.5%, n = 246) and year 2018 (log-log least-square regression R² = 0.87, SMAPE = 16%, n = 245). During these two decades, lake expansions on the

inner Tibetan Plateau triggered several cases of lake capture and basin coalition (Liu et al., 2021). These ongoing basin reorganizations explain some of the basin area discrepancies. Despite such uncertainty, our delineated endorheic basin areas in TopoCat agree well with the benchmark basin areas in both 2000 and 2018 except a few outliers such as Pangong Lake. TopoCat shows slightly overestimated endorheic basin areas in a few cases. In these cases, small endorheic basins were merged with the adjacent basins during the dummy depression removal process, thus increasing the sizes of the adjacent endorheic basins.

## 5 Database availability and versioning

The Lake-TopoCat database and its possible future versions are available in both shapefile and geodatabase formats through Zenodo (https://doi.org/10.5281/zenodo.7916729; Sikder et al., 2023) under the Creative Commons Attribution 4.0 International license. While the method to develop TopoCat is generic and adaptive, the version illustrated in this paper was generated for the global lake mask HydroLAKES v1.0 using primarily MERIT Hydro v1.0.1. As our intention is to make TopoCat an input-data-flexible and recursively improved database, future versions may be produced using other global or regional lake masks and hydrography datasets. To avoid version ambiguity, the specific version illustrated in this paper should be referred to as “HydroLAKES-TopoCat v1.1”, where the first part of the name “HydroLAKES” indicates the applied lake mask and where the second part “TopoCat v1.1” indicates the version of the algorithm and other specifics to construct lake topology and catchments. We will follow the same naming convention for future versioning. Similarly, “Lake-TopoCat” should be used to refer to the generic or collective version of our current and future lake topology and catchment data products.

## 6 Applications and limitations

We developed the initial version of Lake-TopoCat (HydroLAKES-TopoCat v1.1), which contains the topological orders, drainage paths, and unit catchments for 1.4 million lakes and reservoirs larger than 10 ha in the world (HydroLAKES v1.0). Among them, the inter-lake reaches are one of the key features. About 3 million inter-lake reaches were delineated, ~80% of which are shorter than 10 km. These reaches stretch ~10 × 10⁶ km, which is about 4.6 times longer than global drainage networks with rivers wider than 30 m (Altenau et al., 2021; Allen and Pavelsky, 2018). The detailed inter-lake drainage networks in TopoCat allow us to decipher lake drainage topology, which may help model how aquatic species possibly migrate upstream or downstream through the connecting reaches as a response to climatic or ecological changes. This lake
<figure><img src="image_20.png" /><figcaption>**Figure 14.** Skill of TopoCat endorheic basin delineation on the Tibetan Plateau. (a) Comparison of Lake-TopoCat endorheic basin boundaries with the benchmark data (Liu21 endorheic basins 2018). Different colors illustrate the critical success index (CSI) values, calculated based on the area of the benchmark basin overlapped by the corresponding TopoCat basin (see Eq. 1). (b) Cross-comparison between TopoCat basin areas versus those in Liu21 endorheic basins (2000 and 2018, respectively).</figcaption></figure>

topology is required to integrate medium- and small-sized lakes and reservoirs, which are largely lacking in the existing routing models, into the global river networks. Such integrated lake–river networks will be essential to route water storage variability in lakes and reservoirs, such as observed by satellite altimeters including the Surface Water and Ocean Topography (SWOT) mission (Biancamaria et al., 2016), into hydrological models and thus benefit the estimates of river discharge, fluvial sediment transport, and carbon flux cycling. For instance, Abril and Borges (2019) studied carbon leaks from flooded land, and they argued that a global topology of fluvial and lacustrine wetlands is necessary to improve the quantification of hydrological carbon fluxes in flooded ecosystems. This argument also applies to lakes and reservoirs, which are important for carbon sequestration.

The spatiotemporal variability of lake–river connections is important to understand surface water dynamics, especially in the Arctic and arid regions, where lake–river connections are intermittent and complex. In Arctic environments, the fluvial and glaciofluvial processes heavily depend on lake–river connectivity. For example, Dolan et al. (2021) found that about 10% of the lakes in Colville Delta in Alaska have variable temporal connectivity, and the lake–river connectivity is highly correlated with river discharge. This connectivity influences ice breakup timing, where ice in highly connected lakes breaks up earlier than in less connected lakes. The inter-lake drainage networks of TopoCat offer all possible hydrologic connections among global lakes, which can serve as a priori information to monitor the intermittency and dynamics in fluvial–lacustrine connectivity. The lake topology is essential to identify upstream lake abundance as well (Gardner et al., 2019), which is useful for understanding the

density of lakes in a region and their accumulative downstream impacts. From TopoCat, we found that ~82% of the global lakes have a drainage proximity to other lakes of less than 10 km, and the percentage increases to ~95% if the proximity increases to 100 km. Lake-TopoCat data can be used to analyze geomorphic scaling principles that have been shown to be related to lake formation processes and even lake water residence time (Seekell et al., 2021).

Our delineated lake catchments reveal that nearly 57% of the Earth's landmass (excluding Antarctica) drains into a lake larger than 10 ha. These unit lake catchments will allow us to understand lake water mass balance and water quality changes by explicitly considering the geophysical and anthropogenic characteristics of the lake drainage basins and thus the contributions of the basins. The amount of crop-land in a lake catchment can affect the amount of nutrient-driven primary production in lentic ecosystems (Balmer and Downing, 2011). Similarly, the amount of wetland coverage in a lake's catchment controls the amount of dissolved organic material transported downstream, cascading to primarily production in lakes, and ultimately carbon dioxide emissions from these waterbodies (Borges et al., 2022; Maberly et al., 2013). The geomorphology (e.g., gradient, altitude) of a lake's catchment may also play a role in controlling lake biogeochemistry and greenhouse gas emissions of lakes (Casas-Ruiz et al., 2020; Lapierre and Giorgio, 2012).

A few recent studies (e.g., GLCP and LakeATLAS) used the lake drainage basins to compute basin-averaged river area, temperature, precipitation, population, and many other hydrological–limnological variables (e.g., Lehner et al., 2022; Meyer et al., 2020). However, lake catchment perimeters are not directly available from these studies and often do
not accurately represent the actual drainage areas of the lakes (Meyer et al., 2020). Together with other thematic datasets, the high-resolution (i.e., 3 arcsec) lake catchment boundaries in TopoCat provide users with more accuracy and flexibility to customize the calculations of basin-level limnological and socioeconomic attributes. Furthermore, the lake drainage topology along with lake basin boundaries will help users better understand lake water quality issues, such as tracking temperature and nutrient flows from the upstream lakes and catchments that affect downstream lake algae contributions and water color (e.g., Hou et al., 2022; Yang et al., 2022).

Despite the improvements over existing datasets identified above and a variety of potential applications of TopoCat, we would like to acknowledge a few notable limitations. First, TopoCat aims to construct a priori drainage topology and connectivity among global lakes. Thus, the “headwater” extent in our delineated inter-lake reach network was determined by the presence of the most upstream lakes rather than river channels. This is partially due to the ambiguity of what and where should be considered a headwater river and the variability of drainage density across the continents (Lin et al., 2021). However, our lake-based drainage network can be potentially nested into a broader and more detailed drainage network to enable a more complete global lake–river integration.

Second, although our algorithm includes a mechanism to consider possible lake multifurcation, and this mechanism design did successfully capture some of the multifurcation cases (see the example in Fig. 3), we acknowledge that our lake multifurcation attribute has not been fully validated. Our identified multifurcation was based on MERIT Hydro. Although MERIT Hydro incorporated long-term maximum water extents to ensure drainage connectivity and to calibrate drainage directions, topographic errors still remain, particularly in floodplains and wetlands with large drainage complexity and uncertainties. Some of the major multifurcation lakes, such as lakes Vesijako (61.38° N, 25° E), Lummene (61.48° N, 25.05° E), and Diefenbaker (51.03° N, 106.84° W), were unfortunately not captured by our automated algorithm. These errors may be due in part to errors in lake mask polygons from HydroLAKES. HydroLAKES does not always represent precise lake boundaries, and some polygons span multiple lakes and even infringe on other rivers and drainage watersheds (leading to artifacts or pseudo bifurcation issues). These errors have not been fully considered in the current TopoCat version and may require future manual QA/QC and validation. A higher-resolution, remote-sensing-based, and quality-controlled lake mask, such as the SWOT prior lake database (Sheng et al., 2016), may also help improve the representation of lake multifurcation.

Third, the accuracy of the delineated basin boundaries in Lake-TopoCat depends on the quality of MERIT-Hydro, which applied an automatic approach to remove dummy depressions from the DEM (Yamazaki et al., 2019). This automated DEM processing occasionally lumped actual small en-

dorheic basins to their adjacent exorheic or endorheic basins, which could lead to incomplete lake boundaries and misclassification of the lake drainage types. For instance, Pangong Lake (33.72° N, 78.9° E) on the Tibetan Plateau is an endorheic lake, but in MERIT Hydro the Pangong basin was incorrectly merged with the Indus Basin.

Last, the current version of Lake-TopoCat emphasizes lake drainage topology and catchments, rather than explicit lake–river integration and segmentation. Therefore, only lake outlets were considered at this moment. While the position of the lake inlet is subject to upstream and downstream shifts as the lake area expands and shrinks, the outlet position is often more stable, particularly in reservoirs and perennial lakes. Located at the most downstream boundary of a lake, the outlet also determines the full catchment domain governing the water mass and quality of the lake, and it is more straightforward to use outlets to configure lake drainage topology and derive unit lake catchments. For these reasons, we consider the delineation of lake outlets to be both more meaningful and practical. However, we admit that the addition of lake inlets is necessary for a more complete segmentation and coupling between lentic and lotic units along the drainage networks; therefore, the inclusion of both lake inlets and outlets will be considered in our future lake–river harmonization efforts.

**Author contributions.** MSS contributed to data curation, formal analysis, investigation, methodology, programming, quality assurance, quality control, validation, visualization, and writing – original draft preparation and review and editing. JW contributed to conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, programming, project administration, quality assurance, quality control, supervision, validation, visualization, and writing – original draft preparation and review and editing. GHA contributed to data curation, methodology, validation, and writing – review and editing. YS contributed to data curation, methodology, and writing – review and editing. DY contributed to data curation, methodology, and writing – review and editing. CS contributed to data curation, methodology, validation, and writing – review and editing. MD contributed to quality control, and writing – review and editing. JFC contributed to methodology, and writing – review and editing. TMP contributed to methodology, and writing – review and editing.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Authors of this paper claim no responsibility or liability for any consequences related to the use, citation, or dissemination of Lake-TopoCat.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
**Acknowledgements.** This work was supported by a NASA Surface Water and Ocean Topography (SWOT) Science Team grant. The authors thank Kai Liu at Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, for providing the validation data for endorheic basin boundaries on the Tibetan Plateau. The authors are also grateful to Bernhard Lehner at McGill University, Michael F. Meyer at USGS, Richard Iestyn Woolway at Bangor University, Daisuke Tokuda at The University of Tokyo, and Alberto V. Borges at University of Liège for their constructive insights into the construction and applications of Lake-TopoCat.

**Financial support.** This research has been supported by the National Aeronautics and Space Administration, NASA Headquarters (grant no. 80NSSC20K1143).

**Review statement.** This paper was edited by Kang Yang and reviewed by two anonymous referees.

## References

Abbott, B. W., Bishop, K., Zarnetske, J. P., Minaudo, C., Chapin III, F. S., Krause, S., Hannah, D. M., Conner, L., Ellison, D., Godsey, S. E., Plont, S., Marçais, J., Kolbe, T., Huebner, A., Frei, R. J., Hampton, T., Gu, S., Buhman, M., Sayedi, S. S., Ursache, O., Chapin, M., Henderson, K. D., and Pinay, G.: Human domination of the global water cycle absent from depictions and perceptions, Nat. Geosci., 12, 533–540, https://doi.org/10.1038/s41561-019-0374-y, 2019.

Abril, G. and Borges, A. V.: Ideas and perspectives: Carbon leaks from flooded land: do we need to replumb the inland water active pipe?, Biogeosciences, 16, 769–784, https://doi.org/10.5194/bg-16-769-2019, 2019.

Allen, G. H. and Pavelsky, T. M.: Global extent of rivers and streams, Science, 361, 585–588, https://doi.org/10.1126/science.aat0636, 2018.

Altenau, E. H., Pavelsky, T. M., Durand M. T., Yang X., Frasson, R. P. M., and Bendezu L.: The Surface Water and Ocean Topography (SWOT) Mission River Database (SWORD): A Global River Network for Satellite Data Products, Water Resour. Res., 57, e2021WR030054, https://doi.org/10.1029/2021WR030054, 2021.

Amatulli, G., Garcia Marquez, J., Sethi, T., Kiesel, J., Grigoropoulou, A., Üblacker, M. M., Shen, L. Q., and Domisch, S.: Hydrography90m: a new high-resolution global hydrographic dataset, Earth Syst. Sci. Data, 14, 4525–4550, https://doi.org/10.5194/essd-14-4525-2022, 2022.

Balmer, M. B. and Downing, J. A.: Carbon dioxide concentrations in eutrophic lakes: undersaturation implies atmospheric uptake, Inland Waters, 1, 125–132, https://doi.org/10.5268/IW-1.2.366, 2011.

Biancamaria, S., Lettenmaier, D. P., and Pavelsky, T. M.: The SWOT Mission and Its Capabilities for Land Hydrology, Surv. Geophys., 37, 307–337, https://doi.org/10.1007/s10712-015-9346-y, 2016.

Borges, A. V., Deirmendjian, L., Bouillon, S., Okello, W., Lambert, T., Roland, F. A. E., Razanamahandry, V. F., Voarintsoa, N. R. G., Darchambeau, F., Kimirei, I. A., Descy, J.-P., Allen, G. H., and Morana, C.: Greenhouse gas emissions from African lakes are no longer a blind spot, Science Advances, 8, eabi8716, https://doi.org/10.1126/sciadv.abi8716, 2022.

Bowling, L. C. and Lettenmaier, D. P.: Modeling the Effects of Lakes and Wetlands on the Water Balance of Arctic Environment, J. Hydrometeorol., 11, 276–295, https://doi.org/10.1175/2009JHM1084.1, 2010.

Cardille, J. A., Carpenter, S. R., Coe, M. T., Foley, J. A., Hanson, P. C., Turner, M. G., and Vano, J. A.: Carbon and water cycling in lake-rich landscapes: Landscape connections, lake hydrology, and biogeochemistry, J. Geophys. Res., 112, G02031, https://doi.org/10.1029/2006JG000200, 2007.

Casas-Ruiz, J. P., Hutchins, R. H. S., and del Giorgio, P. A.: Total Aquatic Carbon Emissions Across the Boreal Biome of Québec Driven by Watershed Slope, J. Geophy. Res.-Biogeo., 126, e2020JG005863, https://doi.org/10.1029/2020JG005863, 2020.

Chen, T., Song, C., Fan, C., Cheng, J., Duan, X., Wang, L., Liu, K., Deng, S., and Che Y.: A comprehensive data set of physical and human-dimensional attributes for China’s lake basins, Scientific Data, 9, 519, https://doi.org/10.1038/s41597-022-01649-z, 2022.

Chen, W., Liu, Y., Zhang, G., Yang, K., Zhou, T., Wang, J., and Shum, C. K.: What controls lake contraction and then expansion in Tibetan Plateau’s endorheic basin over the past half century?, Geophys. Res. Lett., 49, e2022GL101200, https://doi.org/10.1029/2022GL101200, 2022.

Cheruvelil, K. S., Soranno, P. A., McCullough, I. M., Webster, K. E., Rodriguez, L. K., Smith, N. J.: LAGOS-US LOCUS v1.0: Data module of location, identifiers, and physical characteristics of lakes and their watersheds in the conterminous U. S., Limnol. Oceanogr. Lett., 6, 270–292, https://doi.org/10.1002/olol.10203, 2021.

Dolan, W. Pavelsky, T. M., and Yang, X.: Functional Lake-to-Channel Connectivity Impacts Lake Ice in the Colville Delta, Alaska, J. Geophy. Res.-Earth, 126, e2021JF006362, https://doi.org/10.1029/2021JF006362, 2021.

Earth Resources Observation and Science (EROS) Center: HYDRO1k, a global hydrologic database derived from 1996 GTOPO30 data, USGS EROS Center, Sioux Falls, South Dakota [data set], https://doi.org/10.5066/F77P8WN0, 2018.

Farr, T. G., Rosen, P. A., Caro, E., Crippen, R., Duren, R., Hensley, S., Kobrick, M., Paller, M., Rodriguez, E., Roth, L., Seal, D., Shaffer, S., Shimada, J., Umland, J., Werner, M., Oskin, M., Burbank, D., and Alsdorf, D.: The Shuttle Radar Topography Mission, Rev. Geophys., 45, RG2004, https://doi.org/10.1029/2005RG000183, 2007.

Fergus, C. E., Lapierre, J. F., Oliver, S. K., Skaff, N. K., Cheruvelil, K. S., Webster, K., Scott, C., and Soranno, P.: The freshwater landscape: lake, wetland, and stream abundance and connectivity at macroscales, Ecosphere, 8, e01911, https://doi.org/10.1002/ecs2.1911, 2017.

Gardner, J. R., Pavelsky, T. M., and Doyle, M. W.: The Abundance, Size, and Spacing of Lakes and Reservoirs Connected to River Networks, Geophys. Res. Lett., 46, 2592–2601, https://doi.org/10.1029/2018GL080841, 2019.

Grimaldi, S., Nardi, F., Di Benedetto, F., Istanbulluoglu, E., Bras, R. L.: A physically-based method for removing pits in digital elevation models, Adv. Water Resour., 30, 2151–2158, https://doi.org/10.1016/j.advwatres.2006.11.016, 2007.
Han, M., Mai, J., Tolson, B. A., Craig, J. R., Gaborit, É., Liu, H., and Lee, K.: Subwatershed-based lake and river routing products for hydrologic and land surface models applied over Canada, Can. Water Resour. J., 45, 237–251, https://doi.org/10.1080/07011784.2020.1772116, 2020.

Han, M., Shen, H., Tolson, B. A., Craig, J. R., Mai, J., Lin, S. G. M., Basu, N. B., and Awol, F. S.: BasinMaker 3.0: A GIS toolbox for distributed watershed delineation of complex lake-river routing networks, Environ. Modell. Softw., 164, 105688, https://doi.org/10.1016/j.envsoft.2023.105688, 2023.

Hill, R. A., Weber, M. H., Debbout, R. M., Leibowitz, S. G., and Olsen, A. R.: The Lake-Catchment (LakeCat) Dataset: characterizing landscape features for lake basins within the conterminous USA, Freshw. Sci., 37, 208–221, https://doi.org/10.1086/697966, 2018.

Hou, X., Feng, L., Dai, Y., Hu, C., Gibson, L., Tang, J., Lee, Z., Wang, Y., Cai, X., Liu, J., Zheng, Y., and Zheng, C.: Global mapping reveals increase in lacustrine algal blooms over the past decade, Nat. Geosci., 15, 130–134, https://doi.org/10.1038/s41561-021-00887-x, 2022.

Huziy, O. and Sushama, L.: Impact of lake–river connectivity and interflow on the Canadian RCM simulated regional climate and hydrology for Northeast Canada, Clim. Dynam., 48, 709–725, https://doi.org/10.1007/s00382-016-3104-9, 2017.

Khzaei, B., Read, L. K., Casali, M., Sampson, K. M., and Yates, D. N.: GLOBathy, the global lakes bathymetry dataset, Scientific Data, 9, 36, https://doi.org/10.1038/s41597-022-01132-9, 2022.

Lapierre, J.-F., and del Giorgio, P. A.: Geographical and environmental drivers of regional differences in the lake pCO2 versus DOC relationship across northern landscapes, J. Geophy. Res.-Biogeo., 117, G03015, https://doi.org/10.1029/2012JG001945, 2012.

Lehner, B. and 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, https://doi.org/10.1002/hyp.9740, 2013.

Lehner, B., Verdin, K., and Jarvis, A.: New global hydrography derived from spaceborne elevation data, Eos T. Am. Geophys. Un., 89, 93–94, https://doi.org/10.1029/2008EO100001, 2008.

Lehner, B., Messager, M. L., Korver, M. C., and Linke, S.: Global hydro-environmental lake characteristics at high spatial resolution, Scientific Data, 9, 351, https://doi.org/10.1038/s41597-022-01425-z, 2022.

Leibowitz, S. G., Mushet, D. M., and Newton, W. E.: Intermittent Surface Water Connectivity: Fill and Spill Vs. Fill and Merge Dynamics, Wetlands, 36, 323–342, https://doi.org/10.1007/s13157-016-0830-z, 2016.

Leibowitz, S. G., Wigington Jr., P. J., Schofield, K. A., Alexander, L. C., Vanderhoof, M. K., and Golden, H. E.: Connectivity of streams and wetlands to downstream waters: an integrated systems framework, J. Am. Water Resour. As., 54, 298–322, https://doi.org/10.1111/1752-1688.12631, 2018.

Lin, P., Yang, Z. L., Cai, X., and David, C. H.: Development and evaluation of a physically-based lake level model for water resource management: A case study for Lake Buchanan, Texas, Journal of Hydrology: Regional Studies, 4, 661–674, https://doi.org/10.1016/j.ejrh.2015.08.005, 2015.

Lin, P., Pan, M., Wood, E. F., and Allen, G. H.: A new vector-based global river network dataset accounting for variable drainage

density, Scientific Data, 8, 28, https://doi.org/10.1038/s41597-021-00819-9, 2021.

Liu, K., Song, C., Ke, L., Jiang, L., and Ma, R.: Automatic watershed delineation in the Tibetan endorheic basin: A lake-oriented approach based on digital elevation models, Geomorphology, 358, 107127, https://doi.org/10.1016/j.geomorph.2020.107127, 2020.

Liu, K., Ke, L., Wang, J., Jiang, L., Richards, K. S., Sheng, Y., Zhu, Y., Fan, C., Zhan, P., Luo, S., Cheng, J., Chen, T., Ma, R., Liang, Q., Madson, A., and Song, C.: Ongoing Drainage Reorganization Driven by Rapid Lake Growths on the Tibetan Plateau, Geophys. Res. Lett., 48, e2021GL095795, https://doi.org/10.1029/2021GL095795, 2021.

Liu, S., Kuhn, C., Amatulli, G., and Raymond P. A.: The importance of hydrology in routing terrestrial carbon to the atmosphere via global streams and rivers, P. Natl. Acad. Sci. USA, 119, e2106322119, https://doi.org/10.1073/pnas.2106322119, 2022.

Maberly, S., Barker, P., Stott, A., and De Ville, M. M.: Catchment productivity controls CO2 emissions from lakes, Nat. Clim. Change, 3, 391–394, https://doi.org/10.1038/nclimate1748, 2013.

Martin, S. L. and Soranno, P. A.: Lake landscape position: Relationships to hydrologic connectivity and landscape features, Limnol. Oceanogr., 51, 801–814, https://doi.org/10.4319/lo.2006.51.2.0801, 2006.

McDonnell, J. J., Spence, C., Karran, D. J., Meerveld, H. J., and Harman C. J.: Fill-and-Spill: A Process Description of Runoff Generation at the Scale of the Beholder, Water Resour. Res., 57, e2020WR027514, https://doi.org/10.1029/2020WR027514, 2021.

McKay, L., Bondelid, T., Dewald, T., Johnston, J., Moore, R., and Reah, A.: NHDPlus Version 2: user guide, US Geological Survey, Reston, Virginia, https://nhdplus.com/NHDPlus/ (last access: 10 December 2022), 2012.

Mendonça, R., Müller, R. A., Clow, D., Verpoorter, C., Raymond, P., Tranvik, L. J., and Sobek, S.: Organic carbon burial in global lakes and reservoirs, Nat. Commun., 8, 1694, https://doi.org/10.1038/s41467-017-01789-6, 2017.

Messager, M. L., Lehner, B., Grill, G., Nedeva, I., and Schmitt, O.: Estimating the volume and age of water stored in global lakes using a geo-statistical approach, Nat. Commun., 7, 13603, https://doi.org/10.1038/ncomms13603, 2016.

Meyer, M. F., Labou, S. G., Cramer, A. N., Brousil M. R., and Luff, B. T.: The global lake area, climate, and population dataset, Scientific Data, 7, 174, https://doi.org/10.1038/s41597-020-0517-4, 2020.

Oki, T. and Kanae, S.: Global Hydrological Cycles and World Water Resources, Science, 313, 1068–1072, https://doi.org/10.1126/science.1128845, 2006.

Pekel, J.-F., Cottam, A., Gorelick, N., and Belward, A. S.: High-resolution mapping of global surface water and its long-term changes, Nature, 540, 418–422, https://doi.org/10.1038/nature20584, 2016.

Pi, X., Luo, Q., Feng, L., Xu, Y., Tang, J., Liang, X., Ma, E., Cheng, R., Fensholt, R., Brandt, M., Cai, X., Gibson, L., Liu, J., Zheng, C., Li, W., and Bryan, B. A.: Mapping global lake dynamics reveals the emerging roles of small lakes, Nat. Commun., 13, 5777, https://doi.org/10.1038/s41467-022-33239-3, 2022.
Schmadel, N. M., Harvey, J. W., Alexander, R. B., Schwarz, G. E., Moore, R. B., Eng, E., Gomez-Velez, J. D., Boyer, E. W., and Scott, D.: Thresholds of lake and reservoir connectivity in river networks control nitrogen removal, Nat. Commun., 9, 2779, https://doi.org/10.1038/s41467-018-05156-x, 2018.

Seekell, D., Cael, B., Lindmark, E., and Byström, P.: The Fractal Scaling Relationship for River Inlets to Lakes, Geophys. Res. Lett., 48, e2021GL093366, https://doi.org/10.1029/2021GL093366, 2021.

Sheng, Y., Song, C., Wang, J., Lyons, E. A., Knox, B. R., Cox, J. S., and Gao, F.: Representative lake water extent mapping at continental scales using multi-temporal Landsat-8 imagery, Remote Sens. Environ., 185, 129–41, https://doi.org/10.1016/j.rse.2015.12.041, 2016.

Sikder, M. S., Wang, J., Allen, G. H., Sheng, Y., Yamazaki, D., Song, C., Ding, M., Crétaux, J.-F., and Pavelsky, T. M.: Lake-TopoCat: A global Lake drainage Topology and Catchment database (v1.1), Zenodo [data set], https://doi.org/10.5281/zenodo.7916729, 2023.

Stieglitz, M., Shaman, J., McNamara, J., Engel, V., Shanley, J., and Kling, G. W.: An approach to understanding hydrologic connectivity on the hillslope and the implications for nutrient transport, Global Biogeochem. Cy., 17, 1105, https://doi.org/10.1029/2003GB002041, 2003.

Strahler, A. N.: Quantitative analysis of watershed geomorphology, Eos T. Am. Geophys. Un., 38, 913–920, https://doi.org/10.1029/TR038i006p00913, 1957.

Tan, Z., Wang, X., Chen, B., Liu, X., and Zhang Q.: Surface water connectivity of seasonal isolated lakes in a dynamic lake-floodplain system, J. Hydrol., 579, 124154, https://doi.org/10.1016/j.jhydrol.2019.124154, 2019.

Tokuda, D., Kim, H., Yamazaki, D., and Oki, T.: Development of a coupled simulation framework representing the lake and river continuum of mass and energy (TCHOIR v1.0), Geosci. Model Dev., 14, 5669–5693, https://doi.org/10.5194/gmd-14-5669-2021, 2021.

Tranvik, L. J., Downing, J. A., Cotner, J. B., Loiselle, S. A., Striegl, R. G., Ballatore, T. J., Dillon, P., Finlay, K., Fortino, K., Knoll, L. B., Kortelainen, P. L., Kutser, T., Larsen, S., Laurion, I., Leech, D. M., McCallister, S. L., McKnight, D. M., Melack, J. M., Overholt, E., Porter, J. A., Prairie, Y., Renwick, W. H., Roland, F., Sherman, B. S., Schindler, D. W., Sobek, S., Tremblay, A., Vanni, M. J., Verschoor, A. M., von Wachenfeldt, E., and Weyhenmeyer, G. A.: Lakes and reservoirs as regulators of carbon cycling and climate, Limnol. Oceanogr., 54, 2298–2314, https://doi.org/10.4319/lo.2009.54.6_part_2.2298, 2009.

United States Geological Survey (USGS): National Hydrography Dataset (NHD), US Geological Survey, Reston, Virginia, https://nhd.usgs.gov/ (last access: 10 December 2022), 2001.

Vanderhoof, M. K., Alexander L. C., and Todd M. J.: Temporal and spatial patterns of wetland extent influence variability of surface water connectivity in the Prairie Pothole Region, United States, Landscape Ecol., 31, 805–824, https://doi.org/10.1007/s10980-015-0290-5, 2016.

Vanderhoof, M. K., Christensen, J. R., and Alexander, L. C.: Patterns and drivers for wetland connections in the Prairie Pothole Region, United States, Wetl. Ecol. Manag., 25, 275–297, https://doi.org/10.1007/s11273-016-9516-9, 2017.

Wang, J.: Endorheic water, in: The International Encyclopedia of Geography, edited by: Richardson, D., Castree, N., Goodchild, M. F., Kobayashi, A., Liu, W., and Marston, R. A., The American Association of Geographers (AAG), John Wiley & Sons, Ltd, Malden, MA, USA, https://doi.org/10.1002/9781118786352.wbieg2001, 2020.

Wang, J., Song, C., Reager, J. T., Yao, F., Famiglietti, J. S., Sheng, Y., MacDonald, G. M., Brun, F., Müller Schmied, H., Marston, R. A., and Wada, Y.: Recent global decline in endorheic basin water storages, Nat. Geosci., 11, 926–932, https://doi.org/10.1038/s41561-018-0265-7, 2018.

Woolway, R. I., Kraemer, B. M., Lenters, J. D., Merchant, C. J., O'Reilly, C. M., and Sharma, S.: Global lake responses to climate change, Nature Reviews Earth & Environment, 1, 388–403, https://doi.org/10.1038/s43017-020-0067-5, 2020.

Wu, Q. and Lane, C. R.: Delineating wetland catchments and modeling hydrologic connectivity using lidar data and aerial imagery, Hydrol. Earth Syst. Sci., 21, 3579–3595, https://doi.org/10.5194/hess-21-3579-2017, 2017.

Yamazaki, D., Ikeshima, D., Tawatari, R., Yamaguchi, T., O'Loughlin, F., Neal, J. C., Sampson, C. C., Kanae, S., and Bates, P. D.: A high-accuracy map of global terrain elevations, Geophys. Res. Lett., 44, 5844–5853, https://doi.org/10.1002/2017GL072874, 2017.

Yamazaki, D., Ikeshima, D., Sosa, J., Bates, P. D., Allen, G. H., and Pavelsky, T. M.: MERIT Hydro: a high-resolution global hydrography map based on latest topography dataset, Water Resour. Res., 55, 5053–5073, https://doi.org/10.1029/2019WR024873, 2019.

Yan, D., Li, M., Bi, W., Weng, B., Qin, T., Wang, J., and Do, P.: A data set of inland lake catchment boundaries for the Qiangtang Plateau, Scientific Data, 6, 62, https://doi.org/10.1038/s41597-019-0066-x, 2019.

Yang, X., O'Reilly, C. M., Gardner, J. R., Ross, M. R. V., Topp, S. N., Wang, J., and Pavelsky, T. M.: The color of Earth's lakes, Geophys. Res. Lett., 49, e2022GL098925, https://doi.org/10.1029/2022GL098925, 2022.

Yao, F., Wang, J., Yang, K., Wang, C., Walter, B., and Crétaux, J.-F.: Lake storage variation on the endorheic Tibetan Plateau and its attribution to climate change since the new millennium, Environ. Res. Lett., 13, 064011, https://doi.org/10.1088/1748-9326/aab5d3, 2018.

Zhao, G., Li, Y., Zhou, L., and Gao, H.: Evaporative water loss of 1.42 million global lakes, Nat. Commun., 13, 3686, https://doi.org/10.1038/s41467-022-31125-6, 2022.