The WBM Model Family

The Water Balance Model (WBM) and its descendants predict spatially and temporally varying components of the hydrological cycle and multi constituent water quality variables. Capabilities of the current version of the model include predictions of river discharge, water temperature, dissolved nitrogen, the impacts of irrigation, and the distortion of hydrographs through the operation of reservoirs. Other important features include freeze-thaw dynamics, snowmelt runoff, surface runoff due to impervious surfaces, and a series of physically-based evapotranspiration functions dependent on air temperature, vapor pressure, wind speed, and solar radiation. The WBM model family has been applied for a number of research projects at scales varying from local watersheds (Wollheim et al. 2008b, Stewart 2009) to the global system of rivers. Applications include quantifying regional impacts of irrigation on climate (Douglas et al., 2006), distribution of land use and hydrologic vulnerability (Douglas et al., 2007), magnitude of fluvial sediment trapping by reservoirs (Vorsmarty et al. 2003, Syvitzki et al., 2005), distribution of irrigation water requirements (Wisser et al. 2008) and its impact on the water cycle (Wisser et al. 2010, Fekete et al., 2010), the significance of small reservoirs as a means of rainwater harvesting and storage for irrigation(Wisser et al. 2009); global scale modeling of isotopic signatures in the hydrologic cycle (Fekete et al., 2006), and global estimate of N removal by aquatic systems (Wollheim et al. 2008b).


The original version of WBM (Vörösmarty et al. 1998) was one of the first macroscale hydrological models and forms the basis of all its descendents.

WBMplus Model

WBMplus (Wisser et al. 2010) extends WBM by explicitly accounting for the effects of irrigation and reservoirs on the hydrological cycle. Furhtermore, it implements an improved snow melt routine, a daily time step, and a Muskingum-Cunge flood routing scheme. WBMplus is implemented in the Framework for Aquatic Modeling of the Earth System (FrAMES, Wollheim et al. 2008a). FrAMES is a grid-based, spatially distributed watershed modeling tool that allows application of biogeophysical processes over a wide range of resolutions. The modularity of FrAMES lends itself to staging and testing of new geospatial data sets as well as computer model functions depicting individual hydrologic or biogeochemical processes.

WBMplus Irrigation

WBMplus uses a tiling approach to partition each grid cell into an irrigated and non-irrigated fraction. The irrigated part can be subdivided into any number of parts representing different crops with different physiological properties and hence different water requirements (see Figure). The estimated irrigation water demand for the entire grid cell is met by abstracting water from small reservoirs (if available), the runoff pool, representing groundwater, or river corridor discharge from rivers flowing in the same grid cell. If the demand is higher than those three sources combined, irrigation water is still applied with water taken from a virtual pool, representing abstractions from fossil groundwater that is not connected to the hydrological cycle (see Figure).

The estimated irrigation water demand for the entire grid cell can be met my abstracting water from small reservoirs (if available), the runoff pool, representing groundwater, or river corridor discharge from rivers flowing the same grid cell. If the demand is higher than those three sources combined irrigation water is still applied with water taken from a virtual pool, representing abstactions from fossil groundwater that is not connected to the hydrological cycle (see Figure).

WBMplus Reservoirs

The current version of the model computes water release from large reservoirs (defined at a point along the river network) as a function of inflow to the reservoir, mean annual inflow, current storage, and maximum capacity. The flexible and modular structure of the software framework allows for quick and easy implementation of other reservoir operation policies if more detailed data is available.

WBMplus Regional

The WBMplus Regional model was developed specifically for watershed-scale analyses and incorporates urban land features that are less important at larger spatial scales. This model extends the WBMplus version to combine the previous functionality of water withdrawals, and reservoir routing with hydrological processes associated with impervious surfaces, lawns, and wetland storage. In this update, a fraction of runoff from impervious surfaces is routed directly to receiving waters during each storm event and is based on previous studies of effective impervious surfaces in the basin (Pellerin et al. 2007). Lawns and impervious surfaces decrease the available water capacity in soils based on their spatial coverage and a priori estimates for lawn rooting depths. A wetland detention pool captures and slowly releases surface runoff based on a calibrated parameter. Total runoff from each grid cell is the sum of storm, surface, and groundwater release.


PWBM (Rawlins et al., 2003) was one the first macroscale hydrological models that explictely accounted for the effects of seasonally of permanently frozen soil on components of the hydrological cycle.

HDTM 1.0

The HydroDynamic Model (HDTM 1.0) is the youngest member in the WBM model family. It couples WBMplus with the Geophysical Institute Permafrost Lab permafrost model (GIPL 2.0, Marchenko et al., 2008) to numerically simulate permafrost and active layer parameters and components of the hydrological cycle taking into account the hydraulic properties of frozen soil. To accomplish this, the bucket soil moisture representation in WBM was replaced with a multi layer soil moisture routine that simulates the vertical water transport in the soil column. Hydraulic and thermal properties of the soil layers are estimated from the composition of the soil (fractions of sand, silt, clay, and organic matter), thereby taking into account the thermal properties of peatlands.


Water Temperature Model

A water temperature model has been coupled to the WBM/WTM to predict water temperatures under current and future climate conditions. The model mixes temperatures from surface and groundwater runoff. As water is routed, temperature is exchanged with the atmosphere using the equilibrium temperature approach from Dingman (1974). Water temperatures approach an equilibrium (defined by air temperature and solar radiation) at a rate controlled by river surface area, discharge, and exchange coefficients governed by cloud cover and wind speed. Resulting model outputs allow development of national and global scale maps of water temperature under various climate scenarios.

Nitrogen Removal Model

A nitrogen (N) module has been combined with the WBM/WTM to predict N exports from watersheds, in-channel denitrification, and resulting N concentrations throughout the river network (Wollheim et al. 2008a). The model accounts for both point and nonpoint sources of N inputs to the river network. Nonpoint sources are a function of land use type and runoff conditions. In-channel denitrification is a function of discharge, channel characteristics, and biological uptake parameters.

Key publications (WBMplus):

  1. Wisser, D., S. Frolking, E. M. Douglas, B. M. Fekete, C. J. Vörösmarty, and A. H. Schumann (2008), Global irrigation water demand: Variability and uncertainties arising from agricultural and climate data sets, Geophysical Research Letters, 35, L24408, doi:10.1029/2008GL035296. link to article
  2. Wisser, D., S. Frolking, M. E. Douglas, B. M. Fekete, A. H. Schumann and C. J. Vörösmarty (2009), Blue and green water: The significance of local water resources captured in small reservoirs for crop production- a global scale analysis, Journal of Hydrology, in press, doi: 10.1016/j.jhydrol.2009.07.032. link to article
  3. Wisser, D., B. M. Fekete, C. J. Vörösmarty, and A. H. Schumann (2010): Reconstructing 20th century global hydrography: a contribution to the Global Terrestrial Network- Hydrology (GTN-H), Hydrolology and Earth System Science, 14,1-24. link to pdf
  4. Fekete, B. M, D. Wisser, C. Kroeze, A. Bouwman, Emilio Mayorga, W. Wollheim, C. J. Vörösmarty, Millennium Ecosystem Assessment Scenario drivers (1970-2050): Climate and hydrological alterations, Global Biogeochemical Cycles
  5. Stewart, Robert J. (2009), Separation of River Network Scale Nitrogen Removal Among Main Channel and Two Transient Storage Compartments, MSc. Thesis. University of New Hampshire

Key publications (WBM/WTM):

  1. Vörösmarty, C.J., C.A. Federer and A. Schloss (1998). Potential evaporation functions compared on U.S. watersheds: Implications for global-scale water balance and terrestrial ecosystem modeling. Journal of Hydrology 207: 147-69. link to article
  2. Vörösmarty, C.J., Green, P., Salisbury, J., and Lammers, R. 2000. Global water resources: Vulnerability from climate change and population growth. Science 289 (5477) pp. 284-288.DOI: 10.1126/science.289.5477.284
  3. link to article
  4. Wollheim, W.M., B.J. Peterson, C.J. Vorosmarty, C.S. Hopkinson, and S.A. Thomas. 2008a. Dynamics of N removal over annual time scales in a suburban river network. Journal of Geophysical Research - Biogeosciences. G03038, doi:10.1029/2007JG000660 link to article
  5. Wollheim, W.M., C.J. Vörösmarty, A.F. Bouwman, P. Green, J.A. Harrison, M. Meybeck, B.J. Peterson, S.P. Seitzinger, and J.P. Syvitski 2008. A spatially distributed framework for aquatic modeling of the Earth system (FrAMES). Global Biogeochemical Cycles 22, GB2026, doi:10.1029/2007GB002963.
  6. link to article
  7. Vörösmarty, C. J., M. Meybeck, B. Fekete, K. Sharma, P. A. Green, and J. Syvitski, 2003. Anthropogenic sediment retention: major global impact from registered river impoundments, Global and Planetary Change 39, 169-190.
    link to article
  8. Douglas, E. M., D. Niyogi, S. Frolking, J. B. Yeluripati, R. A. Pielke, Sr., N. Niyogi, C. J. Vörösmarty, and U. C. Mohanty (2006). Changes in moisture and energy fluxes due to agricultural land use and irrigation in the Indian Monsoon Belt. Geophysical Research Letters, 33, L14403.doi:10.1029/2006GL026550. link to article
  9. Vörösmarty, C.J., E.M. Douglas, P.A. Green, and C. Revenga (2005). Geospatial indicators of emerging water stress: An application to Africa. Ambio. 34: 230-236. link to article
  10. Fekete, B. M., J. J. Gibson, P. Aggarwal, and C. J. Vörösmarty, 2006. Application of isotope tracers in continental scale hydrological modeling, Journal of Hydrology 330, 444-456. link to article
  11. Syvitski, J. P. M., C. J. Vörösmarty, A. J. Kettner, and P. Green, 2005. Impact of humans on the flux of terrestrial sediment to the global coastal ocean, Science 308(5720), 376-380. link to article
  12. Douglas,E. M., S. Wood, K. Sebastian and C. J. Vörösmarty. Policy implications of a pan-tropic assessment of the simultaneous hydrological and biodiversity impacts of deforestation. Water Resources Management, 21 (1), 211-232, 2007. link to article

Key publications GIPL:

  1. Marchenko, S., V. Romanovsky, and G. Tipenko, Numerical Modeling of Spatial Permafrost Dynamics in Alaska, In Proceedings of the Ninth International Conference on Permafrost, June 29-July 3, Fairbanks, Alaska, 2008, Vol. 2, pp. 1125-1130, 2008.

Key publications (PWBM):

  1. Rawlins, M. A. , R.B. Lammers, S. Frolking, B. Fekete, C.J. Vörösmarty, 2003. Simulating Pan-Arctic Runoff with a Macro-Scale Terrestrial Water Balance Model. Hydrological Processes, 17, pp. 2521-2539. link to article
  2. Rawlins, M.A., McDonald, K.C., Frolking, S., Lammers, R.B., Fahnestock, M., Kimball, J.S., and Vörösmarty, C.J. 2005. Remote sensing of snow thaw at the pan-Arctic scale using the SeaWinds scatterometer. Journal of Hydrology 312(1-4):294-311. link to article
  3. Rawlins, M., Frolking, S., Lammers, R.B., and Vörösmarty, C.J. 2006. Effects of uncertainty in climate inputs on simulated evapotranspiration and runoff in the Western Arctic. Earth Interactions 10 doi: 10.1175/EI182.1. link to article