i CRWR Online Report 06-12 WEAP Hydrology Model Applied: The Rio Conchos Basin by Charlotte C. Amato, M.S., P.E. Daene C. McKinney, Ph. D., PE Eusebio Ingol-Blanco, M.S. and Rebecca L. Teasley, M.S. December 2006 CENTER FOR RESEARCH IN WATER RESOURCES Bureau of Engineering Research ? The University of Texas at Austin J.J. Pickle Research Campus ? Austin, TX 78712-4497 This document is available online via World Wide Web at http://www.ce.utexas.edu/centers/crwr/reports/online.html ii Contents Section Page 1 Introduction ........................................................................................................................... 1 1.1 Background .......................................................................................................................... 1 1.2 Objective .............................................................................................................................. 2 2 WEAP Model Preparation ................................................................................................... 3 2.1 Model Structure .............................................................................................................. 3 2.2 Catchment Area .............................................................................................................. 4 2.3 Soil and Land Use Groups .............................................................................................. 8 2.4 River Reaches ............................................................................................................... 14 2.5 Catchment Runoff to River Connectivity ..................................................................... 16 2.6 Deep Water Capacity .................................................................................................... 17 2.7 Deep Conductivity ........................................................................................................ 18 2.8 Initial Z2........................................................................................................................ 19 2.9 Soil Water (Root Zone) Capacity ................................................................................. 19 2.10 Root Zone Conductivity ................................................................................................ 22 2.11 Preferred Flow Direction .............................................................................................. 22 2.12 Initial Z1........................................................................................................................ 23 2.13 Crop Coefficient, Kc ..................................................................................................... 23 2.14 Leaf Area Index ............................................................................................................ 24 2.15 Precipitation .................................................................................................................. 25 2.16 Temperature, Wind and Humidity ................................................................................ 27 2.17 Latitude ......................................................................................................................... 28 2.18 Melting Point, Freezing Point and Initial Snow............................................................ 29 3 HEC-HMS Data .................................................................................................................. 30 4. Calibration Process for 1980 .............................................................................................. 32 4.1. Root Zone Water Capacity............................................................................................ 33 4.2. Root Zone Conductivity ................................................................................................ 33 4.3. Initial Root Zone Water Capacity ................................................................................. 33 4.4. Lower Deep Water Capacity ......................................................................................... 34 4.5. Lower Deep Conductivity ............................................................................................. 34 4.6. Initial Lower Layer Capacity ........................................................................................ 35 5 Results .................................................................................................................................. 36 5.1. Annual Streamflow Comparison ................................................................................... 37 5.2. Monthly Streamflow Comparison................................................................................. 38 6. Recommendations ............................................................................................................... 45 References .................................................................................................................................... 46 Appendix 1. Soil Land Use Intersects by Sub-basin ............................................................... 47 Appendix 2. CRWR Geodatabase Reach Lengths .................................................................. 67 1 1 Introduction 1.1 Background The Rio Grande headwaters in Colorado and flows south through New Mexico to Paso del Norte, the point where New Mexico, Texas and the Mexican state of Chihuahua meet. Paso del Norte is an exceptionally large bi-national metropolitan area with three large and rapidly growing cities: Las Cruces, New Mexico, El Paso, Texas and Ciudad Juarez, Chihuahua. Flowing southeast from Paso del Norte the Rio Grande, known as the Rio Bravo in Mexico, forms the Texas- Mexico border all the way to the Gulf of Mexico. In both the United States and Mexico the Rio Grande/Rio Bravo region is experiencing population growth well above the national average of the respective nations. Due the population growth, this primarily desert river basin is further stressed by an increasing agricultural demand. These factors culminate to make the basin one of the most water stressed basins in the world, with less than 500 cubic meters of water available per capita per year. The bi-national aspect and water stressed status of the basin make it the subject of many studies, each with the goal of improving resource management and sustainable solutions to meet the increasing water supply demands. The Center for Research in Water Resources (CRWR) is a member of the RioGrande/Bravo Physical Assessment Project consortium of U.S. and Mexican universities, governmental and non-governmental agencies. The consortium is working to create a basin wide hydrologic planning model to provide improvement in water resource management. The Water Evaluation and Planning (WEAP) system developed by the Stockholm Environmental Institute is a powerful modeling tool. WEAP is river basin simulation software which includes opportunities for scenarios evaluation as well as water balance and allocation calculations. A WEAP model for the Rio Bravo basin has recently been developed (Danner et al., 2006). The missing component is the hydrologic element of the model that allows for predictions of water flows and availability from precipitation sequences. The objective of this project is to develop a hydrologic model of the Rio Conchos basin using WEAP with the ultimate goal of determining the feasibility and practicality of a WEAP hydrologic model for the entire Rio Bravo basin. The confluence of the Rio Conchos and Rio Bravo is near Ojinaga, Chihuahua and Presidio, Texas and is significant because the Rio Conchos provides approximately two thirds of the total annual water volume to the lower Rio Bravo. The Rio Conchos basin is primarily located in the Mexican state of Chihuahua with the southern edge spilling into Durango. Figure 1 depicts the Rio Conchos basin in beige and the remainder of the Rio Bravo basin in grey. 2 Figure 1. Rio Bravo Basin 1.2 Objective As stated above, the objective of this project is to explore the hydrologic capabilities of WEAP through the development of a Rio Conchos basin model with the intent of evaluating the practicality of incorporating a WEAP hydrology model for the entire Rio Bravo basin into the Physical Assessment Project. In creating the hydrologic model, data from several sources was compiled and pre-processed for use by WEAP. The WEAP model structure, data sources and parameter input techniques employed are discussed in Section 2 of this report. The most significant sources include: ? Mexican Institute of Water Technology (IMTA) (Martinez et al., 2005) ? Soil Water Assessment Tool (World Bank, 2006) The hydrologic capabilities of WEAP are evaluated by comparing the flows simulated by WEAP with the flows simulated by the HEC-HMS model prepared by IMTA (Martinez et al., 2005). Flow results for six locations within the Rio Conchos basin are compared in Section 4. 3 The practical utilization of WEAP as a large, basin scale hydrologic model to assist in water resources planning and management will be evaluated by comparing the general performance of the WEAP model with HEC-HMS model. 2 WEAP Model Preparation 2.1 Model Structure WEAP supports the use of three hydrologic modeling methods: the Rainfall Runoff Method FAO (Food and Agriculture Organization of the United Nations), the Irrigation Demands Only of the FAO, and the Rainfall Runoff Soil Moisture Method. The goal of this project is to create a hydrologic model that can be calibrated in the future, therefore the Rainfall Runoff Soil Moisture Method was chosen because it offers the most comprehensive analysis by allowing for the characterization of land use and/or soil type impacts to hydrological processes (Sieber, 2005). The Rainfall Runoff Soil Moisture Method, or Soil Moisture Method, is a one-dimensional, two soil layer algorithm for calculating evapotranspiration, surface runoff, sub-surface runoff and deep percolation for a defined land area unit. A conceptual diagram of the equations incorporated into the Soil Moisture Method water balance calculations are shown in Figure 2. Figure 2. Soil Moisture Method Model (Source: Sieber, 2005) 4 Using the Soil Moisture Method to more accurately describe the hydrologic response of the basin has the implication that more detailed hydrologic and climatic parameters are required for the model. Consequently, the parameters and data are often difficult to define with certainty. The basic input parameters are listed in Table 1 along with the sensitivities identified for each parameter which are a result of the work of Jantzen et al. (2006). WEAP imposes a model structure in terms input parameter resolution, meaning WEAP forces certain parameters to describe the entire catchment and others to describe smaller land unit areas such as the soil classification or land use category. In the remainder of Section 2 each parameter along with the respective data source and implication of WEAP?s model structure is discussed in greater detail. Table 1. Input Parameters and Sensitivity 2.2 Catchment Area A fundamental parameter of any hydrologic model is the catchment area. The Bi-National Rio Grande/Bravo Geodatabase contains delineated sub-basins for all of the Rio Bravo basins including the Rio Conchos (Pati?o-Gomez and McKinney, 2005). Thirty-eight sub-basins were delineated for the Rio Conchos; these are illustrated in Figure 3. P ar am e t e r U nit s R e s ol ut ion Se ns it ivit y L and U s e A r e a s q km C a t c h m e n t H i gh D e e p W a t e r C a pa c i t y mm C a t c h m e n t H i gh D e e p C o n duc t i v i t y m m /da y C a t c h m e n t M o de r a t e I n i t i a l Z 2 n o un i t C a t c h m e n t N o I n f l ue n c e So i l W a t e r C a pa c i t y mm So i l M o de r a t e R o o t Z o n e C o n duc t i v i t y m m /da y So i l M o de r a t e Pr e f e r e d F l o w D i r e c t i o n n o un i t So i l M o de r a t e I n i t i a l Z 1 n o un i t So i l N o I n f l ue n c e C r o p C o e f f i c i e n t , K c n o un i t L a n d U s e H i gh L e a f A r e a I n de x n o un i t L a n d U s e H i gh C l im at e Pr e c i pi t a t i o n m m /da y C a t c h m e n t H i gh T e m pe r a t ur e C C a t c h m e n t M o de r a t e W i n d m /s C a t c h m e n t L o w H um i di t y % C a t c h m e n t L o w M e l t i n g Po i n t C C a t c h m e n t N o t e v a l ua t e d F r e e z i n g Po i n t C C a t c h m e n t N o t e v a l ua t e d L a t i t ude de gr e e C a t c h m e n t N o t e v a l ua t e d I n i t i a l Sn o w mm C a t c h m e n t N o t e v a l ua t e d 5 Figure 3. Bi-National Rio Grande/Bravo Geodatabase Rio Conchos Sub-basins Using different control points 20 sub-basins were delineated by IMTA for use in their HEC- HMS model (Martinez, 2005). One of the tasks of this project is to compare the WEAP model results to the results of IMTA?s HEC-HMS model; hence, it is desirable to use the IMTA sub- basin configuration. Sub-basin areas shown in the IMTA report, however, do not appear to be delineated consistently with the divides of the river systems observed in the hydro-edge shapefile. Figure 4 shows the Rio Conchos sub-basins as they appear in the IMTA report. Using the WRAPHydro process (Pati?o-Gomez and McKinney, 2005) the sub-basin areas of the Rio Conchos basin were re-delineated using IMTA?s basin outlet control points, yielding the results shown in Figure 5. Significant differences occur only between the La Boquilla (sub-basin 19) and Llanitos (12) sub-basins, highlighted by the red boxes. Table 2 summarizes the area for each of the twenty sub-basins as determined by IMTA and the CRWR. 6 Figure 4. IMTA (Martinez, 2005) Rio Conchos Sub-basins Figure 5. CRWR Rio Conchos Sub-basins 7 Table 2. Comparison of IMTA and CRWR Sub-basins Basin Name IMTA ID IMTA Area CRWR Area sq km sq km Peguis 1 7544.99 7999.30 Sacramento 2 1067.03 1042.61 Las Burras 3 11090.21 11309.47 Luis L. Leon 4 5059.43 5085.51 FCO. Madero 5 1209.10 1211.35 Villalba 6 9327.44 9556.86 Conchos 7 1137.42 1114.39 Jimenez 8 4392.87 4422.96 Chuviscar 9 98.26 106.09 El Rejon 10 154.53 146.85 Chihuahua 11 401.77 399.99 Llanitos 12 1483.88 1829.93 Pico del Aguila 13 658.33 647.61 San Antonio 14 843.80 821.16 San Gabriel 15 267.95 305.85 Puente FFCC 16 1251.12 1270.66 Parral 17 347.95 363.79 Colina 18 244.97 259.06 La Boquilla 19 19054.13 18931.98 Ojinaga 20 1004.28 983.47 Conchos Basin Area 66639.45 67808.88 WEAP does not support a geo-referenced map within the program but will allow shapefiles to be imported as a background map for the WEAP schematic. The sub-basin shapefile shown in Figure 5 was imported into WEAP as a vector layer so that the WEAP sub-basins created could be placed in a manner that is visually consistent with the geography of the Rio Conchos basin. The actual placement of the catchments is arbitrary. The WEAP schematic displays in GCS_WGS_1984 coordinates with a degree bound frame while the CRWR geodatabase is in GCS_North_American_1983 coordinates, NAD_1983_Albers projection (Pati?o-Gomez and McKinney, 2005). The shapefile was reprojected from the CRWR geodatabase to the coordinate system used by WEAP. Figure 6 shows the CRWR sub-basins as they are displayed as a background map layer in WEAP with the catchments added. 8 Figure 6. CRWR Rio Conchos Sub-basins Displayed in WEAP 2.3 Soil and Land Use Groups The twenty sub-basins were sub-divided again by soil groups and land use categories. The land use and soil coverages employed by IMTA (Martinez 2005) are applied in the WEAP model rather than those available from the SWAT model (World Bank, 2006). Again, this relates to the task of evaluating the flows simulated by WEAP by comparing results with IMTA?s HEC-HMS model. Figure 7a and b compare the soil coverages from the SWAT and IMTA data sets, respectively. The SWAT soil classification and land use categories are per the Food and Agricultural Organization (FAO) guidelines. A future model version may find it helpful to work under an internationally referenced system like that provided by the FAO. No work was done to investigate and compare the accuracy or quality of the two datasets. In receiving the soil coverage from IMTA an error occurred in that data for the Ojinaga sub- basin (No. 20) was omitted. For this reason soil areas were approximated from what overlapping data was received for the basin. Figure 7c shows the soil coverage that was received from IMTA in October of 2006. 9 Figure 7a. SWAT Soil Coverage (World Bank 2006) Figure 7b. IMTA Soil Coverage (Martinez, 2005) 10 Figure 7c. IMTA Soil Coverage Applied to WEAP Figure 7. Soil Coverage Datasets To reduce the number of soil categories applied to each sub-basin, with the intent of thereby reducing computation time, four hydrologic soil group classification, e.g. A, B, C, or D, were used instead of the soil series classification. Normally this would decrease the resolution of parameter inputs such as hydraulic conductivity and soil capacities. However, most of the data applied in the WEAP model is extracted from IMTA?s HEC-HMS report, which reports hydraulic parameters calibrated on a sub-basin scale rather than a soil group scale. For this reason the consolidation of soil areas will not affect the results of the model because the limitation is inherent with the resolution of the input parameter data. Figure 8a and b compare the land use coverages from the SWAT and IMTA data sets, respectively. The IMTA shapefile data set was applied to the WEAP model to determine percentages of land use areas per basin. HEC-HMS however, does not utilize parameters such as crop coefficient and leaf area index which apply to land use categories and these are required in WEAP?s Soil Moisture Method. These parameter values were taken from the SWAT dataset and from literature, respectively. The significant difference between the two land use coverages is the number of categories and type of categories defined. Table 3 defines the relationship between the categories as applied to the WEAP model. In cases where multiple SWAT land use categories are shown to correlate to one IMTA land use category, the multiple values from the SWAT data set were averaged and applied to the single IMTA category in WEAP. Values were estimated for categories with no correlation. 11 Figure 8a. SWAT Land Use Coverage (World Bank 2006) Figure 8b. IMTA Soil Coverage Applied to WEAP Figure 8. Land Use Coverage Datasets 12 Table 3. SWAT-IMTA Land Use Category Relationship as Applied in WEAP The area of intersection, or overlap, of each land use category and soil hydrologic group within each sub-basin was determined as a percentage of the total sub-basin area using Arc Toolbox. In Arc Toolbox the Analysis/Intersect tool was used to intersect the land use map shapefile and the soil map shapefile. Then the tool was used again to intersect the land use-soil group intersection with the sub-basin shapefile. This final shapefile intersect allowed the area of all of the land use- soil group categories to be determined as shown below in Table 4 for the Peguis sub-basin. Appendix 1. Soil Land Use Intersects by Sub-basin contains a similar table for each of the twenty sub-basins. 85 U r b a n A r e a s U r b a n A r e a 30 W a t e r B o di e s W a t e r B o di e s I r r i ga t e d A gr i c ul t ur e ( de l t a ) I r r i ga t e d A gr i c ul t ur e ( v a l l e y ) 50 N a t ur a l l y I r r i ga t e d A r e a s Suppl e m e n t a l I r r i ga t i o n 20 F o r r e s t G r a s s e s L o w O pe n F o r r e s t O a k F o r e s t Pi n e F o r r e s t 70 H i gh G r a s s e s a n d Sm a l l B r us h C h a pa r r a l M i c r o ph y l l o us Sc r ub l a n ds Sc r ub l a n d w i t h R o s e t t e d V e ge t a t i o n T h o r n s c r ub l a n d T a m a ul i pa n Sub m o n t a n e Sc r ub l a n d 80 G r a z i n g Pa s t ur e s C ul t i v a t e d G r a s s l a n d 60 Sm a l l Pa s t ur e G r a s s e s N a t ur a l G r a s s l a n d 90 W e t l a n d V e ge t a t i o n - 95 W i t h o ut A ppa r e n t V e ge t a t i o n - I M T A L a n d U s e C o de I M T A L a n d U s e C a t e go r y SW A T L a n d U s e C a t e go r y 75 O t h e r V e ge t a t i o n 40 I r r i ga t e d A r e a s 10 F o r r e s t 13 Table 4. Peguis Sub-basin Soil- Land Use Intersect P e g u i s S u b b a s i n N o . 1 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 10809 0 .0 1 0 .0 2 20 0 0 .0 0 0 .0 0 30 38164 0 .0 4 0 .0 5 40 6380071 6 .3 8 8 .9 1 50 6502437 6 .5 0 9 .0 8 60 20395521 2 0 .4 0 2 8 .4 9 70 27555676 2 7 .5 6 3 8 .4 9 75 5846984 5 .8 5 8 .1 7 80 2563815 2 .5 6 3 .5 8 85 576152 0 .5 8 0 .8 0 90 1721146 1 .7 2 2 .4 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 7 1 .5 9 % o f T o ta l B a si n A r e a 0 .9 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 28325 0 .0 3 0 .0 0 20 1509 0 .0 0 0 .0 0 30 210334 0 .2 1 0 .0 1 40 718386 0 .7 2 0 .0 3 50 5966293 5 .9 7 0 .2 8 60 379643859 3 7 9 .6 4 1 7 .5 1 70 1699345239 1 6 9 9 .3 5 7 8 .3 6 75 8675007 8 .6 8 0 .4 0 80 63771996 6 3 .7 7 2 .9 4 85 83173 0 .0 8 0 .0 0 90 4760341 4 .7 6 0 .2 2 95 5466388 5 .4 7 0 .2 5 T o ta l S o i l A r e a ( sq k m ) 2 1 6 8 .6 7 % o f T o ta l B a si n A r e a 2 8 .1 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 497 0 .0 0 0 .0 0 20 30915 0 .0 3 0 .0 0 30 61411 0 .0 6 0 .0 0 40 1953737 1 .9 5 0 .0 6 50 18742950 1 8 .7 4 0 .5 9 60 1081764821 1 0 8 1 .7 6 3 3 .8 2 70 1677846763 1 6 7 7 .8 5 5 2 .4 5 75 117288454 1 1 7 .2 9 3 .6 7 80 287986809 2 8 7 .9 9 9 .0 0 85 935075 0 .9 4 0 .0 3 90 3920024 3 .9 2 0 .1 2 95 8285730 8 .2 9 0 .2 6 T o ta l S o i l A r e a ( sq k m ) 3 1 9 8 .8 2 % o f T o ta l B a si n A r e a 4 1 .5 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 784935 0 .7 8 0 .0 3 20 0 0 .0 0 0 .0 0 30 779358 0 .7 8 0 .0 3 40 268422 0 .2 7 0 .0 1 50 1145767 1 .1 5 0 .0 5 60 196367414 1 9 6 .3 7 8 .6 8 70 1896450935 1 8 9 6 .4 5 8 3 .8 1 75 43258657 4 3 .2 6 1 .9 1 80 119299233 1 1 9 .3 0 5 .2 7 85 152992 0 .1 5 0 .0 1 90 348263 0 .3 5 0 .0 2 95 3877069 3 .8 8 0 .1 7 T o ta l S o i l A r e a ( sq k m ) 2 2 6 2 .7 3 % o f T o ta l B a si n A r e a 2 9 .3 8 14 2.4 River Reaches River reaches cannot be input directly to WEAP from the geodatabase because the WEAP map serves only as a conceptual schematic. The most significant rivers were added to the WEAP schematic by first exporting and projecting the rivers from the CRWR geodatabase to the WEAP display coordinate system and then adding the projected shapefile to WEAP as a vector layer. The rivers were then drawn in WEAP by tracing over the vector layer. Figure 9a shows the main rivers in the Rio Conchos basin that were imported as a background vector layer and Figure 9b shows the selected river reach segments that were drawn in WEAP. Figure 9a. Geodatabase Main Rivers Figure 9b. WEAP River Reaches Figure 9. Rio Conchos CRWR Bi-National Geodatabase Hydroedges Figure 10 compares the river reaches from the CRWR geodatabase with those drawn in WEAP; the geodatabase reaches are shown in red and the WEAP reaches are in blue. Inaccuracies of the required tracing are evident by the red lines visible adjacent to the blue lines. The river reaches drawn in WEAP can be edited to more precise river reach lengths as measured from the CRWR geodatabase by right clicking on the river reach and selecting Edit Data and then Reach Length. WEAP uses reach length data when modeling the groundwater and surface water interactions. Groundwater was not modeled in this project so the reach lengths were not adjusted in this model. Although the reach length data was not used for this project, the values were determined and appear in Appendix 2 for use in future work. 15 Figure 10. WEAP River Reaches The river reaches were drawn into WEAP in segments relating to the sub-basin boundaries for reasons related to catchment connectivity. Relating reach lengths to the sub-basin boundaries allows for the runoff from a catchment to link to the river system in a physically appropriate location. For example the Rio Conchos was traced as one segment in each sub-basin and named accordingly, e.g. Rio_Conchos_20 represents the Rio Conchos through the Ojinaga sub-basin which is identified as sub-basin number 20. The reach names are also shown in Figure 10. There are two exceptions to the naming convention. The first exception is River_10, which drains the El_Rejon sub-basin, sub-basin 10, through the Chuviscar sub-basin, sub-basin 9, to the confluence with the Arroyo_Sacramento_9. This reach is generically named because it did not have a name field in the CRWR geodatabase. The second exception is the WEAP river reach in the La Boquilla sub-basin, sub-basin 19, which is named the RioBalleza_Conchos_19, because it represents the Rio Balleza to the confluence with the Rio Conchos and the Rio Conchos from the confluence to the sub-basin outlet. The river was traced in such a manner so that the Llanitos sub-basin, sub-basin 12, could be entered as headwater to the river reach segment and so that the La Boquilla sub-basin could be connected to the river in a more physically appropriate location, this is further discussed subsequently. 16 2.5 Catchment Runoff to River Connectivity Once the river system is traced and named the catchment data is edited to direct the runoff from the sub-basin to the appropriate river reach. For example the Rio Florido headwaters are in the Puente FFC sub-basin, sub-basin 16, and flows through the San Gabriel sub-basin, sub-basin 15, the San Antonio sub-basin, sub-basin 14, the Pico de Aquila sub-basin, sub-basin 13, the Jimenez sub-basin, sub-basin 18, and part of the Las Burras sub-basin, sub-basin 3, to the confluence with the Rio Conchos. Figure 11 shows the upper Rio Florido portion of the WEAP schematic. The solid blue lines are rivers and the dashed blue lines are the conceptual connection of the catchment to the river reach. Figure 11. WEAP Schematic of the Rio Florido Connectivity The schematic shows the connectivity either connecting directly to the head of a reach or as a perpendicular line to the closest point on the identified river reach. The actual surface flow contribution, however, connects to the center of the selected river reach. For example there is no river reach in the Puente FFC sub-basin, because runoff from the sub-basin is entered as head flow for the segment through the downstream San Gabriel sub-basin, Rio_Florido_15. Runoff from the San Gabriel sub-basin enters the Rio_Florido_15 at the center of the river segment. The runoff contribution point is dictated by right clicking on the catchment, selecting General Info, and then selecting the river reach from the pull down menu and selecting the head flow option if desired. Figure 12 shows the WEAP window settings for the above example. 17 Figure 12a. Sub-basin 16 Runoff Settings Figure 12b. Sub-basin 15 Runoff Settings Figure 12. Runoff Settings 2.6 Deep Water Capacity Deep Water Capacity is the effective water holding capacity, in millimeters, of the deep soil layer, or the second bucket in the Soil Moisture Method. WEAP applies this parameter to an entire catchment so that the parameter cannot be characterized by land use or soil area. IMTA (Martinez, 2005) determined the storage capacity and percolation rate for each sub-basin for a three-soil layer, or three-bucket model. Figure 13 is a conceptual diagram of the soil moisture method applied to HEC-HMS by IMTA. Figure 13. IMTA Soil Moisture Method (Martinez, 2005) 18 The first bucket in the WEAP Soil Moisture Method relates to the ?Soil Profile? layer in IMTA?s model. ?Ground Water 1? is the second layer in IMTA?s model that relates to the second bucket in the WEAP Soil Moisture Method. Data for the third bucket of IMTA?s model was not utilized because the Soil Moisture Method in WEAP only supports two layers. Table 5 lists the values for Deep Water Capacity found in the IMTA report and those applied in the WEAP model. Table 5. Deep Water Capacity Values 2.7 Deep Conductivity The Deep Conductivity parameter represents the conductivity rate of the second bucket, in millimeters per day. As Figure 2 shows, Deep Conductivity controls the transmission of base flow. WEAP applies a single value of Deep Conductivity to the entire catchment. IMTA determined percolation rates for each sub-basin, which are applied as conductivity values by setting the Preferred Flow Direction Values to 0, which indicates 100% vertical flow as percolation implies. The Preferred Flow Direction is discussed in a subsequent section. Table 6 lists the values IMTA determined for percolation rates and those applied to Deep Conductivity in the WEAP model. I M T A Gr o un d W a t e r 1 W E A P S t o r a ge C a pac i t y ( m m ) De e p W a t e r C a pac i t y ( m m ) P e gui s 1 4 4 S a c r a m e n t o 2 5 5 L a s B ur r a s 3 1 1 L ui s L . L e o n 4 25 25 F C O. M a der o 5 10 10 Vi l l a l b a 6 8 8 C o n c h o s 7 3 3 J i m i n e z 8 3 3 C h uv i s c a r 9 3 3 E l R e j o n 10 10 10 C h i h uah ua 11 5 5 L l a n i t o s 12 1 1 P i c o del A gui l a 13 5 5 S a n A n t o n i o 14 20 20 S a n Ga b r i e l 15 25 25 P uen t e F F C C 16 25 25 P a r r a l 17 8 8 C o l i n a 18 15 15 L a B o qui l l a 19 15 15 Oj i n a ga 20 4 4 I M T A I DB a s i n Na m e 19 Table 6. Deep Conductivity Values I M T A Gr o un d W a t e r 1 W E A P P e r c o l a t i o n M a x R a t e ( m m /h ) De e p C o n duct i v i t y ( m m /d ) P e gui s 1 4 96 S a c r a m e n t o 2 20 480 L a s B ur r a s 3 5 120 L ui s L . L e o n 4 15 360 F C O. M a der o 5 25 600 Vi l l a l b a 6 25 600 C o n c h o s 7 13 312 J i m i n e z 8 5 120 C h uv i s c a r 9 13 312 E l R e j o n 10 20 480 C h i h uah ua 11 20 480 L l a n i t o s 12 8 192 P i c o del A gui l a 13 22 528 S a n A n t o n i o 14 35 840 S a n Ga b r i e l 15 10 240 P uen t e F F C C 16 10 240 P a r r a l 17 30 720 C o l i n a 18 10 240 L a B o qui l l a 19 45 1080 Oj i n a ga 20 4 96 I M T A I DB a s i n Na m e 2.8 Initial Z2 The ?Initial Z2? parameter is the relative storage given as a percentage of the total effective storage of the Deep Water Capacity at the beginning of a simulation. WEAP, like Deep Water Capacity, forces Initial Z2 to be constant for each basin. A value of 50 percent was assigned to every sub-basin. Refinement of these three deep soil layer parameters would require finer catchment delineation, referring to Table 1, calibration was determined to be sensitive to Deep Capacity and Conductivity but not Initial Z2 (Jantzen, 2006). Future work may investigate potential value, if any, of such a course. 2.9 Soil Water (Root Zone) Capacity Soil Water or Root Zone Capacity is the effective water holding capacity, in millimeters, of the first bucket in the Soil Moisture Method. The WEAP model structure allows this parameter to characterize the soils groups within a sub-basin. Typically in WEAP, values of Soil Water Capacity are applied to the land use groups delineated within each sub-basin. However, as previously stated IMTA (Martinez, 2005) determined storage capacity and percolation rate on a sub-basin scale so the value for each sub-basin was entered multiple times for each land use category in the sub-basin. Table 7 lists the storage capacity values determined by IMTA and those applied in WEAP. 20 Table 7. Soil Capacity Values WEAP employs a method known as ?Key Assumption?, which allows parameter values that will be applied frequently to be coded in once as a Key Assumption and then referenced throughout the model. The Soil Water Capacity values were coded using the Key Assumption function to assign the sub-basin value to each land use within each soil group. The Root Zone Capacity Key Assumption setup and application are shown in Figure 14 and Figure 15. I M T A S o i l P r o f i l e W E A P S t o r a ge C a pac i t y ( m m ) R o o t Z o n e C a pac i t y ( m m ) P e gui s 1 20 20 S a c r a m e n t o 2 20 20 L a s B ur r a s 3 10 10 L ui s L . L e o n 4 15 15 F C O. M a der o 5 15 15 Vi l l a l b a 6 9 9 C o n c h o s 7 10 10 J i m i n e z 8 25 25 C h uv i s c a r 9 10 10 E l R e j o n 10 5 5 C h i h uah ua 11 10 10 L l a n i t o s 12 6 6 P i c o del A gui l a 13 10 10 S a n A n t o n i o 14 10 10 S a n Ga b r i e l 15 8 8 P uen t e F F C C 16 8 8 P a r r a l 17 15 15 C o l i n a 18 10 10 L a B o qui l l a 19 5 5 Oj i n a ga 20 20 20 I M T A I DB a s i n Na m e 21 Figure 14. Root Zone Capacity Key Assumptions Setup Figure 15. Root Zone Capacity Key Assumptions Applied 22 2.10 Root Zone Conductivity Root Zone Conductivity or soil conductivity is the conductivity in the first bucket. Conductivity rate typically varies among soil and land use classifications. As with the second bucket IMTA determined percolation rates for the first bucket for each sub-basin. Again the percolation values are applied as conductivity values by setting the Preferred Flow Direction Values to 0. Table 8 lists the values for percolation in the first bucket determined by IMTA and those applied to Root Zone Conductivity in the WEAP model. A Key Assumption for the Root Zone Conductivity of each basin was created. The land use branches within each sub-basin and soil group then referenced the Key Assumption. Table 8. Root Zone Conductivity Basin Values 2.11 Preferred Flow Direction The Preferred Flow Direction parameter is used to partition flow out of the root zone layer to the lower soil layer or groundwater. Preferred flow direction can vary by land use classification and ranges from 0 to 1. A preferred flow direction of 1 indicates 100% horizontal flow direction while 0 indicates 100% vertical flow direction. A Key Assumption for the Preferred Flow Direction of each land use category was created with a value equal to zero to effectively apply the available percolation data from IMTA to hydraulic conductivity input, refer to Figure 2. Again, each land use category branch then references the Key Assumption for Preferred Flow Direction. I M T A S o i l P r o f i l e W E A P P e r c o l a t i o n M a x R a t e ( m m /h ) S o i l C o n duct i v i t y ( m m /d ) P e gui s 1 6. 00 144 S a c r a m e n t o 2 25. 00 600 L a s B ur r a s 3 2. 00 48 L ui s L . L e o n 4 8. 00 192 F C O. M a der o 5 25. 00 600 Vi l l a l b a 6 23. 00 552 C o n c h o s 7 1. 25 30 J i m i n e z 8 26. 00 624 C h uv i s c a r 9 1. 25 30 E l R e j o n 10 10. 00 240 C h i h uah ua 11 1. 00 24 L l a n i t o s 12 2. 00 48 P i c o del A gui l a 13 5. 00 120 S a n A n t o n i o 14 25. 00 600 S a n Ga b r i e l 15 0. 28 7 P uen t e F F C C 16 0. 28 7 P a r r a l 17 15. 00 360 C o l i n a 18 15. 00 360 L a B o qui l l a 19 5. 00 120 Oj i n a ga 20 6. 00 144 I M T A I DB a s i n Na m e 23 2.12 Initial Z1 The Initial Z1 parameter is the relative storage given as a percentage of the total effective storage of the Root Zone Water Capacity at the beginning of a simulation. Therefore, like Root Zone Water Capacity this parameter typically varies with the land use, however, because the Root Zone Capacity parameter is coded by basin Initial Z1 was done in the same way. A Key Assumption for the Initial Z1 value of each sub-basin was created with a value equal to 20 percent. Each land use category within the sub-basins then references the Key Assumption. 2.13 Crop Coefficient, Kc The crop coefficient, Kc, parameter represents the effects of vegetative evapotranspiration and soil evaporation, for this reason the parameter varies by land class type. The parameter was created to study the required soil moisture to maximize crop biomass production; hence, Kc is typically used to calculate the required evapotranspiration using the equation: (Evapotranspiration)required = Kc * (Evapotranspiration)reference The Rio Bravo study (World Bank, 2006) determined actual and potential evapotranspiration, which relates to the required and potential evapotranspiration, respectively. These values were determined for fifteen land use categories which differ from the land use categories determined by IMTA and used in the WEAP model. Section 2.3 discusses these differences. Table 9 lists World Banks data and the crop coefficient values applied in WEAP. A Key Assumption for the Kc value of each land use category was created and the land use branches within each sub-basin reference the Key Assumptions. Table 9. Crop Coefficient Values ET a c t ET pot ( m m ) ( m m ) U r b a n A r e a 803 1048 0. 77 85 U r b a n A r e a s 0. 77 W a t e r B o di e s 1578 1578 1. 00 30 W a t e r B o di e s 1. 00 I r r i ga t e d A gr i c ul t ur e ( de l t a ) 1202 1346 0. 89 I r r i ga t e d A gr i c ul t ur e ( v a l l e y ) 898 1040 0. 86 S upp l e m e n t a l I r r i ga t i o n 1242 1298 0. 96 50 N a t ur r a l l y I r r i ga t e d A r e a s 0. 96 L o w O pe n F o r r e s t 483 1272 0. 38 20 F o r r e s t G r a s s e s 0. 38 O a k F o r e s t 538 1747 0. 31 P i n e F o r r e s t 487 1272 0. 38 C h a ppa r r a l 481 1424 0. 34 70 H i gh G r a s s e s a n d S m a l l B r us h 0. 34 M i c r o ph y l l o us S c r ub l a n ds 237 501 0. 47 S c r ub l a n d w i t h R o s e t t e d V e ge t a t i o n 263 616 0. 43 T h o r n s c r ub l a n d T a m a ul i pa n 583 1254 0. 46 S ub m o n t a n e S c r ub l a n d 711 1605 0. 44 C ul t i v a t e d G r a s s l a n d 516 1129 0. 46 80 G r a z i n g P a s t ur e s 0. 46 N a t ur a l G r a s s l a n d 342 642 0. 53 60 S m a l l P a s t ur e G r a s s e s 0. 53 - - - - 90 W e t l a n d V e ge t a t i o n 0. 9 - - - - 95 W i t h o ut A ppa r e n t V e ge t a t i o n 0. 3 0. 35F o r r e s t10 75 O t h e r V e ge t a t i o n 0. 45 K c U s e d i n W e a p 40 I r r i ga t e d A r e a s 0. 88 L a n d U s e I M T A L a n d U s e C o de I M T A L a n d U s e C a t e go r yKc 24 2.14 Leaf Area Index Leaf Area Index (LAI) is a parameter that varies by land use and is used to control the surface runoff response. Runoff tends to decrease with higher values of LAI. LAI is not an easily determined parameter and there are many discrepancies in definitions and values, which apply to the same land use. Scurlock et al. (2001) compiled estimates of LAI published between 1932 and 2000 and produced a table of LAI values for fifteen categories of vegetative land use. The LAI values produced by Scurlock et al. are shown in Table 10. The values as applied to the WEAP model are listed in Table 11. Referring to Table 1, simulated flow results from WEAP are highly sensitive to Kc and LAI parameters. Future work should include an investigation into the importance of these values and the cost benefit for applying site specific LAI versus literature values and importance of seasonal Kc and LAI values Table 10. LAI Values Scurlock et al. , 2001 25 Table 11. LAI Values 2.15 Precipitation Precipitation data was obtained from the IMTA DSS file as spatially weighted daily values for each sub-basin. IMTA (Martinez, 2005) applied the Theissen polygon method to sixty-five climate stations, shown in Figure 16, to determine the incremental precipitation per day for each sub-basin for the period of 1980 to 1999. The differences between the sub-basins as IMTA delineated them and as they were delineated for this study (Figure 5) transfer to the application of the spatially weighted data obtained from IMTA, however, the effect is expected to be minor is not investigated or addressed in this report. Figure 16. IMTA Precipitation Stations (Martinez, 2005) Sc u r l o c k e t a l . , 2 0 0 1 L a n d U s e C a t e g o r y F o r e s t A v e r a g e 10 F o r r e s t 5 . 1 8 F o r e s t A v e r a g e ( B o D B L a n d B O E N L ) 20 F o r r e s t G r a s s e s 3 . 0 7 30 W a t e r B o d i e s 0 . 1 0 C r o p s 40 I r r i g a t e d A r e a s 4 . 2 2 C r o p s 50 N a t u r r a l l y I r r i g a t e d A r e a s 4 . 2 2 G r a s s l a n d 60 Sm a l l Pa s t u r e G r a s s e s 2 . 5 0 Sh r u b 70 H i g h G r a s s e s a n d Sm a l l B r u s h 2 . 0 8 Sh r u b 75 O t h e r V e g e t a t i o n 2 . 0 8 G r a s s l a n d 80 G r a z i n g Pa s t u r e s 2 . 5 0 85 U r b a n A r e a s 8 . 0 0 W e t l a n d s 90 W e t l a n d V e g e t a t i o n 6 . 3 4 D e s e r t 95 W i t h o u t A p p a r e n t V e g e t a t i o n 1 . 3 1 I M T A L a n d U s e C o d e I M T A L a n d U s e C a t e g o r y L A I U s e d i n W e a p 26 The precipitation data extracted from the DSS file was formatted in Excel and saved as a CSV file, which can be read by WEAP as a daily time series expression. An excerpt of the Excel file is shown in Figure 17 and Figure 18 shows how the time series for each sub-basin are entered in WEAP. Figure 17. Excerpt of Precipitation Data CSV File Figure 18. WEAP Precipitation Time Series Expression # W E A P CO L UM N NU M B E R 3 4 5 # IM T A B A S IN NUM B E R 1 2 3 # RE P O RT B A S IN NA M E P E G UIS S A CR A M E NT O L A S B UR RA S # HM S B A S IN NA M E P E G UIS S A CR A M E NT O B UR RA S # G A G E G A G E G A G E # MM MM MM # Y e a r Day Co u n t M o n th Day P E R- CU M P E R- CU M P E R- CU M 1980 1 1 1 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 2 1 2 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 3 1 3 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 4 1 4 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 5 1 5 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 6 1 6 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 7 1 7 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 8 1 8 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 9 1 9 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 10 1 10 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 11 1 11 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 12 1 12 0 .5 0 0 0 .0 0 0 0 .2 0 0 1980 13 1 13 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 14 1 14 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 15 1 15 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 16 1 16 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 17 1 17 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 18 1 18 0 .0 0 0 0 .0 0 0 0 .0 0 0 1980 19 1 19 0 .0 0 0 0 .0 0 0 0 .0 0 0 27 2.16 Temperature, Wind and Humidity Temperature data is entered in degrees Celsius. Humidity is the relative humidity entered as a percentage and Wind values are entered in meters per second. Ideally, each of the parameters should be entered as time series data following the chosen time step of the model, days in this case. However, only averaged monthly data in raster format was available for the Rio Bravo basin from the SWAT data set. The Zonal Statistics tool in Arc Toolbox was used to determine a single average monthly value for each sub-basin for each of the three parameters. The monthly values were repeated for every day of the corresponding month and read into WEAP as a time series expression from a CSV file. Figure 19. Excerpt of Temperature Data CSV File # W E A P C o l um n 3 4 5 # SW A T D A T A # I M T A B A SI N N U M B E R 1 2 3 # B A SI N N A M E PE G U I S SA C R A M E N T O L A S B U R R A S # Y e a r M o n t h T e m p T e m p T e m p # 1 9. 37 8. 8 10 # 2 11. 57 10. 46 12. 24 # 3 14. 9 13. 52 15. 63 # 4 19. 87 17. 56 20. 35 # 5 24. 22 21. 3 23. 84 # 6 28. 16 24. 88 26. 87 # 7 28. 32 24. 42 25. 95 # 8 27. 4 23. 5 25. 29 # 9 25 21. 38 23. 06 # 10 20. 33 17. 46 19. 32 # 11 13. 92 12. 4 13. 98 # 12 9. 86 9 10. 51 # Y e a r D a y C o un t M o n t h D a y F l o w ( M m ^ 3)F l o w ( M m ^ 3) F l o w ( M m ^ 3) 1980 1 1 1 9. 37 8. 8 10 1980 2 1 2 9. 37 8. 8 10 1980 3 1 3 9. 37 8. 8 10 1980 4 1 4 9. 37 8. 8 10 1980 5 1 5 9. 37 8. 8 10 1980 6 1 6 9. 37 8. 8 10 1980 7 1 7 9. 37 8. 8 10 28 Figure 20. WEAP Temperature Data for the Peguis Sub-basin 2.17 Latitude WEAP uses latitude in calculating water temperature in water quality models and losses in the Soil Moisture Method. The latitude at the centroid of each catchment was estimated from Figure 21. The estimated values are summarized in Table 12. 29 Figure 21. Latitude and Longitude (Martinez 2005) Table 12. Sub-basin Latitude Values Basin Name IMTA ID Latitude DD Peguis 1 29.50 Sacramento 2 29.00 Las Burras 3 28.00 Luis L. Leon 4 29.00 FCO. Madero 5 28.50 Villalba 6 28.00 Conchos 7 28.00 Jimenez 8 27.50 Chuviscar 9 27.00 El Rejon 10 28.50 Chihuahua 11 28.50 Llanitos 12 26.50 Pico del Aguila 13 26.50 San Antonio 14 26.50 San Gabriel 15 26.50 Puente FFCC 16 26.50 Parral 17 26.50 Colina 18 27.00 La Boquilla 19 29.50 Ojinaga 20 29.50 2.18 Melting Point, Freezing Point and Initial Snow The remaining three climatic parameters are Melting Point, Freezing Point and Initial Snow value. Melting Point is the threshold for snow melt in degrees Celsius. A value of 0 was used for 30 all basins. Freezing Point is the threshold for snow accumulation in degrees Celsius. A value of 0 was used for all basins. The Initial Snow value is the snow accumulation at the beginning of the first month of the simulation. Zero was chosen as the initial value for each of the twenty catchments. 3 HEC-HMS Data The HEC-HMS model created by IMTA is a tool against which to compare and evaluate the results of the WEAP model. Six locations were selected as comparison points between the WEAP and HEC-HMS simulated flows. These particular points were selected because they are also naturalized and historical flow gage locations. Although, these flow data are not addressed in this project, the points chosen for comparison are significant to future work to allow for congruency in locations at which to evaluate and assess the model. Physical locations of the gages are shown in Figure 22. The virtual location of the gages within the WEAP model structure are listed in Table 13. Figure 22. Gauges Locations 31 Table 13. Rio Conchos Gages Gage Name Gage Location HydroID Rio Conchos at Ojinaga Bottom of Reach Rio_Conchos_20 2020200051 Rio Conchos at El Granero Top of Reach Rio_Conchos_1 2020200004 Rio Conchos at Las Burras Bottom of Reach Rio_Conchos_3 2020200003 Rio San Pedro at Villalba Top of Reach Rio_San_Pedro_5 2020200001 Rio Conchos at Presa La Boquilla Top of Reach Rio_Conchos_18 2020200005 Rio Florido at Cd. Jimenez Top of Reach Rio_Florido_3 2020200002 Many problems were encountered while using the HEC-HMS model received from IMTA including errors for missing input data and incomplete run configurations to correlate to the gage locations. Model runs for La Boquilla, Jimenez and Villalba were completed but with uncertainty. Run configurations for the remaining three gage locations were not completed due to either parameter errors or seemingly missing model components. For these reasons a request was made to IMTA for the direct results of the model; subsequently IMTA provided a DSS file with the hourly computed flow rate for all run configurations for the year 1980. Table 15 contains the run configuration names within the DSS file and the correlation to the gage locations. There were large discrepancies between the computed volumes determined from the model for La Boquilla, Jimenez and Villalba and the values from the DSS file. The source of this discrepancy is unknown at this time possible causes are user error or the use of an incomplete model. The values obtained from each source are listed below: Table 14. Comparison of HEC-HMS Model and DSS Flow Values Gage Name HEC-HMS Flow (Mm3) DSS Flow from IMTA (Mm3) Rio Conchos at Ojinaga - 2,740 Rio Conchos at El Granero - 2,620 Rio Conchos at Las Burras - 2,409 Rio San Pedro at Villalba 4,171 864 Rio Florido at Cd. Jimenez 2,064 262 Rio Conchos at Presa La Boquilla 2,064 923 32 Table 15. Sub-basins Contributing to Gages Gage Name Contributing Sub-basins Rio Conchos at Ojinaga All 23 runs Rio Conchos at El Granero All Runs Except C-PEGUIS OJINAGA Rio Conchos at Las Burras BURRAS ALTA C-LAS BURRAS C-FCO I MAD C-VILLALBA C-CONCHOS C- JIMENEZ C-LLANITOS C- PICO DE AGUILA C- SAN ANTONIO C- SAN GABRIEL C-PUENTE FFC C- PARRAL C.-COLINA C. BOQUILLA BOQUILLA1 BOQUILLA2 Rio San Pedro at Villalba C-VILLALBA Rio Florido at Cd. Jimenez C- JIMENEZ C- PICO DE AGUILA C- SAN ANTONIO C- SAN GABRIEL C-PUENTE FFC Rio Conchos at Presa La Boquilla C-LLANITOS C. BOQUILLA BOQUILLA1 BOQUILLA2 4. Calibration Process for 1980 For this initial calibration, naturalized flows from TCEQ (Brandes, 2003) for the period from January to December 1980 were used to calibrate the WEAP model. The calibration involved both quantitative and qualitative evaluation of the hydrologic response of each tributary in each 33 sub-basin. After that, parameters were adjusted to reproduce the naturalized monthly and annual stream flow. To this end, the soil moisture method in WEAP model was used and the relevant parameters are described below. 4.1. Root Zone Water Capacity Initially, values of Root Zone Capacity above 1000 mm were assumed according to the land use; however, the results did not reproduce the trend of naturalized monthly flows in each stream gage considered for the analysis. For this reason, the values of root zone water capacity were reduced; from 160 mm to 400 mm in the Puente FFC, Ojinaga, and Peguis sub-basins, respectively. Values found for each sub-basin can be seen Table 24. 4.2. Root Zone Conductivity Root Zone Conductivity is a very important parameter in the calibration process which controls the transmission of flow to the lower soil layer and the interflow. The inter flow depends of the preferred flow direction; for our study case, we have assumed it is equal cero or vertical flow which means that there is not inter flow. The flow volume of each catchment from the upper layer to the lower layer is estimated with a simple expression of its relative storage. ? ?? ?Ni ip zf l o w d i rp r e fkAtV 1 211 *).1(*)( According to the values found for the Villalba sub-basin (see table 24), and substituting in the expression above, we have an average volume of percolation of 9.057 Mm3/day. If this parameter is reduced, the stream flow is increased; therefore, the transmission of flow volume toward the lower layer is also reduced. 4.3. Initial Root Zone Water Capacity An Initial Root Zone Water Capacity at the beginning of a simulation was assumed for each sub- basin. This parameter varies from 5 to 30 % in some sub-basins. Surface runoff is directly correlated with the initial storage, 1z ; if 1z is increased, the runoff as well. The values for this parameter are shown in Table 16. 34 Table 16. WEAP Upper Layer Soil Parameters Calibrated for the Conchos River Basin Drainage Bucket 1 Sub-basin Area Root Zone Root Zone Initial Km2 Capacity Conductivity Z1 Mm mm/day % Peguis 7999.2972 400 120 10 Sacramento 1042.6059 250 70 10 Las Burras 11309.4666 350 80 10 Luis L. Leon 5085.5131 350 120 5 FCO. I Madero 1211.3488 260 24 20 Villalba 9556.8624 253 13 27 Conchos 1114.3944 280 14 30 Jimenez 4422.9591 257 28 10 Chuviscar 106.0884 250 75 10 El Rejon 146.8494 250 70 10 Chihuahua 399.9897 250 80 10 Llanitos 1829.9295 270 4 30 Pico de Aguila 647.6067 250 45 5 San Antonio 821.1609 210 24 10 San Gabriel 305.8525 210 24 10 Puente FFCC 1270.6609 160 14 10 Parral 363.7890 250 12 10 Colina 259.0569 280 28 30 La Boquilla 18931.9788 315 4 30 Ojinaga 983.4705 400 100 10 4.4. Lower Deep Water Capacity Values assumed for the Lower Deep Water Capacity are shown in Table 17. It is likely the high values found in some sub-basins shows the existence of deep aquifers. Initially, values between 100 mm to 300 mm as deep water capacity were proposed for the sub-basins. This range of values gave a high stream flow per year (relative to the naturalized flows), and a great amount of base flow was generated from September to December. For this reason, upper values 2800 mm were assumed and evaluated. 4.5. Lower Deep Conductivity The Lower Deep Conductivity controls the transmission of base flow in each sub-basin. This parameter can be estimated with the following expression: 35 )()( 1 222?? ?Ni if zkAtB where iA is the area of the land use cover fraction, i, 2k is the conductivity rate of the lower layer at full saturation ( 00.12?z ) in mm/month, and 2z is the relative storage given as a percentage of the effective storage of the lower soil layer. From the expression mentioned above, initial hydraulic conductivity was estimated of the following way: 222 )/( z ABk if? The baseflow ( fB ) can be estimated with different methods depending of hydrologic behavior of basin in study. To this end, considering the limitation of information, it is possible to make rough calculations of the baseflow using the straight line method. For example, for the Villalba sub-basin with a drainage area of 9,557 km2, the baseflow was estimated to be 1.5 Mm3/month for 1980, and the average initial store 2z was assumed to be 4% for all fractions i, the resulting hydraulic conductivity is m onthmmk /982 ? , or equal to 3.2 mm/day. In this case, in order to have more accurate results from the hydrologic simulation, this value was adjusted to 60 mm/month. The 2k parameter found for each sub-basin can be seen in Table 17. 4.6. Initial Lower Layer Capacity Different values of Lower Layer Initial Storage were assumed in the hydrologic simulation. At the beginning of the simulation, percentages around 40 ? 50 % were assumed; nevertheless, in many cases, the baseflow was more than 50 % with regard to the stream flow. For example, in the Villalba Sub-basin, the baseflow was more than 70%; for this reason, small values of 2z were assumed varying from 4% to 10%. The initial storage 2z assumed for each sub-basin is shown in Table 17. 36 Table 17. WEAP Lower Layer Soil Parameters Calibrated for the Rio Conchos Basin Drainage Bucket 2 Sub-basin Area Deep Water Deep Water Conductivity Initial km2 Capacity Z2 Mm mm/day mm/month % Peguis 7999.2972 4500 1.0 30.0 5 Sacramento 1042.6059 3500 2.0 60.0 5 Las Burras 11309.4666 4500 3.0 90.0 10 Luis L. Leon 5085.5131 30000 1.0 30.0 5 FCO. I Madero 1211.3488 3500 2.0 60.0 6 Villalba 9556.8624 2800 2.0 60.0 4 Conchos 1114.3944 3700 3.0 90.0 10 Jimenez 4422.9591 5000 0.5 15.0 4 Chuviscar 106.0884 3500 2.0 60.0 5 El Rejon 146.8494 3500 2.0 60.0 10 Chihuahua 399.9897 3500 2.0 60.0 5 Llanitos 1829.9295 25000 2.0 60.0 10 Pico de Aguila 647.6067 4500 0.5 15.0 5 San Antonio 821.1609 4000 0.5 15.0 5 San Gabriel 305.8525 4000 0.5 15.0 5 Puente FFCC 1270.6609 4000 0.5 15.0 5 Parral 363.7890 4000 3.0 90.0 5 Colina 259.0569 4000 2.0 60.0 10 La Boquilla 18931.9788 30000 3.5 105.0 10 Ojinaga 983.4705 4500 1.0 30.0 5 5 Results This section compares the simulated flows from the IMTA HEC-HMS model and the WEAP model for the calendar year 1980. This specific year (1980) was selected because it is the first year within in the period of record, 1980 to 1999, for which IMTA utilized climatic data to develop their HEC-HMS model. The year is also within the period of record, 1980 to 1985, against which IMTA calibrated the values of Soil Capacity and Soil Conductivity. Before comparing the flows of the two models the total annual precipitation was determined for each basin and the average annual precipitation was then determined for each of the gage locations. The average annual precipitation was then multiplied by the drainage area to determine to the maximum possible flow volume at each gage site. This serves as a good reference point to evaluate the rough accuracy or plausibility of the simulated flow output. Table 18 contains the total annual precipitation of each sub-basin. 37 Table 18. Sub-basin Annual Precipitation for 1980 Basin Name IMTA ID 1980 Annual Precipitation Total (mm) Area (sq km) Peguis 1 350.9 7999.30 Sacramento 2 410.5 1042.61 Las Burras 3 332.4 11309.47 Luis L. Leon 4 388.9 5085.51 FCO. Madero 5 415.8 1211.35 Villalba 6 447.2 9556.86 Conchos 7 357.7 1114.39 Jimenez 8 400.8 4422.96 Chuviscar 9 424.4 106.09 El Rejon 10 441.6 146.85 Chihuahua 11 379.5 399.99 Llanitos 12 705.9 1829.93 Pico del Aguila 13 447.5 647.61 San Antonio 14 451.5 821.16 San Gabriel 15 441.1 305.85 Puente FFCC 16 414.5 1270.66 Parral 17 442.6 363.79 Colina 18 400.5 259.06 La Boquilla 19 514.935 18931.98 Ojinaga 20 419.196 983.47 5.1. Annual Streamflow Comparison Table 19 shows the annual stream flows simulated by WEAP in each selected station; likewise, the maximum possible flows as resulted to the integration of the drainage area and average precipitation upstream of each gage station. The largest difference between the naturalized flow and WEAP simulated stream flow is at the Ojinaga station where the model tends to overestimate the stream flow by about 16%. Most likely, this behavior is because for the 1980 calibration year, the naturalized flow estimated for the Ojinaga station is smaller than those located upstream of the basin; on the other hand, the model tends to over-estimate the flow at the outlet of the basin. 38 Table 19. Annual TCEQ Naturalized and WEAP Simulated Flows for the Rio Conchos Basin Gage Name Drainage Area (km2) Ave. Precipitation (mm) Stream flows (Mm3) Ratio WEAP/TCEQ Ratio WEAP/HMS TCEQ WEAP HMS Rio Conchos at Ojinaga 67,809 429 2,029 2,362 2,740 1.16 0.86 Rio Conchos at El Granero 58,826 434 2,192 2,299 2,620 1.05 0.88 Rio Conchos at Las Burras 52,045 444 2,220 2,230 2,409 1.00 0.93 Rio San Pedro at Villalba 9,557 447 341 340 864 1.00 0.39 Rio Florido at Cd. Jiminez 7,468 431 123 123 262 1.00 0.47 Rio Conchos at La Boquilla 20,762 610 1,446 1,439 923 1.00 1.56 Figures Figure 23, Figure 24, Figure 25, Figure 26, Figure 27, and Figure 28 show the TCEQ naturalized flows and the WEAP simulated flows for the selected gage stations. Figure 23 is a plot of the monthly mean observed and simulated stream flow for the Rio San Pedro at Villalba (period from Jan to Dec 1980). In this stream gage, the model simulates less baseflow in the period January to April (33% on average of the naturalized flow). However, the simulated flows are more accurate in the summer and fall seasons; for example, the WEAP flow in August represents 95 % of the naturalized flow and 104% in September; which means that in this month, the model only overestimates the monthly peak flow in 4%. It is possible to increase the base flow in the first months of year, by using a larger 2z or 2k . However, if 2z is increased from 4% to 10%, the monthly stream flows in November and December increase to more than 80% and 100 % respectively. 5.2. Monthly Streamflow Comparison Figure 24 shows the monthly WEAP simulated flow for the Rio Conchos at La Boquilla. The most notable difficulty of the hydrologic simulation is to predict the spring and summer flow (Feb ?Jun) when the simulated values are around 45 % of the naturalized flows, except in January when the simulation is almost perfect. However, the simulation is more accurate for the fall and winter flows; even though they are overestimated. For example, from August to October when the most important flows occur the ratio the simulated to naturalized flow is 1.09; which means that the model only overestimated the stream flow by 9% on average during this period. Figure 25 shows the WEAP simulated and TCEQ naturalized flows for the Rio Florido at Cd. Jimenez. Despite of its perfect approximation of the annual flow, there are important differences in the monthly flows. From January to June, simulated flows are smaller than the naturalized flows, while peak flows are close to the naturalized flows (at least in September). It is possible to improve this performance, increasing the hydraulic conductivity, 1k of the upper soil layer to reduce the stream flows in August and September as well as 1z in order to increase the stream flows in April, May, June, and November, in all sub-basin located upstream of the Cd. Jimenez stream gage. Similar behavior is presented in Las Burras, El Granero, and Ojinaga gage stations whose tendencies can be seen in Figure 26, Figure 27, and Figure 28. However, the differences between observed and simulated flows are more noticeable in the last two stations located 39 downstream of the Rio Conchos at Las Burras, because the natural flows are reduced in those points. The rainfall contributions from the Luis Leon, Peguis, and Ojinaga sub-basins are quite small (< 1%) and most of this water is lost to evaporation and seepage along the river. In general, the model overestimates the peak flows produced in September by around 15 %. Figure 29, Figure 30, Figure 31, Figure 32, Figure 33, and Figure 34 are plots showing the daily stream flows simulated by WEAP in the selected points. Sim ulat ed and Nat ural stream fl ow s in t he Rio San P edro at Villal ba 0 50 100 150 200 250 M onth s Strea m fl ow in M m 3 N a t. Flo w S im u la te d W E A P Nat. F l ow 5.200 4.450 3.596 2.648 1.647 1.345 2.398 44.155 233.339 27.286 5.321 9.314 S i m ul ate d W E A P 1.142 1.346 1.391 1.367 1.419 1.441 2.538 42.099 242.464 24.046 11.790 9.101 E ne F eb Mar A pr May J un J ul A ug S ep Oct Nov Dec T o ta l y e a r : W E A P = 3 4 0 .1 4 N a t. Flo w = 3 4 0 .7 0 R a tio = W E A P /Na t. Flo w = 0 .9 9 8 Figure 23. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio San Pedro at Villalba 40 Sim ualt ed and Nat ural stream fl ow s in t he Rio Conchos at P res a La B oquil la 0 100 200 300 400 500 600 700 800 Mon ths Stre am fl ow i n Mm3 N a t. Flo w S im u la te d W E A P N a t. Flo w 1 9 .7 6 0 3 8 .4 3 0 5 3 .1 8 0 5 6 .7 1 0 4 7 .5 8 0 6 6 .0 5 0 7 7 .6 7 0 1 6 5 .5 4 0 7 1 3 .8 1 0 1 6 8 .6 0 0 9 .9 3 0 2 8 .8 2 0 S im u la te d W E A P 2 2 .4 4 6 2 1 .0 5 7 2 2 .6 1 7 2 1 .1 8 0 2 1 .8 9 1 2 2 .5 1 2 3 8 .3 9 1 1 8 9 .0 1 9 7 8 3 .5 0 0 1 7 8 .2 6 3 7 4 .2 1 4 4 3 .9 7 2 E n e Fe b M a r A p r M a y Ju n Ju l A g o S e p O ct N o v Dic T o ta l y e a r : W E A P = 1 4 3 9 . 0 6 2 N a t. Flo w = 1 4 4 6 .0 8 R a tio = W E A P /Na t. Flo w = 0 .9 9 5 Figure 24. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio Conchos at La Boquilla Sim ulat ed and Nat ural stream fl ow s in t he Rio F lorido at CD Gim enez 0 10 20 30 40 50 60 70 80 M onths S t rea m f l ow i n Mm 3 N a t. Flo w S im u la te d W E A P Nat. F l ow 0.000 0.000 0.052 1.497 2.396 3.685 3.156 15.067 67.362 27.124 1.821 0.730 S i m ul ate d W E A P 0.286 0.345 0.245 0.236 0.260 0.408 3.062 26.600 77.525 11.184 1.917 1.069 E ne F eb Ma r c h A pr May J un J ul A ug S ep Oct Nov Di c T o ta l y e a r : W E A P = 1 2 3 .1 3 6 N a t. Flo w = 1 2 2 .8 9 R a tio = W E A P /Na t. Flo w =1 . 0 0 Figure 25. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio Conchos at Jimenez 41 Sim ulat ed and Nat ural stream fl ow s in t he rio Conchos at Las Burr as 0 200 400 600 800 1000 1200 1400 Mon ths S t ream f l ow i n M m3 N a t. Flo w S im u la te d W E A P N a t. Flo w 3 7 .8 6 0 6 5 9 .9 8 9 6 4 5 .6 3 8 8 4 4 .4 0 8 3 6 4 .5 3 0 6 8 4 .1 0 7 7 8 6 .8 5 7 5 3 0 3 .3 9 3 5 1 0 4 3 .7 3 2 9 2 .0 5 4 4 7 1 .4 8 5 8 5 .8 1 5 4 S im u la te d W E A P 3 6 .4 2 1 5 3 4 .4 7 4 3 3 7 .2 7 1 3 3 5 .3 6 5 9 3 6 .5 9 3 3 3 7 .3 6 0 6 5 7 .9 9 7 2 2 9 6 .5 8 5 4 1 2 1 5 .9 1 4 1 2 4 5 .8 2 7 9 1 1 5 .1 0 4 1 8 1 .1 2 5 3 E n e Fe b M a r A p r M a y Ju n Ju l A u g S e p O ct N o v D e c T o ta l y e a r : W E A P = 2 2 3 0 . 0 4 N a t. Flo w = 2 2 1 9 .8 7 R a tio = W E A P /Na t. Flo w =1 . 0 0 4 Figure 26. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio Conchos at Las Burras Sim ulat ed and Nat ural stream fl ow s in t he Rio Conchos at El Graner o 0 200 400 600 800 1000 1200 1400 M onth s S t ream Flow i n M m3 N a t. Flo w S im u la te d W E A P N a t. Flo w 3 6 .3 4 6 5 5 .3 6 6 5 7 .2 4 0 5 0 .2 2 1 7 3 .3 2 0 8 6 .1 1 6 8 3 .5 9 6 3 1 6 .4 8 4 1 0 2 2 .7 7 5 2 8 4 .4 7 4 5 7 .4 3 4 6 9 .1 1 9 S im u la te d W E A P 3 7 .1 8 2 3 5 .2 2 0 3 8 .0 5 6 3 6 .1 8 5 3 7 .3 9 7 3 8 .3 9 5 5 9 .1 2 3 3 2 3 .8 0 0 1 2 4 1 .8 8 7 2 4 8 .9 7 4 1 1 8 .6 6 0 8 3 .7 6 4 E n e Fe b M a r A p r M a y Ju n Ju l A u g S e p O ct N o v D e c T o ta l y e a r : W E A P = 2 2 9 8 .6 4 N a t. Flo w = 2 1 9 2 .4 9 R a tio = W E A P /Na t. Flo w = 1 .0 4 Figure 27. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio Conchos at El Granero 42 Sim ulat ed and Nat ural stream fl ow s in t he Rio Conchos at Oj inaga 0 200 400 600 800 1000 1200 1400 M onth s S t ream f l ow i n M m3 N a t. Flo w S im u la te d W E A P Nat. F l ow 38.437 51.580 46.447 40.870 62.910 73.315 56.817 366.781 891.026 264.565 67.680 68.682 S i m ul ate d W E A P 37.981 35.981 38.904 37.032 38.272 39.333 60.723 346.393 1263.838 254.399 122.531 86.837 E ne F eb Mar A pr May J un J ul A ug S ep Oct Nov Dec T o ta l y e a r : W E A P = 2 3 6 2 .2 2 N a t. Flo w = 2 0 1 9 .1 0 8 R a tio = W E A P /Na t. Flo w = 1 .1 7 Figure 28. Monthly WEAP simulated and TCEQ naturalized streamflow in the Rio Conchos at Ojinaga Figure 29. Daily WEAP simulated streamflow in the Rio San Pedro at Villalba 43 Figure 30. Daily WEAP simulated streamflow in the Rio Conchos at La Boquilla Figure 31. Daily WEAP simulated streamflow in the Rio Florido at Jimenez 44 Figure 32. Daily WEAP simulated streamflow in the Rio Conchos at Las Burras Figure 33. Daily WEAP simulated streamflow in the Rio Conchos at El Granero 45 Figure 34. Daily WEAP simulated streamflow in the Rio Conchos at Ojinaga 6. Recommendations The results show a good approximation to both annual and monthly flows. However, the hydrologic model needs to be validated for a period of time larger, making adjustments as necessary in some soil parameters in order to improve the accuracy in the hydrologic response of the whole basin, considering that the hydraulic conductivity and initial storage increase when the available water in the soil layers increases. In other words, the hydraulic conductivity should be larger in a period with significant rainfall than a period with little to no rainfall. The integration of WEAP?s hydrologic flow predication capabilities into the existing WEAP model of the Rio Bravo basin creates a powerful tool for regional planners. There remains a great deal of work to bring this idea to fruition. The model could be useful in generating inflows to the basin under various sequences of future precipitation. These inflows could be used in the WEAP water management model (Danner et al., 2006) to assess the result for basin stakeholders on different scenarios of basin operation. 46 References Brandes, R.J. Company, Water Availability Modeling for the Rio Grande Basin: Naturalized Streamflow Data, Texas Commission on Environmental Quality, Austin, Texas, October 2003 Constance L. Danner, B.S., Daene C. McKinney, PhD., PE, and Rebecca L. Teasley, M.S.: Documentation and Testing of the WEAP Model for the Rio Grande/Bravo Basin, August 2006 (http://www.crwr.utexas.edu/reports/2006/rpt06-8.shtml) Jantzen, T., B. Klezendorf, J. Middleton, and J. Smith. WEAP Hydrology Modeling Applied: The Upper Rio Florido Rive Basin Martinez, J. and R.M. Zermeno and A.G. Lopez. Estudio Para la Gestion Integrada del Agua en la Cuenca del Rio Bravo, Instituto Mexicano De Tecnologia Del Agua. Report No. 2002- C-01-0569.A3, 2005 (in Spanish). Pati?o-Gomez, C. and D.C. McKinney. GIS for Large-Scale Watershed Observational Data Model, The University of Texas at Austin, Center for Research in Water Resources Online Report No. 2005-05, 2005 (http://www.crwr.utexas.edu/reports/2005/rpt05- 5.shtml) Scurlock, J. M. O., G. P. Asner, and S. T. Gower. 2001. Global Leaf Area Index Data from Field Measurements, 1932-2000. Data set. Available on-line [http://www.daac.ornl.gov] from the Oak Ridge National Laboratory Distributed Active Archive Center, Oak Ridge, Tennessee, U.S.A. Sieber, J., Swartz, C., and Huber-Lee, A. (2005). Water Evaluation and Planning System User Guide for WEAP 21. Stockholm Environment Institute, Tellus Institute, Boston, Massachusetts. Texas Commission on Environmental Quality. Water Availability Modeling for the Rio Grande Basin: Naturalized Streamflow Data: Final Report. Prepared by R.J. Brandes Company. October 2003. World Bank. The Hydrological Flow Path and Options for Sustainable Water Resources Management in the Overexploited Rio Bravo Basin: A Preliminary Analysis from Remote Sensing and Hydrological Modeling: Final Draft. Washington DC, March, 2006. 47 Appendix 1. Soil Land Use Intersects by Sub-basin P e g u i s S u b b a s i n N o . 1 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 10809 0 .0 1 0 .0 2 20 0 0 .0 0 0 .0 0 30 38164 0 .0 4 0 .0 5 40 6380071 6 .3 8 8 .9 1 50 6502437 6 .5 0 9 .0 8 60 20395521 2 0 .4 0 2 8 .4 9 70 27555676 2 7 .5 6 3 8 .4 9 75 5846984 5 .8 5 8 .1 7 80 2563815 2 .5 6 3 .5 8 85 576152 0 .5 8 0 .8 0 90 1721146 1 .7 2 2 .4 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 7 1 .5 9 % o f T o ta l B a si n A r e a 0 .9 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 28325 0 .0 3 0 .0 0 20 1509 0 .0 0 0 .0 0 30 210334 0 .2 1 0 .0 1 40 718386 0 .7 2 0 .0 3 50 5966293 5 .9 7 0 .2 8 60 379643859 3 7 9 .6 4 1 7 .5 1 70 1699345239 1 6 9 9 .3 5 7 8 .3 6 75 8675007 8 .6 8 0 .4 0 80 63771996 6 3 .7 7 2 .9 4 85 83173 0 .0 8 0 .0 0 90 4760341 4 .7 6 0 .2 2 95 5466388 5 .4 7 0 .2 5 T o ta l S o i l A r e a ( sq k m ) 2 1 6 8 .6 7 % o f T o ta l B a si n A r e a 2 8 .1 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 497 0 .0 0 0 .0 0 20 30915 0 .0 3 0 .0 0 30 61411 0 .0 6 0 .0 0 40 1953737 1 .9 5 0 .0 6 50 18742950 1 8 .7 4 0 .5 9 60 1081764821 1 0 8 1 .7 6 3 3 .8 2 70 1677846763 1 6 7 7 .8 5 5 2 .4 5 75 117288454 1 1 7 .2 9 3 .6 7 80 287986809 2 8 7 .9 9 9 .0 0 85 935075 0 .9 4 0 .0 3 90 3920024 3 .9 2 0 .1 2 95 8285730 8 .2 9 0 .2 6 T o ta l S o i l A r e a ( sq k m ) 3 1 9 8 .8 2 % o f T o ta l B a si n A r e a 4 1 .5 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 784935 0 .7 8 0 .0 3 20 0 0 .0 0 0 .0 0 30 779358 0 .7 8 0 .0 3 40 268422 0 .2 7 0 .0 1 50 1145767 1 .1 5 0 .0 5 60 196367414 1 9 6 .3 7 8 .6 8 70 1896450935 1 8 9 6 .4 5 8 3 .8 1 75 43258657 4 3 .2 6 1 .9 1 80 119299233 1 1 9 .3 0 5 .2 7 85 152992 0 .1 5 0 .0 1 90 348263 0 .3 5 0 .0 2 95 3877069 3 .8 8 0 .1 7 T o ta l S o i l A r e a ( sq k m ) 2 2 6 2 .7 3 % o f T o ta l B a si n A r e a 2 9 .3 8 48 S a c r a m e n t o S u b b a s i n N o . 2 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 2080698 2 .0 8 3 .7 5 60 8281386 8 .2 8 1 4 .9 3 70 10034072 1 0 .0 3 1 8 .0 9 75 0 0 .0 0 0 .0 0 80 34897129 3 4 .9 0 6 2 .9 2 85 136785 0 .1 4 0 .2 5 90 0 0 .0 0 0 .0 0 95 28815 0 .0 3 0 .0 5 T o ta l S o i l A r e a ( sq k m ) 5 5 .4 6 % o f T o ta l B a si n A r e a 5 .5 5 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 421382 0 .4 2 0 .1 1 20 9461908 9 .4 6 2 .3 6 30 0 0 .0 0 0 .0 0 40 864396 0 .8 6 0 .2 2 50 106013400 1 0 6 .0 1 2 6 .4 6 60 57268101 5 7 .2 7 1 4 .2 9 70 27067093 2 7 .0 7 6 .7 6 75 0 0 .0 0 0 .0 0 80 177389566 1 7 7 .3 9 4 4 .2 8 85 22133664 2 2 .1 3 5 .5 2 90 0 0 .0 0 0 .0 0 95 3602 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 4 0 0 .6 2 % o f T o ta l B a si n A r e a 4 0 .1 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 58019667 5 8 .0 2 1 0 .6 9 20 256652163 2 5 6 .6 5 4 7 .2 7 30 54019 0 .0 5 0 .0 1 40 3601 0 .0 0 0 .0 0 50 3521854 3 .5 2 0 .6 5 60 9638093 9 .6 4 1 .7 8 70 20659107 2 0 .6 6 3 .8 0 75 0 0 .0 0 0 .0 0 80 194298060 1 9 4 .3 0 3 5 .7 9 85 18005 0 .0 2 0 .0 0 90 0 0 .0 0 0 .0 0 95 86430 0 .0 9 0 .0 2 T o ta l S o i l A r e a ( sq k m ) 5 4 2 .9 5 % o f T o ta l B a si n A r e a 5 4 .3 5 T o ta l B a si n A r e a ( sq k m ) 9 9 9 .0 3 49 L a s B u r r a s S u b b a s i n N o . 3 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 1348679 1 .3 5 0 .2 2 20 3491398 3 .4 9 0 .5 6 30 268862 0 .2 7 0 .0 4 40 51729070 5 1 .7 3 8 .2 5 50 86410404 8 6 .4 1 1 3 .7 8 60 82011311 8 2 .0 1 1 3 .0 8 70 116617155 1 1 6 .6 2 1 8 .6 0 75 34740718 3 4 .7 4 5 .5 4 80 239524731 2 3 9 .5 2 3 8 .2 1 85 9928058 9 .9 3 1 .5 8 90 0 0 .0 0 0 .0 0 95 816475 0 .8 2 0 .1 3 T o ta l S o i l A r e a ( sq k m ) 6 2 6 .8 9 % o f T o ta l B a si n A r e a 5 .7 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 10808 0 .0 1 0 .0 0 20 57645 0 .0 6 0 .0 0 30 33211 0 .0 3 0 .0 0 40 146417057 1 4 6 .4 2 6 .8 1 50 152207712 1 5 2 .2 1 7 .0 8 60 630808464 6 3 0 .8 1 2 9 .3 4 70 489318169 4 8 9 .3 2 2 2 .7 6 75 44466833 4 4 .4 7 2 .0 7 80 666539542 6 6 6 .5 4 3 1 .0 0 85 9861356 9 .8 6 0 .4 6 90 10522179 1 0 .5 2 0 .4 9 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 1 5 0 .2 4 % o f T o ta l B a si n A r e a 1 9 .5 4 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 100871 0 .1 0 0 .0 0 20 8485561 8 .4 9 0 .1 3 30 3383403 3 .3 8 0 .0 5 40 520324527 5 2 0 .3 2 7 .7 9 50 837594285 8 3 7 .5 9 1 2 .5 3 60 2302947107 2 3 0 2 .9 5 3 4 .4 6 70 790753723 7 9 0 .7 5 1 1 .8 3 75 618490117 6 1 8 .4 9 9 .2 5 80 1451828717 1 4 5 1 .8 3 2 1 .7 2 85 24897806 2 4 .9 0 0 .3 7 90 123750614 1 2 3 .7 5 1 .8 5 95 887566 0 .8 9 0 .0 1 T o ta l S o i l A r e a ( sq k m ) 6 6 8 3 .4 4 % o f T o ta l B a si n A r e a 6 0 .7 5 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 8725885 8 .7 3 0 .5 7 20 7897203 7 .9 0 0 .5 1 30 5232815 5 .2 3 0 .3 4 40 123439312 1 2 3 .4 4 8 .0 1 50 98788610 9 8 .7 9 6 .4 1 60 177172330 1 7 7 .1 7 1 1 .4 9 70 346895520 3 4 6 .9 0 2 2 .5 0 75 87087745 8 7 .0 9 5 .6 5 80 661684417 6 6 1 .6 8 4 2 .9 2 85 21056387 2 1 .0 6 1 .3 7 90 3291647 3 .2 9 0 .2 1 95 429987 0 .4 3 0 .0 3 T o ta l S o i l A r e a ( sq k m ) 1 5 4 1 .7 0 % o f T o ta l B a si n A r e a 1 4 .0 1 T o ta l B a si n A r e a ( sq k m ) 1 1 0 0 2 .2 8 50 L u i s L L e o n S u b b a s i n N o . 4 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 244063 0 .2 4 3 .0 8 50 978266 0 .9 8 1 2 .3 5 60 1685046 1 .6 9 2 1 .2 7 70 3233952 3 .2 3 4 0 .8 1 75 0 0 .0 0 0 .0 0 80 1782518 1 .7 8 2 2 .5 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 7 .9 2 % o f T o ta l B a si n A r e a 0 .1 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 22933 0 .0 2 0 .0 0 20 183851 0 .1 8 0 .0 2 30 481725 0 .4 8 0 .0 5 40 15393824 1 5 .3 9 1 .6 9 50 45813616 4 5 .8 1 5 .0 3 60 207607031 2 0 7 .6 1 2 2 .7 8 70 388365033 3 8 8 .3 7 4 2 .6 1 75 7949771 7 .9 5 0 .8 7 80 234443639 2 3 4 .4 4 2 5 .7 2 85 9379959 9 .3 8 1 .0 3 90 7206 0 .0 1 0 .0 0 95 1717651 1 .7 2 0 .1 9 T o ta l S o i l A r e a ( sq k m ) 9 1 1 .3 7 % o f T o ta l B a si n A r e a 1 8 .0 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 513492 0 .5 1 0 .0 2 20 132648 0 .1 3 0 .0 1 30 460195 0 .4 6 0 .0 2 40 47912375 4 7 .9 1 1 .8 9 50 175992744 1 7 5 .9 9 6 .9 3 60 1116136960 1 1 1 6 .1 4 4 3 .9 3 70 671332564 6 7 1 .3 3 2 6 .4 3 75 58332673 5 8 .3 3 2 .3 0 80 348657461 3 4 8 .6 6 1 3 .7 2 85 111326273 1 1 1 .3 3 4 .3 8 90 52882 0 .0 5 0 .0 0 95 9619594 9 .6 2 0 .3 8 T o ta l S o i l A r e a ( sq k m ) 2 5 4 0 .4 7 % o f T o ta l B a si n A r e a 5 0 .2 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 5980885 5 .9 8 0 .3 7 20 10580280 1 0 .5 8 0 .6 6 30 13338357 1 3 .3 4 0 .8 3 40 2970647 2 .9 7 0 .1 9 50 14201108 1 4 .2 0 0 .8 9 60 223999441 2 2 4 .0 0 1 4 .0 1 70 937009929 9 3 7 .0 1 5 8 .5 9 75 3783948 3 .7 8 0 .2 4 80 368206205 3 6 8 .2 1 2 3 .0 2 85 16899482 1 6 .9 0 1 .0 6 90 414581 0 .4 1 0 .0 3 95 1867838 1 .8 7 0 .1 2 T o ta l S o i l A r e a ( sq k m ) 1 5 9 9 .2 5 % o f T o ta l B a si n A r e a 3 1 .6 1 T o ta l B a si n A r e a ( sq k m ) 5 0 5 9 .0 1 51 F c o . I M a d e r o S u b b a s i n N o . 5 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 14409 0 .0 1 0 .0 1 20 3602 0 .0 0 0 .0 0 30 718311 0 .7 2 0 .2 9 40 617966 0 .6 2 0 .2 5 50 96833 0 .1 0 0 .0 4 60 96709578 9 6 .7 1 3 8 .8 6 70 127171804 1 2 7 .1 7 5 1 .1 1 75 3362445 3 .3 6 1 .3 5 80 19442604 1 9 .4 4 7 .8 1 85 0 0 .0 0 0 .0 0 90 320622 0 .3 2 0 .1 3 95 384658 0 .3 8 0 .1 5 T o ta l S o i l A r e a ( sq k m ) 2 4 8 .8 4 % o f T o ta l B a si n A r e a 2 0 .5 5 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 10806 0 .0 1 0 .0 1 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 529854 0 .5 3 0 .5 7 60 8404140 8 .4 0 9 .0 1 70 64451280 6 4 .4 5 6 9 .1 0 75 244893 0 .2 4 0 .2 6 80 19180655 1 9 .1 8 2 0 .5 7 85 0 0 .0 0 0 .0 0 90 199476 0 .2 0 0 .2 1 95 246604 0 .2 5 0 .2 6 T o ta l S o i l A r e a ( sq k m ) 9 3 .2 7 % o f T o ta l B a si n A r e a 7 .7 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 21615 0 .0 2 0 .0 1 20 177022 0 .1 8 0 .0 5 30 850377 0 .8 5 0 .2 3 40 454043 0 .4 5 0 .1 2 50 5502779 5 .5 0 1 .4 6 60 66640433 6 6 .6 4 1 7 .6 8 70 85946828 8 5 .9 5 2 2 .8 0 75 6576606 6 .5 8 1 .7 4 80 209290902 2 0 9 .2 9 5 5 .5 2 85 0 0 .0 0 0 .0 0 90 1296603 1 .3 0 0 .3 4 95 218890 0 .2 2 0 .0 6 T o ta l S o i l A r e a ( sq k m ) 3 7 6 .9 8 % o f T o ta l B a si n A r e a 3 1 .1 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 9180988 9 .1 8 1 .8 7 20 9229651 9 .2 3 1 .8 8 30 16087745 1 6 .0 9 3 .2 7 40 930897 0 .9 3 0 .1 9 50 2272540 2 .2 7 0 .4 6 60 26940121 2 6 .9 4 5 .4 7 70 137512567 1 3 7 .5 1 2 7 .9 4 75 5368029 5 .3 7 1 .0 9 80 284000170 2 8 4 .0 0 5 7 .7 1 85 479061 0 .4 8 0 .1 0 90 56554 0 .0 6 0 .0 1 95 36020 0 .0 4 0 .0 1 T o ta l S o i l A r e a ( sq k m ) 4 9 2 .0 9 % o f T o ta l B a si n A r e a 4 0 .6 3 T o ta l B a si n A r e a ( sq k m ) 1 2 1 1 .1 8 52 V i l l a l b a S u b b a s i n N o . 6 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 198132641 1 9 8 .1 3 8 .5 6 20 413894292 4 1 3 .8 9 1 7 .8 7 30 187225 0 .1 9 0 .0 1 40 8933646 8 .9 3 0 .3 9 50 533238747 5 3 3 .2 4 2 3 .0 3 60 98264833 9 8 .2 6 4 .2 4 70 107702593 1 0 7 .7 0 4 .6 5 75 0 0 .0 0 0 .0 0 80 952504833 9 5 2 .5 0 4 1 .1 3 85 1555012 1 .5 6 0 .0 7 90 0 0 .0 0 0 .0 0 95 1228024 1 .2 3 0 .0 5 T o ta l S o i l A r e a ( sq k m ) 2 3 1 5 .6 4 % o f T o ta l B a si n A r e a 2 4 .7 8 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 125748168 1 2 5 .7 5 4 .6 5 20 161637049 1 6 1 .6 4 5 .9 8 30 0 0 .0 0 0 .0 0 40 5899623 5 .9 0 0 .2 2 50 555046757 5 5 5 .0 5 2 0 .5 4 60 182512530 1 8 2 .5 1 6 .7 5 70 291867882 2 9 1 .8 7 1 0 .8 0 75 375438 0 .3 8 0 .0 1 80 1375486854 1 3 7 5 .4 9 5 0 .9 0 85 439389 0 .4 4 0 .0 2 90 381 0 .0 0 0 .0 0 95 3179484 3 .1 8 0 .1 2 T o ta l S o i l A r e a ( sq k m ) 2 7 0 2 .1 9 % o f T o ta l B a si n A r e a 2 8 .9 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 362443991 3 6 2 .4 4 1 2 .5 6 20 443870114 4 4 3 .8 7 1 5 .3 8 30 7201 0 .0 1 0 .0 0 40 7068322 7 .0 7 0 .2 4 50 269218659 2 6 9 .2 2 9 .3 3 60 220214230 2 2 0 .2 1 7 .6 3 70 126818601 1 2 6 .8 2 4 .3 9 75 0 0 .0 0 0 .0 0 80 1452923855 1 4 5 2 .9 2 5 0 .3 3 85 3010849 3 .0 1 0 .1 0 90 3221 0 .0 0 0 .0 0 95 965096 0 .9 7 0 .0 3 T o ta l S o i l A r e a ( sq k m ) 2 8 8 6 .5 4 % o f T o ta l B a si n A r e a 3 0 .8 9 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 335729645 3 3 5 .7 3 2 3 .3 1 20 370882733 3 7 0 .8 8 2 5 .7 5 30 0 0 .0 0 0 .0 0 40 21603 0 .0 2 0 .0 0 50 54475310 5 4 .4 8 3 .7 8 60 28195313 2 8 .2 0 1 .9 6 70 156117713 1 5 6 .1 2 1 0 .8 4 75 0 0 .0 0 0 .0 0 80 492312880 4 9 2 .3 1 3 4 .1 8 85 2211562 2 .2 1 0 .1 5 90 10906 0 .0 1 0 .0 0 95 317715 0 .3 2 0 .0 2 T o ta l S o i l A r e a ( sq k m ) 1 4 4 0 .2 8 % o f T o ta l B a si n A r e a 1 5 .4 1 T o ta l B a si n A r e a ( sq k m ) 9 3 4 4 .6 5 53 C o n c h o s S u b b a s i n N o . 7 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 0 .0 0 20 50436 0 .0 5 0 .0 2 30 0 0 .0 0 0 .0 0 40 408759 0 .4 1 0 .1 8 50 309305 0 .3 1 0 .1 4 60 83971705 8 3 .9 7 3 6 .9 8 70 35919750 3 5 .9 2 1 5 .8 2 75 43267506 4 3 .2 7 1 9 .0 6 80 57232239 5 7 .2 3 2 5 .2 1 85 0 0 .0 0 0 .0 0 90 158360 0 .1 6 0 .0 7 95 5743406 5 .7 4 2 .5 3 T o ta l S o i l A r e a ( sq k m ) 2 2 7 .0 6 % o f T o ta l B a si n A r e a 2 0 .3 8 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 110478 0 .1 1 0 .0 2 40 43823515 4 3 .8 2 8 .9 5 50 17907504 1 7 .9 1 3 .6 6 60 237194242 2 3 7 .1 9 4 8 .4 7 70 70536060 7 0 .5 4 1 4 .4 1 75 57198139 5 7 .2 0 1 1 .6 9 80 47225575 4 7 .2 3 9 .6 5 85 2708745 2 .7 1 0 .5 5 90 12673210 1 2 .6 7 2 .5 9 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 4 8 9 .3 8 % o f T o ta l B a si n A r e a 4 3 .9 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 347434 0 .3 5 0 .1 2 50 1253316 1 .2 5 0 .4 4 60 72240662 7 2 .2 4 2 5 .4 1 70 6163864 6 .1 6 2 .1 7 75 69507528 6 9 .5 1 2 4 .4 5 80 124531733 1 2 4 .5 3 4 3 .8 0 85 0 0 .0 0 0 .0 0 90 8414032 8 .4 1 2 .9 6 95 1852499 1 .8 5 0 .6 5 T o ta l S o i l A r e a ( sq k m ) 2 8 4 .3 1 % o f T o ta l B a si n A r e a 2 5 .5 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 0 0 .0 0 0 .0 0 20 36025 0 .0 4 0 .0 3 30 0 0 .0 0 0 .0 0 40 9770 0 .0 1 0 .0 1 50 133881 0 .1 3 0 .1 2 60 23126719 2 3 .1 3 2 0 .3 8 70 31533083 3 1 .5 3 2 7 .7 9 75 3644145 3 .6 4 3 .2 1 80 52251006 5 2 .2 5 4 6 .0 6 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 2717052 2 .7 2 2 .3 9 T o ta l S o i l A r e a ( sq k m ) 1 1 3 .4 5 % o f T o ta l B a si n A r e a 1 0 .1 8 T o ta l B a si n A r e a ( sq k m ) 1 1 1 4 .2 0 54 J i m e n e z S u b b a s i n N o . 8 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 34913902 3 4 .9 1 3 .7 6 20 28441541 2 8 .4 4 3 .0 7 30 1389805 1 .3 9 0 .1 5 40 4202684 4 .2 0 0 .4 5 50 97692805 9 7 .6 9 1 0 .5 3 60 84047669 8 4 .0 5 9 .0 6 70 1924469 1 .9 2 0 .2 1 75 0 0 .0 0 0 .0 0 80 666241918 6 6 6 .2 4 7 1 .8 4 85 5153187 5 .1 5 0 .5 6 90 3341768 3 .3 4 0 .3 6 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 9 2 7 .3 5 % o f T o ta l B a si n A r e a 2 0 .9 8 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 23815323 2 3 .8 2 6 .2 1 20 48552171 4 8 .5 5 1 2 .6 7 30 0 0 .0 0 0 .0 0 40 3189980 3 .1 9 0 .8 3 50 29959563 2 9 .9 6 7 .8 2 60 43308066 4 3 .3 1 1 1 .3 0 70 21440609 2 1 .4 4 5 .5 9 75 0 0 .0 0 0 .0 0 80 212842931 2 1 2 .8 4 5 5 .5 4 85 113220 0 .1 1 0 .0 3 90 0 0 .0 0 0 .0 0 95 32421 0 .0 3 0 .0 1 T o ta l S o i l A r e a ( sq k m ) 3 8 3 .2 5 % o f T o ta l B a si n A r e a 8 .6 7 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 15956187 1 5 .9 6 0 .6 0 20 31508941 3 1 .5 1 1 .1 8 30 1505934 1 .5 1 0 .0 6 40 63100775 6 3 .1 0 2 .3 5 50 464100157 4 6 4 .1 0 1 7 .3 2 60 558066867 5 5 8 .0 7 2 0 .8 2 70 191146341 1 9 1 .1 5 7 .1 3 75 20821027 2 0 .8 2 0 .7 8 80 1317856274 1 3 1 7 .8 6 4 9 .1 7 85 2615272 2 .6 2 0 .1 0 90 13393245 1 3 .3 9 0 .5 0 95 3602 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 6 8 0 .0 7 % o f T o ta l B a si n A r e a 6 0 .6 2 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 48816092 4 8 .8 2 1 1 .3 4 20 33443498 3 3 .4 4 7 .7 7 30 0 0 .0 0 0 .0 0 40 78064788 7 8 .0 6 1 8 .1 4 50 68821328 6 8 .8 2 1 5 .9 9 60 54500399 5 4 .5 0 1 2 .6 6 70 11771107 1 1 .7 7 2 .7 3 75 0 0 .0 0 0 .0 0 80 128987283 1 2 8 .9 9 2 9 .9 7 85 1247547 1 .2 5 0 .2 9 90 4748616 4 .7 5 1 .1 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 4 3 0 .4 0 % o f T o ta l B a si n A r e a 9 .7 4 T o ta l B a si n A r e a ( sq k m ) 4 4 2 1 .0 8 55 C h u v i s c a r S u b b a s i n N o . 9 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 64830 0 .0 6 0 .2 0 20 0 0 .0 0 0 .0 0 30 374542 0 .3 7 1 .1 4 40 0 0 .0 0 0 .0 0 50 4194101 4 .1 9 1 2 .7 5 60 5450720 5 .4 5 1 6 .5 6 70 10403415 1 0 .4 0 3 1 .6 1 75 0 0 .0 0 0 .0 0 80 6810254 6 .8 1 2 0 .7 0 85 5608765 5 .6 1 1 7 .0 4 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 3 2 .9 1 % o f T o ta l B a si n A r e a 3 1 .6 9 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 0 0 .0 0 0 .0 0 20 1085320 1 .0 9 1 .5 3 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 362367 0 .3 6 0 .5 1 60 6543030 6 .5 4 9 .2 3 70 48892701 4 8 .8 9 6 8 .9 4 75 0 0 .0 0 0 .0 0 80 6040955 6 .0 4 8 .5 2 85 7997881 8 .0 0 1 1 .2 8 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 7 0 .9 2 % o f T o ta l B a si n A r e a 6 8 .3 1 T o ta l B a si n A r e a ( sq k m ) 1 0 3 .8 3 56 E l R e j o n S u b b a s i n N o . 1 0 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 6387713 6 .3 9 3 8 .3 9 60 3037632 3 .0 4 1 8 .2 6 70 3075528 3 .0 8 1 8 .4 8 75 0 0 .0 0 0 .0 0 80 4138956 4 .1 4 2 4 .8 7 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 6 .6 4 % o f T o ta l B a si n A r e a 1 1 .3 4 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 6056313 6 .0 6 4 .6 5 20 41631691 4 1 .6 3 3 1 .9 9 30 0 0 .0 0 0 .0 0 40 28813 0 .0 3 0 .0 2 50 2486618 2 .4 9 1 .9 1 60 4600997 4 .6 0 3 .5 4 70 11264658 1 1 .2 6 8 .6 6 75 0 0 .0 0 0 .0 0 80 64079278 6 4 .0 8 4 9 .2 4 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 3 0 .1 5 % o f T o ta l B a si n A r e a 8 8 .6 6 T o ta l B a si n A r e a ( sq k m ) 1 4 6 .7 9 57 C h i h u a h u a S u b b a s i n N o . 1 1 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 5980028 5 .9 8 2 0 .0 5 20 2063619 2 .0 6 6 .9 2 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 15126355 1 5 .1 3 5 0 .7 1 60 234106 0 .2 3 0 .7 8 70 314262 0 .3 1 1 .0 5 75 0 0 .0 0 0 .0 0 80 6111300 6 .1 1 2 0 .4 9 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 9 .8 3 % o f T o ta l B a si n A r e a 7 .4 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 241982 0 .2 4 3 2 .5 7 20 25687 0 .0 3 3 .4 6 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 300345 0 .3 0 4 0 .4 2 60 0 0 .0 0 0 .0 0 70 51689 0 .0 5 6 .9 6 75 0 0 .0 0 0 .0 0 80 123307 0 .1 2 1 6 .6 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 0 .7 4 % o f T o ta l B a si n A r e a 0 .1 9 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 43220 0 .0 4 0 .0 6 20 1825414 1 .8 3 2 .3 3 30 551099 0 .5 5 0 .7 0 40 0 0 .0 0 0 .0 0 50 11020715 1 1 .0 2 1 4 .0 6 60 3614801 3 .6 1 4 .6 1 70 32032113 3 2 .0 3 4 0 .8 5 75 0 0 .0 0 0 .0 0 80 29321473 2 9 .3 2 3 7 .4 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 7 8 .4 1 % o f T o ta l B a si n A r e a 1 9 .6 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 25436490 2 5 .4 4 8 .7 5 20 79428062 7 9 .4 3 2 7 .3 1 30 0 0 .0 0 0 .0 0 40 32412 0 .0 3 0 .0 1 50 5477565 5 .4 8 1 .8 8 60 1361756 1 .3 6 0 .4 7 70 44840350 4 4 .8 4 1 5 .4 2 75 0 0 .0 0 0 .0 0 80 134267868 1 3 4 .2 7 4 6 .1 6 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 18008 0 .0 2 0 .0 1 T o ta l S o i l A r e a ( sq k m ) 2 9 0 .8 6 % o f T o ta l B a si n A r e a 7 2 .7 4 T o ta l B a si n A r e a ( sq k m ) 3 9 9 .8 4 58 L l a n i t o s S u b b a s i n N o . 1 2 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 207646613 2 0 7 .6 5 6 6 .4 6 20 46749497 4 6 .7 5 1 4 .9 6 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 622991 0 .6 2 0 .2 0 60 567584 0 .5 7 0 .1 8 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 55468865 5 5 .4 7 1 7 .7 5 85 1378106 1 .3 8 0 .4 4 90 0 0 .0 0 0 .0 0 95 7203 0 .0 1 0 .0 0 T o ta l S o i l A r e a ( sq km ) 3 1 2 .4 4 % o f T o ta l B a si n A r e a 1 8 .8 5 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 424400774 4 2 4 .4 0 4 7 .5 9 20 289603320 2 8 9 .6 0 3 2 .4 7 30 0 0 .0 0 0 .0 0 40 115251 0 .1 2 0 .0 1 50 4628554 4 .6 3 0 .5 2 60 9684430 9 .6 8 1 .0 9 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 162060957 1 6 2 .0 6 1 8 .1 7 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 1311011 1 .3 1 0 .1 5 T o ta l S o i l A r e a ( sq km ) 8 9 1 .8 0 % o f T o ta l B a si n A r e a 5 3 .8 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 6519778 6 .5 2 8 .6 8 20 15990376 1 5 .9 9 2 1 .2 9 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 9375068 9 .3 8 1 2 .4 8 60 1530698 1 .5 3 2 .0 4 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 41693464 4 1 .6 9 5 5 .5 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 14406 0 .0 1 0 .0 2 T o ta l S o i l A r e a ( sq km ) 7 5 .1 2 % o f T o ta l B a si n A r e a 4 .5 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 260274965 2 6 0 .2 7 6 8 .8 5 20 71044864 7 1 .0 4 1 8 .7 9 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 3075931 3 .0 8 0 .8 1 60 375746 0 .3 8 0 .1 0 70 0 0 .0 0 0 .0 0 75 64830 0 .0 6 0 .0 2 80 41570696 4 1 .5 7 1 1 .0 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 1615693 1 .6 2 0 .4 3 T o ta l S o i l A r e a ( sq km ) 3 7 8 .0 2 % o f T o ta l B a si n A r e a 2 2 .8 1 T o ta l B a si n A r e a ( sq km ) 1 6 5 7 .3 9 59 P i c o d e l A g u i l a S u b b a s i n N o . 1 3 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 563035 0 .5 6 0 .2 9 20 3633432 3 .6 3 1 .9 0 30 3499174 3 .5 0 1 .8 3 40 248823 0 .2 5 0 .1 3 50 3411540 3 .4 1 1 .7 9 60 6280169 6 .2 8 3 .2 9 70 10814 0 .0 1 0 .0 1 75 0 0 .0 0 0 .0 0 80 173088051 1 7 3 .0 9 9 0 .5 7 85 0 0 .0 0 0 .0 0 90 365534 0 .3 7 0 .1 9 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 9 1 .1 0 % o f T o ta l B a si n A r e a 2 9 .6 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 1943288 1 .9 4 1 .0 7 20 4702495 4 .7 0 2 .5 9 30 0 0 .0 0 0 .0 0 40 1456622 1 .4 6 0 .8 0 50 30130067 3 0 .1 3 1 6 .6 1 60 22652882 2 2 .6 5 1 2 .4 9 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 118993787 1 1 8 .9 9 6 5 .6 1 85 555938 0 .5 6 0 .3 1 90 936755 0 .9 4 0 .5 2 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 8 1 .3 7 % o f T o ta l B a si n A r e a 2 8 .1 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 35215145 3 5 .2 2 1 3 .9 1 20 11339547 1 1 .3 4 4 .4 8 30 183745 0 .1 8 0 .0 7 40 270212 0 .2 7 0 .1 1 50 20267463 2 0 .2 7 8 .0 0 60 8301945 8 .3 0 3 .2 8 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 177661692 1 7 7 .6 6 7 0 .1 6 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 5 3 .2 4 % o f T o ta l B a si n A r e a 3 9 .2 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 0 0 .0 0 0 .0 0 20 525530 0 .5 3 2 .6 5 30 0 0 .0 0 0 .0 0 40 1396 0 .0 0 0 .0 1 50 2943179 2 .9 4 1 4 .8 4 60 3456129 3 .4 6 1 7 .4 3 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 12888235 1 2 .8 9 6 5 .0 0 85 13305 0 .0 1 0 .0 7 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 9 .8 3 % o f T o ta l B a si n A r e a 3 .0 7 T o ta l B a si n A r e a ( sq k m ) 6 4 5 .5 4 60 S a n A n t o n i o S u b b a s i n N o . 1 4 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 12245098 1 2 .2 5 3 .3 7 20 10223355 1 0 .2 2 2 .8 1 30 0 0 .0 0 0 .0 0 40 6803211 6 .8 0 1 .8 7 50 70427757 7 0 .4 3 1 9 .3 6 60 93638097 9 3 .6 4 2 5 .7 3 70 24514 0 .0 2 0 .0 1 75 0 0 .0 0 0 .0 0 80 168993744 1 6 8 .9 9 4 6 .4 4 85 529094 0 .5 3 0 .1 5 90 974165 0 .9 7 0 .2 7 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 3 6 3 .8 6 % o f T o ta l B a si n A r e a 4 4 .5 8 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 1006218 1 .0 1 5 .1 5 20 1137574 1 .1 4 5 .8 2 30 0 0 .0 0 0 .0 0 40 59303 0 .0 6 0 .3 0 50 601006 0 .6 0 3 .0 8 60 3761115 3 .7 6 1 9 .2 5 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 12589580 1 2 .5 9 6 4 .4 5 85 0 0 .0 0 0 .0 0 90 379312 0 .3 8 1 .9 4 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 9 .5 3 % o f T o ta l B a si n A r e a 2 .3 9 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 9386 0 .0 1 0 .0 1 20 7938507 7 .9 4 4 .5 5 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 66525125 6 6 .5 3 3 8 .1 4 60 14968165 1 4 .9 7 8 .5 8 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 84987486 8 4 .9 9 4 8 .7 2 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 7 4 .4 3 % o f T o ta l B a si n A r e a 2 1 .3 7 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 7055845 7 .0 6 2 .7 3 20 34584464 3 4 .5 8 1 3 .3 8 30 0 0 .0 0 0 .0 0 40 5054908 5 .0 5 1 .9 6 50 23854670 2 3 .8 5 9 .2 3 60 30056910 3 0 .0 6 1 1 .6 3 70 225807 0 .2 3 0 .0 9 75 0 0 .0 0 0 .0 0 80 151562386 1 5 1 .5 6 5 8 .6 5 85 2558464 2 .5 6 0 .9 9 90 3472480 3 .4 7 1 .3 4 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 5 8 .4 3 % o f T o ta l B a si n A r e a 3 1 .6 6 T o ta l B a si n A r e a ( sq k m ) 8 1 6 .2 5 61 S a n G a b r i e l S u b b a s i n N o . 1 5 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 9792047 9 .7 9 6 .3 1 20 2553014 2 .5 5 1 .6 4 30 2113306 2 .1 1 1 .3 6 40 2898077 2 .9 0 1 .8 7 50 23577971 2 3 .5 8 1 5 .1 8 60 49325485 4 9 .3 3 3 1 .7 6 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 64658424 6 4 .6 6 4 1 .6 4 85 0 0 .0 0 0 .0 0 90 376775 0 .3 8 0 .2 4 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 5 5 .3 0 % o f T o ta l B a si n A r e a 5 0 .9 7 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 16709344 1 6 .7 1 2 9 .9 0 20 7096729 7 .1 0 1 2 .7 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 268279 0 .2 7 0 .4 8 60 1166782 1 .1 7 2 .0 9 70 41548 0 .0 4 0 .0 7 75 203977 0 .2 0 0 .3 7 80 30392781 3 0 .3 9 5 4 .3 9 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 5 5 .8 8 % o f T o ta l B a si n A r e a 1 8 .3 4 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 7575164 7 .5 8 8 .1 0 20 33131067 3 3 .1 3 3 5 .4 3 30 7361479 7 .3 6 7 .8 7 40 2260 0 .0 0 0 .0 0 50 3276839 3 .2 8 3 .5 0 60 7467932 7 .4 7 7 .9 9 70 650342 0 .6 5 0 .7 0 75 0 0 .0 0 0 .0 0 80 34040505 3 4 .0 4 3 6 .4 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 9 3 .5 1 % o f T o ta l B a si n A r e a 3 0 .6 9 T o ta l B a si n A r e a ( sq k m ) 3 0 4 .6 8 62 P u e n t e F F C C S u b b a s i n N o . 1 6 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 13674219 1 3 .6 7 9 .3 6 20 5234650 5 .2 3 3 .5 8 30 0 0 .0 0 0 .0 0 40 1159997 1 .1 6 0 .7 9 50 53005090 5 3 .0 1 3 6 .2 8 60 15190917 1 5 .1 9 1 0 .4 0 70 889732 0 .8 9 0 .6 1 75 0 0 .0 0 0 .0 0 80 55854892 5 5 .8 5 3 8 .2 3 85 1018642 1 .0 2 0 .7 0 90 64841 0 .0 6 0 .0 4 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 4 6 .0 9 % o f T o ta l B a si n A r e a 1 2 .2 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 0 0 .0 0 0 .0 0 60 0 0 .0 0 0 .0 0 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 366379 0 .3 7 1 0 0 .0 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 0 .3 7 % o f T o ta l B a si n A r e a 0 .0 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 91085210 9 1 .0 9 2 1 .3 0 20 133778292 1 3 3 .7 8 3 1 .2 9 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 60150210 6 0 .1 5 1 4 .0 7 60 8993159 8 .9 9 2 .1 0 70 204982 0 .2 0 0 .0 5 75 558517 0 .5 6 0 .1 3 80 132792471 1 3 2 .7 9 3 1 .0 6 85 10807 0 .0 1 0 .0 0 90 0 0 .0 0 0 .0 0 95 3552 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 4 2 7 .5 8 % o f T o ta l B a si n A r e a 3 5 .8 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 341348993 3 4 1 .3 5 5 5 .0 3 20 206538030 2 0 6 .5 4 3 3 .3 0 30 0 0 .0 0 0 .0 0 40 37 0 .0 0 0 .0 0 50 3998262 4 .0 0 0 .6 4 60 2448281 2 .4 5 0 .3 9 70 2417668 2 .4 2 0 .3 9 75 21612 0 .0 2 0 .0 0 80 60283452 6 0 .2 8 9 .7 2 85 231431 0 .2 3 0 .0 4 90 0 0 .0 0 0 .0 0 95 3004692 3 .0 0 0 .4 8 T o ta l S o i l A r e a ( sq k m ) 6 2 0 .2 9 % o f T o ta l B a si n A r e a 5 1 .9 4 T o ta l B a si n A r e a ( sq k m ) 1 1 9 4 .3 3 63 P a r r a l S u b b a s i n N o . 1 7 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 516654 0 .5 2 0 .6 7 20 4846086 4 .8 5 6 .2 7 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 19084748 1 9 .0 8 2 4 .6 8 60 2770671 2 .7 7 3 .5 8 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 49733890 4 9 .7 3 6 4 .3 3 85 364463 0 .3 6 0 .4 7 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 7 7 .3 2 2 1 .2 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 18939225 1 8 .9 4 1 1 .5 7 20 38996986 3 9 .0 0 2 3 .8 2 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 4084334 4 .0 8 2 .4 9 60 2998509 3 .0 0 1 .8 3 70 319330 0 .3 2 0 .2 0 75 0 0 .0 0 0 .0 0 80 95929175 9 5 .9 3 5 8 .6 0 85 2445131 2 .4 5 1 .4 9 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 6 3 .7 1 % o f T o ta l B a si n A r e a 4 5 .0 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 0 0 .0 0 0 .0 0 20 356330 0 .3 6 5 7 .6 4 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 0 0 .0 0 0 .0 0 60 0 0 .0 0 0 .0 0 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 261833 0 .2 6 4 2 .3 6 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 0 .6 2 % o f T o ta l B a si n A r e a 0 .1 7 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 35647474 3 5 .6 5 2 9 .1 9 20 7732552 7 .7 3 6 .3 3 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 19077829 1 9 .0 8 1 5 .6 2 60 3005443 3 .0 1 2 .4 6 70 280 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 56054692 5 6 .0 5 4 5 .9 1 85 587348 0 .5 9 0 .4 8 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 1 2 2 .1 1 % o f T o ta l B a si n A r e a 3 3 .5 7 T o ta l B a si n A r e a ( sq k m ) 3 6 3 .7 5 64 C o l i n a S u b b a s i n N o . 1 8 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 0 0 .0 0 0 .0 0 60 47971579 4 7 .9 7 3 9 .9 3 70 20500895 2 0 .5 0 1 7 .0 6 75 20308646 2 0 .3 1 1 6 .9 0 80 26849590 2 6 .8 5 2 2 .3 5 85 0 0 .0 0 0 .0 0 90 1253886 1 .2 5 1 .0 4 95 3249572 3 .2 5 2 .7 0 T o ta l S o i l A r e a ( sq k m ) 1 2 0 .1 3 % o f T o ta l B a si n A r e a 4 6 .9 4 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 121537 0 .1 2 0 .2 5 40 2379373 2 .3 8 4 .9 7 50 615402 0 .6 2 1 .2 9 60 15650886 1 5 .6 5 3 2 .6 9 70 23176041 2 3 .1 8 4 8 .4 1 75 0 0 .0 0 0 .0 0 80 5893782 5 .8 9 1 2 .3 1 85 0 0 .0 0 0 .0 0 90 36028 0 .0 4 0 .0 8 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 4 7 .8 7 % o f T o ta l B a si n A r e a 1 8 .7 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 258856 0 .2 6 0 .4 0 20 0 0 .0 0 0 .0 0 30 5161004 5 .1 6 8 .0 6 40 4142957 4 .1 4 6 .4 7 50 3509971 3 .5 1 5 .4 8 60 25459693 2 5 .4 6 3 9 .7 8 70 4716003 4 .7 2 7 .3 7 75 0 0 .0 0 0 .0 0 80 20750688 2 0 .7 5 3 2 .4 2 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 6 4 .0 0 % o f T o ta l B a si n A r e a 2 5 .0 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 274349 0 .2 7 1 .1 5 20 0 0 .0 0 0 .0 0 30 2650974 2 .6 5 1 1 .0 8 40 2098721 2 .1 0 8 .7 7 50 1447651 1 .4 5 6 .0 5 60 7382191 7 .3 8 3 0 .8 6 70 9806207 9 .8 1 4 0 .9 9 75 0 0 .0 0 0 .0 0 80 261624 0 .2 6 1 .0 9 85 0 0 .0 0 0 .0 0 90 3603 0 .0 0 0 .0 2 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 2 3 .9 3 % o f T o ta l B a si n A r e a 9 .3 5 T o ta l B a si n A r e a ( sq k m ) 2 5 5 .9 3 65 L a B o q u i l l a S u b b a s i n N o . 1 9 S o i l - L a n d U s e I n t e r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 1170628229 1 1 7 0 .6 3 2 4 .8 6 20 746690894 7 4 6 .6 9 1 5 .8 6 30 400181 0 .4 0 0 .0 1 40 10911772 1 0 .9 1 0 .2 3 50 619673766 6 1 9 .6 7 1 3 .1 6 60 422781624 4 2 2 .7 8 8 .9 8 70 85146360 8 5 .1 5 1 .8 1 75 317979 0 .3 2 0 .0 1 80 1633173132 1 6 3 3 .1 7 3 4 .6 8 85 1313984 1 .3 1 0 .0 3 90 50402 0 .0 5 0 .0 0 95 17931306 1 7 .9 3 0 .3 8 T o ta l S o i l A r e a ( sq k m ) 4 7 0 9 .0 2 % o f T o ta l B a si n A r e a 2 6 .5 6 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 2076083278 2 0 7 6 .0 8 3 0 .1 8 20 1610168768 1 6 1 0 .1 7 2 3 .4 1 30 0 0 .0 0 0 .0 0 40 9992579 9 .9 9 0 .1 5 50 497251893 4 9 7 .2 5 7 .2 3 60 257303211 2 5 7 .3 0 3 .7 4 70 42635586 4 2 .6 4 0 .6 2 75 523486 0 .5 2 0 .0 1 80 2380528535 2 3 8 0 .5 3 3 4 .6 0 85 1871345 1 .8 7 0 .0 3 90 1260822 1 .2 6 0 .0 2 95 1808601 1 .8 1 0 .0 3 T o ta l S o i l A r e a ( sq k m ) 6 8 7 9 .4 3 % o f T o ta l B a si n A r e a 3 8 .8 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 494974896 4 9 4 .9 7 1 2 .7 9 20 471496196 4 7 1 .5 0 1 2 .1 8 30 17042455 1 7 .0 4 0 .4 4 40 4733926 4 .7 3 0 .1 2 50 343967487 3 4 3 .9 7 8 .8 9 60 446845917 4 4 6 .8 5 1 1 .5 4 70 63865488 6 3 .8 7 1 .6 5 75 1594508 1 .5 9 0 .0 4 80 2005201600 2 0 0 5 .2 0 5 1 .8 0 85 5061537 5 .0 6 0 .1 3 90 0 0 .0 0 0 .0 0 95 16256688 1 6 .2 6 0 .4 2 T o ta l S o i l A r e a ( sq k m ) 3 8 7 1 .0 4 % o f T o ta l B a si n A r e a 2 1 .8 3 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 343710053 3 4 3 .7 1 1 5 .1 3 20 395442448 3 9 5 .4 4 1 7 .4 1 30 32836869 3 2 .8 4 1 .4 5 40 3428887 3 .4 3 0 .1 5 50 108937067 1 0 8 .9 4 4 .8 0 60 392471158 3 9 2 .4 7 1 7 .2 8 70 70823994 7 0 .8 2 3 .1 2 75 2142447 2 .1 4 0 .0 9 80 909521959 9 0 9 .5 2 4 0 .0 4 85 284184 0 .2 8 0 .0 1 90 57591 0 .0 6 0 .0 0 95 11918114 1 1 .9 2 0 .5 2 T o ta l S o i l A r e a ( sq k m ) 2 2 7 1 .5 7 % o f T o ta l B a si n A r e a 1 2 .8 1 T o ta l B a si n A r e a ( sq k m ) 1 7 7 3 1 .0 6 66 O j i n a g a Su b b a s i n N o . 2 0 So i l - L a n d U s e I n te r s e c t S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p A 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p B 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 0 0 .0 0 0 .0 0 60 16361 0 .0 2 5 9 .0 7 70 11338 0 .0 1 4 0 .9 3 75 0 0 .0 0 0 .0 0 80 0 0 .0 0 0 .0 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 0 .0 3 % o f T o ta l B a si n A r e a 6 0 .0 1 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p C 10 0 0 .0 0 # DIV /0 ! 20 0 0 .0 0 # DIV /0 ! 30 0 0 .0 0 # DIV /0 ! 40 0 0 .0 0 # DIV /0 ! 50 0 0 .0 0 # DIV /0 ! 60 0 0 .0 0 # DIV /0 ! 70 0 0 .0 0 # DIV /0 ! 75 0 0 .0 0 # DIV /0 ! 80 0 0 .0 0 # DIV /0 ! 85 0 0 .0 0 # DIV /0 ! 90 0 0 .0 0 # DIV /0 ! 95 0 0 .0 0 # DIV /0 ! T o ta l S o i l A r e a ( sq k m ) 0 .0 0 % o f T o ta l B a si n A r e a 0 .0 0 S o il G r o u p L a n d Us e A r e a ( s q m ) A r e a ( s q k m ) % o f S o il G r o u p D 10 0 0 .0 0 0 .0 0 20 0 0 .0 0 0 .0 0 30 0 0 .0 0 0 .0 0 40 0 0 .0 0 0 .0 0 50 0 0 .0 0 0 .0 0 60 18458 0 .0 2 1 0 0 .0 0 70 0 0 .0 0 0 .0 0 75 0 0 .0 0 0 .0 0 80 0 0 .0 0 0 .0 0 85 0 0 .0 0 0 .0 0 90 0 0 .0 0 0 .0 0 95 0 0 .0 0 0 .0 0 T o ta l S o i l A r e a ( sq k m ) 0 .0 2 % o f T o ta l B a si n A r e a 3 9 .9 9 T o ta l B a si n A r e a ( sq k m ) 0 .0 5 67 Appendix 2. CRWR Geodatabase Reach Lengths Basin Name IMTA Basin ID Reach Name Reach Length Notes (km) Peguis 1 Rio_Conchos_1 142.32 DS Reach Rio_Conchos_20 Sacramento 2 - - Headwater to Rio_Sacramento_4 Las Burras 3 Rio_Conchos_3 135.41 DS Reach Rio_Conchos_4 Rio_Florido_3 89.05 DS Reach Rio_Conchos_3 Arroyo_El_Parral_3 148.77 DS Reach Rio_Florido_3 Rio_San_Pedro_3 38.65 DS Reach Rio_Conchos_3 Luis L. Leon 4 Rio_Conchos_4 88.74 DS Reach Rio_Conchos_1 Rio_Sacramento_4 103.49 DS Reach Rio_Conchos_4 Arroyo_Sacramento_4 5.55 DS Reach Rio_Sacramento_4 FCO. Madero 5 Rio_San_Pedro_5 29.15 Headwater to Rio_San_Pedro_3 Villalba 6 - - Headwater to Rio_San_Pedro_5 Conchos 7 Rio_Conchos_7 25.94 DS Reach: Rio_Conchos_3 Jimenez 8 Rio_Florido_8 88.08 DS Reach: Rio_FLorido_3 Chuviscar 9 Arroyo_Sacramento_9 140.41 DS Reach: Arroyo_Sacramento_4 River_10 3.12 DS Reach Arroyo_Sacramento_9 El Rejon 10 - - Headwater to River_10 Chihuahua 11 - - Headwater to Arroyo_Sacramento_9 Llanitos 12 - - Headwater to RioBalleza_Conchos_19 Pico del Aguila 13 Rio_Florido_13 24.05 DS Reach Rio_Florido_8 San Antonio 14 Rio_Florido_14 28.33 DS Reach Rio_Florido_13 San Gabriel 15 Rio_Florido_15 12.54 DS Reach Rio_Florido_14 Puente FFCC 16 - - Headwater to Rio_Florido_15 Parral 17 - - Headwater to Arroyo_El_Parral_3 Colina 18 Rio_Conchos_18 18210.59 DS Reach Rio_Conchos_7 La Boquilla 19 RioBalleza_Conchos_19 259.86 DS Reach Rio_Conchos_18 Ojinaga 20 Rio_Conchos_20 45.90 Outfall