1. Introduction

RHESSys has a hierarchical structure with three main modules that allow hydrological, microclimate, and ecological processes to be simulated separately in different layers and to reflect the multi-scale feature of watersheds^[1][2] [1,2]. By parametrizing regional eco-hydrological processes through a series of coupled physical mechanism models at different levels, RHESSys can model the interactions between hydrological, climatic and ecosystem processes in a watershed, and thus, can simulate the regional multi-components cycle of nitrogen, carbon, and water^[2] [2]. In particular, RHESSys has been applied for various research fields due to the characteristics of the hierarchical structure and coupled multiple physical processes. First, the hierarchical structure defined by multiple processes improves the simulation efficiency as the multiple processes are operated individually at multiple spatio-temporal scales. In addition, RHESSys incorporates a plant physiological process that simulates the carbon and nitrogen cycling of vegetation and soil to reflect the nonlinear ecosystem response^[2] [2]. Moreover, RHESSys has a flexible structure to be coupled with other models (e.g., WMFire, phenology models)^[3][4][5] [3–5] and consequently is adapted for various research fields. In summary, RHESSys has been designed to assess the interactions between vegetation and water for ecohydrological research and sustainable management of watersheds by simulating regional carbon, nitrogen, and water cycles and distributions. 

Over the last 30 years, RHESSys has been continuously advanced in model structures and algorithms, and applied for various basins to support local water resource management ^[6][7][8][6–8]. For example, Zabalza–Martinez et al. [8] applied RHESSys to simulate hydrologic responses to climate and land-use change scenarios for a basin controlled by the Boadella–Darnius Dam in Spain and suggested water resource management strategies for the reservoir and corresponding sub-basins. Peng et al.^[9] [9] evaluated the impacts of soil and water conservation measures on the runoff of the Jinghe basin in China for not only filling gaps of local assessment of the effectiveness of soil and water conservation measures but also supporting watershed management. Martin et al.^[6] [6] also simulated the water yield of the Yadkin-Pee basin in North Carolina, USA to evaluate the impacts of climate change and human activities for reasonable water resource management. As many studies have shown the advantages and limitations of RHESSys, a systematic evaluation of the application progress of RHESSys can provide useful and scientific information for the in-depth understanding of the strengths and weaknesses of the RHESSys model. Furthermore, such a systematic review may help in improving models and optimizing the application of RHESSys^[2] [2].

This study aims at reviewing the progress of RHESSys-based research. Firstly, we introduced the basic structure, principles, and development history for the RHESSys model. By a systematic review of relevant literature, the progress of RHESSys-based research was summarized including the calibration approaches, verification methods, uncertainty analysis, and applications. The ultimate objective of this paper is to reconsider the structure, principles, main research topics, and future trend of RHESSys for further support of the improvement and even broader application of RHESSys.

2. The Basic Structure and Development History of RHESSys

Since Band et al. released the initial version of RHESSys in 1993^[1][10] [1,10], RHESSys has become a matured and popular ecohydrological model (Figure 1) over the last 30 years. Initially, RHESSys was designed by explicitly coupling the Forest Biogeochemical Cycles (FOREST-BGC) canopy model^[11] [11] with a Mountain Climate Simulator (MT-CLIM)^[10] [10], and advanced by coupling with a topography based hydrological model (TOPMODEL)^[12] [12] for the hydrologic process. This version of RHESSys is capable of simulating water, carbon, and the nitrogen cycle in a forest-dominated basin. The Forest-BGC model can simulate vegetation growth, nutrient, and water cycle of the forest ecosystem while MT-CLIM mainly conducts interpolating meteorological variables at a climate station to target points. TOPMODEL is a physical-based quasi-distributed hydrological model. In the coupled RHESSys, the simple soil–water module in the FOREST-BGC was replaced by the vertical infiltration and soil flow process in TOPMODEL [2]. In the updated version of RHESSys, the Forest-BGC was replaced by Biome-BGC to simulate the eco-hydrological processes of multiple ecosystems, while it was coupled with a soil–plant nutrient cycling model (CENTURY_NGAS)^[13] [13] to optimize the simulation process of the nitrogen cycle, especially nitrification and denitrification. For an advance in hydrological processes, RHESSys incorporated an explicit hydrologic routing model (DHSVM) that could account for non-grid-based patches and nonexponential transmissivity profiles. More details on each module and algorithms of RHESSys have been described by model developers^[1][2] [1,2].^{[14][15][16][17][18][19]}

In addition, RHESSys provides a useful tool to simulate the surface processes of watersheds by coupling with other models (Figure 1). For instance, the impacts of fire on the ecohydrological process were evaluated by loosely coupling a fire spread model (WMfire) with RHESSys [3,4]. Moreover, some studies examined the effects of phenology changes on the watershed runoff and evapotranspiration process by coupling a phenology model^[5] [5].

Figure 1. The development history of RHESSys.

RHESSys describes a basin as an object containment hierarchy of Basin, Zone, Hillslope, Patch, and Canopy strata, which allows different hydrological and ecological processes to be modeled at different scales^[20] [14] (Figure 2). The patch represents the smallest unit that has similar soil moisture and land cover. The soil representation is a relatively simple bilayer generalization, i.e., unsaturated and saturated layers. Vertical soil moisture processing and soil biogeochemical cycles are modeled at this level considering snowpack and litter stores. Patches can be derived by multiple layers of land use, soil moisture distribution, or topographic map. Canopy strata describes the vertical process above the ground at the same resolution and partition with the patch, which mainly refers to the physiological processes for plants such as respiration and photosynthesis. In short, the subsurface process is modeled in the patch while the aboveground process is modeled in the canopy strata. The hillslope defines horizontal water movement and redistribution between patches to produce streamflow in a sub-catchment that drains into a stream reach. The hillslope is often derived by GIS-based terrain-partitioning algorithms. The zone defines a region that is usually partitioned by the distribution of climate stations or elevation bands. The zone contains meteorological variables and uses extrapolation methods to characterize spatial variation in these variables. The basin defines a spatial boundary for a catchment, which generally refers to the entire watershed simulated by the model. The basin typically aggregates the net flux of water, carbon, and nitrogen across the whole study area^[2][20] [2,14].

Figure 2. The hierarchical structure of RHESSys (Edited based on [15,16]).

Table 1 illustrates structures, key processes, and applications for various ecohydrological models. It is prominent that RHESSys has a hierarchical structure (Figure 2) to better reflect the multi-scale characteristics of ecohydrological processes in a watershed while other ecohydrological models have ‘basin-grid’ and 'basin-subbasin-grid' structures. For instance, Gorelick et al.^[16] [17] found that RHESSys can handle mixed and heterogeneous land cover at a fine spatial resolution and is suitable for more detailed ecohydrological modeling in small catchments. However, SWAT (Soil and Water Assessment Tool) is spatially lumped at the subbasin level and applies to model large basins with spatially well-segregated landscapes. Moreover, RHESSys realized the bidirectional coupling of the eco-hydrological processes, depicting not only the effects of soil water processes on the plant physiological processes, but also the impacts of vegetation growth on hydrological processes. In contrast, most earlier ecohydrological models (i.e., SWAT and TOPOG_IRM) simplified the complicated vegetation–water interactions. Morán–Tejeda et al.^[17] [18] compared the performance of RHESSys and SWAT with the same input and application areas. The results suggested that RHESSys was more sensitive to land cover and vegetation change while SWAT produced larger changes under climate change. The major underlying cause was that SWAT uses empirical functions of potential evapotranspiration to calculate evapotranspiration, but RHESSys estimates evapotranspiration in a more process-based way, as a complex representation of canopy transpiration controlled by rooting depth, stomatal conductance, etc. Therefore, RHESSys has the advantage for watershed simulations that focus on land cover or vegetation–water interactions. In addition, RHESSy has a flexible structure to further dynamically couple with other models such as phenology, fire, and land-use models, leading to a wider range of applications to support the water resources management under assorted conditions^[2][3][4][5] [2–5].

Table 1. Characteristics of RHESSys and ecohydrological models used in literature.

3. Future Perspectives of RHESSys

3.1. Key Challenges

The application of RHESSys has been hindered by the complexity and availability of parameters, and the requirement of detailed data. As a physical process-based model, RHESSys simulates the ecohydrological processes at the expense of involving substantial parameters and data support, which impedes the application of RHESSys and induces more uncertainties^[23][24] [43,71]. Son et al. ^[25][31] found difficulties in determining snow and soil parameters for dissolved carbon simulations, which may cause a bias in model outputs. Martin et al.^[6] [6] pointed out the lack of detailed rainfall intensity and urban drainage data, which may induce underestimating the peak streamflow. Although empirical parameters have been provided for some biomes^[2][26] [2,69], the model users often need to modify the ecological parameters to cope with localized vegetations^[27] [76].

The calibration approaches incorporated into RHESSys also need to be improved. Currently, RHESSys employs the Monte Carlo method that optimizes the parameters by randomly sampling paired parameter groups and picking out a group with the best performance. This approach usually requires tremendous computational resources. Furthermore, the automatic optimization is available only for soil-related and water quality-related parameters in RHESSys. Although Reyes et al.^[28] [77] applied a Latin super-square sampling method to optimize carbon-allocation parameters, most studies have calibrated vegetation-related parameters manually. Thus, the automated and more systematic and efficient calibration methods for vegetation-related parameters are needed^[29][30] [38,66].

As RHESSys has initially been designed for small-scale basins^[2][9] [2,9], application for a large-scale basin is a challenge. Since large watersheds have strong spatial heterogeneity in the ecosystem and usually include data-sparse areas, it is often challenging to obtain sufficient data and parameters to be used for RHESSys, resulting in high uncertainties. Besides, the applications for a large-scale basin require high-performance hardware to perform huge computational tasks. Therefore, the current applications often simplify vegetation types and degrade the spatial resolution to reduce the complexity of modeling for large basins^[9][26] [9,26].

RHESSys has been applied for a variety of fields such as climate change and land management. However, the current studies have mainly focused on natural systems in watersheds while few studies have taken into account socio-economic systems. Few studies have paid attention to linking the outputs of the carbon, nitrogen, and water simulations from RHESSys to the regional water supply–demand balance, ecosystem services, or human well-being. RHESSys has a lot of potentials to quantify ecosystem services more accurately and to be a useful tool for studies on watershed sustainability. However, there are few relevant studies in the literature^[25][31] [25,78].

3.2. Future Directions

RHESSys needs to enhance the abilities of data collection for improving the simulation capabilities of regional water resources and land management [79]. A single data source may often result in the overfitting phenomenon of parameters, leading to unreliable simulations and predictions [73]. In this sense, multi-source data composition and data assimilation methods can be an alternative to be employed into RHESSys. For instance, Sakas et al.^[32] [37] effectively improved the calibration efficiency by incorporating multi-sources data such as remote-sensing products, ground observations, and field measurements. Hanan et al. [34] also incorporated remote sensing data to set a goal state in the spin-up process, which enhanced the reliability and accuracy of model outputs. Moreover, a more comprehensive parameter library can be built from various applications in parameter localization and advanced remote-sensing technologies.

A number of calibration methods have been developed and applied for hydrologic models such as simulated annealing (SA)^[33] [80], genetic algorithm (GA)^[34] [81] and shuffled complex evolution method (SCE-UA)^[35] [82]. Therefore, RHESSys needs to incorporate the most suitable calibration methods in the future to improve calibration efficiency.

In addition, RHESSys is necessary to be further adapted to large-scale basins. Over the last years, ecohydrological simulations at a large-scale have received attention more and more as climate and land-use change have intertwined with ecohydrology^[36] [83]. As the current version of RHESSys may not be suitable for simulating a large-scale basin, the model structure and some mechanisms need to be modified to adapt to a large-scale basin. Moreover, a parallel computation module is also very useful to reduce computational burdens.

As human activities have intensively increased in watersheds, previous RHESSys-related studies have investigated the impacts of human activities on the ecohydrological process, such as urban planning, agricultural irrigation, soil and water conservation, and reservoir construction^[8][9][37] [8,9,44]. It is necessary to project the impacts of human activities on sustainability in a watershed. Consequently, RHESSys provides a useful tool to quantify watershed ecosystem services and to assess regional sustainability, resulting in promoting sustainable development for watersheds.