Personal tools

Mohid SWAT

From MohidWiki

Revision as of 18:36, 7 February 2011 by PedroChambel (talk | contribs)
Jump to: navigation, search

Introduction

Two modified versions of the SWAT model were developed, based on the two SWAT releases (SWAT2000 and SWAT2005). The functionalities described herein are equal in both versions. So we refer to this modified versions as SWAT-MOHID.

Outputs

This modified version mainly produce new outputs of SWAT in MOHID format. The aim of these developments is not to replace the use of SWAT graphical user interfaces for input or output, but to complement it. This also means that the developed outputs do not include all the SWAT variables but instead, they were developed on a need basis either for easier exploring, results analysis, model coupling, etc. This development benefits from many current developments and improvements of MOHID tools made by its users and will continue to benefit in the future.

Time series output

A first modification was the implementation of time series outputs. The modified SWAT source code enables the user to output time series for two geometric entities found in SWAT: (i) hydrological response units (HRU) and (ii) or sub-basins/reaches. The output of the time series is done with a frequency which is independent of the “normal” SWAT output and the format of the time series is done in MOHID format. This allows for example the use of Mohid Statistics Analyser which calculates statistics and obtains graphs of comparison between two time series data. Time series written for selected HRU’s are organized in four categoiries: (i) meteorological information, (ii) plant growth, (iii) nutrient concentration and (iv) erosion. Time series written from selected sub-basins/reaches are nutrient concentration, sediments concentrations and the water flows into the reaches. Time series generated by SWAT-MOHID:

Results per sub-basin

  • RCH – results from the SWAT river – Flows, properties concentrations, temperature, etc
  • DIS – results per sub-basin – contribution of each sub-basin to the reach) – This files were created for coupling SWAT with Mohid River Network. Units are in m3/s and mg/l

Results per HRU

  • PLT – results related with plant growth – LAI, Biomass, Growth Stress, etc
  • MTO – results related with climatic model variables
  • USL – results at HRU level related with erosion
  • NIT – results related with the nitrogen cycle in the soil – plant uptake of N, N remove by leaching, etc
  • PHS – results related with the phosphorus cycle in the soil

Global results

  • LDU – average results related per land use of all watershed
  • GLB – results related with overall watershed
    • GlobalWater.glb - results related with the global water
    • GlobalNutrient.glb – Results related with global sediment budgets
    • Average_Basin_kg_ha.glb – Results related with global model pathways

Animated map results (HDF)

MOHID’s HDF class is a class on the top of the Hierarchical Data Format (HDF library). HDF is a general purpose library and file format for storing scientific data, developed and maintained by the United States National Center for Supercomputing Applications. The main functionality of this library is to store matrix data in a structured way. The inclusion of output result files in HDF5 format is currently under development. This allows the visualization of watershed properties (modeled by SWAT) in animated maps using MOHID animation tools like MOHID GIS (Figure 1). These tools allow the production of animated graphic files showing the spatial and time evolution of the modeled properties. The user decides on the time step output. For each output a vector with the values of soil water content of each sub-basin is stored.

Figure 1 SWAT Soil Water content results and visualization in MOHID GIS using the HDF result files produced by SWAT-MOHID

Code

You can get the source code from MOHID download area in www.mohid.com.

SWAT-MOHID uses the following MOHID modules:

Coupling SWAT and MRN

Why couple SWAT and MRN

Common problems for the application of basin models in the temporary catchments are related to:

  1. Periods without runoff (which results in numerical problems for most models)
  2. Extreme first flush effects with the beginning of the rain period (sub-hourly time steps required in simulations)
  3. Quality of the water (sediments, solutes) is frequently poor described in the models and highly variable in time.

In order to respond to such demands a physical based model, the TempQsim STREAM model, was developed (Galvão et al., 2005). This model is currently maintained by MOHID group under the name MRN.

What is MRN

MRN computes water, sediments and properties transport in a river network. The model is written in FORTRAN 95 and follows an object oriented programming philosophy with a finite volume approach. The different processes occurring in the river are programmed in different modules. This model has been calibrated for the Vène watershed (France) with a special focus on the transport of particulates for the first significant flood events (Obermann, 2007). Fluid flow in this model is governed by conservation equations for mass, momentum, energy and any additional constituents and the numerical algorithm is based on the finite volume approach. Following this strategy it is easier to build conservative transport models and coupling between modules is also simpler because it is based on fluxes. Object oriented programming was also used, which facilitates model coupling.

How does the coupling work

An interface between this model and SWAT was developed in order to simulate agriculture in the catchment using MRN for simulating the river network and the corresponding sediment transport and biochemistry (Chambel-Leitão et al., 2006). SWAT source code was slightly changed so time series of flow/properties are produced for each sub-basin. A “watsub” theme with the location of the outlets of each sub-basin (created by SWAT ESRI ® ArcView extension) is read by MOHID GIS and dynamically construct links between the time series locations of SWAT and discharge nodes of the MRN. MRN runs using surface runoff, lateral flow and groundwater flow from SWAT as input discharges. Property concentrations can be considered constant (user supplied) or can be an output of modified SWAT.

Learn to setup MRN project

You have to create a MNR project similar to the SWAT project in terms of threshold area needed to form a stream. This way you garantee you are using the same river network. How to create a MOHID River Project Step-by-Step

Learn to setup a SWAT MOHID Project

Step 1

Step 2

Go to SWAT Project folder

There go to scenarios\default (i.e. C:\AVSWATX\mulargia1\scenarios\default)

In the default folder create a folder named RES.

In RES folder create a folder for each simulation (name it what better describes it)

For each simulation folder create a folder called Output, a folder called res and copy scenarios\default\txtinout folder into it

Folder structure for simulation

Step 3

Get the latest SWAT-MOHID package (SWAT-MOHID V2) from www.mohid.com (downloads). If you are not registered jet, please do it to download the software. Copy the ascii file named “nomfich.dat” (which is inside the SWAT-MOHID package in the folder ..\SWAT Package\txtinout\) in to the folder ..\scenarios\default\res\[SimulationName]\txtinout and edit the ascii file to 'ROOT_SRT : ..\Output\

Step 4

Copy the executable SWAT.EXE to the ..\scenarios\default\res\[SimulationName]\txtinout folder (you can find this in ..\SWAT Package\txtinout\) Use SWAT2000 or SWAT2005 depending on the version of your SWAT project. Do not run SWAT2005 in a SWAT2000 and vice versa

Step 5

Copy all the ascii files from ..\SWAT Package\res\ into folder ..\scenarios\default\res\[SimulationName]\res. These files include outputs.dat, ts_global.dat, ts_subbasin.dat, etc…

Edit this ascii file for the subbasins and HRU that you want to output. In most of the cases it would be needed an output file for each branch of the network, if you want to couple with Mohid River Network. If you just want to explore SWAT results, just oyt put in Subbasins and HRU you want to analyse.

In file TS_SubBasin.dat the keyword LOCALIZATION_I is the number of the Subbasin to which is desired the output:

In which subbasin (river station) do you want to produce timeserie outputs and which names will be given

The TS_HRU.dat the keyword LOCALIZATION_I is the number of HRU (to get this number go to SWAT file ..\Watershed\text\HRULandUseSoilsReport.txt is the numbers on the right) the keyword LOCALIZATION_J is the number of the Subbasin. This two numbers define the HRU that will return the output:

In which HRU do you want to produce timeserie outputs

Step 6

Run the executable you copied to txtinout

Step 7

Install the latest Mohid package 4.9 available at Mohid Downloads

Step 8

Vizualize output files from folder ..\scenarios\default\res using Mohid Time Series Editor

References

Demissie, G (2010) Integrated Sediment Transport Modeling Using SWAT-MOHID Models in Lake Tana Subbasin, Upper Blue Nile, Ethiopia. MSc UNESCO-IHE. April 2010. [abstract] ([1])

Chambel-Leitão P., F. Braunschweig, L. Fernandes, R. Neves, P. Galvão. (2007) Integration of MOHID model and tools with SWAT model, submitted to the Proceedings of the, 4th International SWAT Conference, July 2-6 2007. ([2])

Obermann, M (2007). Nutrient dynamics in temporary waters of Mediterranean catchments. PhD thesis. University of Hannover. March 2007. ([3])

Chambel-Leitao P, Braunschweig F, Obermann M, Trancoso R, Galvao P, Lo Porto, A (2006). Coupling SWAT and TempQsim Mohid River Network. Geophysical Research Abstracts, EGU Meeting, Wien 25-29 April (Austria). ([4])

See Also