Codicil to the geophysical expression of selected mineral deposit models
Explore the source record for details and available documents.
SEARCH · USGS Science
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
This appendix describes how a-priori rupture models were developed for the northern California Type-A faults. As described in the main body of this report, and in Appendix G, “a-priori” models represent an initial estimate of the rate of single and multi-segment surface ruptures on each fault. Whether or not a given model is moment balanced (i.e., satisfies section slip-rate data) depends on assumptions made regarding the average slip on each segment in each rupture (which in turn depends on the chosen magnitude-area relationship). Therefore, for a given set of assumptions, or branch on the logic tree, the methodology of the present Working Group (WGCEP-2007) is to find a final model that is as close as possible to the a-priori model, in the least squares sense, but that also satisfies slip rate and perhaps other data. This is analogous the WGCEP- 2002 approach of effectively voting on the relative rate of each possible rupture, and then finding the closest moment-balance model (under a more limiting set of assumptions than adopted by the present WGCEP, as described in detail in Appendix G). The 2002 Working Group Report (WCCEP, 2003, referred to here as WGCEP-2002), created segmented earthquake rupture forecast models for all faults in the region, including some that had been designated as Type B faults in the NSHMP, 1996, and one that had not previously been considered. The 2002 National Seismic Hazard Maps used the values from WGCEP-2002 for all the faults in the region, essentially treating all the listed faults as Type A faults. As discussed in Appendix A, the current WGCEP found that there are a number of faults with little or no data on slip-per-event, or dates of previous earthquakes. As a result, the WGCEP recommends that faults with minimal available earthquake recurrence data: the Greenville, Mount Diablo, San Gregorio, Monte Vista-Shannon and Concord-Green Valley be modeled as Type B faults to be consistent with similarly poorly-known faults statewide. As a result, the modified segmented models discussed here only concern the San Andreas, Hayward-Rodgers Creek, and Calaveras faults.
This data report presents long-term oceanographic observations made in western Massachusetts Bay at two locations: (1) 42 deg 22.6' N., 70 deg 47.0' W. (Site A, 33 m water depth) from December 1989 through December 2002 (figure 1), and (2) 42 deg 9.8' N., 70 deg 38.4' W. (Site B, 21 m water depth) from October 1997 through December 2002. Site A is approximately 1 km south of the new ocean outfall that began discharging treated sewage effluent from the Boston metropolitan area into Massachusetts Bay on September 6, 2000. These long-term oceanographic observations have been collected by the U.S. Geological Survey (USGS) in partnership with the Massachusetts Water Resources Authority (MWRA) and with logistical support from the U.S. Coast Guard (USCG - http://www.uscg.mil). This report presents time series data through December 2002, updating a similar report that presented data through December 2000 (Butman and others, 2002). In addition, the Statistics and Mean Flow sections include some new plots and tables and the format of the report has been streamlined by combining yearly figures into single .pdfs. Figure 1 (PDF format) The long-term measurements are planned to continue at least through 2005. The long-term oceanographic observations at Sites A and B are part of a USGS study designed to understand the transport and long-term fate of sediments and associated contaminants in the Massachusetts bays. (See http://woodshole.er.usgs.gov/project-pages/bostonharbor/ and Butman and Bothner, 1997.) The long-term observations document seasonal and inter-annual changes in currents, hydrography, and suspended-matter concentration in western Massachusetts Bay, and the importance of infrequent catastrophic events, such as major storms or hurricanes, in sediment resuspension and transport. They also provide observations for testing numerical models of circulation. This data report presents a description of the field program and instrumentation, an overview of the data through summary plots and statistics, and the data in NetCDF and ASCII format for the period December 1989 through December 2002 for Site A and October 1997 through December 2002 for Site B. The objective of this report is to make the data available in digital form and to provide summary plots and statistics to facilitate browsing of the long-term data set.
Epithermal gold-silver deposits are vein, stockwork, disseminated, and replacement deposits that are mined primarily for their gold and silver contents; some deposits also contain substantial resources of lead, zinc, copper, and (or) mercury. These deposits form in the uppermost parts of the crust, at depths less than about 1,500 meters below the water table, and at temperatures below about 300 °C. Most epithermal gold-silver deposits are genetically related to hydrothermal systems associated with subaerial volcanism and intrusion of calc-alkaline magmas along convergent plate margins. These deposits formed throughout most of geologic time, although most known deposits are Cenozoic, which reflects preferential preservation of these shallowly formed deposits in tectonically unstable regions. Epithermal gold-silver deposits range in size from tens of thousands to greater than 1 billion metric tons of ore and have gold contents of 0.1 to greater than 30 grams per metric ton and silver contents of less than 1 to several thousand grams per metric ton. Historically, these deposits have been an important source of gold and silver and are estimated to contain about 8 percent of global gold. The wide range of tonnage-grade characteristics makes epithermal gold-silver deposits an attractive target for small and large exploration and mining companies. This report constitutes a new descriptive model for epithermal gold-silver deposits. It summarizes characteristics of known deposits, including their geological, geophysical, geochemical, and geoenvironmental aspects. Models concerning the genesis of epithermal gold-silver deposits are discussed. The application of descriptive and genetic aspects of the model to mineral exploration and resource assessment of undiscovered deposits is described. Finally, areas where additional research is needed to better understand the genesis of these deposits are identified. An extensive summary table outlining the characteristics of about 100 epithermal gold-silver deposits is included as an appendix; this summary table includes most of the world’s largest epithermal gold-silver deposits, and many smaller, well-studied deposits.
This report presents a new three-dimensional geologic framework model (GFM) of the northern Great Plains region, encompassing parts of Montana, North Dakota, South Dakota, and Wyoming. The model provides a regionally consistent, geographic information system (GIS)-ready representation of Phanerozoic sedimentary strata, major fault systems, and Precambrian basement geometry across two sedimentary basins and adjacent uplifts. More than 300,000 geologic and geophysical data inputs were synthesized to model 41 stratigraphic horizons and 47 faults, yielding an internally coherent, sealed-volume interpretation of the subsurface. The modeling workflow developed for this study demonstrates an efficient and scalable approach for constructing basin-to regional-scale GFMs in geologically complex and data-variable settings. Although model fidelity varies with data density and quality, the resulting geometry is broadly consistent with 1:500,000-scale geologic mapping and highlights areas where additional geologic study is most needed. The three-dimensional GFM provides a foundational framework to support groundwater, energy, and mineral resource assessments, and offers a transferable methodology for potential future U.S. Geological Survey efforts to build large-area subsurface models in underexplored regions of the United States.
Estuaries are biogeochemical hot spots because they receive large inputs of nutrients and organic carbon from land and oceans to support high rates of metabolism and primary production. We synthesize published rates of annual phytoplankton primary production (APPP) in marine ecosystems influenced by connectivity to land – estuaries, bays, lagoons, fjords and inland seas. Review of the scientific literature produced a compilation of 1148 values of APPP derived from monthly incubation assays to measure carbon assimilation or oxygen production. The median value of median APPP measurements in 131 ecosystems is 185 and the mean is 252 g C m −2 yr −1 , but the range is large: from −105 (net pelagic production in the Scheldt Estuary) to 1890 g C m −2 yr −1 (net phytoplankton production in Tamagawa Estuary). APPP varies up to 10-fold within ecosystems and 5-fold from year to year (but we only found eight APPP series longer than a decade so our knowledge of decadal-scale variability is limited). We use studies of individual places to build a conceptual model that integrates the mechanisms generating this large variability: nutrient supply, light limitation by turbidity, grazing by consumers, and physical processes (river inflow, ocean exchange, and inputs of heat, light and wind energy). We consider method as another source of variability because the compilation includes values derived from widely differing protocols. A simulation model shows that different methods reported in the literature can yield up to 3-fold variability depending on incubation protocols and methods for integrating measured rates over time and depth. Although attempts have been made to upscale measures of estuarine-coastal APPP, the empirical record is inadequate for yielding reliable global estimates. The record is deficient in three ways. First, it is highly biased by the large number of measurements made in northern Europe (particularly the Baltic region) and North America. Of the 1148 reported values of APPP, 958 come from sites between 30 and 60° N; we found only 36 for sites south of 20° N. Second, of the 131 ecosystems where APPP has been reported, 37% are based on measurements at only one location during 1 year. The accuracy of these values is unknown but probably low, given the large interannual and spatial variability within ecosystems. Finally, global assessments are confounded by measurements that are not intercomparable because they were made with different methods. Phytoplankton primary production along the continental margins is tightly linked to variability of water quality, biogeochemical processes including ocean–atmosphere CO 2 exchange, and production at higher trophic levels including species we harvest as food. The empirical record has deficiencies that preclude reliable global assessment of this key Earth system process. We face two grand challenges to resolve these deficiencies: (1) organize and fund an international effort to use a common method and measure APPP regularly across a network of coastal sites that are globally representative and sustained over time, and (2) integrate data into a unifying model to explain the wide range of variability across ecosystems and to project responses of APPP to regional manifestations of global change as it continues to unfold.
A three-dimensional hydrodynamic model was developed as part of a cooperative U.S. Geological Survey/Massachusetts Water Resources Authority program to study contaminated sediment accumulation and transport in Massachusetts Bay. This report details the development of the model and assesses how well the model represents observed currents and water properties in the bay. It also summarizes circulation and comparative effluent dilution simulations from existing and future Boston sewage outfalls over a three-year period from October 1, 1989 to December 31, 1992. The ECOM-si model, a semi-implicit version of the Blumberg and Mellor (1987) Estuarine, Coastal and Ocean Model, is shown to reproduce many of the important hydrodynamical features of Massachusetts Bay: the seasonal evolution of the pycnocline, the mean flow pattern, and the strength of sub-tidal current fluctuations. Throughout the simulation period, during both vertically well-mixed and stratified conditions, the seasonal statistics of observed currents are well-represented by the model. The model is therefore appropriate for studying the average dilution of sewage effluent and other continuously discharged substances over seasonal time scales. The ability of the model to reproduce individual flow events varies with season and location within the bay. Flow events during unstratified conditions in western Massachusetts Bay are particularly well-represented, indicating that the model is appropriate for studying processes such as the transport of suspended material from the future outfall site due to winter storms. Individual flow events during stratified conditions and in the offshore Stellwagen Bank region, however, are less well-represented due to small length scales (caused by upwelling and river discharge events) coupled with insufficient data to specify open boundary forcing from the Gulf of Maine. Thus while the model might be used to answer issues such as the frequency with which Gulf of Maine river plumes visit the new outfall site, attempting to predict whether a particular plume would visit the outfall site could be problematic. Comparative simulations of effluent discharged from the existing and future Boston outfalls show that the region of relatively high effluent concentrations (1 part effluent to 200 parts sea water) is significantly smaller with the future outfall and is limited to Western Massachusetts Bay during both unstratified and stratified seasons. The region of even higher concentration (1 part effluent to 50 parts sea water) that covers much of Boston Harbor with the existing outfall is non-existent in the future outfall simulation. Additional simulations of chlorination plant failure predict that the offshore location of the future outfall will lead to dramatically lower levels of pathogens at area beaches.
Uranium in the central mining area of Marysvale, Utah occurs in hydrothermal veins cutting granitic and volcanic rocks in the eastern source area of the Mount Belknap Volcanics. A preliminary model for the origin of the veins envisages deposition in near-surface fractures above an unexposed pluton that may host a porphyry-type ore deposit. This model is based on the work in progress by the U.S. Geological Survey, and embodies all presently available data from field mapping, literature study, fluid inclusion studies, and diverse geochemical and isotopic studies. Recent advances in uranium geochemistry have been particularly helpful. The work is not yet complete so this model should be considered as a progress report suggesting possible targets for exploration and testing; such testing would in turn lead to refinements in the model and to a clearer understanding of vein-type uranium deposits in general.
The history of remote sensing and development of different sensors for environmental and natural resources mapping and data acquisition is reviewed and reported. Application examples in urban studies, hydrological modeling such as land-cover and floodplain mapping, fractional vegetation cover and impervious surface area mapping, surface energy flux and micro-topography correlation studies is discussed. The review also discusses the use of remotely sensed-based rainfall and potential evapotranspiration for estimating crop water requirement satisfaction index and hence provides early warning information for growers. The review is not an exhaustive application of the remote sensing techniques rather a summary of some important applications in environmental studies and modeling.
Potential changes in climate could alter interactions between environmental and societal systems and adversely affect the availability of water resources in many coastal communities. Changes in streamflow patterns in conjunction with sea-level rise may change the salinity-intrusion dynamics of coastal rivers. Several municipal water-supply intakes are located along the Georgia and South Carolina coast that are proximal to the present day saltwater-freshwater interface of tidal rivers. Increases in the extent of salinity intrusion resulting from climate change could threaten the availability of freshwater supplies in the vicinity of these intakes. To effectively manage these supplies, water-resource managers need estimates of potential changes in the frequency, duration, and magnitude of salinity intrusion near their water-supply intakes that may occur as a result of climate change. This study examines potential effects of climate change, including altered streamflow and sea-level rise, on the dynamics of saltwater intrusion near municipal water-supply intakes in two coastal areas. One area consists of the Atlantic Intracoastal Waterway (AIW) and the Waccamaw River near Myrtle Beach along the Grand Strand of the South Carolina Coast, and the second area is on or near the lower Savannah River near Savannah, Georgia. The study evaluated how future sea-level rise and a reduction in streamflows can potentially affect salinity intrusion and threaten municipal water supplies and the biodiversity of freshwater tidal marshes in these two areas. Salinity intrusion occurs as a result of the interaction between three principal forces—streamflow, mean coastal water levels, and tidal range. To analyze and simulate salinity dynamics at critical coastal gaging stations near four municipal water-supply intakes, various data-mining techniques, including artificial neural network (ANN) models, were used to evaluate hourly streamflow, salinity, and coastal water-level data collected over a period exceeding 10 years. The ANN models were trained (calibrated) to learn the specific interactions that cause salinity intrusions, and resulting models were able to accurately simulate historical salinity dynamics in both study areas. Changes in sea level and streamflow quantity and timing can be simulated by the salinity intrusion models to evaluate various climate-change scenarios. The salinity intrusion models for the study areas are deployed in a decision support system to facilitate the use of the models for management decisions by coastal water-resource managers. The report describes the use of the salinity-intrusion models decision support system to evaluate salinity-intrusion dynamics for various climate-change scenarios, including incremental increases in sea level in combination with incremental decreases in streamflow. Operation of municipal water-treatment plants is problematic when the specific-conductance values for source water are greater than 1,000 to 2,000 microsiemens per centimeter (µS/cm). High specific-conductance values contribute to taste problems that require treatment. Data from a gage downstream from a municipal water intake indicate specific conductance exceeded 1,000 µS/cm about 5.4 percent of the time over the 14-year period from August 1995 to August 2008. Simulations of specific conductance at this gaging station that incorporates sea-level rises resulted in a doubling of the exceedances to 11.0 percent for a 1-foot increase and 17.6 percent for a 2-foot increase. The frequency of intrusion of water with specific conductance values of 1,000 µS/cm was less sensitive to incremental reductions in streamflow than to incremental increases in sea level. Simulations of conditions associated with a 10-percent reduction in streamflow, in combination with a 1-foot rise in sea level, increased the percentage of time specific conductance exceeded 1,000 µS/cm at this site from 11.0 to 13.3 percent, and a 20-percent reduction in streamflow increased the percentage of time to 16.6 percent. Precipitation and temperature data from a global circulation model were used, after scale adjustments, as input to a watershed model of the Yadkin-Pee Dee River basin, which flows into the Waccamaw River and Atlantic Intracoastal Waterway study area in South Carolina. The simulated streamflow for historical conditions and projected climate change in the future was used as input for the ANN model in decision support system. Results of simulations incorporating climate-change projections for alterations in streamflow indicate an increase in the frequency of salinity-intrusion events and a shift in the seasonal occurrence of the intrusion events from the summer to the fall.
This report documents a review of the hydraulic and sediment-transport models developed by the City of Salisbury, Alcoa Power Generating, Inc., and the Federal Energy Regulatory Commission to address issues of flooding and sedimentation in the vicinity of Salisbury's water-supply intake 19.4 miles upstream from High Rock Dam. The objective of the review was to determine if the modeling results submitted by Salisbury clearly demonstrate that the presence of High Rock Dam has led to an increase in water levels at Salisbury facilities or, conversely, if the documents of Alcoa Power Generating, Inc., demonstrate that High Rock Dam has not had an effect on water levels at Salisbury facilities. No new data were collected as a part of the review, and the models developed by involved parties were not tested during the review. Some historical discharge-measurement notes and previously published reports were checked as part of the review. The one-dimensional hydraulic modeling results submitted by Alcoa Power Generating, Inc., did not assess the effects of changes in bathymetry on changes in flood levels at Salisbury's facilities because pre-impoundment conditions were not simulated. Hydraulic modeling performed by consultants for the City of Salisbury seems to indicate that both the presence of the dam in the absence of any post-impoundment sedimentation and changes in bathymetry between pre-impoundment and 1997 conditions have resulted in increased water levels relative to pre-impoundment conditions at Salisbury facilities on the Yadkin River for a fairly wide range of flows. The degree to which the dam and the changes in bathymetry have affected flood levels at the Salisbury facilities relative to pre-impoundment conditions is open to discussion because of uncertainty in topographic/bathymetric data and the absence of calibration and sensitivity testing of the hydraulic models. None of the three hydraulic models appears to have been calibrated to or tested against measurements, and no sensitivity testing was reported. Sediment-transport modeling results submitted by the City of Salisbury were calibrated, well documented, and provide a good understanding of the expected growth of the sediment delta in the upper end of High Rock Lake. Simulations made using this model seem to have demonstrated that the presence of the dam and the growth of the delta have resulted in increases in water-surface elevations at Salisbury's facilities over a range of flows and that these increases are expected to increase through time if current conditions remain unchanged.
The Miocene Columbia River Basalt Group and younger sedimentary deposits of lacustrine, fluvial, eolian, and cataclysmic-flood origins compose the aquifer system of the Quincy Basin in eastern Washington. Irrigation return flow and canal leakage from the Columbia Basin Project have caused groundwater levels to rise substantially in some areas. Water resource managers are considering extraction of additional stored groundwater to supply increasing demand. To help address these concerns, the transient groundwater model of the Quincy Basin documented in this report was developed to quantify the changes in groundwater flow and storage. The model based on the U.S. Geological Survey modular three-dimensional finite-difference numerical code MODFLOW uses a 1-kilometer finite-difference grid and is constrained by logs from 698 wells in the study area. Five model layers represent two sedimentary hydrogeologic units and underlying basalt formations. Head-dependent flux boundaries represent the Columbia River and other streams, lakes and reservoirs, underflow to and (or) from adjacent areas, and discharge to agricultural drains and springs. Specified flux boundaries represent recharge from precipitation and anthropogenic sources, including irrigation return flow and leakage from water-distribution canals and discharge through groundwater withdrawal wells. Transient conditions were simulated from 1920 to 2013 using annual stress periods. The model was calibrated with the parameter-estimation code PEST to a total of 4,064 water levels measured in 710 wells. Increased recharge since predevelopment resulted in an 11.5 million acre-feet increase in storage in the Quincy Groundwater Management Subarea of the Quincy Basin. Four groundwater-management scenarios were formulated with input from project stakeholders and were simulated using the calibrated model to provide representative examples of how the model could be used to evaluate the effect on groundwater levels as a result of potential changes in recharge, groundwater withdrawals, or increased flow in Crab Creek. Decreased recharge and increased groundwater withdrawals both resulted in declines in groundwater levels over 2013 conditions, whereas increasing the flow in Crab Creek resulted in increased groundwater levels over 2013 conditions.
This report documents the addition of shared node Local Grid Refinement (LGR) to MODFLOW-2005, the U.S. Geological Survey modular, transient, three-dimensional, finite-difference ground-water flow model. LGR provides the capability to simulate ground-water flow using one block-shaped higher-resolution local grid (a child model) within a coarser-grid parent model. LGR accomplishes this by iteratively coupling two separate MODFLOW-2005 models such that heads and fluxes are balanced across the shared interfacing boundary. LGR can be used in two-and three-dimensional, steady-state and transient simulations and for simulations of confined and unconfined ground-water systems. Traditional one-way coupled telescopic mesh refinement (TMR) methods can have large, often undetected, inconsistencies in heads and fluxes across the interface between two model grids. The iteratively coupled shared-node method of LGR provides a more rigorous coupling in which the solution accuracy is controlled by convergence criteria defined by the user. In realistic problems, this can result in substantially more accurate solutions and require an increase in computer processing time. The rigorous coupling enables sensitivity analysis, parameter estimation, and uncertainty analysis that reflects conditions in both model grids. This report describes the method used by LGR, evaluates LGR accuracy and performance for two- and three-dimensional test cases, provides input instructions, and lists selected input and output files for an example problem. It also presents the Boundary Flow and Head (BFH) Package, which allows the child and parent models to be simulated independently using the boundary conditions obtained through the iterative process of LGR.
The comparison between a spatially distributed crop index and reported yield was evaluated. The crop index was generated by a crop water balance model that simulates the reduction of crop yield (as a percentage of the potential) due to water deficit. This model is currently operational as a monitoring and forecasting tool for region-wide food security analyses in drought prone countries in Sub-Saharan Africa. The most important input parameters of the model are Africa-wide satellite-based rainfall estimate images and spatially distributed potential evapotranspiration images calculated from 6-hourly numerical weather model data. District level crop yield data from Ethiopia were used to evaluate the performance of the model. Historical sorghum yield data from 1996-1999 were used to evaluate the performance of a seasonal water requirement satisfaction index (WRSI) for sorghum. WRSI values and reported district yield data were significantly correlated. The WRSI model was particularly successful in capturing the response of the crop during a relatively dry year. In districts that never experienced water deficit during the study period, it was possible to infer the magnitude of yield variability that was caused by factors other than water supply.
This report includes a digital, three-dimensional (3-D) faulted hydrostratigraphic model we have constructed to represent the geologic framework of the Edwards aquifer system in the area of San Antonio, northern Bexar County, Texas. The model is based on mapped geologic relationships that reflect the complex structures of the Balcones fault zone, detailed lithologic descriptions and interpretations of about 40 principal wells (and qualified data from numerous other wells), and a conceptual model of the gross geometry of the Edwards Group units derived from prior interpretations of depositional environments and paleogeography. The EarthVision 3-D software included in this publication is provided for use on four platforms: Windows 2000/XP, SGI, SUN, and Linux. Please see the Readme1.txt file for information about starting and using the software on each platform. The 3-D model and viewer software can be downloaded from this Web site, and they are also available on the CD-ROM version of this report, which is available as a sales item from Information Services, U.S. Geological Survey (1-888-ASK-USGS). Earth Vision software is the intellectual property of Dynamic Graphics, Inc., and is used here with their permission. Other than those rights granted to the USGS and third parties to assist them in using this product (any venue), Dynamic Graphics, Inc., reserves all rights in their software. Any appearances of trade, product, or firm names, or private logos, in or on this product do not, and are not intended to, imply endorsement by the U.S. Government.
This report provides an overview of existing hydrologic information describing the quality, quantity, and extent of the major surface-water and groundwater resources in the Cheyenne and Arapaho Tribal jurisdictional area, west-central Oklahoma. Hydrologic information is provided for five major river systems (Cimarron River, North Canadian River, Canadian River, Washita River, and North Fork Red River), two reservoirs (Foss Reservoir and Canton Lake), and eight aquifers consisting of the alluvial aquifers associated with each of the five major river systems and three major bedrock aquifers (Ogallala aquifer, Elk City aquifer, and Rush Springs aquifer). Types of information provided about rivers and reservoirs for the Cheyenne and Arapaho Tribal jurisdictional area include diversion sites and amounts of water allocated and diverted for permitted uses in 2015; treated wastewater discharge sites and amounts discharged in 2015; and characteristics describing water-quality field properties, major ions, nutrients, and selected trace elements. Major ions, nutrients, and selected trace elements are compared to secondary maximum contaminant levels and maximum contaminant levels for finished drinking water. Additionally, statistics are provided describing daily, monthly, and annual streamflow characteristics at 12 U.S. Geological Survey streamgages. Streamflow statistics include the magnitudes and frequencies of floods, base-flow characteristics, and long-term streamflow trends. Types of information provided about the aquifers include amounts of water allocated and pumped for permitted uses in 2015; characteristics of groundwater describing water-quality field properties, major ions, nitrate (measured as nitrogen), and selected trace elements with comparisons to secondary maximum contaminant levels and maximum contaminant levels for finished drinking water; groundwater levels and long-term changes in water levels; and ranges of hydraulic conductivity, aquifer recharge, specific yield, transmissivity, and well yields from reports and groundwater-flow models. Surface water is used primarily for irrigation and mining and other nonconsumptive uses in the Cheyenne and Arapaho Tribal jurisdictional area, except from the Washita and North Fork Red Rivers, where water is treated for use as a public-water supply. Large concentrations of dissolved solids are the primary limiting factor affecting the use of surface water. Median concentrations of dissolved solids in surface water range from less than 1,000 milligrams per liter (mg/L) in samples from the North Canadian River to greater than 9,000 mg/L in samples from the Cimarron River. Large dissolved solids concentrations are correlated with hard water. Median hardness as calcium carbonate concentrations in surface water ranges from 427 mg/L in samples from Canton Lake to 1,000 mg/L in samples from the Washita River. In 2015, groundwater was used at more than twice the rate of surface water in the Cheyenne and Arapaho Tribal jurisdictional area. Although the alluvial aquifers are considered reliably good sources of water in the Cheyenne and Arapaho Tribal jurisdictional area, concentrations of nitrate (measured as nitrogen) exceed the maximum contaminant level of 10 mg/L established by the U.S. Environmental Protection Agency for finished drinking water in parts of all of the alluvial aquifers. Water from the three major bedrock aquifers is used for irrigation, mining, public-water supply, and other uses; however, large concentrations of dissolved solids, nitrate (measured as nitrogen), and naturally occurring trace elements such as arsenic and uranium may limit the use of groundwater as a source of public-water supply in some areas. As of 2015, the depletion of groundwater from the major aquifers in west-central Oklahoma is a minor concern to the Oklahoma Water Resources Board. Groundwater levels and other hydrologic information show that recharge rates exceed the rates of water pumped from aquifers, except in areas that may be affected locally by groundwater depletions.
This report is one in a series resulting from the U.S. Geological Survey's Snake River Plain Rasa (Regional Aquifer System Analysis) study that was initiation in October 1979. Purposes of the RASA study were to (1) refine knowledge of the regional ground-water flow system, (2) determine effects of conjunctive use of ground water and surface water, and (3) describe water chemistry. This purpose of this report is to describe flow characteristics of the Snake River and tributaries and to develop a water budget for the Snake River Plain. Data and interpretation in this report will be used in development of ground-water flow models.
Streamflow and water quality vary spatially and temporally in the Rio Grande from Del Norte, Colorado, to El Paso, Texas. The variations in streamflow and in concentrations of selected waterquality constituents—dissolved solids, dissolved nitrite plus nitrate as nitrogen, total phosphorus, and suspended sediment—are described in this report. A multivariate linear regression model, ESTIMATOR2000, was used to estimate loads for selected constituents. Streamflow decreases in the downstream direction throughout most of the basin because outflows (due to agricultural use, leakage to ground water, and evapotranspiration) are greater than inflows. Streamflow increases between Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, to Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, because ground-water and tributary inflow are greater than outflow. Concentrations of dissolved solids, dissolved nitrite plus nitrate, total phosphorus, and suspended sediment generally increase in the downstream direction. Concentrations of dissolved solids, dissolved nitrite plus nitrate, and total phosphorus decrease between Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, and Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, because of dilution by tributary inflow. Concentrations of dissolved nitrite plus nitrate, total phosphorus, and suspended sediment decrease between Rio Grande Floodway at San Marcial, New Mexico, and Rio Grande below Leasburg Dam, near Leasburg, New Mexico, because of reservoir effects (nutrient uptake and settling of sediment). Several instances of decreasing streamflow and increasing loads indicate the presence of inflows with large constituent concentrations (relative to those of the Rio Grande immediately upstream from that inflow); this occurs (1) between Rio Grande near Del Norte, Colorado, and Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, for dissolved solids, (2) between Rio Grande at Otowi Bridge, near San Ildefonso, New Mexico, and Rio Grande Floodway at San Marcial, NewMexico, for all constituents, and (3) between Rio Grande below Leasburg Dam, near Leasburg, New Mexico, and Rio Grande at El Paso, Texas, for all constituents. Streamflow increases along every reach of the Rio Grande between the streamflow-gaging station Rio Grande above the mouth of Trinchera Creek, near Lasauses, Colorado, and the station Rio Grande at Otowi Bridge, near San Ildefonso, NewMexico. These increases in streamflow result in increases in the loads of dissolved solids, total phosphorus, and suspended sediment regardless of changes in concentrations.