Improving the representation of the prairie pothole dynamics in land surface models

2020 ◽  
Author(s):  
Mohamed I. Ahmed ◽  
Amin Elshorbagy ◽  
Alain Pietroniro
Keyword(s):  
Surface Runoff ◽  
Land Surface ◽  
Flow Simulation ◽  
Hydrologic Model ◽  
Surface Model ◽  
Prairie Pothole ◽  
Runoff Generation ◽  
Dem Resolution ◽  
Physically Based ◽  
Streamflow Simulation

<p>The hydrography of the prairie basins is complicated by the existence of numerous land depressions, known as prairie potholes, which can retain a substantial amount of surface runoff. Consequently, the runoff production in the prairies follows a fill, spill, and merging mechanism, which results in a dynamic contributing area that makes the streamflow simulation challenging. Existing approaches to represent the potholes’ dynamics, in different hydrological models, use either a lumped or a series of reservoirs that contribute flow after exceeding a certain storage threshold. These approaches are simplified and do not represent the actual dynamics of the potholes nor their spatial water extents. Consequently, these approaches may not be useful in capturing the potholes’ complexities and may not be able to accurately simulate the complex prairie streamflow. This study advances towards more accurate and physically-based streamflow simulation in the prairies by implanting a physically-based runoff generation algorithm (Prairie Region Inundation MApping, PRIMA model) within the MESH land surface model, and is referred to as MESH-PRIMA. PRIMA is a recently developed hydrological routing model that can simulate the lateral movement of water over prairie landscape using topographic data provided via DEMs. In MESH-PRIMA, MESH handles the vertical water balance calculations, whereas PRIMA routes the water and determines the amount of water storage and surface runoff. The streamflow simulations of MESH-PRIMA (using different DEM resolution as a topographic input) and MESH with its existing conceptual pothole dynamics algorithm are tested on a number of pothole-dominated watersheds within Saskatchewan, Canada, and compared against observed flows. MESH-PRIMA provides improved streamflow and peak flow simulation, compared to that of MESH with its conceptual pothole algorithm, based on the metrics evaluated for the simulations. MESH-PRIMA shows potential for simulating the actual pothole water extents when compared against water areas obtained from remote sensing data. The use of different DEM resolution changes the resulting pothole water extent, especially for the small potholes as they are not detected in the coarse DEM. MESH-PRIMA can be considered as a hydraulic-hydrologic model that can be used for better understanding and accurate representation of the complex prairie hydrology.</p>

scholarly journals Exploring the Sensitivity of Photosynthesis and Stomatal Resistance Parameters in a Land Surface Model

2017 ◽  
Vol 18 (3) ◽  
pp. 897-915 ◽  
Author(s):  
Jennifer L. Jefferson ◽  
Reed M. Maxwell ◽  
Paul G. Constantine
Keyword(s):  
Land Surface ◽  
Stomatal Resistance ◽  
Hydrologic Model ◽  
Surface Model ◽  
Land Surface Models ◽  
Common Land Model ◽  
Surface Models ◽  
Input Parameters ◽  
Rate Of Photosynthesis ◽  
Quantities Of Interest

Abstract Land surface models, like the Common Land Model component of the ParFlow integrated hydrologic model (PF-CLM), are used to estimate transpiration from vegetated surfaces. Transpiration rates quantify how much water moves from the subsurface through the plant and into the atmosphere. This rate is controlled by the stomatal resistance term in land surface models. The Ball–Berry stomatal resistance parameterization relies, in part, on the rate of photosynthesis, and together these equations require the specification of 20 input parameters. Here, the active subspace method is applied to 2100 year-long PF-CLM simulations, forced by atmospheric data from California, Colorado, and Oklahoma, to identify which input parameters are important and how they relate to three quantities of interest: transpiration, stomatal resistance from the sunlit portion of the canopy, and stomatal resistance from the shaded portion. The slope (mp) and intercept (bp) parameters associated with the Ball–Berry parameterization are consistently important for all locations, along with five parameters associated with ribulose bisphosphate carboxylase/oxygenase (RuBisCO)- and light-limited rates of photosynthesis [CO2 Michaelis–Menten constant at 25°C (kc25), maximum ratio of oxygenation to carboxylation (ocr), quantum efficiency at 25°C (qe25), maximum rate of carboxylation at 25°C (vcmx25), and multiplier in the denominator of the equation used to compute the light-limited rate of photosynthesis (wj1)]. The importance of these input parameters, quantified by the active variable weight, and the relationship between the input parameters and quantities of interest vary seasonally and diurnally. Input parameter values influence transpiration rates most during midday, summertime hours when fluxes are large. This research informs model users about which photosynthesis and stomatal resistance parameters should be more carefully selected. Quantifying sensitivities associated with the stomatal resistance term is necessary to better understand transpiration estimates from land surface models.

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

2021 ◽  
Author(s):  
Aristeidis Koutroulis ◽  
Manolis Grillakis ◽  
Camilla Mathison ◽  
Eleanor Burke
Keyword(s):  
Spatial Resolution ◽  
Land Surface ◽  
River Flow ◽  
Land Surface Model ◽  
Flow Simulation ◽  
Surface Model ◽  
Routing Scheme ◽  
Basin Scale ◽  
River Flow Simulation ◽  
River Routing

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

scholarly journals Continental Runoff into the Oceans (1950–2008)

2015 ◽  
Vol 16 (4) ◽  
pp. 1502-1520 ◽  
Author(s):  
Elizabeth A. Clark ◽  
Justin Sheffield ◽  
Michelle T. H. van Vliet ◽  
Bart Nijssen ◽  
Dennis P. Lettenmaier
Keyword(s):  
South America ◽  
Land Surface ◽  
Water Cycle ◽  
Land Surface Model ◽  
Hydrologic Model ◽  
Surface Model ◽  
Infiltration Capacity ◽  
Continental Scale ◽  
Wide Range ◽  
Continental Runoff

Abstract A common term in the continental and oceanic components of the global water cycle is freshwater discharge to the oceans. Many estimates of the annual average global discharge have been made over the past 100 yr with a surprisingly wide range. As more observations have become available and continental-scale land surface model simulations of runoff have improved, these past estimates are cast in a somewhat different light. In this paper, a combination of observations from 839 river gauging stations near the outlets of large river basins is used in combination with simulated runoff fields from two implementations of the Variable Infiltration Capacity land surface model to estimate continental runoff into the world’s oceans from 1950 to 2008. The gauges used account for ~58% of continental areas draining to the ocean worldwide, excluding Greenland and Antarctica. This study estimates that flows to the world’s oceans globally are 44 200 (±2660) km3 yr−1 (9% from Africa, 37% from Eurasia, 30% from South America, 16% from North America, and 8% from Australia–Oceania). These estimates are generally higher than previous estimates, with the largest differences in South America and Australia–Oceania. Given that roughly 42% of ocean-draining continental areas are ungauged, it is not surprising that estimates are sensitive to the land surface and hydrologic model (LSM) used, even with a correction applied to adjust for model bias. The results show that more and better in situ streamflow measurements would be most useful in reducing uncertainties, in particular in the southern tip of South America, the islands of Oceania, and central Africa.

Distributed Hydrologic Modeling in Northwest Mexico Reveals the Links between Runoff Mechanisms and Evapotranspiration

2012 ◽  
Vol 13 (3) ◽  
pp. 785-807 ◽  
Author(s):  
Agustín Robles-Morua ◽  
Enrique R. Vivoni ◽  
Alex S. Mayer
Keyword(s):  
North American ◽  
Land Surface ◽  
Temporal Variations ◽  
Hydrologic Model ◽  
Runoff Generation ◽  
Limited Region ◽  
The North ◽  
Wide Range ◽  
Infiltration Excess ◽  
The Impact

Abstract A distributed hydrologic model is used to evaluate how runoff mechanisms—including infiltration excess (RI), saturation excess (RS), and groundwater exfiltration (RG)—influence the generation of streamflow and evapotranspiration (ET) in a mountainous region under the influence of the North American monsoon (NAM). The study site, the upper Sonora River basin (~9350 km2) in Mexico, is characterized by a wide range of terrain, soil, and ecosystem conditions obtained from best available data sources. Three meteorological scenarios are compared to explore the impact of spatial and temporal variations of meteorological characteristics on land surface processes and to identify the value of North American Land Data Assimilation System (NLDAS) forcing products in the NAM region. The following scenarios are considered for a 1-yr period: 1) a sparse network of ground-based stations, 2) raw forcing products from NLDAS, and 3) NLDAS products adjusted using available station data. These scenarios are discussed in light of spatial distributions of precipitation, streamflow, and runoff mechanisms during annual, seasonal, and monthly periods. This study identified that the mode of runoff generation impacts seasonal relations between ET and soil moisture in the water-limited region. In addition, ET rates at annual and seasonal scales were related to the runoff mechanism proportions, with an increase in ET when RS was dominant and a decrease in ET when RI was more important. The partitioning of runoff mechanisms also helps explain the monthly progression of runoff ratios in these seasonally wet hydrologic systems. Understanding the complex interplay between seasonal responses of runoff mechanisms and evapotranspiration can yield information that is of interest to hydrologists and water managers.

Improvement of the simulation of the water and energy cycle using Multiscale Parameter Regionalization (MPR)

2020 ◽  
Author(s):  
Stephan Thober ◽  
Matthias Kelbling ◽  
Florian Pappenberger ◽  
Christel Prudhomme ◽  
Gianpaolo Balsamo ◽  
...  
Keyword(s):  
Parameter Estimation ◽  
High Resolution ◽  
Parameter Space ◽  
Land Surface ◽  
Hydrologic Model ◽  
Estimation Methods ◽  
Surface Model ◽  
Model Parameters ◽  
Environmental Models ◽  
Energy Cycle

<p>The representation of the water and energy cycle in environmental models is closely linked to the parameter values used in the process parametrizations. The dimension of the parameter space in spatially distributed environmental models corresponds to the number of grid cells multiplied by the number of parameters per grid cell. For large-scale simulations on national and continental scales, the dimensionality of the parameter space is too high for efficient parameter estimation using inverse estimation methods. A regularization of the parameter space is necessary to reduce its dimensionality. The Multiscale Parameter Regionalization (MPR) is one approach to achieve this.</p><p>MPR translates local geophysical properties into model parameters. It consists of two steps: 1) local high-resolution geophysical data sets (e.g. soil maps) are translated into model parameters using a transfer function. 2) the high-resolution model parameters are scaled to the model resolution using suitable upscaling operators (e.g., harmonic mean). The MPR technique was introduced into the mesoscale hydrologic model (mHM, Samaniego et al. 2010, Kumar et al. 2013) and it is key factor for its success on transferring parameters across scales and locations.  </p><p>In this study, we apply MPR to vegetation and soil parameters in the land surface model HTESSEL. This model is the land-surface component of the European Centre for Medium-Range Weather Forecasting seasonal forecasting system. About 100 hard-coded parameters have been extracted to allow for a comprehensive sensitivity analysis and parameter estimation.</p><p>We analyze simulated evaporation and runoff fluxes by HTESSEL using parameters estimated by MPR in comparison to a default HTESSEL setup over Europe. The magnitude of simulated long-term fluxes deviates the most (up to 10% and 20% for evapotranspiration and runoff, respectively) in regions with a large subgrid variability in geophysical attributes (e.g., soil texture). The choice of transfer functions and upscaling operators influences the magnitude of these differences and governs model performance assessed after calibration against observations (e.g. streamflow).</p><p><strong>References:</strong></p><p>Samaniego L., et al.  <strong>https://doi.org/10.1029/2008WR007327</strong></p><p>Kumar, R., et al.  <strong>https://doi.org/10.1029/2012WR012195</strong></p>

scholarly journals Introducing Hysteresis in Snow Depletion Curves to Improve the Water Budget of a Land Surface Model in an Alpine Catchment

2014 ◽  
Vol 15 (2) ◽  
pp. 631-649 ◽  
Author(s):  
Claire Magand ◽  
Agnès Ducharne ◽  
Nicolas Le Moine ◽  
Simon Gascoin
Keyword(s):  
Snow Cover ◽  
Land Surface ◽  
Water Budget ◽  
Land Surface Model ◽  
Surface Model ◽  
Depletion Curve ◽  
Physically Based ◽  
Moderate Resolution Imaging Spectroradiometer ◽  
The Difference ◽  
Snow Model

Abstract The Durance watershed (14 000 km2), located in the French Alps, generates 10% of French hydropower and provides drinking water to 3 million people. The Catchment land surface model (CLSM), a distributed land surface model (LSM) with a multilayer, physically based snow model, has been applied in the upstream part of this watershed, where snowfall accounts for 50% of the precipitation. The CLSM subdivides the upper Durance watershed, where elevations range from 800 to 4000 m within 3580 km2, into elementary catchments with an average area of 500 km2. The authors first show the difference between the dynamics of the accumulation and ablation of the snow cover using Moderate Resolution Imaging Spectroradiometer (MODIS) images and snow-depth measurements. The extent of snow cover increases faster during accumulation than during ablation because melting occurs at preferential locations. This difference corresponds to the presence of a hysteresis in the snow-cover depletion curve of these catchments, and the CLSM was adapted by implementing such a hysteresis in the snow-cover depletion curve of the model. Different simulations were performed to assess the influence of the parameterizations on the water budget and the evolution of the extent of the snow cover. Using six gauging stations, the authors demonstrate that introducing a hysteresis in the snow-cover depletion curve improves melting dynamics. They conclude that their adaptation of the CLSM contributes to a better representation of snowpack dynamics in an LSM that enables mountainous catchments to be modeled for impact studies such as those of climate change.

scholarly journals Assessment of Precipitation Error Propagation in Multi-Model Global Water Resources Reanalysis

2018 ◽  
Author(s):  
Md Abul Ehsan Bhuiyan ◽  
Efthymios I. Nikolopoulos ◽  
Emmanouil N. Anagnostou ◽  
Clement Albergel ◽  
Emanuel Dutra ◽  
...  
Keyword(s):  
Water Resources ◽  
Surface Runoff ◽  
Land Surface ◽  
Runoff Generation ◽  
Surface Exchange ◽  
Hydrologic Modelling ◽  
Study Region ◽  
Precipitation Characteristics ◽  
Model Structures ◽  
Global Water

Abstract. This study focuses on the Iberian Peninsula and investigates the propagation of precipitation uncertainty, and its interaction with hydrologic modelling, in global water resources reanalysis. Analysis is based on ensemble hydrologic simulations for a period spanning 11 years (2000–2010). To simulate the hydrological variables of surface runoff, subsurface runoff, and evapotranspiration, we used four land surface models—JULES (Joint UK Land Environment Simulator), ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems), SURFEX (Surface Externalisée), and HTESSEL (Hydrology-Tiled ECMWF Scheme for Surface Exchange over Land)—and one global hydrological model, WaterGAP3 (Water–Global Assessment and Prognosis). Simulations were carried out for five precipitation products—CMORPH, PERSIANN, 3B42 (V7), ECMWF reanalysis, and a machine learning-based blended product. As reference, we used a ground-based observation-driven precipitation dataset, named SAFRAN, available at 5 km/1  h resolution. We present relative performances of hydrologic variables for the different multi-model/multi-forcing scenarios. Overall, results reveal the complexity of the interaction between precipitation characteristics and different modelling schemes and show that uncertainties in the model simulations are attributed to both uncertainty in precipitation forcing and the model structure. Surface runoff is strongly sensitive to precipitation uncertainty and the degree of sensitivity depends significantly on the runoff generation scheme of each model examined. Evapotranspiration fluxes are comparatively less sensitive for this study region. Finally, our results suggest that there is no single model/forcing combination that can outperform all others consistently for all variables examined and thus reinforce the fact that there are significant benefits in exploring different model structures as part of the overall modelling approaches used for water resources applications.

scholarly journals Identifying uncertainties in simulated streamflow from hydrologic model components for climate change impact assessments

2019 ◽  
Author(s):  
Dongmei Feng ◽  
Edward Beighley
Keyword(s):  
Climate Change ◽  
General Circulation ◽  
Bayesian Model Averaging ◽  
Hydrologic Model ◽  
Process Models ◽  
Runoff Generation ◽  
Modeling Framework ◽  
Total Uncertainty ◽  
Physically Based ◽  
Model Components

Abstract. Assessing the impacts of climate change on hydrologic systems is critical for developing adaptation and mitigation strategies for water resource management, risk control and ecosystem conservation practices. Such assessments are commonly accomplished using outputs from a hydrologic model forced with future precipitation and temperature projections. The algorithms used in the hydrologic model components (e.g., runoff generation) can introduce significant uncertainties in the simulated hydrologic variables, yet the identification and quantification of such uncertainties is rarely studied. Here, a modeling framework is developed that integrates multiple runoff generation algorithms with a routing model and associated parameter optimizations. This framework is able to identify uncertainties from both hydrologic model components and climate forcings as well as associated parameterization. Three fundamentally different runoff generation approaches: runoff coefficient method (RCM, conceptual), variable infiltration capacity (VIC, physically-based, infiltration excess) and simple-TOPMODEL (STP, physically-based, saturation excess), are coupled with Hillslope River Routing model to simulate streamflow. A case study conducted in Santa Barbara County, California, reveals that the median changes are 1–10 % increases in mean annual discharge (Qm) and 10–40 % increases in annual maximum daily discharge (Qp) and 100-yr flood discharge (Q100). The Bayesian Model Averaging analysis indicates that the probability of increase in streamflow can be up to 85 %. However, the simulated discharge uncertainties are large (i.e., 230 % for Qm and 330 % for Qp and Q100) with general circulation models (GCMs) and emission scenarios accounting for more than half of the total uncertainty. Hydrologic process models contribute 10–30 % of the total uncertainty, while uncertainty due to hydrologic model parameterization is almost negligible (

scholarly journals Evaluation of the JULES land surface model in simulating catchment hydrology in Southern Africa

2013 ◽  
Vol 10 (8) ◽  
pp. 11093-11128 ◽  
Author(s):  
N. C. MacKellar ◽  
S. J. Dadson ◽  
M. New ◽  
P. Wolski
Keyword(s):  
Southern Africa ◽  
Land Surface ◽  
River Discharge ◽  
Water Resource Management ◽  
River Flow ◽  
Land Surface Model ◽  
Dry Season ◽  
Surface Model ◽  
Runoff Generation ◽  
Subsurface Runoff

Abstract. Land surface models (LSMs) are advanced tools which can be used to estimate energy, water and biogeochemical exchanges at regional scales. The inclusion of a river flow routing module in an LSM allows for the simulation of river discharge from a catchment and offers an approach to evaluate the response of the system to variations in climate and land-use, which can provide useful information for regional water resource management. This study offers insight into some of the pragmatic considerations of applying an LSM over a regional domain in Southern Africa. The objectives are to identify key parameter sensitivities and investigate differences between two runoff production schemes in physically contrasted catchments. The Joint UK Land Environment Simulator (JULES) LSM was configured for a domain covering Southern Africa at a 0.5° resolution. The model was forced with meteorological input from the WATCH Forcing Data for the period 1981–2001 and sensitivity to various model configurations and parameter settings were tested. Both the PDM and TOPMODEL sub-grid scale runoff generation schemes were tested for parameter sensitivities, with the evaluation focussing on simulated river discharge in sub-catchments of the Orange, Okavango and Zambezi rivers. It was found that three catchments respond differently to the model configurations and there is no single runoff parameterization scheme or parameter values that yield optimal results across all catchments. The PDM scheme performs well in the upper Orange catchment, but poorly in the Okavango and Zambezi, whereas TOPMODEL grossly underestimates discharge in the upper Orange and shows marked improvement over PDM for the Okavango and Zambezi. A major shortcoming of PDM is that it does not realistically represent subsurface runoff in the deep, porous soils typical of the Okavango and Zambezi headwaters. The dry-season discharge in these catchments is therefore not replicated by PDM. TOPMODEL, however, simulates a more realistic seasonal cycle of subsurface runoff and hence improved dry-season flow.

Sign in / Sign up

Export Citation Format

Share Document