scholarly journals Application of a distributed physically-based hydrological model to a medium size catchment

2000 ◽  
Vol 4 (1) ◽  
pp. 47-63 ◽  
Author(s):  
L. Feyen ◽  
R. Vázquez ◽  
K. Christiaens ◽  
O. Sels ◽  
J. Feyen

Abstract. Physically based distributed models are rarely calibrated and validated thoroughly because of lack of data. In practice, validation is limited to comparison of simulated and predicted discharges in a catchment, or of simulated and observed piezometric levels in some calibrated wells. Rarely, internal noncalibrated wells or discharge stations are included in model evaluation. In this study, the fully distributed physically based MIKE SHE model was applied to the 600-km2 catchment of the Grote and the Kleine Gete, Belgium. Firstly, the MIKE SHE model was calibrated against both daily discharge measurements and observed water levels and then validated using a simple split-sample test. The observed discharges were simulated successfully in both the calibration and the validation period, while results for the piezometric levels differed considerably among the wells. In addition, a multi-site validation test for 2 internal discharge stations and 6 observation wells showed inferior results for the discharge stations and comparable results for the water table wells. As in the calibration and the split-sample test validation, water table fluctuations were predicted well in some wells, but with little agreement in others. This may be due to scale effects and to the poor quality of the data in certain areas of the catchment. Mainly, the lack of data made it difficult to simulate time series of internal catchment variables with acceptable accuracy so that even the calibrated and validated model could not provide reliable predictions of the water table over the entire catchment. Keywords: integral hydrological modelling; distributed code; MIKE-SHE; model performance; model calibration; model validation

2010 ◽  
Vol 18 (4) ◽  
pp. 30-40 ◽  
Author(s):  
M. Tegelhoffová

Analysis of the development of a hydrological balance for future decades in the Senianska depression in the Eastern Slovak lowlandThe goal of the article was to analyze the hydrological balance for future decades in a pilot area in the Eastern Slovak lowland. The aim was to set up the physically-based Mike SHE hydrological model for the modeling hydrological balance in the selected wetland ecosystem in the Eastern Slovak Lowland. The pilot area - the Senianska depression is located near the village of Senne, between the Laborec and Uh Rivers. Specifically, it is a traditional landscape of meadows, marshes, cultivated soil, small water control structures and forests. To get a complete model set up for simulating elements of the hydrologic balance in the pilot area, it was necessary to devise a model for a larger area, which includes the pilot area - the Senianska depression. Therefore, both the Mike SHE model was set up for the Laborec River basin (a model domain of 500 × 500 m) and the Čierna voda River basin (a model domain of 100 × 100 m), for the simulation period of 1981-2007, is order to get the boundary conditions (overland flow depth, water levels, discharges and groundwater table) for the model of the pilot area. The Mike SHE model constructed for the pilot area - the Senianska depression (a model domain of 1 × 1 m) -was used to simulate the elements of the hydrological balance for the existing conditions during the simulation period of 1983-2007 and for climate scenarios for the simulation period of 1983-2100. The results of the simulated elements of the hydrological balance for the existing conditions were used for a comparison of the evolution of the hydrologic conditions in the past, for identifying wet and flooded areas and for identifying the spatial distribution of the actual evapotranspiration in the pilot area. The built-up model with setting values was used for modeling the hydrological balance in changed conditions - climate change.


2019 ◽  
Vol 23 (5) ◽  
pp. 2245-2260 ◽  
Author(s):  
Mohammad Bizhanimanzar ◽  
Robert Leconte ◽  
Mathieu Nuth

Abstract. We present a new conceptual scheme of the interaction between unsaturated and saturated zones of the MOBIDIC (MOdello Bilancio Idrologico DIstributo e Continuo) hydrological model which is applicable to shallow water table conditions. First, MODFLOW was coupled to MOBIDIC as the physically based alternative to the conceptual groundwater component of the MOBIDIC–MODFLOW. Then, assuming a hydrostatic equilibrium moisture profile in the unsaturated zone, a dynamic specific yield that is dependent on the water table level was added to MOBIDIC–MODFLOW, and calculation of the groundwater recharge in MOBIDIC was revisited using a power-type equation based on the infiltration rate, soil moisture deficit, and a calibration parameter linked to the initial water table depth, soil type, and rainfall intensity. Using the water table fluctuation (WTF) method for a homogeneous soil column, the parameter of the proposed groundwater recharge equation was determined for four soil types, i.e. sand, loamy sand, sandy loam, and loam under a pulse of rain with different intensities. The fidelity of the introduced modifications in MOBIDIC–MODFLOW was assessed by comparison of the simulated water tables against those of MIKE SHE, a physically based integrated hydrological modelling system simulating surface and groundwater flow, in two numerical experiments: a two-dimensional case of a hypothetical watershed in a vertical plane (constant slope) under a 1 cm d−1 uniform rainfall rate and a quasi-real three-dimensional watershed under 1 month of a measured daily rainfall hyetograph. The comparative analysis confirmed that the simplified approach can mimic simple and complex groundwater systems with an acceptable level of accuracy. In addition, the computational efficiency of the proposed approach (MIKE SHE took 180 times longer to solve the three-dimensional case than the MOBIDIC–MODFLOW framework) demonstrates its applicability to real catchment case studies.


2018 ◽  
Author(s):  
Mohammad Bizhanimanzar ◽  
Robert Leconte ◽  
Mathieu Nuth

Abstract. We present a new conceptual scheme of the interaction between unsaturated and saturated zones of the MOBIDIC (Modello Bilancio Idrologico DIstributo e Continuo) hydrological model which is applicable to shallow water table conditions. First, a hydrostatic equilibrium moisture profile was assumed for simulating changes in water table levels. This resulted in a water table based expression of specific yield, which was included in the coupled MOBIDIC-MODFLOW modelling framework for capturing shallow water tables fluctuations. Second, the groundwater recharge was defined using a power type equation based on infiltration rate, soil moisture deficit and a calibration parameter linked to initial water table level, soil type and rainfall intensity. Using the Water Table Fluctuation (WTF) method, the water table rise for a homogeneous soil column under a pulse of rain with different intensities (up to 30 mm/day) the parameter of the proposed groundwater recharge equation was determined for four soil types i.e., sand, loamy sand, sandy loam and loam. The simulated water table levels were compared against those simulated by MIKE-SHE, a physically based integrated hydrological modelling system simulating surface and groundwater flow. Two numerical experiments were carried out: a two-dimensional case of a hypothetical watershed in a vertical plane (constant slope) under a 1 cm/day uniform rainfall rate, and a quasi-real three dimensional watershed under one month of measured daily rainfall hyetograph. The comparative analysis confirmed that the simplified approach can mimic simple and complex groundwater systems with an acceptable level of accuracy. In addition, the computational efficiency of the proposed approach (MIKE-SHE took 180 times longer to solve the 3D case than the MOBIDIC-MODFLOW framework) demonstrates its applicability to real catchment case studies.


2012 ◽  
Vol 9 (10) ◽  
pp. 12061-12102 ◽  
Author(s):  
A. M. Ireson ◽  
A. P. Butler

Abstract. A framework for the rigorous quantification of the timing and magnitude of groundwater recharge is proposed. This involves developing a physically based model for the flow processes in the unsaturated and saturated zones that is consistent with the conceptualisation of the system, and with field observations. Subsequently, the essential behaviour of this model is emulated using a simpler model that can be applied within operational groundwater models. We take a UK Chalk aquifer as a case study. Flow processes are simulated convincingly using a dual permeability, equivalent continuum, Richards' equation model, applied to a 2-D hillslope transect along which four monitoring wells recorded water levels in the unconfined aquifer. A simple conventional recharge model that has been widely used was calibrated to reproduce the water table response simulated by the physically based model. The performance in reproducing the water table was surprisingly good, given the known discrepancies between the actual processes and the model representation. However, comparisons of recharge fluxes simulated by each model highlighted problems with the recharge processes in the simple model. Specifically, artificial bypass flow events during the summer were compensating for recharge that should have come from slow, continual drainage of the unsaturated zone. Such a model may still be useful for assessment of groundwater resources on a monthly basis, under non-extreme climatic conditions. However, under extreme wet or dry conditions, or under a changed climate the predictive capacity of such models is likely to be inadequate.


Water ◽  
2020 ◽  
Vol 12 (2) ◽  
pp. 363 ◽  
Author(s):  
Mohammad Bizhanimanzar ◽  
Robert Leconte ◽  
Mathieu Nuth

This paper presents a comparative analysis of the use of an externally linked (MOBIDIC-MODFLOW) and a physically based (MIKE SHE) surface water-groundwater model to capture the integrated hydrologic responses of the Thomas Brook catchment, in Canada. The main objective of the study is to investigate the effect of simplification in representation of the hydrological processes in MOBIDIC-MODFLOW on its simulation accuracy. To this aim, MOBIDIC and MODFLOW were coupled in order to sequentially exchange the groundwater recharge and baseflow discharges within each computation time step. Using identical sets of hydrogeological properties for the two models, the coefficients of the gravity and capillary reservoirs in MOBIDIC were calibrated so as to closely predict the hydrological budget of the catchment simulated with MIKE SHE. The simulated results show that the two models can closely replicate the observed water table responses at two monitoring wells. However, in very shallow water table locations, the instantaneous response of the water table was not precisely captured in MOBIDIC-MODFLOW. Additionally, the simplified conceptualization of the unsaturated flow in MOBIDIC-MODFLOW resulted in overestimated groundwater recharge during spring and underestimation during summer. Moreover, the computational efficiency of MOBIDIC-MODFLOW, as compared to MIKE SHE, along with less required input data, confirms its potential for regional scale groundwater-surface water interaction modelling applications.


Author(s):  
Sandeep Samantaray ◽  
Abinash Sahoo

Accurate prediction of water table depth over long-term in arid agricultural areas are very much important for maintaining environmental sustainability. Because of intricate and diverse hydrogeological features, boundary conditions, and human activities researchers face enormous difficulties for predicting water table depth. A virtual study on forecast of water table depth using various neural networks is employed in this paper. Hybrid neural network approach like Adaptive Neuro Fuzzy Inference System (ANFIS), Recurrent Neural Network (RNN), Radial Basis Function Neural Network (RBFN) is employed here to appraisal water levels as a function of average temperature, precipitation, humidity, evapotranspiration and infiltration loss data. Coefficient of determination (R2), Root mean square error (RMSE), and Mean square error (MSE) are used to evaluate performance of model development. While ANFIS algorithm is used, Gbell function gives best value of performance for model development. Whole outcomes establish that, ANFIS accomplishes finest as related to RNN and RBFN for predicting water table depth in watershed.


2013 ◽  
Vol 17 (7) ◽  
pp. 2981-3004 ◽  
Author(s):  
E. Vannametee ◽  
D. Karssenberg ◽  
M. R. Hendriks ◽  
M. F. P. Bierkens

Abstract. This paper presents an evaluation of the closure relation for Hortonian runoff, proposed in Vannametee et al. (2012), that incorporates a scaling component to explicitly account for the process heterogeneity and scale effects in runoff generation for the real-world case studies. We applied the closure relation, which was embedded in an event-based lumped rainfall–runoff model, to a 15 km2 catchment in the French Alps. The catchment was disaggregated into a number of landform units, referred to as Geomorphologic Response Units (GRUs), to each of which the closure relation was applied. The scaling component in the closure relation was identified using the empirical relations between rainstorm characteristics, geometry, and local-scale measurable properties of the GRUs. Evaluation of the closure relation performance against the observed discharge shows that the hydrograph and discharge volume were quite satisfactorily simulated even without calibration. Performance of the closure relation can be mainly attributed to the use of scaling component, as it is shown that our closure relation outperforms a benchmark closure relation that lacks this scaling component. The discharge prediction is significantly improved when the closure relation is calibrated against the observed discharge, resulting in local-scale GRU-properties optimal for the predictions. Calibration was done by changing one local-scale observable, i.e. hydraulic conductivity (Ks), using a single pre-factor for the entire catchment. It is shown that the calibrated Ks values are somewhat comparable to the observed Ks values at a local scale in the study catchment. These results suggest that, in the absence of discharge observations, reasonable estimates of catchment-scale runoff responses can possibly be achieved with the observations at the sub-GRU (i.e. plot) scale. Our study provides a platform for the future development of low-dimensional, semi-distributed, physically based discharge models in ungauged catchments.


Author(s):  
Francesco Balduzzi ◽  
Alessandro Bianchini ◽  
Riccardo Maleci ◽  
Giovanni Ferrara ◽  
Lorenzo Ferrari

Darrieus wind turbines are experiencing a renewed interest in the wind energy scenario, in particular whenever small and medium-size installations are considered. In these contexts, the average wind speeds are generally quite low due to scale effects and therefore the most exploited design choices for the turbines are the H-shape configuration, as the entire blade can take advantage of the maximum rotational radius, and high chord to radius ratios, in order to ensure suitable Reynolds numbers on the airfoils. By doing so, the aerodynamic effects induced by the motion of the airfoils in a curved flowpath become more evident and the airfoils themselves have to be designed to compensate these phenomena if conventional design tools based on the BEM theory are used. In this study, fully unsteady 2D simulations were exploited to analyze a three-bladed H-Darrieus wind turbine in order to define the real flow structure and its effects on the turbine performance; in detail, the influence of both the virtual camber and the virtual incidence were investigated. CFD results were supported by experimental data collected on full-scale models reproducing two different airfoil mountings. Finally, the proper design criteria to compensate these phenomena are proposed and their benefits on a conventional simulation with a BEM approach are discussed.


2003 ◽  
Vol 28 ◽  
Author(s):  
Barendra Purkait

The Ganga-Brahmaputra river system together forms one of the largest deltas in the world comprising some 59570 sq km. The waterpower resources of the Brahmaputra have been presumed to be the fourth biggest in the world being 19.83 x 103 m3s1. The entire lower portion of the Brahmaputra consists of a vast network of distributary channels, which are dry in the cold season but are inundated during monsoon. The catchment area of the entire river is about 580,000 sq km, out of which 195,000 sq km lies in India. The maximum discharge as measured at Pandu in 1962 was of the order of 72800 m3 s-1 while the minimum was 1750 m3 s-1 in 1968. The drainage pattern in the valley is of antecedent type while the yazoo drainage pattern is most significant over the composite flood plain to the south of the Brahmaputra. The Brahmaputra valley is covered by Recent alluvium throughout its stretch except a few isolated sedimentary hills in the upper Assam, inselbergs/bornhardt of gneissic hills in the Darrang, Kamrup and Goalpara districts and a few inlying patches of Older Alluvium in the Darrang and Goalpara districts. The basin is very unstable. The present configuration of the basin is the result of uplift and subsidence of the Precambrian crystalline landmasses. Four geotectonic provinces can be delineated in the N-E India through which the Brahmaputra flows. These are bounded by major tectonic lineaments such as the basement E-W trending Dauki fault, a NE-SW trending structural feature of imbricate thrusts known as 'belt of Schuppen' and the NW-SE trending Mishmi thrust. Hydrogeologically, the Brahmaputra basin can be divided into two distinct categories, viz(a) dissected alluvial plain and (b) the inselberg zone. The first category is rep resented in the flood plain extending from the south of Sub-Himalayan piedmont fan zone in the north to right upto the main rock promontory of Garo Hills and Shillong Plateau. The inselberg zone is characterized by fractured, jointed and weathered ancient crystalline rocks with interhill narrow valley plains, consisting of thin to occasionally thick piles of assorted sediments. From the subsurface lithological data, two broad groups of aquifers are identified. These are i) shallow water table and ii) deeper water table or confined ones, separated by a system of aquicludes. The shallow aquifer materials, in general, consist of white to greyish white, fine grained micaceous sand and the thickness ranges from 1.2 to 10.3 m. The sand and clay ratio varies from 1: 2.5 to 1:26. The bedrock occurs at depth ranges of 30.4 to 39.5 m. The materials of the deeper aquifers comprise grey to greyish white, fine to medium grained sand. The sand and clay ratio varies from 1:2 to 1:7. The effective size of the aquifer materials varies from 0.125 to 0.062 mm with uniformity co-efficient around 4.00, porosity 38 to 42%, co-efficient of permeability 304 to 390 galls per day/0.3m2. The ground water is mildly alkaline with pH value 6.5 to 8.5, chloride 10 to 40 ppm, bi-carbonate 50 to 350 ppm, iron content ranges from a fraction of a ppm to 50 ppm. Total dissolved solids are low, hardness as CaCo3 50 to 300 ppm, specific conductance at 25 °C 150 to 650 mhos/cm. The yield from shallow aquifers is 1440 litres to 33750 litres/hour and for deeper aquifers ~ 1700 litres/hour at a drawdown of 13.41 m, specific capacity 21 litres/minute. The temperatures of ground water are 23°-25° C during winter, 24°-26° C during pre-monsoon and 27°- 28° C during peak monsoon. The general hydraulic gradient in the north bank is 1:800 whereas in the south bank it is 1: 300-400 The Tertiary sediments yield a range of water from 200 to 300 l.p.m whereas the yield from the Older Alluvium is 500 to 700 1.p.m. The estimated transmissibility and co-efficient of storage is of the order of ~ 800 1.p.m/ m and 8.2 x 10-3 respectively. Depths to water levels range from 5.3 to 10m below land surface (b.l.s). In the Younger or Newer Alluvium, ground water occurs both under water table and confined conditions. Depths to water levels vary from ground level to 10 m b.l.s. Depth to water ranges from 6 m b.l.s. to 2 m above land surface. The yield of the deep tubewells ranges from 2 to 4 kl/minute for a drawdown of 3 m to 6 m. The transmissibility of the aquifers varies from 69 to 1600 l.p.m/m and the storage co-fficient is of the order of 3.52 x 10-2.


2018 ◽  
Author(s):  
Tim Eckhardt ◽  
Christian Knoblauch ◽  
Lars Kutzbach ◽  
Gillian Simpson ◽  
Evgeny Abakumov ◽  
...  

Abstract. Arctic tundra ecosystems are currently facing rates of amplified climate change. This is critical as these ecosystems store significant amounts of carbon in their soils, which can be mineralized to CO2 and CH4 and released to the atmosphere. To understand how the CO2 net ecosystem exchange (NEE) fluxes will react to changing climatic conditions, it is necessary to understand the individual responses of the physiological processes contributing to CO2 NEE. Therefore, this study aimed: (i) to partition NEE fluxes at the soil-plant-atmosphere interface in an arctic tundra ecosystem; and (ii) to identify the main environmental drivers of these fluxes. Hereby, the NEE fluxes were partitioned into gross primary productivity (GPP) and ecosystem respiration (Reco) and further into autotrophic (RA) and heterotrophic respiration (RH). The study examined flux data collected during the growing season in 2015 using closed chamber measurements in a polygonal tundra landscape in the Lena River Delta, northeastern Siberia. The measured fluxes on the microscale (1 m–10 m) were used to model the NEE, GPP, Reco, RH, RA and net ecosystem production (NPP) over the growing season. Here, for the first time, the differing response of in situ measured RA and RH fluxes from permafrost-affected soils to hydrological conditions have been examined. It was shown that low RA fluxes are associated to a high water table, most likely due to the submersion of mosses, while an effect of water table fluctuations on RH fluxes was not observed. Furthermore, this work found the polygonal tundra in the Lena River Delta to be a sink for atmospheric CO2 during the growing season. Spatial heterogeneity was apparent with the net CO2 uptake at a wet, depressed polygon center being more than twice as high as that measured at a drier polygon rim. In addition to higher GPP fluxes, the differences in NEE between the two microsites were caused by lower Reco fluxes at the center compared to the rim. Here, the contrasting hydrological conditions caused the CO2 flux differences between the microsites, where high water levels lad to lower decomposition rates due to anoxic conditions.


Sign in / Sign up

Export Citation Format

Share Document