Geology ReportsSearch

Geology topics

Edward R. Banta

Publications and source records attributed to Edward R. Banta.

At least 19 recordsLinked to original sources

Documentation for the MODFLOW 6 Groundwater Flow Model

This report documents the Groundwater Flow (GWF) Model for a new version of MODFLOW called MODFLOW 6. The GWF Model for MODFLOW 6 is based on a generalized control-volume finite-difference approach in which a cell can be hydraulically connected to any number of surrounding cells. Users can define the model grid using one of three discretization packages, including (1) a structured discretization package for defining regular MODFLOW grids consisting of layers, rows, and columns, (2) a discretization by ver­tices package for defining layered unstructured grids consisting of layers and cells, and (3) a general unstruc­tured discretization package for defining flexible grids comprised of cells and their connection properties. For layered grids, a new capability is available for removing thin cells and vertically connecting cells overlying and underlying the thin cells. For complex problems involving water-table conditions, an optional Newton-Raphson formulation, based on the formulations in MODFLOW-NWT and MODFLOW-USG, can be acti­vated. Use of the Newton-Raphson formulation will often improve model convergence and allow solutions to be obtained for difficult problems that cannot be solved using the traditional wetting and drying approach. The GWF Model is divided into “packages,” as was done in previous MODFLOW versions. A package is the part of the model that deals with a single aspect of simulation. Packages included with the GWF Model include those related to internal calculations of groundwater flow (discretization, initial conditions, hydraulic conduc­tance, and storage), stress packages (constant heads, wells, recharge, rivers, general head boundaries, drains, and evapotranspiration), and advanced stress packages (streamflow routing, lakes, multi-aquifer wells, and unsaturated zone flow). An additional package is also available for moving water available in one package into the individual features of the advanced stress packages. The GWF Model also has packages for obtaining and controlling output from the model. This report includes detailed explanations of physical and mathematical concepts on which the GWF Model and its packages are based. Like its predecessors, MODFLOW 6 is based on a highly modular structure; however, this structure has been extended into an object-oriented framework. The framework includes a robust and generalized numeri­cal solution object, which can be used to solve many different types of models. The numerical solution object has several different matrix preconditioning options as well as several methods for solving the linear system of equations. In this new framework, the GWF Model itself is an object as are each of the GWF Model packages. A benefit of the object-oriented structure is that multiple objects of the same type can be used in a single sim­ulation. Thus, a single forward run with MODFLOW 6 may contain multiple GWF Models. GWF Models can be hydraulically connected using GWF-GWF Exchange objects. Connecting GWF models in different ways permits the user to utilize a local grid refinement strategy consisting of parent and child models or to couple adjacent GWF Models. An advantage of the approach implemented in MODFLOW 6 is that multiple models and their exchanges can be incorporated into a single numerical solution object. With this design, models can be tightly coupled at the matrix level.

Techniques and Methods

Documentation for the MODFLOW 6 framework

MODFLOW is a popular open-source groundwater flow model distributed by the U.S. Geological Survey. Growing interest in surface and groundwater interactions, local refinement with nested and unstructured grids, karst groundwater flow, solute transport, and saltwater intrusion, has led to the development of numerous MODFLOW versions. Often times, there are incompatibilities between these different MODFLOW versions. The report describes a new MODFLOW framework called MODFLOW 6 that is designed to support multiple models and multiple types of models. The framework is written in Fortran using a modular object-oriented design. The primary framework components include the simulation (or main program), Timing Module, Solutions, Models, Exchanges, and Utilities. The first version of the framework focuses on numerical solutions, numerical models, and numerical exchanges. This focus on numerical models allows multiple numerical models to be tightly coupled at the matrix level.

Techniques and Methods

Two graphical user interfaces for managing and analyzing MODFLOW groundwater-model scenarios

Scenario Manager and Scenario Analyzer are graphical user interfaces that facilitate the use of calibrated, MODFLOW-based groundwater models for investigating possible responses to proposed stresses on a groundwater system. Scenario Manager allows a user, starting with a calibrated model, to design and run model scenarios by adding or modifying stresses simulated by the model. Scenario Analyzer facilitates the process of extracting data from model output and preparing such display elements as maps, charts, and tables. Both programs are designed for users who are familiar with the science on which groundwater modeling is based but who may not have a groundwater modeler’s expertise in building and calibrating a groundwater model from start to finish. With Scenario Manager, the user can manipulate model input to simulate withdrawal or injection wells, time-variant specified hydraulic heads, recharge, and such surface-water features as rivers and canals. Input for stresses to be simulated comes from user-provided geographic information system files and time-series data files. A Scenario Manager project can contain multiple scenarios and is self-documenting. Scenario Analyzer can be used to analyze output from any MODFLOW-based model; it is not limited to use with scenarios generated by Scenario Manager. Model-simulated values of hydraulic head, drawdown, solute concentration, and cell-by-cell flow rates can be presented in display elements. Map data can be represented as lines of equal value (contours) or as a gradated color fill. Charts and tables display time-series data obtained from output generated by a transient-state model run or from user-provided text files of time-series data. A display element can be based entirely on output of a single model run, or, to facilitate comparison of results of multiple scenarios, an element can be based on output from multiple model runs. Scenario Analyzer can export display elements and supporting metadata as a Portable Document Format file.

Techniques and Methods

GWM-VI: groundwater management with parallel processing for multiple MODFLOW versions

Groundwater Management–Version Independent (GWM–VI) is a new version of the Groundwater Management Process of MODFLOW. The Groundwater Management Process couples groundwater-flow simulation with a capability to optimize stresses on the simulated aquifer based on an objective function and constraints imposed on stresses and aquifer state. GWM–VI extends prior versions of Groundwater Management in two significant ways—(1) it can be used with any version of MODFLOW that meets certain requirements on input and output, and (2) it is structured to allow parallel processing of the repeated runs of the MODFLOW model that are required to solve the optimization problem. GWM–VI uses the same input structure for files that describe the management problem as that used by prior versions of Groundwater Management. GWM–VI requires only minor changes to the input files used by the MODFLOW model. GWM–VI uses the Joint Universal Parameter IdenTification and Evaluation of Reliability Application Programming Interface (JUPITER-API) to implement both version independence and parallel processing. GWM–VI communicates with the MODFLOW model by manipulating certain input files and interpreting results from the MODFLOW listing file and binary output files. Nearly all capabilities of prior versions of Groundwater Management are available in GWM–VI. GWM–VI has been tested with MODFLOW-2005, MODFLOW-NWT (a Newton formulation for MODFLOW-2005), MF2005-FMP2 (the Farm Process for MODFLOW-2005), SEAWAT, and CFP (Conduit Flow Process for MODFLOW-2005). This report provides sample problems that demonstrate a range of applications of GWM–VI and the directory structure and input information required to use the parallel-processing capability.

Techniques and Methods

Demonstration optimization analyses of pumping from selected Arapahoe aquifer municipal wells in the west-central Denver Basin, Colorado, 2010–2109

Declining water levels caused by withdrawals of water from wells in the west-central part of the Denver Basin bedrock-aquifer system have raised concerns with respect to the ability of the aquifer system to sustain production. The Arapahoe aquifer in particular is heavily used in this area. Two optimization analyses were conducted to demonstrate approaches that could be used to evaluate possible future pumping scenarios intended to prolong the productivity of the aquifer and to delay excessive loss of saturated thickness. These analyses were designed as demonstrations only, and were not intended as a comprehensive optimization study. Optimization analyses were based on a groundwater-flow model of the Denver Basin developed as part of a recently published U.S. Geological Survey groundwater-availability study. For each analysis an optimization problem was set up to maximize total withdrawal rate, subject to withdrawal-rate and hydraulic-head constraints, for 119 selected municipal water-supply wells located in 96 model cells. The optimization analyses were based on 50- and 100-year simulations of groundwater withdrawals. The optimized total withdrawal rate for all selected wells for a 50-year simulation time was about 58.8 cubic feet per second. For an analysis in which the simulation time and head-constraint time were extended to 100 years, the optimized total withdrawal rate for all selected wells was about 53.0 cubic feet per second, demonstrating that a reduction in withdrawal rate of about 10 percent may extend the time before the hydraulic-head constraints are violated by 50 years, provided that pumping rates are optimally distributed. Analysis of simulation results showed that initially, the pumping produces water primarily by release of water from storage in the Arapahoe aquifer. However, because confining layers between the Denver and Arapahoe aquifers are thin, in less than 5 years, most of the water removed by managed-flows pumping likely would be supplied by depleting overlying hydrogeologic units, substantially increasing the rate of decline of hydraulic heads in parts of the overlying Denver aquifer.

Colorado

MODFLOW-CDSS, a version of MODFLOW-2005 with modifications for Colorado Decision Support Systems

MODFLOW-CDSS is a three-dimensional, finite-difference groundwater-flow model based on MODFLOW-2005, with two modifications. The first modification is the introduction of a Partition Stress Boundaries capability, which enables the user to partition a selected subset of MODFLOW's stress-boundary packages, with each partition defined by a separate input file. Volumetric water-budget components of each partition are tracked and listed separately in the volumetric water-budget tables. The second modification enables the user to specify that execution of a simulation should continue despite failure of the solver to satisfy convergence criteria. This modification is particularly intended to be used in conjunction with automated model-analysis software; its use is not recommended for other purposes.

Open-File Report

ModelMate - A graphical user interface for model analysis

ModelMate is a graphical user interface designed to facilitate use of model-analysis programs with models. This initial version of ModelMate supports one model-analysis program, UCODE_2005, and one model software program, MODFLOW-2005. ModelMate can be used to prepare input files for UCODE_2005, run UCODE_2005, and display analysis results. A link to the GW_Chart graphing program facilitates visual interpretation of results. ModelMate includes capabilities for organizing directories used with the parallel-processing capabilities of UCODE_2005 and for maintaining files in those directories to be identical to a set of files in a master directory. ModelMate can be used on its own or in conjunction with ModelMuse, a graphical user interface for MODFLOW-2005 and PHAST.

Techniques and Methods

The U.S. Geological Survey Modular Ground-Water Model - PCGN: A Preconditioned Conjugate Gradient Solver with Improved Nonlinear Control

The preconditioned conjugate gradient with improved nonlinear control (PCGN) package provides addi-tional means by which the solution of nonlinear ground-water flow problems can be controlled as compared to existing solver packages for MODFLOW. Picard iteration is used to solve nonlinear ground-water flow equations by iteratively solving a linear approximation of the nonlinear equations. The linear solution is provided by means of the preconditioned conjugate gradient algorithm where preconditioning is provided by the modi-fied incomplete Cholesky algorithm. The incomplete Cholesky scheme incorporates two levels of fill, 0 and 1, in which the pivots can be modified so that the row sums of the preconditioning matrix and the original matrix are approximately equal. A relaxation factor is used to implement the modified pivots, which determines the degree of modification allowed. The effects of fill level and degree of pivot modification are briefly explored by means of a synthetic, heterogeneous finite-difference matrix; results are reported in the final section of this report. The preconditioned conjugate gradient method is coupled with Picard iteration so as to efficiently solve the nonlinear equations associated with many ground-water flow problems. The description of this coupling of the linear solver with Picard iteration is a primary concern of this document.

Open-File Report

User Guide for HUFPrint, A Tabulation and Visualization Utility for the Hydrogeologic-Unit Flow (HUF) Package of MODFLOW

This report documents HUFPrint, a computer program that extracts and displays information about model structure and hydraulic properties from the input data for a model built using the Hydrogeologic-Unit Flow (HUF) Package of the U.S. Geological Survey's MODFLOW program for modeling ground-water flow. HUFPrint reads the HUF Package and other MODFLOW input files, processes the data by hydrogeologic unit and by model layer, and generates text and graphics files useful for visualizing the data or for further processing. For hydrogeologic units, HUFPrint outputs such hydraulic properties as horizontal hydraulic conductivity along rows, horizontal hydraulic conductivity along columns, horizontal anisotropy, vertical hydraulic conductivity or anisotropy, specific storage, specific yield, and hydraulic-conductivity depth-dependence coefficient. For model layers, HUFPrint outputs such effective hydraulic properties as horizontal hydraulic conductivity along rows, horizontal hydraulic conductivity along columns, horizontal anisotropy, specific storage, primary direction of anisotropy, and vertical conductance. Text files tabulating hydraulic properties by hydrogeologic unit, by model layer, or in a specified vertical section may be generated. Graphics showing two-dimensional cross sections and one-dimensional vertical sections at specified locations also may be generated. HUFPrint reads input files designed for MODFLOW-2000 or MODFLOW-2005.

Techniques and Methods

User Guide and Documentation for Five MODFLOW Ground-Water Modeling Utility Programs

This report documents five utility programs designed for use in conjunction with ground-water flow models developed with the U.S. Geological Survey's MODFLOW ground-water modeling program. One program extracts calculated flow values from one model for use as input to another model. The other four programs extract model input or output arrays from one model and make them available in a form that can be used to generate an ArcGIS raster data set. The resulting raster data sets may be useful for visual display of the data or for further geographic data processing. The utility program GRID2GRIDFLOW reads a MODFLOW binary output file of cell-by-cell flow terms for one (source) model grid and converts the flow values to input flow values for a different (target) model grid. The spatial and temporal discretization of the two models may differ. The four other utilities extract selected 2-dimensional data arrays in MODFLOW input and output files and write them to text files that can be imported into an ArcGIS geographic information system raster format. These four utilities require that the model cells be square and aligned with the projected coordinate system in which the model grid is defined. The four raster-conversion utilities are * CBC2RASTER, which extracts selected stress-package flow data from a MODFLOW binary output file of cell-by-cell flows; * DIS2RASTER, which extracts cell-elevation data from a MODFLOW Discretization file; * MFBIN2RASTER, which extracts array data from a MODFLOW binary output file of head or drawdown; and * MULT2RASTER, which extracts array data from a MODFLOW Multiplier file.

Techniques and Methods

UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API

This report documents the computer codes UCODE_2005 and six post-processors. Together the codes can be used with existing process models to perform sensitivity analysis, data needs assessment, calibration, prediction, and uncertainty analysis. Any process model or set of models can be used; the only requirements are that models have numerical (ASCII or text only) input and output files, that the numbers in these files have sufficient significant digits, that all required models can be run from a single batch file or script, and that simulated values are continuous functions of the parameter values. Process models can include pre-processors and post-processors as well as one or more models related to the processes of interest (physical, chemical, and so on), making UCODE_2005 extremely powerful. An estimated parameter can be a quantity that appears in the input files of the process model(s), or a quantity used in an equation that produces a value that appears in the input files. In the latter situation, the equation is user-defined. UCODE_2005 can compare observations and simulated equivalents. The simulated equivalents can be any simulated value written in the process-model output files or can be calculated from simulated values with user-defined equations. The quantities can be model results, or dependent variables. For example, for ground-water models they can be heads, flows, concentrations, and so on. Prior, or direct, information on estimated parameters also can be considered. Statistics are calculated to quantify the comparison of observations and simulated equivalents, including a weighted least-squares objective function. In addition, UCODE_2005 can be used fruitfully in model calibration through its sensitivity analysis capabilities and its ability to estimate parameter values that result in the best possible fit to the observations. Parameters are estimated using nonlinear regression: a weighted least-squares objective function is minimized with respect to the parameter values using a modified Gauss-Newton method or a double-dogleg technique. Sensitivities needed for the method can be read from files produced by process models that can calculate sensitivities, such as MODFLOW-2000, or can be calculated by UCODE_2005 using a more general, but less accurate, forward- or central-difference perturbation technique. Problems resulting from inaccurate sensitivities and solutions related to the perturbation techniques are discussed in the report. Statistics are calculated and printed for use in (1) diagnosing inadequate data and identifying parameters that probably cannot be estimated; (2) evaluating estimated parameter values; and (3) evaluating how well the model represents the simulated processes. Results from UCODE_2005 and codes RESIDUAL_ANALYSIS and RESIDUAL_ANALYSIS_ADV can be used to evaluate how accurately the model represents the processes it simulates. Results from LINEAR_UNCERTAINTY can be used to quantify the uncertainty of model simulated values if the model is sufficiently linear. Results from MODEL_LINEARITY and MODEL_LINEARITY_ADV can be used to evaluate model linearity and, thereby, the accuracy of the LINEAR_UNCERTAINTY results. UCODE_2005 can also be used to calculate nonlinear confidence and predictions intervals, which quantify the uncertainty of model simulated values when the model is not linear. CORFAC_PLUS can be used to produce factors that allow intervals to account for model intrinsic nonlinearity and small-scale variations in system characteristics that are not explicitly accounted for in the model or the observation weighting. The six post-processing programs are independent of UCODE_2005 and can use the results of other programs that produce the required data-exchange files. UCODE_2005 and the other six codes are intended for use on any computer operating system. The programs consist of algorithms programmed in Fortran 90/95, which efficiently performs numerical calculations. The model runs required to obtain perturbation sensitivities can be performed using multiple processors. The programs are constructed in a modular fashion using JUPITER API conventions and modules. For example, the data-exchange files and input blocks are JUPITER API conventions and many of those used by UCODE_2005 are read or written by JUPITER API modules. UCODE-2005 includes capabilities likely to be required by many applications (programs) constructed using the JUPITER API, and can be used as a starting point for such programs.

Techniques and Methods

A new streamflow-routing (SFR1) package to simulate stream-aquifer interaction with MODFLOW-2000

The increasing concern for water and its quality require improved methods to evaluate the interaction between streams and aquifers and the strong influence that streams can have on the flow and transport of contaminants through many aquifers. For this reason, a new Streamflow-Routing (SFR1) Package was written for use with the U.S. Geological Survey's MODFLOW-2000 ground-water flow model. The SFR1 Package is linked to the Lake (LAK3) Package, and both have been integrated with the Ground-Water Transport (GWT) Process of MODFLOW-2000 (MODFLOW-GWT). SFR1 replaces the previous Stream (STR1) Package, with the most important difference being that stream depth is computed at the midpoint of each reach instead of at the beginning of each reach, as was done in the original Stream Package. This approach allows for the addition and subtraction of water from runoff, precipitation, and evapotranspiration within each reach. Because the SFR1 Package computes stream depth differently than that for the original package, a different name was used to distinguish it from the original Stream (STR1) Package. The SFR1 Package has five options for simulating stream depth and four options for computing diversions from a stream. The options for computing stream depth are: a specified value; Manning's equation (using a wide rectangular channel or an eight-point cross section); a power equation; or a table of values that relate flow to depth and width. Each stream segment can have a different option. Outflow from lakes can be computed using the same options. Because the wetted perimeter is computed for the eight-point cross section and width is computed for the power equation and table of values, the streambed conductance term no longer needs to be calculated externally whenever the area of streambed changes as a function of flow. The concentration of solute is computed in a stream network when MODFLOW-GWT is used in conjunction with the SFR1 Package. The concentration of a solute in a stream reach is based on a mass-balance approach and accounts for exchanges with (inputs from or losses to) ground-water systems. Two test examples are used to illustrate some of the capabilities of the SFR1 Package. The first test simulation was designed to illustrate how pumping of ground water from an aquifer connected to streams can affect streamflow, depth, width, and streambed conductance using the different options. The second test simulation was designed to illustrate solute transport through interconnected lakes, streams, and aquifers. Because of the need to examine time series results from the model simulations, the Gage Package first described in the LAK3 documentation was revised to include time series results of selected variables (streamflows, stream depth and width, streambed conductance, solute concentrations, and solute loads) for specified stream reaches. The mass-balance or continuity approach for routing flow and solutes through a stream network may not be applicable for all interactions between streams and aquifers. The SFR1 Package is best suited for modeling long-term changes (months to hundreds of years) in ground-water flow and solute concentrations using averaged flows in streams. The Package is not recommended for modeling the transient exchange of water between streams and aquifers when the objective is to examine short-term (minutes to days) effects caused by rapidly changing streamflows.

Open-File Report

MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model; documentation of packages for simulating evapotranspiration with a segmented function (ETS1) and drains with return flow (DRT1)

Two new packages for the U.S. Geological Survey modular finite-difference ground-water-flow model MODFLOW-2000 are documented. The new packages provide flexibility in simulating evapotranspiration and drain features not provided by the MODFLOW-2000 Evapotranspiration (EVT) and Drain (DRN) Packages. The report describes conceptualization of the packages, input instructions, listings and explanations of the source code, and example simulations. The new Evapotranspiration Segments (ETS1) Package allows simulation of evapotranspiration with a user-defined relation between evapotranspiration rate and hydraulic head. This capability provides a degree of flexibility not supported by the EVT Package, which has been available in MODFLOW since its initial release. In the ETS1 Package, the relation of evapotranspiration rate to hydraulic head is conceptualized as a segmented line between an evaporation surface, defined as the elevation where the evapotranspiration rate reaches a maximum, and an elevation located at an extinction depth below the evaporation surface, where the evapotranspiration rate reaches zero. The user supplies input to define as many intermediate segment endpoints as desired to define the relation of evapotranspiration rate to head between these two elevations. The EVT Package, in contrast, simulates evapotranspiration with a single linear function. The new Drain Return (DRT1) Package can be used to simulate the return flow of water discharged from a drain feature back into the ground-water system. The DRN Package, which has been available in MODFLOW since its initial release, does not have the capability to simulate return of flow. If the return-flow option of the DRT1 Package is selected, for each cell designated as a drain-return cell, the user has the option of specifying a proportion of the water simulated as leaving the ground-water system through the drain feature that is to be simulated as returning simultaneously to one other cell in the model.

Open-File Report

MODFLOW-2000, The U.S. Geological Survey modular ground-water model: User guide to modularization concepts and the ground-water flow process

MODFLOW is a computer program that numerically solves the three-dimensional ground-water flow equation for a porous medium by using a finite-difference method. Although MODFLOW was designed to be easily enhanced, the design was oriented toward additions to the ground-water flow equation. Frequently there is a need to solve additional equations; for example, transport equations and equations for estimating parameter values that produce the closest match between model-calculated heads and flows and measured values. This report documents a new version of MODFLOW, called MODFLOW-2000, which is designed to accommodate the solution of equations in addition to the ground-water flow equation. This report is a user's manual. It contains an overview of the old and added design concepts, documents one new package, and contains input instructions for using the model to solve the ground-water flow equation.

Open-File Report

Ground Water Atlas of the United States: Segment 2, Arizona, Colorado, New Mexico, Utah

This chapter of the Ground Water Atlas of the United States describes the aquifers in Arizona, Colorado, New Mexico, and Utah. These four States, which comprise Segment 2 of this Atlas, are located in the Southwestern United States and extend from the rolling grasslands of the Great Plains on the east across the Rocky Mountains and Continental Divide to the desert basins of the Southwest. The 425,000-square-mile area ranges in altitude from about 14,400 feet above sea level in the Rocky Mountains of Colorado to about 100 feet near the lower Colorado River in southwestern Arizona. All the ground water in Segment 2 ultimately is derived from infiltration of precipitation, which varies considerably with the altitude and topography of the area. The Great Plains Physiographic Province of the Central United States extends into eastern Colorado and New Mexico (fig. 1), where flat to rolling prairie (fig. 2) with scattered hills and bluffs gradually rises westward to 5,000 to 7,000 feet above sea level and abruptly gives way to the frontal ranges of the Rocky Mountains in the Southern Rocky Mountain and Basin and Range Physiographic Provinces. West of the frontal ranges in Colorado and northern New Mexico are additional and higher mountain ranges generally oriented north-south but with many spurs and extensions oriented in other directions. The many ranges of the Rocky Mountains are separated by valleys and high mountain parks (fig. 3). Colorado contains about three-fourths of the Nation's land area above 10,000 feet and has 53 mountain peaks higher than 14,000 feet. Most of these high peaks are located near the Continental Divide (fig. 1), which extends approximately north-south through central Colorado and western New Mexico. The altitude of the divide decreases in southern New Mexico to less than 4,500 feet in a few areas. Farther westward, the mountains are less prevalent and are interspersed with broad structural basins. These basins and the broad valleys of the middle Colorado River and its tributaries form the irregular intermontane topography of the Colorado Plateaus Physiographic Province (fig. 4). Plateaus and high mesas are formed where the surface has been dissected by rugged canyons carved by the Colorado River and its tributaries (fig. 5). The largest of these canyons-the Grand Canyon-extends about 220 miles southwestward from the mouth of the Little Colorado River in Arizona and ranges from 4 to 18 miles in width and from 2,700 to 5,700 feet in depth below the rim. Small mountain ranges and intervening broad desert valleys of the Basin and Range Physiographic Province are prevalent to the west and south of the Colorado Plateaus in western Utah, southern Arizona, and southern New Mexico (fig. 1). These mountain ranges generally protrude 3,000 to 6,000 feet above the surrounding valley floor (fig. 6) and commonly extend from 20 to 50 miles in a north or northwesterly direction. Small parts of the Middle Rocky Mountains and Wyoming Basin Physiographic Provinces extend into northwestern Colorado and northeastern Utah (fig. 1). The topography, geology, and hydrology of the two areas are described in Chapter 1 of this Atlas. Four of the Nation's major river systems have headwaters in the mountainous areas of Segment 2. The South Platte River of the Missouri River system drains the eastern slope of northern Colorado; the Arkansas River and its tributary, the Canadian River, drain southeastern Colorado and northeastern New Mexico; the Rio Grande and its tributary, the Pecos River, drain south-central Colorado and central New Mexico; and the Colorado River and its tributaries drain Arizona, eastern Utah, Northwestern New Mexico, and western Colorado (fig. 1). Western Utah is drained by numerous streams that terminate in local desert basins, the Great Salt Lake, or other local lakes and reservoirs. Because the Great Salt Lake lies in the Great Basin, which is the largest closed basin in North America, it has no outlet to the sea. The salinity of the lake water is about 20 percent or about 6 times the salinity of seawater. Most of Segment 2 is sparsely populated. The average population density of counties is less than 8 persons per square mile in about 65 percent of the four-State area (fig. 7). Population densities range from less than 0.5 person per square mile in a few rural counties to more than 4 ,000 persons per square mile in populous urban areas. The 1990 population of the four States was about 10 million; almost 70 percent of this population was in Arizona and Colorado. Most land in Segment 2 is undeveloped forest grassland, or desert shrubland, much of which is used for livestock grazing. Land used for production of commercial crops primarily is in eastern Colorado and eastern New Mexico.

Arizona, Colorado, New Mexico, Utah