Run SWAT Model - ArcView Interface - Revision history

Davidbrito at 14:56, 2 November 2011

2011-11-02T14:56:51Z

Davidbrito: /* In Menu >Edit Input > Point Sources Discharges */

2011-11-02T14:46:06Z

‎In Menu >Edit Input > Point Sources Discharges

Davidbrito: /* Run SWAT with common calibrations */

2011-11-02T14:45:48Z

‎Run SWAT with common calibrations

Davidbrito: /* Using a modified SWAT executable to produce outputs in MOHID format */

2011-11-02T14:25:38Z

‎Using a modified SWAT executable to produce outputs in MOHID format

Davidbrito: /* Run SWAT with Point Sources */

2011-11-02T14:18:35Z

‎Run SWAT with Point Sources

Davidbrito: /* Run SWAT with Point Sources */

2011-11-02T14:17:51Z

‎Run SWAT with Point Sources

Davidbrito: /* Run SWAT with common calibrations */

2011-11-02T14:17:27Z

‎Run SWAT with common calibrations

Davidbrito: /* Run SWAT with common calibrations */

2011-11-02T14:17:09Z

‎Run SWAT with common calibrations

Davidbrito: Created page with "==SWAT Project and first simulation with default parameters== ===Create the SWAT Project=== #Run ArcView #If a Window appears "Welcome to ArcView GIS" asking to create a new view..."

2011-11-02T14:08:54Z

Created page with "==SWAT Project and first simulation with default parameters== ===Create the SWAT Project=== #Run ArcView #If a Window appears "Welcome to ArcView GIS" asking to create a new view..."

New page

==SWAT Project and first simulation with default parameters==
===Create the SWAT Project===
#Run ArcView
#If a Window appears "Welcome to ArcView GIS" asking to create a new view, click "Cancel"
#Go to menu "File"->"Extensions"
#Click on the check box "AVSWATX extendable" and '''NEVER click the check box "Make Default"'''
#[[Image:Swat_ArcViewExtension.png|425px|thumb|center|ArcView SWAT extension]]
#Go to Explorer and create a directory for the project (e.g. Aplica\AVSWAT2005\Portugal\Guadiana\BaciaExemplo\Cenario_teste)
#A new Window appears - SWAT interface. Click on the box "New Project" to create a new project.
#Choose the name for the project in the directory you created (e.g. Bacia_crsr)
#Go back to Explorer and to your project directory (the just created by the project) and copy the "data" file from \\Datacenter\Albufeiras_Eutrofizacao\ComoFazer\1.PreparaInputs\3.SWAT\ExemploAplicacao\data
#A window appears "SWAT - ArcView: DataSets". Choose the data directory just created and named "data".

===Create the project options text file===
Copy from \\Datacenter\Albufeiras_Eutrofizacao\ComoFazer\1.PreparaInputs\3.SWAT\ExemploAplicacao\SimulationProperties_xxx.txt the text file where you will write all the options that you make in this interface. This file is not used by the model, it is only a memo and will be very useful when you have a lot of old simulations and have to understand what options have you chosen and difference betwen them.

Place the file in the ..\watershed\text folder that meanwhile was created by SWAT.

===Watershed Delineation===

A window appears "Watershed delineation".

[[Image:Swat_WatershedDelin.png|425px|thumb|center|ArcView SWAT watershed delineation]]

====In Dem Setup====
'''Remind to write down relevant information to your memo file. Follow the example file that you copied previously (project options text file)'''

#Click the folder icon to open the DEM (digital elevation model), choose to open from disk and select the DEM in the folder data\dem.
#If you do not have yet the DEM go to [[Get_NASA_DEM_for_SWAT|Get NASA DEM for SWAT]]
#Click Properties button
#Click Projection button
#Message "Undefined Projection". Click "OK"
#Message "Do you want to modify your projection?". Click "YES"
#Click "Custom Projection" button
#"Projection Properties" window. Select radio buttom "Custom"
#On "Projection" choose "Transverse Mercator"
#On "Spheroid" choose "Internation 1909"
#Write the coordinate settings in the next picture so that the final it should look like this:
#[[Image:Swat_Projection.png|425px|thumb|center|ArcView SWAT projection]]
#Click "OK" twice.
#Back to "Dem Setup" click "Apply" button
#Message "Do you want the edge cells to flow inside?" Choose "No"
#Message "End of processing". Click "OK"

====In Stream Definition====
'''Remind to write down relevant information to your memo file. Follow the example file that you copied previously (project options text file)'''

#Choose the treshold area for the consideration of a stream line. The value depends on the number of subsbasins generated. As you only will know the exact number of subbasins after completing the Watershed Delineation step, as a first value put in the recommended. Then, in the end of the Watershed Delineation step you should have around 50-80 subbasins; if you have less you shoul decrease the treshold area, if you have more, you should increase.
#Click "Apply"
#When the process is finished, streams will show up in the background

====Outlet and Inlet Definition====
'''Remind to write down relevant information to your memo file. Follow the example file that you copied previously (project options text file)'''

Outlets are places where the model will produce outputs (flow, concentrations), so you should have an outlet in the reservoir wall and in the river monitoring points from INAG (to validate model results)

Inlets are inputs from point sources (e.g. WWTP or industrial sources).
If running with point sources see [[SWAT - Add Point Sources Location]] how to add the location of point sources loads.

#If there is not already an outlet where you whish to have one, add the outlets leaving the radio button on "Outlet" and clicking on "Add". Choose where to add the point. Hint: Use the reservoir shape to know where to locate the downstream watershed outlet (at the reservoir wall so that all drainage area is considered).
#After adding all needed outlets, select the outlet which is the watershed outlet (downstream outlet), clicking first on "Select button" and then pointing at the outlet (or selecting an area where that is the only outlet inside).
#Click "yes" to confirm that you selected one outlet.
#The model generates the watershed and subbasins, then click "Apply" so the model computes the parameters for each subbasin generated.
#Message window saying that the computation was completed. Click "OK"

===Land Use and Soil Definition===
#Add to the ArcView the shapes with the soil data (in data/soil) and land use (in data/uso_solo)
#Go to Menu "AVSWAT" -> "Land Use and Soil Definition". Note that the previous step (watershed delineation) is in this menu and that the next step is unavailable yet.

====In "Land Use Data Layer"====
#Click the folder icon to load the land use shape.
#Choose option: Select shape from watershed view
#Choose option: Format: shape
#Select shape: Corine2000
#Select the field with the land use code: code_00
#Message saying that shape was sucessfully clipped and to reclassify. Click "OK"
#Now we need to translate the land use codes to a land use that SWAT understands. That is done by a table with correspondences. Click "Lookup Table Grid..." button.
#Select: "User Table"
#Select: dbf file
#Browse for the file lucnew_wwater.dbf in data/uso_solo
#Now land uses in swat language appear. Check that all the values in the table have an "Land Use SWAT" code.
#Click Reclassify

====In "Soil Data Layer"====

First you need to copy from SAMSUNG the file "usersoil.dbf" from "C:\AVSWATX\AVSWATDB" to your local folder "C:\AVSWATX\AVSWATDB". This file has the soil names (and hydraulic parameters associated) that you will have in the soil map.

#Load the soil map to the view: Click on the main menu icon "load theme" then select the folder and file ''data/solo/solo_eu_portugal.shp''. This will load the theme in the theme view in ArcView.
#Then, in the "Avswatx"--> "Land Use and Soil Definition" --> "Soil Data Layer" panel, click the folder icon to load the soil shape.
#Choose option: Select shape from watershed view
#Choose option: Format: shape
#Select shape: Solo_eu_Portugal
#Select the field with the land use code: sltxcl
#Message saying that shape was sucessfully clipped and to reclassify. Click "OK"
#Now we need to translate the land use codes to a land use that SWAT understands. That is done by a table with correspondences. Click "Lookup Table Grid..." button.
#Select: dbf file
#Browse for the file soil_crsr.dbf in data/solo
#Now soil types in swat language appear. Check that all the values in the table have an "Name" code.
#Click Reclassify

_____

#Click "Overlay" where model takes out landuse and soil infor for each subbasin.
#Message saying that process is complete. Click "OK".
#This is how it should look like the window after all done:
[[Image:Swat_LandUseSoil.png|650px|thumb|center|ArcView SWAT Land Use and Soil definition]]

===Definition of HRU===
HRU (Hidrological response unit) is a uniform area with the same land use and soil type. It is over these that the model will make the hydraulic and WQ computations. One sub-basin may have one HRU (and the sub-basin will inherit the major land use and major soil type of the sub-basin) which can be a coarser aproximation, or can have multiple HRU's. HRU's do not have georeferenced information, they area areas inside a subbasin with one land use and one soil type.

This is a text phase where you will try to see if one HRU per sub-basin is enough to describe your complete watershed land use and soil distibutions; if not increase the number of HRU.

#First try to have only one HRU per sub-basin (select dominant land use and click ok)
#Click "OK" twice to proceed
#Check on menu "Reports"->"Show List" the original land use and soil type area percentages inside your watershed in "SWAT model Land Use and Soil Distribution" text. Compare these values with the ones generated with your HRU definition (in this case one HRU to one sub-basin in "SWAT model Land Use and Soil Distribution (after treshold application)". The difference between percentages of land use and/or soil shoul not be as higher than 2-3%. If higher, you need to repeat the HRU definition process to have more HRU's.
#If needing more HRU, go back to menu "Window"->"Watershed" and then in menu "AVSWAT"->"HRU distribution" select "Multiple Hydrologic Response Units".
#Two ruler appears one for land use and other to soil. Usullay only land use has high variability comparing to sub-basin sizes (see shapes created by the model) so we only need to move the land use ruler. Try 20% and the click "OK"
#Click "OK" twice to proceed
#Do the same checking again comparing percentages from land use and soil maps with your model definition. If still the difference is not inside 2-3% have to get more HRU and repeat process.

'''What means 20% and 10% Land Use tresholds for HRU definition? It means that SWAT will produce an HRU if inside the sub-basin there is one land use area bigger than 20% the total sub-basin area, meaning that ignores land uses that have areas less than that percentage. The same rational for 10%.'''

===Weather and write default parameters===
# You will need to [[Get_Rain_For_SWAT|Get Rain For SWAT]]
# Create a new (new excel file in office 2003 then save to .dbf) file based on ".\data\weather\LocRain.dbf" and leave the same file name. This file points to the precipitation stations, so substitute (or add if more stations) coordinates and station names. Remind to use a short name for station names and it has to be the same names as the dbf's where precipitation were written in the previous step. Station coordinates may be obtained in SNIRH -->"meteorologia"-->"dados base"-->"Caracteristicas das estações".
# Copy from SAMSUNG file "userwgn.dbf" from "C:\AVSWATX\AVSWATDB" to your local folder "C:\AVSWATX\AVSWATDB".
# Edit the file in ".\data\weather\LocWgn.dbf" with the name and coordinates of the station(s) near to your watershed.

This file points to weather stations that will create climate variables (wind, humidity, etc..) so you should check which is closer to your watershed using the weather stations shapefile in \\datacenter\Albufeiras_Eutrofizacao\Produtos\ProjectosGIS\ShapesBase\EM_WGNSWAT_PT. Station name is in this shape file and coordinates in the file "C:\AVSWATX\AVSWATDB\userwgn.dbf" searching for station name. Please remind to use a station with '''no data fails''' in file userwgn.dbf. <webimage>http://content.screencast.com/users/GRiflet/folders/Jing/media/49b8fe23-0cf1-49ef-8a72-2bb2b93fa851/2009-11-03_1521.png</webimage>
# Go to menu "AVSWAT"->"SWAT View" then menu "Input"->"Weather Stations"
# "Rainfall data" --> "Raingages" --> "Locations Table" [browse] and choose file ".\data\weather\LocRain.dbf"
# "Weather simulation data" --> "Custom database" --> "Locations table" [browse] and choose file ".\data\weather\LocWgn.dbf"
# Click "Ok".
# Substitute file ".\watershed\tables\PCP.dbf" all values -99. by adequate values of precipitation (they should not appear as noted in the first step).
# "input"-->"build all default"-->"yes"-->"continue".
# First window "no"
# Second window "no"
# Third window "yes"
# Ok

===Run SWAT in the Interface===
# "Simulation"-->"RunSwat"
# Click on "Bsn" button...<webimage>http://content.screencast.com/users/GRiflet/folders/Jing/media/ce39e5c4-2fd2-4cd3-8a1e-50f0790e9069/2009-11-03_1303.png</webimage>
# Change the "Potential ET Method" to "Pennman-Monteith method".
# "Setup SWAT Run"
# "Run SWAT". <webimage>http://content.screencast.com/users/GRiflet/folders/Jing/media/0dd9d060-b488-4f9b-acbb-7a8d70e3e04b/2009-11-03_1543.png</webimage> <webimage>http://content.screencast.com/users/GRiflet/folders/Jing/media/1193a02d-4b8f-4e06-8e7d-a27ab5bb9d16/2009-11-03_1554.png</webimage>
#Everything should run fine. Most errors should occur from bad formatting and/or bad data in the dbf input data files. The precipitation file and the LocWgn and LocRain files.
# You may not wait until the end of SWAT simulation and close the run window. This is because SWAT run in its original form creates giant and confuse outputs so there has been created a MOHID version of SWAT. This version runs the original code of SWAT bur produces output in MOHID language (e.g. timeseries that can be read in Timeserie Editor or be used as discharges in MOHID). See next step to know what to do to have your outputs in MOHID way.

==Using a modified SWAT executable to produce outputs in MOHID format==
We are not interested in running in the interface (the previous step was just to see that it runs) because the output files are giant so everytime a change is made in the interface it is needed to run in the folder structure explained in the link below. Each new run with changes has to create a new folder structure and run it using the input files from the interface.
To obtain outputs in a format that suits our needs see the steps in [[Mohid_SWAT|Run SWAT with MOHID outputs]].

==Run SWAT with common calibrations==
The commom calibrations are in groundwater, soil, and management files. The first callibration affects flow and the other two affect water quality.

Just need to change some parameters in the SWAT interface, make a new simulation structure and rerun the model in to that structure so that you will have two simulations (the first, default parameters, and the second with parameters more fitted in to our watershed reality).

===In Menu >Edit Input >SubBasins Data===

By default the first HRU is selected (number, land use, soil are presented and selected)

====1) Groundwater ====
Lets change in all HRU's the hydraulic response of soil to feed the river (groundwater or baseflow)
#Click the .gw
#Change GW_DELAY to 3 days
#Change ALPHA_BF to 1 days
#Click OK
#In >LandUses select radio button to >selected and select all the land uses of your project (old down shift and click to multiple selections)
#In >Soils select radio button to selected and select all the soils of your project (old down shift and click to multiple selections)
It should look like something like this:
[[Image:Swat_GWParameter.png|600px|thumb|center|ArcView SWAT GW parameter change]]
#Click OK to apply it to all the HRU's

====2) Soil Erodibility ====
Lets change in all HRU's the soil erodibility. '''Attention that it has to be applied for each soil type at a time!!''' - just follow instructions

Have to find one HRU for each soil type that exists in the watershed. These HRU's will apply its configuration to all other with the same soil type.

Do the next steps for each soil type (and HRU representing that soil type):
#Click the .sol
#Change USLE_K to 0.3 (originally is 0.1)
#Click OK
#In >LandUses select radio button to >selected and select all the land uses of your project (old down shift and click to multiple selections)
#In >Soils do not do anything! It will apply those settings to all land uses with only that soil type.
It should look like something like this:
[[Image:Swat_SoilParameter.png|650px|thumb|center|ArcView SWAT Soil parameter change]]
#Click OK to apply it to all the HRU's with that soil
#Do not forget to select now other HRU with a differente soil type. Do this steps until all soil types are covered.

====3) Autofertilization ====
Lets put the autofertilization in all HRU's. '''Attention that it has to be applied for each soil type and each land use at a time!!''' - just follow instructions

Autofertilization is disconnected by default we will put the parameters for autofertilization thar are in the manual in the agricultural land uses.

Have to find all the agriculture crops (e.g. AGRC, AGRR, ORCD, RICE...) and select one HRU for each land use (only agricultural) and soil type that exists in the watershed. These HRU's will apply its configuration to all other with the same land use and soil type.
Do the next steps for each land use (agricultural) and soil type (and HRU representing that configuration):
#Click the .mgt
#Double click >Auto fetilization init
#Change >fertilizer to 28-10-10 - it is a standard fertilizer (28%N, 10%P, 10%K)
#Change >AUTO_NSTRS to 0.93 - it is the treshold of the plant stress to start applying fertilizer
#Change >AUTO_NAPP to 50 - it is the maximum amount of fertilizer (N) in one application (kg)
#Change >AUTO_NYR to 300 - it is the maximum amount of fertilizer (N) in one year (kg)
#Change >AUTO_EFF to 1.3 - it is the fertilization efficiency
#Change >AFRT_SURFACE to 0.2 - it is the fraction of fertilizer applied in the first layer of soil [[Image:Swat_MGTParameter.png|650px|thumb|center|ArcView SWAT MGT parameter change]]
#Click OK
#In >LandUses do not do anything! It will apply those settings to all HRUS with this land use.
#In >Soils do not do anything! It will apply those settings to all HRUS with that soil type.
It should look like something like this:
[[Image:Swat_MGTParameter2.png|500px|thumb|center|ArcView SWAT MGT parameter change]]
#Click OK to apply it to all the HRU's with that soil and landuse
#Do not forget to select now other HRU with a different soil type and same crop. Do this steps until all soil types are covered for the first crop type.
#Do not forget to select now other HRU with a different land use (e.g. AGRC, AGRR, ORCD, RICE...).
#Run the simulation with SWAT_MOHID as decribed previously in [[Run_SWAT_Model#Using_SWAT_executable_to_produce_MOHID_outputs]]

==Run SWAT with Point Sources==
If running with point sources:
===In Menu >Edit Input > Point Sources Discharges===
#see [[SWAT - Add Point Sources Loads]] how to add point sources loads.
#Run the simulation with SWAT_MOHID as decribed previously in [[Run_SWAT_Model#Using_SWAT_executable_to_produce_MOHID_outputs]]

===SWAT Project with/without point sources===
A simulation with point sources (defined flow and concentrations) can be easily transformed in one simulation with no point sources if flow and concentrations are set to zero.

Remind that any change made have to run the simulation with SWAT_MOHID as decribed previously in [[Run_SWAT_Model#Using_SWAT_executable_to_produce_MOHID_outputs]]

==Run SWAT with Water Quality processes in the river==

For now water quality in the river is disconnected

@@ Line 228: / Line 228: @@
 #see [[SWAT - Add Point Sources Loads]] how to add point sources loads.
-'''Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT|Run SWAT with MOHID outputs]]'''
+''' Final Step:''' Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT|Run SWAT with MOHID outputs]]
 ===SWAT Project with/without point sources===

← Older revision		Revision as of 14:45, 2 November 2011
Line 221:		Line 221:


−	'''Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT\|Run SWAT with MOHID outputs]]~~'''~~	+	''' Final Step:''' Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT\|Run SWAT with MOHID outputs]]

	==Run SWAT with Point Sources==		==Run SWAT with Point Sources==

@@ Line 150: / Line 150: @@
 # You may not wait until the end of SWAT simulation and close the run window. This is because SWAT run in its original form creates giant and confuse outputs so there has been created a MOHID version of SWAT. This version runs the original code of SWAT bur produces output in MOHID language (e.g. timeseries that can be read in Timeserie Editor or be used as discharges in MOHID). See next step to know what to do to have your outputs in MOHID way.
-==Using a modified SWAT executable to produce outputs in MOHID format==
+===Using a modified SWAT executable to produce outputs in MOHID format===
 We are not interested in running in the interface (the previous step was just to see that it runs) because the output files are giant so everytime a change is made in the interface it is needed to run in the folder structure explained in the link below.
 To obtain outputs in a format that suits our needs see the steps in [[Mohid_SWAT|Run SWAT with MOHID outputs]].
@@ Line 164: / Line 164: @@
 Just need to change some parameters in the SWAT interface, make a new simulation structure and rerun the model in to that structure so that you will have two simulations (the first, default parameters, and the second with parameters more fitted in to our watershed reality).
-===In Menu >Edit Input >SubBasins Data===
+===In the Interface, change parameters===
 By default the first HRU is selected (number, land use, soil are presented and selected)
-====1) Groundwater ====
+====Groundwater ====
 Lets change in all HRU's the hydraulic response of soil to feed the river (groundwater or baseflow)
 #Click the .gw
@@ Line 180: / Line 182: @@
 #Click OK to apply it to all the HRU's
-====2) Soil Erodibility ====
+====Soil Erodibility ====
 Lets change in all HRU's the soil erodibility. '''Attention that it has to be applied for each soil type at a time!!''' - just follow instructions
@@ Line 196: / Line 198: @@
 #Do not forget to select now other HRU with a differente soil type. Do this steps until all soil types are covered.
-====3) Autofertilization ====
+====Autofertilization ====
 Lets put the autofertilization in all HRU's. '''Attention that it has to be applied for each soil type and each land use at a time!!''' - just follow instructions
@@ Line 221: / Line 223: @@
-''' Final Step:''' Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT|Run SWAT with MOHID outputs]]
+===Run the modified executable outside the interface===
 ==Run SWAT with Point Sources==
 If running with point sources:
-===In Menu >Edit Input > Point Sources Discharges===
 #see [[SWAT - Add Point Sources Loads]] how to add point sources loads.
-''' Final Step:''' Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT|Run SWAT with MOHID outputs]]
+===Run the modified executable outside the interface===
 ===SWAT Project with/without point sources===
 A simulation with point sources (defined flow and concentrations) can be easily transformed in one simulation with no point sources if flow and concentrations are set to zero.
 ==Run SWAT with Water Quality processes in the river==
 For now water quality in the river is disconnected

@@ Line 214: / Line 214: @@
 #Do not forget to select now other HRU with a different soil type and same crop. Do this steps until all soil types are covered for the first crop type.
 #Do not forget to select now other HRU with a different land use (e.g. AGRC, AGRR, ORCD, RICE...).
-#Run the simulation outside the interface with SWAT_MOHID as decribed previously in [[Mohid_SWAT|Run SWAT with MOHID outputs]]
 ==Run SWAT with Point Sources==