A dual objective optimization model was developed to restore habitat connectivity while considering water scarcity losses. The model was applied to Utah's Weber River Basin but the model is generalizable to other watersheds.
Software Requirements: 1) GAMS version 24.2.3 or later (http://gams.com, requires license. Use form on Downloads page to ask GAMS for a temporary license) 2) Excel 3) R studio (https://rstudio.com/, if you wish to combine the four csv files in CombineInput_PATH_UP yourself) Directions to Use: 1)Download all materials to a local folder 2) If desired, use the R script in CombineINPUT_PATH_UP folder to combine the 4 csv files into a single Excel file. Move the combiend Excel file back into the main folder. 3) Alternatively, download the already combined Excel file INPUT_PATH_UP.xlsx from Kraft and Null (2018) - https://www.hydroshare.org/resource/889b9ccbb0c7407ea9a5a1b5d2bbb935/ 4) Download and install GAMS version 24.2.3 from. See Older GAMS Distributions at the bottom of https://www.gams.com/download/. 5) Open the GAMS IDE 6) Select File => Project => New Project and create a new project in the folder where you downloaded materials in Step 1 (a GAMS project is a working directory) 7) Open the file Optimization_Model.gms. 8) Select the Run icon or press F9 9) GAMS will solve ~ 1315*1 = 45 scenarios (about 12 hours or more of run time) the program will spit out a bunch of output .gdx and .xlsx files with the names barriers_revomoved.xlsx, Results.xlsx, etc. that are the output 10) To test everything works on just 1 run, comment out line 248 (a LOOP all statement), and uncomment line 250 (A conditional LOOP statement for one run). 11) Postprocessing - more forecoming
NSF, Division of Chemical, Bioengineering, Environmental, and Transport Systems (CBET)
NSF, Division of Chemical, Bioengineering, Environmental, and Transport Systems (CBET) 1653452
CAREER: Robust aquatic habitat representation for water resources decision-making
Kraft, M., Rosenberg, D.E. and Null, S.E., 2019. Prioritizing stream barrier removal to maximize connected aquatic habitat and minimize water scarcity. JAWRA Journal of the American Water Resources Association, 55(2), pp.382-400.
Explanation of Model Input Parameters:
- Parameter distance(j,i) The distance between the two barriers in km ;
- parameter cost(k) the cost of removing the barriers in $;
- Parameter IICnum_month(months) the precalcualted IICnumerator value for the entire basin for each month (unitless) ;
- parameter A(j, months) the habitat or area associated with each barrier. The quality weighted habitat * the penalty ;
- Parameter links(j,i) The topological distance between the barriers (unitless) ;
- Parameter path_up(j,j,k) The barriers located between two barriers as well as the upward path marked (unitless);
- Parameter economic_costs(k, months) The normalized economic losses (between 0 and 1) (unitless/month);
- Parameter dam_costs(k, months) The actual economic losses($) of the barriers ($/month) ;
- parameter objweights(w) The weights to be applied on the objectives (0-1) (unitless) ;
- Parameter habitat(k, months) the monthly habitat without the penalty (km/month) ;
- Parameter Rem_Budget(scenarios) Maximum budget for barrier removal ($/month) ;
- Parameter R(rad_infl) limit the distance from a single reach which the model will consider connectivity(km) ;
- Parameter area(months) the total quality weighted habitat area by month (km/month);
- Parameter trib(i,j) giving the tributary reaches a bonus (unitless) ;
Additional descriptions of variables are available in the GAMS file.
Earth Sciences | Environmental Sciences | Hydrology | Natural Resource Economics | Natural Resources and Conservation | Natural Resources Management and Policy | Water Resource Management
This work is licensed under a Creative Commons Attribution 4.0 License.
Rosenburg, D.E, & Kraft, M. (2020). Optimizing-Stream-Barrier-Removal. Github. https://github.com/MaggiK/Optimizing-Stream-Barrier-Removal