Using evapotranspiration in SWAT-CUP

[Salvi Asefi]

I am working on the impacts of climate change on hydrological balance of Lake Tanganyika in East Africa Great Lakes region. 

I was wondering if anyone could help me with a SWAT-CUP problem that I have been trying to resolve. 

My study region suffers from lack of proper stream discharge and weather data. In the first stage of my study I have used what was available and my calibration result was rather poor so I decided to use evapotranspiration data retrieved from MODIS. 

But in all my attempts I have not been able to get SWAT-CUP to run. 

I have entered information regarding the river discharge data in the Rch section and information regarding ET in HRU section of SWAT-CUP.

 

If I run SWAT-CUP with river discharge alone (Rch section on), SWAT-CUP runs fine and gives me a simulated output. But if I use both Rch and HRU (with the ET information included) it gives me an error saying that it can not open input file SUFI2.out/ET_OUT_16.txt (please see the attaches screen shot of the error). 

I get the same error if I run SWAT-CUP using only ET data (Rch section off, HRU section on).  

I have not been able to fix this problem. So I was wondering if anyone could help me. I have left the swat-cup files under the following folder:  

https://drive.google.com/open?id=1zWYybKuZIkHfJqgBtrzW8bT1bTrbBbc8

I would very much appreciate your help,

Salvi 

[Karim Abbaspour]

well you didn't complete everything. You need to have both variables in the observed.txt file, and in the var_file_name.txt, and also remove REM for the extract_rch file in the extract.bat file so it gets extracted. See attached.

Also your parameter selection is strange. Use EPCO and ESCO either in bsn or hru, not in both. Calibrate CANMX by itself and fix its value and remove from the list. Follow the protocol in the files below and avoid using random parameters to calibrate. 

https://www.sciencedirect.com/science/article/pii/S0022169415001985 

http://www.mdpi.com/2073-4441/10/1/6 

Karim

--
You received this message because you are subscribed to the Google Groups "SWAT-CUP" group.
To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+u...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Salvi Asefi]

Dear Karim, 

Thank you so much indeed. This is working now. 

I will follow your advice in choosing the parameters as well. 

I really appreciate your quick response. 

Best regards,

Salvi 

[Srishti Gaur]

Dear Karim

I am trying to calibrate at HRU (Rch section off, HRU section on), I tried the procedure followed by you but during calibration, new parameter values are not getting updated. Please help me out with this issue.

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

--

[Karim Abbaspour]

Always upload these files with your question:

- par_inf.txt

- extract_rch.def and/or extract_hru.def (which ever being used) 

- observed_rch.txt and/or observed_hru.txt (which ever being used)

- observed.txt

- file.cio

- var_file_rch.txt and/or var_file_hru.txt (which ever being used)

Karim

[Srishti Gaur]

Dear Karim

Here are the required details, please help me out with this issue.

 details.rar

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

new parameter file is not there because you have 10 parameters and only 2 simulations! You should make about 500 simulations.

Also remove r__CANMX.hru    0   20, or calibrate it first by itself and fix its value (set to to the best value).

Karim

[Srishti Gaur]

Dear Karim Sir

I have tried for 500 simulations, but new parameters are not getting generated. For all simulations results are same i.e., parameters are not getting updated.   Should I try for individual HRU, Please suggest me the possible solution?

Thank you

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

I don't know, just upload all the files again to see.

Karim

[Srishti Gaur]

 cali.rar

Dera Karim

here are the files.

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

You are calibrating using Potential ET. This is meaningless to do. PET is the atmospheric demand and has nothing to do with what goes on in the soil. This is why you have the same results 500 times! You should use Actual ET, if you have any, for calibration.

Karim

[Srishti Gaur]

So what do I need to for that?

[Karim Abbaspour]

what?

[Srishti Gaur]

Sorry That was sent due to mistake, I have this actual ET data only. What is variable column number(s) in the swat output file for AET?

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

in output.rch AET is column 12, just count it. 

Also fix this: 

in observed_hru:   ET_OUT_16

in observed.txt         T_OUT_16

They should all be the same.

Karim

[Srishti Gaur]

Do I need to fix, ET_OUT_16 for each HRU?  Actually, before I was using ET_OUT_1 to ET_OUT_142 for 1 to 142 HRUs.

Thank you

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

what do you mean fix ET_OUT_16 for each HRU? You have 29 subbasins and 29 HRUs, or? You have AET for each HRU?

[Srishti Gaur]

I have 29 sub-basins and 142 HRU. Yes, I have data for each HRU. I was asking as you recommand in last mail 

Also fix this: 

in observed_hru:   ET_OUT_16

in observed.txt         T_OUT_16

Why ET_OUT_16??

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

Sender notified by Mailtrack 11/21/18, 10:30:04 PM

[Karim Abbaspour]

I only see 29 HRUs on the output.hru. I must be looking at different files. Pls send me the files again if you need further assistance.

Karim

[Srishti Gaur]

Dear Sir,

Here is the file,

 output.hru

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

ok, you have 142 hrus and AET is column 12.

Karim

[Srishti Gaur]

Dear Karim 

Can you please suggest me some papers to give an idea regarding AET calibration.

Thank you

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Srishti Gaur]

Dear Karim Sir,

Here I am attaching my set-up for AET calibration at HRU scale. I am not getting appropriate values for ranges for set of parameters. I have used MODIS ET data and you could see the results obtained in iteration 34 that patterns for observed and simulated datasets are not matching (see attached figure). Please suggest me: Is it possible to calibrate this AET? Your comments are very valuable to me.

Thank you 

t

 cali.rar

-

Regards:

Srishti Gaur

Research scholar (Land and Water Resources Engineering)

Agricultural and Food Engineering Department

Indian Institute of Technology, Kharagpur

West Bengal-721302

[Karim Abbaspour]

This is one huge file! Please next time delete all swat output files from the main and Backup folders before zipping the project. There should be no output file in the Backup folder. I can just run swat and get these output, which are anyways not needed. 

Also your first and 35th iterations look the same, so what is the point of making so many iterations after a while?! 

Also, you should try to calibrate first for one subbasin before doing a gruesome 142 simultaneous subbasins!

Also, pay attention to value of parameters, ESCO and EPCO cannot be negative or more than 1. Check these parameters first before running them.

Do an initial run, before any calibration, and decide if this model can be calibrated or not. The simulation results say there is something not correct with this setup.

Karim

 

[Rym Hadded]

Dear Srishti Gaur,

please I use the ET for calibrating with the flow_out variable.

there is limited documentation on how to use the both variables.

please could you help me to use both variable?