D 411 Report On JAMS Modelling at The Scale - SMART Transboundary Groundwater Model
D 411 Report On JAMS Modelling at The Scale - SMART Transboundary Groundwater Model
D 411 Report On JAMS Modelling at The Scale - SMART Transboundary Groundwater Model
29.11.2013
Sponsored by the German Federal Ministry of Education and Research (BMBF):
Prof. Dr. Nico Goldscheider & Prof. Dr. Martin Sauter Dr. Roland Mueller &
Dr. Jochen Klinger Dr. Stefan Geyer
The size of the SMART trans-boundary model defines the working frame of the hydrological
model J2000g. In respect to deliverable D601, considering the main groundwater basin along
the western and eastern side of the Lower Jordan Valley, respectively, set up the extension of
the trans-boundary model. The main groundwater aquifers of the region are on the western side
the Judea Group Aquifers and on the Jordanian side the A7-B2 aquifer.
The mentioned input data for JAMS are spatially organized in hydrological response units
(HRU), which lead to a spatially discriminated output. The HRU mesh builds a raster of defined
scale, and therefore, it is necessary to up- or downscale the input data, e.g. satellite scenes,
aerial photographs and field data (Figure 2). To simplify the linkage between the SMART trans-
boundary groundwater model and the hydrological model JAMS, the used HRU mesh is derived
from the finite element grid of the numerical groundwater flow model (Fig. 2 and 3). This mesh
defines the resolution for all input parameters of J2000g. Because the mentioned input data of
J2000g, particularly the characteristics of the relief, the specific soil and the vegetation
parameters are spatially organised in HRUs. In the present case the hydrological response units
are conform to the elements of the trans-boundary model. The descriptions of the used JAMS
attributes are given in the report to deliverable 401.
Fig. 2. The concept of HRU derivation
On the base of the HRU ground mesh an oversimplification is given on same position after the
intersection but the context of the input parameters is conserved.
Fig. 5: The finished intersection between HRU ground mesh and the land use parameter (left site original
data, right site intersected data).
DATA ACQUISITION FOR THE INPUT PARAMETER OF J2000G
To apply J2000g, climatic information with spatially distributed parameters describing the
catchment are needed. Therefore, installing own tipping buckets completed already existing
precipitation and climate stations in Israel and Jordan. The additional tipping buckets were
installed in the West bank, the Ajloun and within the Fuheis demonstration site in Jordan during
December 2011. The tipping buckets are automatically uploading their data via GPRS to the
webpage of Pessl instruments (fieldclimate) (Fig. 6).
Fig. 6: Network of tipping buckets of the UFZ (blue dots existing rain and climate stations; red dots future
station (installation spring 2012)).
The installation of the tipping buckets on the western side was done in close cooperation
between MEKOROT Co. Ltd., the Hydrological Survey of Israel and the UFZ. The cooperation
additionally opened the way to share data. Consequently, we got a much better climate data
pool for the SMART project area. Due to the option of DASIY (WP2), we were able to use
climatic data, generated by the SMART community (e.g. by Ministry of Water and Irrigation,
Jordan) to close data gaps. As a result 102 precipitation stations were available and were used
to set up the J2000g model. The integration of all existing rain and climate stations in the region
provides very reliable information of distribution of precipitation along the western and eastern
part of the Lower Jordan Valley (Fig. 7). The results of regionalisation processing in J2000g of
punctual data of 102 precipitation stations and the annual precipitation map for the region are
shown in Figure 7. It is observable that the interpolations of the 102 stations in J2000g lead to a
well regionalisation and both maps represent each other. However, areas of strong divergence
do exist, particularly in the southern part of the Lower Jordan Valley. Possible reasons for the
discrepancies are non-continuous precipitation data during the modelled period. The
precipitation map shows a good reproduction of the rainfall gradient, with higher values >700
mm/yr at the elevated mountain areas in both flanks and low values of about 100 mm/yr inside
the Jordan Valley.
Fig. 7: Comparison between annual and regionalized precipitation by J2000g for the model area
Parameterisations of physical soil properties based on the soil map classification of Israel (Dan
et al. 1976; Ravikovitch 1969) and National Soil Map of Jordan (Ministry of Agriculture, Jordan,
1994) (Fig. 8). 14 different soil units (soil types) were classified in the project area. To unify the
information of the different soil maps into one soil map (Fig. 2) was a challenge. The average
soil depth and field capacity of each soil unit was estimated by results of field trips and from
literature values (Bashan et al. 1995, Dan 1980, Dan et al. 1983, Ministry of Agriculture Jordan
1994, Mimi et al. 2009, Singer 2007, Smith et al. 1985). The USDA taxonomy 1975 of National
Soil Map of Jordan hindered the parameterisation of the map. Therefore a diploma student
modified the soil map for the Wadi Zarqa catchment (Jordan) and used a morphological
approach to generate a better soil map (Schulz 2011). During the diploma thesis, the approach
was successfully implemented into J2000g. It is to consider that the approach should general
used on the Jordanian site instead of the National Soil Map of Jordan.
Impervious areas 9 10 -
Water 0 10 -
Sand/Rocks 32 70 -
Parameterisations of soil physical properties based on the soil map classification of Israel (Dan
et al. 1976; Ravikovitch 1969) and National Soil Map of Jordan (Ministry of Agriculture, Jordan,
1994) (Fig. 9). 14 main soil classes and 77 subclasses were digitized (Figure 9). These classes
can be distinguished by grain size, porosities, field capacities and thickness.
The thickness of the different soil types range between 30 and 200 cm. The total field capacities
vary between 42 and 743 mm.
The maximum percolation capacities for the geological units were derived from Berndtsson and
Larson (1987). As stated above, to generate the HRU, all input parameters were spatially
integrated, which resulted in a spatially discriminated mesh of about 92,500 HRU-cells.
Fig. 11: Correlation results of surface runoff for the validation areas Wadi al Malih and Wadi Qilt
The validation on the western side was done according to simulation of surface runoff and
groundwater recharge in Wadi al Malih and Wadi Qilt, Westbank. In both catchments the
observed surface runoff events typically react on the heavy precipitation events by a quick
response and with a high flow rate. These flash floods predominantly consist of direct surface
runoff. Therefore the observed and the simulated surface runoff by J2000g were compared. In
Figure 11 it is visible that the results correlate well, but the best fit results in Wadi Qilt (lower
flash flood picture). The circumstance are linked by the model time resolution in month and
mainly by the precipitation frequency within the different catchments. The northern part of the
western catchment side is characterized by a multitude of precipitation events within a month. In
contrast the southern part is affected by single heavy precipitation events during a short time
(rainfall diagram in Figure 11). In regarding to the available precipitation data (mainly information
in month) the less daily information are summarized per month. During that step a falsification of
the rain values was generated because the hydrological model is not able more to distinguish
between rain duration and frequencies. That case lead to an overestimation of the precipitation
events, which affect also the simulation of the surface runoff and groundwater recharge. In the
upper picture in Figure 11 is detectable, that partly flash floods events are computed while they
are not measured. That phenomenon reflects the upper limit of time-steps for calculation of the
hydrological model. In fact to that results future models should have a higher temporally
resolution.
Additionally the dimension of the simulated groundwater recharge is used for the validation of
the trans-boundary model. The produced results correlates well with the results of groundwater
recharge by CMB and GSD (Rödiger and Siebert, 2009). In the combination of all facts it act to
the assumption that the model reflects the main hydrological process well and can be used for
the quantification of the groundwater recharge. The simulation results are summarized in table 2.
MODELING RESULTS
The model simulation period was set up for a time series from 1977-2010. The presented model
approach of J2000g has the ability to reflect the major hydrological processes of the catchment.
The main simulation results are summarized in table 2. Figure 12 shows the final regionalized
and simulated annual groundwater recharge, the surface runoff, the potential and real (actual)
evapotranspiration derived for the individual HRUs for the modelling period. The recharge plot
illustrates the consistent relation between precipitation and groundwater recharge in the study
area. The groundwater recharge fluctuates similar to precipitation. High values of around 250
mm/yr of groundwater recharge are restricted to the elevated parts of the model area.
Fig. 12: Regionalized simulated groundwater recharge, direct runoff, actual and potential evaporation by
J2000g for the trans boundary catchment.
Tab. 2: Results of water budget simulated by J2000g in depend on the catchment scale
Contrastingly, in the Lower Jordan Valley the groundwater recharge decreases to less than 25
mm/yr. All other simulation results (surface runoff, actual and potential evapotranspiration)
reflect also the connection to the precipitation. General a separation between the western and
eastern mountain areas and the Lower Jordan valley is visible.
Fig. 13: Final simulated precipitation (grey column) and groundwater recharge (black line) by J2000g for
the trans boundary catchment (Year 1977 – 2010)
The final simulated mean groundwater recharge for the entire catchment is shown in Figure 13.
The Figure illustrates the rhythm of precipitation only during the wet and very seldom during the
dry seasons in the period 1977-2010. It is detectable that the years 1979 and 1999 were the
driest years, whereas the years 1986, 1992, and 2003 were wettest. It is shown that the moment
of groundwater recharge (black line) is strongly related to the moment of precipitation. The
outcomes reveal only precipitation amounts exceeding 100 mm/month are generating
groundwater recharge.
The spatial setting of the hydrological model allows us to specify groundwater recharge even for
small locations within the model area. Hence, these simulation results are the best available
input-function to toughen up the SMART trans-boundary model to become transient.
REFERENCES
Acherar, M. and Rambal, S. (1992). Comparative water relations of four Mediterranean oak
species. Plant Ecology 99–100, 177–184.
Baldocchi, D.D., Xu, L. and Kiang, N.Y. (2004). How plant functional-type, weather, seasonal
drought, and soil physical properties alter water and energy fluxes of an oak-grass savanna and
an annual grassland. Agr. Forest Meteorol., 123, 13-39,
Breuer, L., Eckhardt, K. and Frede, H.G. (2003). Plant parameter values for models in temperate
climates. Ecological Modelling 169, 237–293
Dan, J., Yaalon, D.H., Koyumdjisky, H. and Raz, Z. (1976). Soils of Israel (with a 1:500.000
map). Volcani Center, Bet-Dagan, Pamphlet No. 159
Dan, J. & Yaalon, D.H. (1980): Origin and distribution of soils and landscapes in the Northern
Negev. Studies in the Geography of Israel 11, 31-56. Israel Exploration Society, Jerusalem.
(Hebrew, with English summary).
Dan, J. and Dagan, B. (1983). Soil Chronosequences in Israel. Cateno, Vol.10, p. 287-319
Dorman, J. L., and Sellers, P. J. (1989). A global climatology of albedo, roughness length and
stomata] resistance for atmospheric general circulation models as represented by the Simple
Biosphere Model, J. Appl. Meteorol. 28, 833- 855.
Heilmeier, H., Schulze, E.D., Fan, J. and Hartung, W. (2007). General relations of stomatal
responses to xylem sap abscisic acid under stress in the rooting zone – A global perspective.
Flora 202, 624–636
Kergoat, L. (1998). A model for hydrological equilibrium of Leaf Area Index on a global scale.
Journal of Hydrology, 212-213, 268-286
Körner, Ch. (1994). Leaf diffusive conductances in the major vegeta- tion types of the globe. In:
Schulze, E.D., Goldwell, M.M. (Eds), Ecophysiology of Photosynthesis, Ecological Studies, vol.
100. Springer, Berlin, pp. 463–490
Kralisch, S. & Krause, P. (2006). JAMS - A framework for natural resource model development
and application. In: Gourbesville, P., Cunge, J., Guinot, V. and Liong, S-Y., ed., 'Proceedings of
the 7th International Conference on Hydroinformatics', pp. 2356-2363
Krause, P. & Hanisch, S. (2009). Simulation and analysis of the impact of projected climate
change on the spatially distributed waterbalance in Thuringia, Germany, Advances in
Geosciences, 21, 33-48
Krause, P., Biskop, S., Helmschrot, J., Flügel, W.-A., Kang, S. and Gao, T. (2010). Hydrological
System Analysis and Modelling of the Nam Co Basin in Tibet, Advances in Geosciences,
Advances in Geosciences, 2010, 27, 29-36
Lambin, E.F. (1999). Monitoring forest degradation in tropical regions by remote sensing: some
methodological issues. Global Ecology and Biogeography 8, 191–198
Ministry of Agriculture, Jordan (1994). National Soil Map and Land Use Project - The Soils of
Jordan, Level 2: Semi-detailed Studies. Vol. 2: Main Report, Ministry of Agriculture, Jordan
Mimi, Z.A. and Assi, A. (2009). Intrinsic vulnerability, hazard and risk mapping for karst aquifers:
A case study. Journal of Hydrology, 364 (2009) 298–310
Niinemets Ü., Cescatti A., Rodeghiero M. & Tosens T. (2005). Leaf internal diffusion
conductance limits photosynthesis more strongly in older leaves of Mediterranean evergreen
broad- leaved species. Plant, Cell & Environment 28, 1552–1566
Ravikovitch, S. (1969). Manual and map of soils of Israel. Magnes Press, Jerusalem, 96 p.
(English summary)
Rödiger, T., & Siebert, C. (2009). Estimation of natural groundwater replenishment, temporal
and spatial variability (pp. 1–23). SMART-WP4-D402.
Rödiger, T., Geyer, S., Mallast, U., Merz, R., Krause, P., Fischer, C. and Siebert, C. (2013).
Multi-response calibration of a conceptual hydrological model in the semiarid catchment of Wadi
al Arab, Jordan. - Journal of Hydrology (accepted)
Schulz, S. (2011). Groundwater recharge estimation using the hydrological model J2000g in the
Zarqa River Catchment, NW-Jordan. Diploma thesis at UFZ- Helmholtz Centre for
Environmental Research
Schulz, S., Siebert, C., Rödiger, T., Al-Raggad, M. and Merz, R. (2013): Application of the water
balance model J2000 to estimate groundwater recharge in a semi-arid environment - a case
study in the Zarqa River catchment, NW-Jordan. Env. Earth Sci. (accepted).
Schulze, E.-D., Kelliher, F.M., Körner, C., Lloyd, J., Leuning, R. (1994). Relationship between
maximum stomatal conductance, ecosystem surface conductance, carbon assimilation rate and
plant nitrogen nutrition: a global ecology scaling exercise. Ann. Rev. Ecol. Syst. 25, 629–660.
Tilch, N., Uhlenbrook, S., Leibundgut, Ch. (2002). Regionalisierungsverfahren zur Ausweisung
von Hydrotopen in von periglazialem Hangschutt geprägten Gebieten. Grundwasser – Zeitschrift
der Fachsektion Hydrogeologie, 4/2002, 206-216.
USDA, (1975). Soil Taxonomy, Agricultural Handbook No.436, United State, Department of
Agriculture
Wilson, M.F. and Henderson-Sellers, A. (1985). A global archive of land cover and soils data for
use in general circulation climate models. J. Clim., 5, 119-143.
Zeng, X. (2001). Global Vegetation Root Distribution for Land Modeling. Notes and
correspondenc