A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics. |
The U.S. Geological Survey, in cooperation with the Texas Water Development Board, evaluated historical channel adjustment and estimated selected hydraulic values at U.S. Geological Survey streamflow-gaging stations in the lower Sabine River Basin in Texas and Louisiana and lower Brazos River Basin in Texas to support geomorphic assessments of the Texas Instream Flow Program. Channel attributes including cross-section geometry, slope, and planform change were evaluated to learn how each river's morphology changed over the years in response to natural and anthropogenic disturbances. Historical and contemporary cross-sectional channel geometries at several gaging stations on each river were compared, planform changes were assessed, and hydraulic values were estimated including mean flow velocity, bed shear stress, Froude numbers, and hydraulic depth. The primary sources of historical channel morphology information were U.S. Geological Survey hard-copy discharge-measurement field notes. Additional analyses were done using computations of selected flow hydraulics, comparisons of historical and contemporary aerial photographs, comparisons of historical and contemporary ground photographs, evaluations of how frequently stage-discharge rating curves were updated, reviews of stage-discharge relations for field measurements, and considerations of bridge and reservoir construction activities. Based on historical cross sections at three gaging stations downstream from Toledo Bend Reservoir, the lower Sabine River is relatively stable, but is subject to substantial temporary scour-and-fill processes during floods. Exceptions to this characterization of relative stability include an episode of channel aggradation at the Sabine River near Bon Wier, Texas, during the 1930s, and about 2 to 3 feet of channel incision at the Sabine River near Burkeville, Texas, since the late 1950s. The Brazos River, at gaging stations downstream from Waco, Texas, has adjusted to a combination of hydrologic, sedimentary, and anthropogenic controls. Since the 1960s, numerous point bars have vertically accreted and vegetation has encroached along the channel margins, which probably promotes channel-bed incision to compensate for a reduction in cross-sectional area. Channel incision was detected at all gaging stations along the Brazos River, and the depth of incision is greatest in the lowermost gaging stations, exemplified by about 5 feet of channel-bed incision between 1993 and 2004 at Richmond, Texas. One notable exception to this pattern of incision was a period of aggradation at U.S. Geological Survey gaging station 08096500 Brazos River at Waco, Texas, during the late 1920s and 1930s, probably associated with upstream dam construction. Lateral channel migration rates along the Brazos River determined from aerial photographs are greatest between Waco and Hempstead, Texas, with numerous bends moving an average of more than 10 feet per year. Migration rates at selected bends downstream from Hempstead were measured as less than 10 feet per year, on average. Two tributaries of the Brazos River, the Little and Navasota Rivers, also were investigated for historical channel adjustment. The Little River near Cameron, Texas (08106500) has incised its channel bed about 12 feet since 1949, and the lower Navasota River shows complex adjustment to bridge construction activities and a channel avulsion.
","language":"English","publisher":"U.S. Geological Survey","publisherLocation":"Reston, VA","doi":"10.3133/sir20095174","collaboration":"Prepared in cooperation with the Texas Water Development Board","usgsCitation":"Heitmuller, F.T., and Greene, L.E., 2009, Historical Channel Adjustment and Estimates of Selected Hydraulic Values in the Lower Sabine River and Lower Brazos River Basins, Texas and Louisiana: U.S. Geological Survey Scientific Investigations Report 2009-5174, xii, 144 p., https://doi.org/10.3133/sir20095174.","productDescription":"xii, 144 p.","onlineOnly":"N","additionalOnlineFiles":"Y","costCenters":[{"id":583,"text":"Texas Water Science Center","active":true,"usgs":true}],"links":[{"id":118473,"rank":0,"type":{"id":24,"text":"Thumbnail"},"url":"https://pubs.usgs.gov/thumbnails/sir_2009_5174.jpg"},{"id":12988,"rank":100,"type":{"id":15,"text":"Index Page"},"url":"https://pubs.usgs.gov/sir/2009/5174/","linkFileType":{"id":5,"text":"html"}}],"geographicExtents":"{ \"type\": \"FeatureCollection\", \"features\": [ { \"type\": \"Feature\", \"properties\": {}, \"geometry\": { \"type\": \"Polygon\", \"coordinates\": [ [ [ -98,27 ], [ -98,33 ], [ -93,33 ], [ -93,27 ], [ -98,27 ] ] ] } } ] }","noUsgsAuthors":false,"publicationStatus":"PW","scienceBaseUri":"4f4e4b25e4b07f02db6aedc3","contributors":{"authors":[{"text":"Heitmuller, Franklin T.","contributorId":67476,"corporation":false,"usgs":true,"family":"Heitmuller","given":"Franklin","email":"","middleInitial":"T.","affiliations":[],"preferred":false,"id":303239,"contributorType":{"id":1,"text":"Authors"},"rank":1},{"text":"Greene, Lauren E.","contributorId":22868,"corporation":false,"usgs":true,"family":"Greene","given":"Lauren","email":"","middleInitial":"E.","affiliations":[],"preferred":false,"id":303238,"contributorType":{"id":1,"text":"Authors"},"rank":2}]}} ,{"id":97816,"text":"sir20095092 - 2009 - Geophysical Characterization of the Quaternary-Cretaceous Contact Using Surface Resistivity Methods in Franklin and Webster Counties, South-Central Nebraska","interactions":[],"lastModifiedDate":"2012-03-08T17:16:30","indexId":"sir20095092","displayToPublicDate":"2009-09-12T00:00:00","publicationYear":"2009","noYear":false,"publicationType":{"id":18,"text":"Report"},"publicationSubtype":{"id":5,"text":"USGS Numbered Series"},"seriesTitle":{"id":334,"text":"Scientific Investigations Report","code":"SIR","onlineIssn":"2328-0328","printIssn":"2328-031X","active":true,"publicationSubtype":{"id":5}},"seriesNumber":"2009-5092","title":"Geophysical Characterization of the Quaternary-Cretaceous Contact Using Surface Resistivity Methods in Franklin and Webster Counties, South-Central Nebraska","docAbstract":"To help manage and understand the Platte River system in Nebraska, the Platte River Cooperative Hydrology Study (COHYST), a group of state and local governmental agencies, developed a regional ground-water model. The southern boundary of this model lies along the Republican River, where an area with insufficient geologic data immediately north of the Republican River led to problems in the conceptualization of the simulated flow system and to potential problems with calibration of the simulation. Geologic descriptions from a group of test holes drilled in south-central Nebraska during 2001 and 2002 indicated a possible hydrologic disconnection between the Quaternary-age alluvial deposits in the uplands and those in the Republican River lowland. This disconnection was observed near a topographic high in the Cretaceous-age Niobrara Formation, which is the local bedrock. In 2003, the U.S. Geological Survey, in cooperation with the COHYST, collected surface geophysical data near these test holes to better define this discontinuity.\r\n\r\nTwo-dimensional imaging methods for direct-current resistivity and capacitively coupled resistivity were used to define the subsurface distribution of resistivity along several county roads near Riverton and Inavale, Nebraska. The relation between the subsurface distribution of resistivity and geology was defined by comparing existing geologic descriptions of test holes to surface-geophysical resistivity data along two profiles and using the information gained from these comparisons to interpret the remaining four profiles. In all of the resistivity profile sections, there was generally a three-layer subsurface interpretation, with a resistor located between two conductors. Further comparison of geologic data with the geophysical data and with surficial features was used to identify a topographic high in the Niobrara Formation near the Franklin Canal which was coincident with a resistivity high. Electrical properties of the Niobrara Formation made accurate interpretation of the resistivity profile sections difficult and less confident because of similar resistivity of this formation and that of the coarser-grained sediment of the Quaternary-age deposits. However, distinct conductive features were identified within the resistivity profile sections that aided in delineating the contact between the resistive Quaternary-age deposits and the resistive Niobrara Formation. Using this information, an interpretive boundary was drawn on the resistivity profile sections to represent the contact between the Quaternary-age alluvial deposits and the Cretaceous-age Niobrara Formation.\r\n\r\nA digital elevation model (DEM) of the top of the Niobrara Formation was constructed using the altitudes from the interpreted contact lines. This DEM showed that the general trend of top of the Niobrara Formation dips to the southeast. At the north edge of the study site, the Niobrara Formation topographic high trends east-west with an altitude range of 559 meters in the west to 543 meters in the east. Based on the land-surface elevation and the Niobrara Formation DEM, the estimated thickness of the Quaternary-age alluvial deposits throughout the study area was mapped and showed a thinning of the Quaternary-age alluvial deposits to the north, approximately where the topographic high of the Niobrara Formation is located. This topographic high in the Niobrara Formation has the potential to act as a barrier to ground-water flow from the uplands alluvial aquifer to the Republican River alluvial aquifer as shown in the resistivity profile sections. The Quaternary-age alluvial deposits in the uplands and those in the Republican River Valley are not fully represented as disconnected because it is possible that there are ground-water flow paths that were not mapped during this study.","language":"ENGLISH","publisher":"U.S. Geological Survey","doi":"10.3133/sir20095092","collaboration":"Prepared in cooperation with the Platte River Cooperative Hydrology Study","usgsCitation":"Teeple, A., Kress, W.H., Cannia, J.C., and Ball, L.B., 2009, Geophysical Characterization of the Quaternary-Cretaceous Contact Using Surface Resistivity Methods in Franklin and Webster Counties, South-Central Nebraska: U.S. Geological Survey Scientific Investigations Report 2009-5092, vi, 35 p., https://doi.org/10.3133/sir20095092.","productDescription":"vi, 35 p.","costCenters":[{"id":464,"text":"Nebraska Water Science Center","active":true,"usgs":true}],"links":[{"id":118632,"rank":0,"type":{"id":24,"text":"Thumbnail"},"url":"https://pubs.usgs.gov/thumbnails/sir_2009_5092.jpg"},{"id":12989,"rank":100,"type":{"id":15,"text":"Index Page"},"url":"https://pubs.usgs.gov/sir/2009/5092/","linkFileType":{"id":5,"text":"html"}}],"noUsgsAuthors":false,"publicationStatus":"PW","scienceBaseUri":"4f4e4ac9e4b07f02db67c47c","contributors":{"authors":[{"text":"Teeple, Andrew 0000-0003-1781-8354 apteeple@usgs.gov","orcid":"https://orcid.org/0000-0003-1781-8354","contributorId":1399,"corporation":false,"usgs":true,"family":"Teeple","given":"Andrew ","email":"apteeple@usgs.gov","affiliations":[{"id":583,"text":"Texas Water Science Center","active":true,"usgs":true}],"preferred":false,"id":303241,"contributorType":{"id":1,"text":"Authors"},"rank":1},{"text":"Kress, Wade H.","contributorId":100475,"corporation":false,"usgs":true,"family":"Kress","given":"Wade","email":"","middleInitial":"H.","affiliations":[],"preferred":false,"id":303243,"contributorType":{"id":1,"text":"Authors"},"rank":2},{"text":"Cannia, James C.","contributorId":94356,"corporation":false,"usgs":true,"family":"Cannia","given":"James","email":"","middleInitial":"C.","affiliations":[],"preferred":false,"id":303242,"contributorType":{"id":1,"text":"Authors"},"rank":3},{"text":"Ball, Lyndsay B. 0000-0002-6356-4693 lbball@usgs.gov","orcid":"https://orcid.org/0000-0002-6356-4693","contributorId":1138,"corporation":false,"usgs":true,"family":"Ball","given":"Lyndsay","email":"lbball@usgs.gov","middleInitial":"B.","affiliations":[{"id":211,"text":"Crustal Geophysics and Geochemistry Science Center","active":true,"usgs":true}],"preferred":true,"id":303240,"contributorType":{"id":1,"text":"Authors"},"rank":4}]}} ,{"id":97813,"text":"ofr20091166 - 2009 - Downsizer - A Graphical User Interface-Based Application for Browsing, Acquiring, and Formatting Time-Series Data for Hydrologic Modeling","interactions":[],"lastModifiedDate":"2012-02-02T00:14:30","indexId":"ofr20091166","displayToPublicDate":"2009-09-11T00:00:00","publicationYear":"2009","noYear":false,"publicationType":{"id":18,"text":"Report"},"publicationSubtype":{"id":5,"text":"USGS Numbered Series"},"seriesTitle":{"id":330,"text":"Open-File Report","code":"OFR","onlineIssn":"2331-1258","printIssn":"0196-1497","active":true,"publicationSubtype":{"id":5}},"seriesNumber":"2009-1166","title":"Downsizer - A Graphical User Interface-Based Application for Browsing, Acquiring, and Formatting Time-Series Data for Hydrologic Modeling","docAbstract":"The U.S. Geological Survey Downsizer is a computer application that selects, downloads, verifies, and formats station-based time-series data for environmental-resource models, particularly the Precipitation-Runoff Modeling System. Downsizer implements the client-server software architecture. The client presents a map-based, graphical user interface that is intuitive to modelers; the server provides streamflow and climate time-series data from over 40,000 measurement stations across the United States. This report is the Downsizer user's manual and provides (1) an overview of the software design, (2) installation instructions, (3) a description of the graphical user interface, (4) a description of selected output files, and (5) troubleshooting information.","language":"ENGLISH","publisher":"U.S. Geological Survey","doi":"10.3133/ofr20091166","usgsCitation":"Ward-Garrison, C., Markstrom, S., and Hay, L.E., 2009, Downsizer - A Graphical User Interface-Based Application for Browsing, Acquiring, and Formatting Time-Series Data for Hydrologic Modeling: U.S. Geological Survey Open-File Report 2009-1166, iv, 27 p., https://doi.org/10.3133/ofr20091166.","productDescription":"iv, 27 p.","onlineOnly":"Y","costCenters":[{"id":595,"text":"U.S. Geological Survey","active":false,"usgs":true}],"links":[{"id":125480,"rank":0,"type":{"id":24,"text":"Thumbnail"},"url":"https://pubs.usgs.gov/thumbnails/ofr_2009_1166.jpg"},{"id":12985,"rank":100,"type":{"id":15,"text":"Index Page"},"url":"https://pubs.usgs.gov/of/2009/1166/","linkFileType":{"id":5,"text":"html"}}],"noUsgsAuthors":false,"publicationStatus":"PW","scienceBaseUri":"4f4e4a54e4b07f02db62c1bd","contributors":{"authors":[{"text":"Ward-Garrison, Christian","contributorId":85682,"corporation":false,"usgs":true,"family":"Ward-Garrison","given":"Christian","affiliations":[],"preferred":false,"id":303233,"contributorType":{"id":1,"text":"Authors"},"rank":1},{"text":"Markstrom, Steven L. 0000-0001-7630-9547 markstro@usgs.gov","orcid":"https://orcid.org/0000-0001-7630-9547","contributorId":1986,"corporation":false,"usgs":true,"family":"Markstrom","given":"Steven L.","email":"markstro@usgs.gov","affiliations":[{"id":5044,"text":"National Research Program - Central Branch","active":true,"usgs":true}],"preferred":false,"id":303232,"contributorType":{"id":1,"text":"Authors"},"rank":2},{"text":"Hay, Lauren E. 0000-0003-3763-4595 lhay@usgs.gov","orcid":"https://orcid.org/0000-0003-3763-4595","contributorId":1287,"corporation":false,"usgs":true,"family":"Hay","given":"Lauren","email":"lhay@usgs.gov","middleInitial":"E.","affiliations":[{"id":5044,"text":"National Research Program - Central Branch","active":true,"usgs":true}],"preferred":true,"id":303231,"contributorType":{"id":1,"text":"Authors"},"rank":3}]}} ,{"id":97811,"text":"ofr20091174 - 2009 - Seismic Velocities and Thicknesses of Alluvial Deposits along Baker Creek in the Great Basin National Park, East-Central Nevada","interactions":[],"lastModifiedDate":"2012-03-08T17:16:27","indexId":"ofr20091174","displayToPublicDate":"2009-09-11T00:00:00","publicationYear":"2009","noYear":false,"publicationType":{"id":18,"text":"Report"},"publicationSubtype":{"id":5,"text":"USGS Numbered Series"},"seriesTitle":{"id":330,"text":"Open-File Report","code":"OFR","onlineIssn":"2331-1258","printIssn":"0196-1497","active":true,"publicationSubtype":{"id":5}},"seriesNumber":"2009-1174","title":"Seismic Velocities and Thicknesses of Alluvial Deposits along Baker Creek in the Great Basin National Park, East-Central Nevada","docAbstract":"To better understand how proposed large-scale water withdrawals in Snake Valley may affect the water resources and hydrologic processes in the Great Basin National Park, the National Park Service needs to have a better understanding of the relations between streamflow and groundwater flow through alluvium and karst topography of the Pole Canyon Limestone. Information that is critical to understanding these relations is the thickness of alluvial deposits that overlay the Pole Canyon Limestone.\r\n\r\nIn mid-April 2009, the U.S. Geological Survey and National Park Service used seismic refraction along three profiles adjacent to Baker Creek to further refine understanding of the local geology. Two refractors and three distinct velocity layers were detected along two of the profiles and a single refractor and two distinct velocity layers were detected along a third profile.\r\n\r\nIn the unsaturated alluvium, average velocity was 2,000 feet per second, thickness ranged from about 7 to 20 feet along two profiles downstream of the Narrows, and thickness was at least 100 feet along a single profile upstream of the Narrows. Saturated alluvium was only present downstream of the Narrows - average velocity was 4,400 feet per second, and thickness ranged from about 40 to 110 feet. The third layer probably represented Pole Canyon Limestone or Tertiary granitic rock units with an average velocity of 12,500 feet per second. Along the upstream and middle profiles (profiles 3 and 1, respectively), the depth to top of the third layer ranged from at least 60 to 110 feet below land surface and is most likely the Pole Canyon Limestone. The third layer at the farthest downstream profile (profile 2) may be a Tertiary granitic rock unit.\r\n\r\nBaker Creek is disconnected from the groundwater system along the upstream profile (profile 3) and streamflow losses infiltrate vertically downward to the Pole Canyon Limestone. Along the downstream and middle profiles (profiles 2 and 1, respectively), the presence of a shallow water table indicates that low permeability Tertiary granitic rock may extend across the Baker Creek Drainage intersecting the Pole Canyon Limestone. The Tertiary granitic rock may be acting as a barrier to groundwater flow within the Pole Canyon Limestone.","language":"ENGLISH","publisher":"U.S. Geological Survey","doi":"10.3133/ofr20091174","collaboration":"Prepared in cooperation with the National Park Service","usgsCitation":"Allander, K.K., and Berger, D.L., 2009, Seismic Velocities and Thicknesses of Alluvial Deposits along Baker Creek in the Great Basin National Park, East-Central Nevada: U.S. Geological Survey Open-File Report 2009-1174, iv, 15 p., https://doi.org/10.3133/ofr20091174.","productDescription":"iv, 15 p.","costCenters":[{"id":465,"text":"Nevada Water Science Center","active":true,"usgs":true}],"links":[{"id":118529,"rank":0,"type":{"id":24,"text":"Thumbnail"},"url":"https://pubs.usgs.gov/thumbnails/ofr_2009_1174.jpg"},{"id":12983,"rank":100,"type":{"id":15,"text":"Index Page"},"url":"https://pubs.usgs.gov/of/2009/1174/","linkFileType":{"id":5,"text":"html"}}],"geographicExtents":"{ \"type\": \"FeatureCollection\", \"features\": [ { \"type\": \"Feature\", \"properties\": {}, \"geometry\": { \"type\": \"Polygon\", \"coordinates\": [ [ [ -114.5,38.666666666666664 ], [ -114.5,39.166666666666664 ], [ -114,39.166666666666664 ], [ -114,38.666666666666664 ], [ -114.5,38.666666666666664 ] ] ] } } ] }","noUsgsAuthors":false,"publicationStatus":"PW","scienceBaseUri":"4f4e4a0de4b07f02db5fcfb7","contributors":{"authors":[{"text":"Allander, Kip K. 0000-0002-3317-298X kalland@usgs.gov","orcid":"https://orcid.org/0000-0002-3317-298X","contributorId":2290,"corporation":false,"usgs":true,"family":"Allander","given":"Kip","email":"kalland@usgs.gov","middleInitial":"K.","affiliations":[{"id":465,"text":"Nevada Water Science Center","active":true,"usgs":true}],"preferred":true,"id":303227,"contributorType":{"id":1,"text":"Authors"},"rank":1},{"text":"Berger, David L. dlberger@usgs.gov","contributorId":1861,"corporation":false,"usgs":true,"family":"Berger","given":"David","email":"dlberger@usgs.gov","middleInitial":"L.","affiliations":[{"id":595,"text":"U.S. Geological Survey","active":false,"usgs":true}],"preferred":false,"id":303226,"contributorType":{"id":1,"text":"Authors"},"rank":2}]}} ,{"id":97802,"text":"tm6A32 - 2009 - The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1","interactions":[],"lastModifiedDate":"2012-03-08T17:16:30","indexId":"tm6A32","displayToPublicDate":"2009-09-05T00:00:00","publicationYear":"2009","noYear":false,"publicationType":{"id":18,"text":"Report"},"publicationSubtype":{"id":5,"text":"USGS Numbered Series"},"seriesTitle":{"id":335,"text":"Techniques and Methods","code":"TM","onlineIssn":"2328-7055","printIssn":"2328-7047","active":true,"publicationSubtype":{"id":5}},"seriesNumber":"6-A32","title":"The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1","docAbstract":"The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'.\r\n\r\nThe first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain.\r\n\r\nAll of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR).\r\n\r\nIn addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas.\r\n\r\nThe example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de","language":"ENGLISH","publisher":"U.S. Geological Survey","doi":"10.3133/tm6A32","usgsCitation":"Schmid, W., and Hanson, R.T., 2009, The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1: U.S. Geological Survey Techniques and Methods 6-A32, x, 103 p., https://doi.org/10.3133/tm6A32.","productDescription":"x, 103 p.","costCenters":[{"id":154,"text":"California Water Science Center","active":true,"usgs":true}],"links":[{"id":118600,"rank":0,"type":{"id":24,"text":"Thumbnail"},"url":"https://pubs.usgs.gov/thumbnails/tm_6_a32.jpg"},{"id":12973,"rank":100,"type":{"id":15,"text":"Index Page"},"url":"https://pubs.usgs.gov/tm/tm6a32/","linkFileType":{"id":5,"text":"html"}}],"noUsgsAuthors":false,"publicationStatus":"PW","scienceBaseUri":"4f4e4ac9e4b07f02db67c861","contributors":{"authors":[{"text":"Schmid, Wolfgang","contributorId":84020,"corporation":false,"usgs":false,"family":"Schmid","given":"Wolfgang","affiliations":[{"id":13040,"text":"Department of Hydrology and Water Resources, University of Arizona","active":true,"usgs":false}],"preferred":false,"id":303209,"contributorType":{"id":1,"text":"Authors"},"rank":1},{"text":"Hanson, R. T.","contributorId":91148,"corporation":false,"usgs":true,"family":"Hanson","given":"R.","email":"","middleInitial":"T.","affiliations":[],"preferred":false,"id":303210,"contributorType":{"id":1,"text":"Authors"},"rank":2}]}} ]}