An enhanced river routing scheme for the closure of global water budget

Author(s):  
Stefano Materia ◽  
Daniele Peano ◽  
Marianna Benassi ◽  
Tomas Lovato ◽  
Silvio Gualdi ◽  
...  

<p>Large-scale river routing schemes are essential to close the hydrological cycle in fully coupled<span> </span>Earth System Models (ESMs). The availability of a realistic water flow is a powerful instrument to evaluate modeled land surface, a crucial component of the global climate whos<span>e </span>properties are often simplified by heavy parameterization, due to lack of process knowledge and<span> </span>validation data.<span> </span>We built up a new concept of river routing model, named HYDROS (HYdro-Dynamic ROuting Scheme), that replaces the present scheme embedded in the CMCC-CM2 global coupled model.<span> </span>The new scheme aims at overcoming one of the current major limitations,<span> </span>that is the use of time-independent flow velocities parameterized as a function of topography.<span> </span>Through the imposition of hydraulic equations, HYDROS defines a time-varying flow velocity<span> </span>associated with the amount of lateral runoff generated by the ESM's land component and the flow through the<span> </span>river system.<span> Compared to the scheme currently in place,</span> HYDROS show improvements in the simulation of mean annual discharge phase, especially for the Arctic rivers and<span> </span>the Amazon. In <span>the Mississippi case</span>,<span> </span><span>an extreme</span> flood episode <span>is</span> better caught by the new representation, indicating that the improved flow velocity better catches the discharge peaks after extreme rainfalls. The new routing model<span> </span>is not able to improve the volumes of simulated river discharge, <span>whose</span> magnitude depends on<span> </span>the ability of the ESM land surface scheme to generate correct surface and sub-surface runoff. Once implemented in coupled mode, HYDROS will guarantee a plausible amount and timing of freshwater discharge into the global ocean, unveiling possible unresolved feedback mechanisms occurring in proximity of river mouths<span>.</span></p>

2020 ◽  
Author(s):  
Nathaniel Chaney ◽  
Noemi Vergopolan ◽  
Colby Fisher

<p>Over the past decade there has been important progress towards modeling the water, energy, and carbon cycles at field scales (10-100 meter) over continental extents. One such approach, named HydroBlocks, accomplishes this task while maintaining computational efficiency via sub-grid hydrologic response units (HRUs); these HRUs are defined via cluster analysis of available field-scale environmental datasets (e.g., elevation). However, until now, there has yet to be complementary advances in river routing schemes that are able to fully harness HydroBlocks’ approach to sub-grid heterogeneity, thus limiting the added value of field-scale resolving land surface models (e.g., riparian zone dynamics, irrigation from surface water, and interactive floodplains). In this presentation, we will introduce a novel large scale river routing scheme that leverages the modeled field-scale heterogeneity in HydroBlocks through more realistic sub-grid stream network topologies, reach-based river routing, and the simulation of floodplain dynamics.</p><p>The primary features of the novel river routing scheme include: 1) each macroscale grid cell is assigned its own river network delineated from field-scale DEMs; 2) similar sub-grid reaches (e.g., Shreve order) are grouped/clustered to ensure computational tractability; 3) the fine-scale inlet/outlet reaches of the macroscale grid cells are linked to assemble the continental river networks; 4) river dynamics are solved at the reach-level via an implicit solution of the Kinematic wave with floodplain dynamics; 5) two way connectivity is established between each cell’s sub-grid HRUs and the river network. The resulting routing scheme is able to effectively represent sub-100 meter-delineated stream networks within Earth system models with relatively minor increases in computation with respect to existing approaches. To illustrate the scheme’s novelty when coupled to the HydroBlocks land surface model, we will present simulation results over the Yellowstone river in the United States between 2002 and 2018. We will show the added value of the scheme when compared to existing approaches with regards to floodplain dynamics, water management, and riparian corridors. Furthermore, we will present results regarding the scheme’s computational tractability to ensure the feasibility of its use within Earth system models. Finally, we will discuss the potential of this approach to enhance flood and drought monitoring tools, numerical weather prediction, and climate models.</p>


2021 ◽  
Author(s):  
Aristeidis Koutroulis ◽  
Manolis Grillakis ◽  
Camilla Mathison ◽  
Eleanor Burke

<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>


2022 ◽  
Author(s):  
Qing-Bin Lu

Abstract Time-series observations of global lower stratospheric temperature (GLST), global land surface air temperature (LSAT), global mean surface temperature (GMST), sea ice extent (SIE) and snow cover extent (SCE), together with observations reported in Paper I, combined with theoretical calculations of GLSTs and GMSTs, have provided strong evidence that ozone depletion and global climate changes are dominantly caused by human-made halogen-containing ozone-depleting substances (ODSs) and greenhouse gases (GHGs) respectively. Both GLST and SCE have become constant since the mid-1990s and GMST/LSAT has reached a peak since the mid-2000s, while regional continued warmings at the Arctic coasts (particularly Russia and Alaska) in winter and spring and at some areas of Antarctica are observed and can be well explained by a sea-ice-loss warming amplification mechanism. The calculated GMSTs by the parameter-free warming theory of halogenated GHGs show an excellent agreement with the observed GMSTs after the natural El Niño southern oscillation (ENSO) and volcanic effects are removed. These results provide a convincing mechanism of global climate change and will make profound changes in our understanding of atmospheric processes. This study also emphasizes the critical importance of continued international efforts in phasing out all anthropogenic halogenated ODSs and GHGs.


2018 ◽  
Author(s):  
Trung Nguyen-Quang ◽  
Jan Polcher ◽  
Agnès Ducharne ◽  
Thomas Arsouze ◽  
Xudong Zhou ◽  
...  

Abstract. This study presents a revised river routing scheme (RRS) for the Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model. The revision is carried out to benefit from the high resolution topography provided the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), processed to a resolution of approximately 1 kilometer. The RRS scheme of the ORCHIDEE uses a unit-to-unit routing concept which allows to preserve as much of the hydrological information of the HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasted size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the revised scheme is evaluated against observations at both monthly and daily timescales. The new RRS captures satisfactorily the seasonal variability of river discharges, although important biases come from the water budget simulated by the ORCHIDEE model. The results highlight that realistic streamflow simulations require accurate precipitation forcing data and a precise river catchment description over a wide range of scales, as permitted by the new RRS. Detailed analyses at the daily timescale show promising performances of this high resolution RRS for replicating river flow variation at various frequencies. Eventually, this RRS is well adapted for further developments in the ORCHIDEE land surface model to assess anthropogenic impacts on river processes (e.g. damming for irrigation operation).


2017 ◽  
Author(s):  
Imme Benedict ◽  
Chiel C. van Heerwaarden ◽  
Albrecht H. Weerts ◽  
Wilco Hazeleger

Abstract. The hydrological cycle of river basins can be simulated by combining global climate models (GCMs) and global hydrological models (GHMs). The spatial resolution of these models is restricted by computational resources and therefore limits the processes and level of detail that can be resolved. To further improve simulations of precipitation and river-runoff on a global scale, we assess and compare the benefits of an increased resolution for a GCM and a GHM. We focus on the Rhine and Mississippi basin. Increasing the resolution of a GCM (1.125° to 0.25°) results in more realistic large-scale circulation patterns over the Rhine and an improved precipitation budget. These improvements with increased resolution are not found for the Mississippi basin, most likely because precipitation is strongly dependent on the representation of still unresolved convective processes. Increasing the resolution of vegetation and orography in the high resolution GHM (from 0.5° to 0.05°) shows no significant differences in discharge for both basins, because the hydrological processes depend highly on other parameter values that are not readily available at high resolution. Therefore, increasing the resolution of the GCM provides the most straightforward route to better results. This route works best for basins driven by large-scale precipitation, such as the Rhine basin. For basins driven by convective processes, such as the Mississippi basin, improvements are expected with even higher resolution convection permitting models.


2012 ◽  
Vol 8 (5) ◽  
pp. 5293-5340 ◽  
Author(s):  
I. Nikolova ◽  
Q. Yin ◽  
A. Berger ◽  
U. K. Singh ◽  
M. P. Karami

Abstract. This paper presents a detailed analysis of the climate of the last interglacial simulated by two climate models of different complexities, LOVECLIM and CCSM3. The simulated surface temperature, hydrological cycle, vegetation and ENSO variability during the last interglacial are analyzed through the comparison with the simulated Pre-Industrial (PI) climate. In both models, the last interglacial period is characterized by a significant warming (cooling) over almost all the continents during boreal summer (winter) leading to a largely increased (reduced) seasonal contrast in the northern (southern) hemisphere. This is mainly due to the much higher (lower) insolation received by the whole Earth in boreal summer (winter) during this interglacial. The arctic is warmer than PI through the whole year, resulting from its much higher summer insolation and its remnant effect in the following fall-winter through the interactions between atmosphere, ocean and sea ice. In the tropical Pacific, the change in the SST annual cycle is suggested to be related to a minor shift towards an El Nino, slightly stronger for MIS-5 than for PI. Intensified African monsoon and vegetation feedback are responsible for the cooling during summer in North Africa and Arabian Peninsula. Over India precipitation maximum is found further west, while in Africa the precipitation maximum migrates further north. Trees and grassland expand north in Sahel/Sahara. A mix of forest and grassland occupies continents and expand deep in the high northern latitudes. Desert areas reduce significantly in Northern Hemisphere, but increase in North Australia. The simulated large-scale climate change during the last interglacial compares reasonably well with proxy data, giving credit to both models and reconstructions. However, discrepancies exist at some regional scales between the two models, indicating the necessity of more in depth analysis of the models and comparisons with proxy data.


2021 ◽  
Author(s):  
Michiel Maertens ◽  
Veerle Vanacker ◽  
Gabriëlle De Lannoy ◽  
Frederike Vincent ◽  
Raul Giménez ◽  
...  

<p>The South-American Dry Chaco is a unique ecoregion as it is one of the largest sedimentary plains in the world hosting the planet’s largest dry forest. The 787.000 km² region covers parts of Argentina, Paraguay, and Bolivia and is characterized by a negative climatic water balance as a consequence of limited rainfall inputs (800 mm/year) and high temperatures (21°C). In combination with the region’s extreme flat topography (slopes < 0.1%) and shallow groundwater tables, saline soils are expected in substantial parts of the region. In addition, it is expected that large-scale deforestation processes disrupt the hydrological cycle resulting in rising groundwater tables and further increase the risk for soil salinization.</p><p>In this study, we identified the regional-scale patterns of subsurface soil salinity in the Dry Chaco.  Field data were obtained during a two-month field campaign in the dry season of 2019. A total of 492 surface- and 142 subsurface-samples were collected along East-West transects to determine soil electric conductivity, pH, bulk density and humidity. Spatial regression techniques were used to reveal the topographic and ecohydrological variables that are associated with subsurface soil salinity over the Dry Chaco. The hydrological information was obtained from a state-of-the-art land surface model with an improved set of satellite-derived vegetation and land cover parameters.</p><p>In the presentation, we will present a subsurface soil salinity map for a part of the Argentinean Dry Chaco and provide relevant insights into the driving mechanisms behind it.</p>


Author(s):  
Jose A. Marengo ◽  
Carlos A. Nobre

The Amazon region is of particular interest because it represents a large source of heat in the tropics and has been shown to have a significant impact on extratropical circulation, and it is Earth’s largest and most intense land-based convective center. During the Southern Hemisphere summer when convection is best developed, the Amazon basin is one of the wettest regions on Earth. Amazonia is of course not isolated from the rest of the world, and a global perspective is needed to understand the nature and causes of climatological anomalies in Amazonia and how they feed back to influence the global climate system. The Amazon River system is the single, largest source of freshwater on Earth. The flow regime of this river system is relatively unimpacted by humans (Vörösmarty et al. 1997 a, b) and is subject to interannual variability in tropical precipitation that ultimately is translated into large variations in downstream hydrographs (Marengo et al. 1998a, Vörösmarty et al. 1996, Richey et al. 1989a, b). The recycling of local evaporation and precipitation by the forest accounts for a sizable portion of the regional water budget (Nobre et al. 1991, Eltahir 1996), and as large areas of the basin are subject to active deforestation there is grave concern about how such land surface disruptions may affect the water cycle in the tropics (see reviews in Lean et al. 1996). Previous studies have emphasized either how large-scale atmospheric circulation or land surface conditions can directly control the seasonal changes in rainfall producing mechanisms. Studies invoking controls of convection and rainfall by large-scale circulation emphasize the relationship between the establishment of upper-tropospheric circulation over Bolivia and moisture transport from the Atlantic ocean for initiation of the wet season and its intensity (see reviews in Marengo et al. 1999). On the other hand, Eltahir and Pal (1996) have shown that Amazon convection is closely related to land surface humidity and temperature, while Fu et al. (1999) indicate that the wet season in the Amazon basin is controlled by both changes in land surface temperature and the sea surface temperature (SST) in the adjacent oceans, depending if the region is north-equatorial or southern Amazonia.


2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).


1999 ◽  
Vol 33 (1) ◽  
pp. 81-84
Author(s):  
Jinro Ukila ◽  
Moloyoshi Ikeda

The Frontier Research System for Global Change—the International Arctic Research Center (Frontier-IARC) is a research program funded by the Frontier Research System for Global Change. The program is jointly run under a cooperative agreement between the Frontier Research System for Global Change and the University of Alaska Fairbanks. The aim of the program is to understand the role of the Arctic region in global climate change. The program concentrates its research effort initially on the areas of air-sea-ice interactions, bio-geochemical processes and the ecosystem. To understand the arctic climate system in the context of global climate change, we focus on mechanisms controlling arctic-subarctic interactions, and identify three key components: the freshwater balance, the energy balance, and the large-scale atmospheric processes. Knowledge of details of these components and their interactions will be gained through long-term monitoring, process studies, and modeling; our focus will be on the latter two categories.


Sign in / Sign up

Export Citation Format

Share Document