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

2020 ◽  
Author(s):  
Nathaniel W. Chaney ◽  
Laura Torres-Rojas ◽  
Noemi Vergopolan ◽  
Colby K. Fisher

Abstract. 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 (

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 ◽  
Vol 14 (11) ◽  
pp. 6813-6832
Author(s):  
Nathaniel W. Chaney ◽  
Laura Torres-Rojas ◽  
Noemi Vergopolan ◽  
Colby K. Fisher

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


2013 ◽  
Vol 6 (1) ◽  
pp. 255-296
Author(s):  
C. Ottlé ◽  
J. Lescure ◽  
F. Maignan ◽  
B. Poulter ◽  
T. Wang ◽  
...  

Abstract. High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover datasets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing earth system models. Earth system models also require specific land cover classification systems based on plant functional types, rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover datasets against one another and with auxiliary data to identify key uncertainties that contribute to variability in Plant Functional Type (PFT) classifications that would introduce errors in earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, continuous vegetation fields (MODIS-VCF) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover dataset, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFTs maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT dataset, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground-truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to better represent the water and carbon cycles in northern latitudes. Updated land cover datasets are critical for improving and maintaining the relevance of earth system models for assessing climate and human impacts on biogeochemistry and biophysics. The new PFT map at 5 km scale is available for download from the PANGAEA website, at: doi:10.1594/PANGAEA.810709.


2018 ◽  
Vol 54 (12) ◽  
Author(s):  
Wondmagegn Yigzaw ◽  
Hong‐Yi Li ◽  
Yonas Demissie ◽  
Mohamad I. Hejazi ◽  
L. Ruby Leung ◽  
...  

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


2018 ◽  
Vol 11 (12) ◽  
pp. 4965-4985 ◽  
Author(s):  
Trung Nguyen-Quang ◽  
Jan Polcher ◽  
Agnès Ducharne ◽  
Thomas Arsouze ◽  
Xudong Zhou ◽  
...  

Abstract. The river routing scheme (RRS) in the Organising Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) land surface model is a valuable tool for closing the water cycle in a coupled environment and for validating the model performance. This study presents a revision of the RRS of the ORCHIDEE model that aims to benefit from the high-resolution topography provided by the Hydrological data and maps based on SHuttle Elevation Derivatives at multiple Scales (HydroSHEDS), which is processed to a resolution of approximately 1 km. Adapting a new algorithm to construct river networks, the new RRS in ORCHIDEE allows for the preservation of as much of the hydrological information from HydroSHEDS as the user requires. The evaluation focuses on 12 rivers of contrasting size and climate which contribute freshwater to the Mediterranean Sea. First, the numerical aspect of the new RRS is investigated, in order to identify the practical configuration offering the best trade-off between computational cost and simulation quality for ensuing validations. Second, the performance of the new scheme is evaluated against observations at both monthly and daily timescales. The new RRS satisfactorily captures the seasonal variability of river discharge, although important biases stem 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 the promising performance of this high-resolution RRS with respect to replicating river flow variation at various frequencies. Furthermore, this RRS may also eventually be 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 ◽  
Vol 21 (1) ◽  
pp. 217-233 ◽  
Author(s):  
Elham Rouholahnejad Freund ◽  
James W. Kirchner

Abstract. Most Earth system models are based on grid-averaged soil columns that do not communicate with one another, and that average over considerable sub-grid heterogeneity in land surface properties, precipitation (P), and potential evapotranspiration (PET). These models also typically ignore topographically driven lateral redistribution of water (either as groundwater or surface flows), both within and between model grid cells. Here, we present a first attempt to quantify the effects of spatial heterogeneity and lateral redistribution on grid-cell-averaged evapotranspiration (ET) as seen from the atmosphere over heterogeneous landscapes. Our approach uses Budyko curves, as a simple model of ET as a function of atmospheric forcing by P and PET. From these Budyko curves, we derive a simple sub-grid closure relation that quantifies how spatial heterogeneity affects average ET as seen from the atmosphere. We show that averaging over sub-grid heterogeneity in P and PET, as typical Earth system models do, leads to overestimations of average ET. For a sample high-relief grid cell in the Himalayas, this overestimation bias is shown to be roughly 12 %; for adjacent lower-relief grid cells, it is substantially smaller. We use a similar approach to derive sub-grid closure relations that quantify how lateral redistribution of water could alter average ET as seen from the atmosphere. We derive expressions for the maximum possible effect of lateral redistribution on average ET, and the amount of lateral redistribution required to achieve this effect, using only estimates of P and PET in possible source and recipient locations as inputs. We show that where the aridity index P/PET increases with altitude, gravitationally driven lateral redistribution will increase average ET (and models that overlook lateral redistribution will underestimate average ET). Conversely, where the aridity index P/PET decreases with altitude, gravitationally driven lateral redistribution will decrease average ET. The effects of both sub-grid heterogeneity and lateral redistribution will be most pronounced where P is inversely correlated with PET across the landscape. Our analysis provides first-order estimates of the magnitudes of these sub-grid effects, as a guide for more detailed modeling and analysis.


2020 ◽  
Author(s):  
Norman Steinert ◽  
Fidel González-Rouco ◽  
Stefan Hagemann ◽  
Philipp de Vrese ◽  
Elena García-Bustamante ◽  
...  

<p>The representation of the thermal and hydrological state in the land model component of Earth System Models is crucial to have a realistic simulation of subsurface processes and the coupling between the atmo-, lito- and biosphere. There is evidence suggesting an inaccurate simulation of subsurface thermodynamics in current-generation Earth System Models, which have Land Surface Models that are too shallow. In simulations with a bottom boundary too close to the surface, the energy propagation and spatio-temporal variability of subsurface temperatures are affected. This potentially restrains the simulation of land-air interactions and subsurface phenomena, e.g. energy/moisture balance and storage capacity, freeze/thaw cycles and permafrost evolution. We introduce modifications for a deeper soil into the JSBACH soil model component of the MPI-ESM for climate projections of the 21st century. Subsurface layers are added progressively to increase the bottom boundary depth from 10m to 1400m. This leads to near-surface cooling of the soil and encourages regional terrestrial energy uptake by one order of magnitude and more. <br>The depth-changes in the soil also have implications for the hydrological regime, in which the moisture between the surface and the bedrock is sensitive to variations in the thermal regime. Additionally, we compare two different global soil parameter datasets that have major implications for the vertical distribution and availability of soil moisture and its exchange with the land surface. The implementation of supercool water and water phase changes in the soil creates a coupling between the soil thermal and hydrological regimes. In both cases of bottom boundary and water depth changes, we explore the sensitivity of JSBACH from the perspective of changes in the soil thermodynamics, energy balance and storage, as well as the effect of including freezing and thawing processes and their influence on the simulation of permafrost areas in the Northern Hemisphere high latitudes. The latter is of particular interest due to their vulnerability to long-term climate change.</p>


2018 ◽  
Vol 22 (6) ◽  
pp. 3311-3330 ◽  
Author(s):  
Nathaniel W. Chaney ◽  
Marjolein H. J. Van Huijgevoort ◽  
Elena Shevliakova ◽  
Sergey Malyshev ◽  
Paul C. D. Milly ◽  
...  

Abstract. The continual growth in the availability, detail, and wealth of environmental data provides an invaluable asset to improve the characterization of land heterogeneity in Earth system models – a persistent challenge in macroscale models. However, due to the nature of these data (volume and complexity) and computational constraints, these data are underused for global applications. As a proof of concept, this study explores how to effectively and efficiently harness these data in Earth system models over a 1/4∘ (∼ 25 km) grid cell in the western foothills of the Sierra Nevada in central California. First, a novel hierarchical multivariate clustering approach (HMC) is introduced that summarizes the high-dimensional environmental data space into hydrologically interconnected representative clusters (i.e., tiles). These tiles and their associated properties are then used to parameterize the sub-grid heterogeneity of the Geophysical Fluid Dynamics Laboratory (GFDL) LM4-HB land model. To assess how this clustering approach impacts the simulated water, energy, and carbon cycles, model experiments are run using a series of different tile configurations assembled using HMC. The results over the test domain show that (1) the observed similarity over the landscape makes it possible to converge on the macroscale response of the fully distributed model with around 300 sub-grid land model tiles; (2) assembling the sub-grid tile configuration from available environmental data can have a large impact on the macroscale states and fluxes of the water, energy, and carbon cycles; for example, the defined subsurface connections between the tiles lead to a dampening of macroscale extremes; (3) connecting the fine-scale grid to the model tiles via HMC enables circumvention of the classic scale discrepancies between the macroscale and field-scale estimates; this has potentially significant implications for the evaluation and application of Earth system models.


Sign in / Sign up

Export Citation Format

Share Document