River network and hydro-geomorphology parametrization for global river routing modelling at 1/12° resolution

Global scale river routing models (RRMs) are commonly used in a variety of studies, including studies on the impact of climate change on extreme flows (floods and droughts), water resources monitoring or large scale flood forecasting. Over the last two decades, the increasing number of observational datasets, mainly from satellite missions, and the increasing computing capacities, have allowed better performances of RRMs, namely by increasing their spatial resolution. The spatial resolution of a RRM corresponds to the spatial resolution of its river network, which provides flow direction of all grid cells. River networks may be derived at various spatial resolution by upscaling high resolution hydrography data. This paper presents a new global scale river network at 1/12° derived from the MERIT-Hydro dataset. The river network is generated automatically using an adaptation of the Hierarchical Dominant River Tracing (DRT) algorithm, and its quality is assessed over the 70 largest basins of the world. Although this new river network may be used for a variety of hydrology-related studies, it is here provided with a set of hydro-geomorphological parameters at the same spatial resolution. These parameters are derived during the generation of the river network and are based on the same high resolution dataset, so that the consistency between the river network and the parameters is ensured. The set of parameters includes a description of river stretches (length, slope, width, roughness, bankfull depth), floodplains (roughness, sub-grid topography) and aquifers (transmissivity, porosity, sub-grid topography). The new river network and parameters are assessed by comparing the performances of two global scale simulations with the CTRIP model, one with the current spatial resolution (1/2°) and the other with the new spatial resolution (1/12°). It is shown that CTRIP at 1/12° overall outperforms CTRIP at 1/2°, demonstrating the added value of the spatial resolution increase. The new river network and the consistent hydro-geomorphology parameters may be useful for the scientific community, especially for hydrology and hydro-geology modelling, water resources monitoring or climate studies.

HydroBlocks v0.2: enabling a field-scale two-way coupling between the land surface and river networks in Earth system models

Over the past decade, there has been appreciable progress towards modeling the water, energy, and carbon cycles at field scales (10–100 m) over continental to global extents in Earth system models (ESMs). One such approach, named HydroBlocks, accomplishes this task while maintaining computational efficiency via Hydrologic Response Units (HRUs), more commonly known as “tiles” in ESMs. In HydroBlocks, these HRUs are learned via a hierarchical clustering approach from available global high-resolution environmental data. However, until now there has yet to be a river routing approach that is able to leverage HydroBlocks' approach to modeling field-scale heterogeneity; bridging this gap will make it possible to more formally include riparian zone dynamics, irrigation from surface water, and interactive floodplains in the model. This paper introduces a novel dynamic river routing scheme in HydroBlocks that is intertwined with the modeled field-scale land surface heterogeneity. Each macroscale polygon (a generalization of the concept of macroscale grid cell) is assigned its own fine-scale river network that is derived from very high resolution (∼ 30 m) digital elevation models (DEMs); the inlet–outlet reaches of a domain's macroscale polygons are then linked to assemble a full domain's river network. The river dynamics are solved at the reach-level via the kinematic wave assumption of the Saint-Venant equations. Finally, a two-way coupling between each HRU and its corresponding fine-scale river reaches is established. To implement and test the novel approach, a 1.0∘ bounding box surrounding the Atmospheric Radiation and Measurement (ARM) Southern Great Plains (SGP) site in northern Oklahoma (United States) is used. The results show (1) the implementation of the two-way coupling between the land surface and the river network leads to appreciable differences in the simulated spatial heterogeneity of the surface energy balance, (2) a limited number of HRUs (∼ 300 per 0.25∘ cell) are required to approximate the fully distributed simulation adequately, and (3) the surface energy balance partitioning is sensitive to the river routing model parameters. The resulting routing scheme provides an effective and efficient path forward to enable a two-way coupling between the high-resolution river networks and state-of-the-art tiling schemes in ESMs.

Parametrization of a lake water dynamics model MLake in the ISBA-CTRIP land surface system (SURFEX v8.1)

Lakes are of fundamental importance in the Earth system as they support essential environmental and economic services, such as freshwater supply. Streamflow variability and temporal evolution are impacted by the presence of lakes in the river network; therefore, any change in the lake state can induce a modification of the regional hydrological regime. Despite the importance of the impact of lakes on hydrological fluxes and the water balance, a representation of the mass budget is generally not included in climate models and global-scale hydrological modeling platforms. The goal of this study is to introduce a new lake mass module, MLake (Mass-Lake model), into the river-routing model CTRIP to resolve the specific mass balance of open-water bodies. Based on the inherent CTRIP parameters, the development of the non-calibrated MLake model was introduced to examine the influence of such hydrological buffer areas on global-scale river-routing performance. In the current study, an offline evaluation was performed for four river networks using a set of state-of-the-art quality atmospheric forcings and a combination of in situ and satellite measurements for river discharge and lake level observations. The results reveal a general improvement in CTRIP-simulated discharge and its variability, while also generating realistic lake level variations. MLake produces more realistic streamflows both in terms of daily and seasonal correlation. Excluding the specific case of Lake Victoria having low performances, the mean skill score of Kling–Gupta efficiency (KGE) is 0.41 while the normalized information contribution (NIC) shows a mean improvement of 0.56 (ranging from 0.15 to 0.94). Streamflow results are spatially scale-dependent, with better scores associated with larger lakes and increased sensitivity to the width of the lake outlet. Regarding lake level variations, results indicate a good agreement between observations and simulations with a mean correlation of 0.56 (ranging from 0.07 to 0.92) which is linked to the capability of the model to retrieve seasonal variations. Discrepancies in the results are mainly explained by the anthropization of the selected lakes, which introduces high-frequency variations in both streamflows and lake levels that degraded the scores. Anthropization effects are prevalent in most of the lakes studied, but they are predominant for Lake Victoria and are the main cause for relatively low statistical scores for the Nile River However, results on the Angara and the Neva rivers also depend on the inherent gap of ISBA-CTRIP process representation, which relies on further development such as the partitioned energy budget between the snow and the canopy over a boreal zone. The study is a first step towards a global coupled land system that will help to qualitatively assess the evolution of future global water resources, leading to improvements in flood risk and drought forecasting.

River parametrisation of the JULES land surface model for improved runoff routing at the global scale.

<p>The JULES land surface model has a wide ranging application in studying different processes of the earth system including hydrological modeling [1]. Our aim is to tune the existing configuration of the global river routing scheme at 0.5<sup>o</sup> spatial resolution [2] and improve river flow simulation performance at finer temporal scales. To do so, we develop a factorial experiment of varying effective river velocity and meander coefficient, components of the Total Runoff Integrating Pathways (TRIP) river routing scheme. We test and adjust best performing configurations at the basin scale based on observations from GRDC 230 stations that exhibiting a variety of hydroclimatic and physiographic conditions. The analysis was focused on watersheds of near-natural conditions [3] to avoid potential influences of human management on river flow. The HydroATLAS database [4] was employed to identify basin scale descriptive hydro-environmental indicators that could be associated with the components of the TRIP. These indicators summarize hydrologic and physiographic characteristics of the drainage area of each flow gauge. For each basin we select the best performing set of TRIP parameters per basin resulting to the optimal efficiency of river flow simulation based on the Nash–Sutcliffe and Kling–Gupta efficiency metrics. We find that better performance is driven predominantly by characteristics related to the stream gradient and terrain slope. These indicators can serve as descriptors for extrapolating the adjustment of TRIP parameters for global land configurations at 0.5<sup>o</sup> spatial resolution using regression models.</p><p> </p><p>[1] Papadimitriou et al 2017, Hydrol. Earth Syst. Sci., 21, 4379–4401</p><p>[2] Falloon et al 2007. Hadley Centre Tech. Note 72, 42 pp.</p><p>[3] Fang Zhao et al 2017 Environ. Res. Lett. 12 075003</p><p>[4] Linke et al 2019, Scientific Data 6: 283.</p>

Parametrization of lakes water dynamics in the ISBA-CTRIP land surface system (SURFEX v8.1)

Lakes are of fundamental importance in the Earth system as they support essential environmental and economic services such as freshwater supply. They also modify the local hydro-meteorological continuum as a lower boundary of the atmosphere. Sentinels of climate change and anthropization, these open water bodies are facing disruptions of their equilibrium generally leading to a notable reduction of their levels worldwide. Stream-flow variability and temporal evolution are impacted by the presence of lakes in the river network, therefore any change in the lake state can induce a modification of the regional hydrological regime. Despite the importance of the impact of lakes on hydrological fluxes and the water balance, a representation of the mass budget is generally not included in climate models and global scale hydrological modeling platforms. The goal of this study is to introduce a new lake mass module, MLake (Mass-Lake model), into the river routing model CTRIP to resolve the specific mass-balance of open water bodies. Based on the inherent CTRIP parameters, the development of the non-calibrated MLake model was introduced to examine the influence of such hydrological buffer areas on the global scale river routing performances. In the current study, an off-line evaluation was performed for four river networks using a set of state-of-the-art quality atmospheric forcings and a combination of in situ and satellite measurements for river discharge and lake level observations. The results reveal a general improvement in CTRIP simulated discharge and its variability, while also generating realistic lake level variations. MLake produces more realistic streamflows both in terms of daily and seasonal correlation. Excluding the specific case of Lake Victoria having low performances, the mean skill score of Kling-Gupta Efficiency (KGE) is 0.41 while the Normalized Information Contribution (NIC) shows a mean improvement of 0.56 (ranging from 0.15 to 0.94). Streamflow results are spatially scale-dependent, with better scores associated with larger lakes, and increased sensitivity to the width of the lake outlet. Regarding lake levels variations, results indicate a good agreement between observations and simulations with a mean correlation of 0.56 (ranging from 0.07 to 0.92) which is linked to the capability of the model to retrieve seasonal variations. Discrepancies in the results are mainly explained by the anthropization of the selected lakes which introduces high-frequency variations in both streamflows and lake levels that degraded the scores. Anthropization effects are prevalent in most of the lakes studied, but they are predominant for Lake Victoria and are the main cause for relatively low statistical scores for this river. However, results on the Angara and the Neva rivers also depend on the inherent gap of ISBA-CTRIP processes representation which relies on further development such as the partitioned energy budget between the snow and the canopy over a Boreal zone. The study is a first step towards a global coupled land system that will help to qualitatively assess the evolution of future global water resources leading to improvements in flood risk and drought forecasting.

Two-way coupling between the sub-grid land surface and river networks in Earth system models

Over the past decade, there has been appreciable progress towards modeling the water, energy, and carbon cycles at field-scales (10–100 m) over continental to global extents. One such approach, named HydroBlocks, accomplishes this task while maintaining computational efficiency via sub-grid tiles, or Hydrologic Response Units (HRUs), learned via a hierarchical clustering approach from available global high-resolution environmental data. However, until now, there has yet to be a macroscale river routing approach that is able to leverage HydroBlocks' approach to sub-grid heterogeneity, thus limiting the added value of field-scale land surface modeling in Earth System Models (e.g., riparian zone dynamics, irrigation from surface water, and interactive floodplains). This paper introduces a novel dynamic river routing scheme in HydroBlocks that is intertwined with the modeled field-scale land surface heterogeneity. The primary features of the routing scheme include: 1) the fine-scale river network of each macroscale grid cell's is derived from very high resolution (