As part of my research, I am modelling runoff discharge on a rural watershed. I would like to activate Land Use Update with SWAT. I am using ArcGIS 10.3.1 with ArcSWAT, along with SWATCUP 5.1.6 for automatic calibration. I used SWAT LUU tool (see https://saraswat-swat.rcac.purdue.edu/swatluu) to prepare inputs lup.dat and relevant files for 3 differents land use scenerios (2004, 2009 and 2017). My simulation goes from 2004 to 2017, with one year warm up (2004).
I noticed that when I run SWAT-CUP, it seems to find lup.dat, since deleting any of the associated *.dat files containing HRU fractions (HRU_FR) values returns an error. However, during the auto-calibration process, *.hru files are constantly edited, but their HRU_FR (first line) does not change. As such, it seems that land use update does not work.
Are there any recommended guidelines or anything I have missed ? In case any information is needed upon this, I would be glad to provided them. Thanks in advance.
Roland.
Hi Roland,
Do you mean that the SWAT_CUP copys the lup.dat but It does not copy the associated files (HRUs fraction files)?
If this is the case, might be you could try to copy the land use change fraction files to the working folder of SWAT_CUP for it.
About the .hru changing constantly but the HRU fraction does not change. I think this is because every time the parameters were change, all of the files will be rewrite by SWAT_CUP. SWAT_CUP does not need to change the HRU_FR in .hru file because SWAT will change it according to the information in the lup.dat file.
I also created an R script to generate input files (TxtInOut) for SWAT (new HRUs emerings due to land use change will also be considered and you can also interpolate the HRU fraction at any time, for example, between 2004,2009,2017). If you are interested in using this script, please let me know.
Dear Tam V. Nguyen , I have tried your R script. However it returns me some errors. First part, creating hru_all_land_use_reports went fine and produced hru_lu_all.txt report file as expected. However the second part, writing TxInOut folder, I ran into some errors. Line 139 startingwith "if (as.vector(nluc_all[6+j+k])[i,1] > 0.0){" reported the following error:
"Error in if (as.vector(nluc_all[6 + j + k])[i, 1] > 0) { : missing value where TRUE/FALSE needed"
It looks like the condition being evaluated falls to NAN. I did not know how to fix it. Thing is, nluc_all table in my case has 6 entries, while 6+j+k seems to evaluate to 7, hence creating the error (I believe).
Would you help me fix that ?
Thanks always for your awesome script.
HI Rohland, could you send me your hru_lu_all.txt file and all of the files in "TxtInOut" folder for all land use scenarios. I could run again and let you know. I will delete your data when it is done (so you don't worry about data privacy)
Dear Tam V. Nguyen
I have just sent you an email containing a link to the relevant files uploaded to my Google Drive.
Thanks for your time.
Hi Rohland, I have checked it, I have updated the code in GitHub and it works now.
- Badges
- Science topic
- Similar topics
- Geoscience
- Cartography
- Mapping