USGS Science⌕ Search

SEARCH · USGS Science

Results for “Modeling Report”

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.

At least 613 records · Page 34Linked to original sources
USGS2009

Groundwater-flow model of the Ozark Plateaus aquifer system, northwestern Arkansas, southeastern Kansas, southwestern Missouri, and northeastern Oklahoma

To assess the effect that increased water use is having on the long-term availability of groundwater within the Ozark Plateaus aquifer system, a groundwater-flow model was developed using MODFLOW 2000 for a model area covering 7,340 square miles for parts of Arkansas, Kansas, Missouri, and Oklahoma. Vertically the model is divided into five units. From top to bottom these units of variable thickness are: the Western Interior Plains confining unit, the Springfield Plateau aquifer, the Ozark confining unit, the Ozark aquifer, and the St. Francois confining unit. Large mined zones contained within the Springfield Plateau aquifer are represented in the model as extensive voids with orders-of-magnitude larger hydraulic conductivity than the adjacent nonmined zones. Water-use data were compiled for the period 1960 to 2006, with the most complete data sets available for the period 1985 to 2006. In 2006, total water use from the Ozark aquifer for Missouri was 87 percent (8,531,520 cubic feet per day) of the total pumped from the Ozark aquifer, with Kansas at 7 percent (727,452 cubic feet per day), and Oklahoma at 6 percent (551,408 cubic feet per day); water use for Arkansas within the model area was minor. Water use in the model from the Springfield Plateau aquifer in 2005 was specified from reported and estimated values as 569,047 cubic feet per day. Calibration of the model was made against average water-level altitudes in the Ozark aquifer for the period 1980 to 1989 and against waterlevel altitudes obtained in 2006 for the Springfield Plateau and Ozark aquifers. Error in simulating water-level altitudes was largest where water-level altitude gradients were largest, particularly near large cones of depression. Groundwater flow within the model area occurs generally from the highlands of the Springfield Plateau in southwestern Missouri toward the west, with localized flow occurring towards rivers and pumping centers including the five largest pumping centers near Joplin, Missouri; Carthage, Missouri; Noel, Missouri; Pittsburg, Kansas; and Miami, Oklahoma. Hypothetical scenarios involving various increases in groundwater-pumping rates were analyzed with the calibrated groundwater-flow model to assess changes in the flow system from 2007 to the year 2057. Pumping rates were increased between 0 and 4 percent per year starting with the 2006 rates for all wells in the model. Sustained pumping at 2006 rates was feasible at the five pumping centers until 2057; however, increases in pumping resulted in dewatering the aquifer and thus pumpage increases were not sustainable in Carthage and Noel for the 1 percent per year pumpage increase and greater hypothetical scenarios, and in Joplin and Miami for the 4 percent per year pumpage increase hypothetical scenarios. Zone-budget analyses were performed to assess the groundwater flow into and out of three zones specified within the Ozark-aquifer layer of the model. The three zones represented the model part of the Ozark aquifer in Kansas (zone 1), Oklahoma (zone 2), and Missouri and Arkansas (zone 3). Groundwater pumping causes substantial reductions in water in storage and induces flow through the Ozark confining unit for all hypothetical scenarios evaluated. Net simulated flow in 2057 from Kansas (zone 1) to Missouri (zone 3) ranges from 74,044 cubic feet per day for 2006 pumping rates (hypothetical scenario 1) to 625,319 cubic feet per day for a 4 percent increase in pumping per year (hypothetical scenario 5). Pumping from wells completed in the Ozark aquifer is the largest component of flow out of zone 3 in Missouri and Arkansas, and varies between 88 to 91 percent of the total flow out of zone 3 for all of the hypothetical scenarios. The largest component of flow into Oklahoma (zone 2) comes from the overlying Ozark confining unit, which is consistently about 45 percent of the total. Flow from the release of water in storage, from general-head boundaries, and from zones 1 and 3 is considerably smaller values that range from 3 to 22 percent of the total flow into zone 2. The largest flow out of the Oklahoma part of the model occurs from pumping from wells and ranges from 52 to 69 percent of the total.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS1991

Use of the Priestley-Taylor evaporation equation for soil water limited conditions in a small forest clearcut

The Priestley-Taylor equation, a simplification of the Penman equation, was used to allow calculations of evapotranspiration under conditions where soil water supply limits evapotranspiration. The Priestley-Taylor coefficient, α, was calculated to incorporate an exponential decrease in evapotranspiration as soil water content decreases. The method is appropriate for use when detailed meteorological measurements are not available. The data required to determine the parameter for the α coefficient are net radiation, soil heat flux, average air temperature, and soil water content. These values can be obtained from measurements or models. The dataset used in this report pertains to a partially vegetated clearcut forest site in southwest Oregon with soil depths ranging from 0.48 to 0.70 m and weathered bedrock below that. Evapotranspiration was estimated using the Bowen ratio method, and the calculated Priestley-Taylor coefficient was fitted to these estimates by nonlinear regression. The calculated Priestley-Taylor coefficient (α′) was found to be approximately 0.9 when the soil was near field capacity (0.225 cm 3 cm −3 ). It was not until soil water content was less than 0.14 cm 3 cm −3 that soil water supply limited evapotranspiration. The soil reached a final residual water content near 0.05 cm 3 cm −3 at the end of the growing season.

Agricultural and Forest Meteorology↗
USGS2022

Greenhouse gas balances in coastal ecosystems: Current challenges in “blue carbon” estimation and significance to national greenhouse gas inventories

Coastal wetlands are defined herein as inundated, vegetated ecosystems with hydrology, and biogeochemistry influenced by sea levels, at timescales of tides to millennia. Coastal wetlands are necessary components of global greenhouse gas estimation and scenario modeling, both for continental and oceanic mass balances. The carbon pools and fluxes on coastal lands, especially those influenced by tidal drivers and sea level rise, are distinct in their magnitude, rates, and uncertainties. We describe herein the pathways taken for a US scale estimation of blue carbon based on annual timesteps and bottom-up modeling, as appropriate for the first effort to include coastal wetlands in the Intergovernmental Panel on Climate Change (IPCC) guidelines for a National Greenhouse Gas Inventory (NGGI). As such, we summarize multiple efforts to reconcile mapping, modeling, and measurement issues and we report the assumptions we made based on data availability. Provided as requested feedback to the IPCC. Subsidiary Body for Scientific and Technological Advice (SBSTA) evaluation of guidance criteria, these analyses synergistically point scientists, practitioners, and policy makers toward the greatest uncertainties to address in future assessments: coastal wetland methane emissions and carbon dioxide emissions associated with the fate of eroded soil. This is a story of what was learned in the 2014–2018 NASA Carbon Monitoring System project (https://carbon.nasa.gov/cgi-bin/cms_projects.pl), how it informs “good practice” (IPCC 2006) in reporting coastal wetland emissions and removals, and where it points scientifically toward data needs at different temporal and spatial scales.

Book chapter↗
USGS2021

Establishment of baseline cytology metrics in nestling American kestrels (Falco sparverius): Immunomodulatory effects of the flame retardant isopropylated triarylphosphate isomers

Avian populations must mount effective immune responses upon exposure to environmental stressors such as avian influenza and xenobiotics. Although multiple immune assays have been tested and applied to various avian species, antibody-mediated immune responses in non-model avian species are not commonly reported due to the lack of commercially available species-specific antibodies. The objectives of the present study were to advance methods for studying wild bird immune responses and to apply these to the evaluation of cytological responses after exposure of American kestrels, Falco sparverius, to a commercial flame retardant mixture containing isopropylated triarylphosphate isomers (ITP). Hatchlings were gavaged daily with safflower oil or 1.5 ug/g bw/day of ITP suspended in safflower oil, then bled on days 9, 17, and 21. The ITP treatment group ( n = 18) and a subset of controls (Poly I:C treatment group; n = 10) were injected on days 9 and 15 with a synthetic analog of viral double-stranded RNA, polyinosinic:polycytidylic acid (Poly I:C), a toll-like receptor ligand and synthetic viral mimic, and responses compared to a sham injected control group (n = 8). The hypotheses tested whether kestrels showed immunological differences among treatment groups, genetic sex, and/or white blood cell (WBC) subpopulation type over time. A flow cytometry (FCM) gating strategy categorized heterophils (H), lymphocytes (L), and monocytes (M) and their proportions, and measured relative fluorescence in response to anti-chicken CD4 binding. Fluorescent cell surfaces and some granular/vacuolar inclusions were visualized by epifluorescence microscopy. A fourth subpopulation with higher levels of granularity than M but less than H became increasingly apparent with time and was gated along with the H subpopulation; its frequency of occurrence was lowest in the ITP group ( P = 0.0023). The percentages of cells differed among treatment groups, days, and sexes ( P = 0.0001). For both sexes, percentages of H and L were higher than M in control and Poly I:C. In the ITP group, L percentages were higher than H and M ( P = 0.0457), and H and L were higher than M on days 9 and 21 ( P = 0.0001). The ratios of H:L and H:WBC, indicators of robust immunity, were also higher on days 9 and 21 than on 17 ( P = 0.0079). For each sex, the highest levels of activity measured by FCM geometric means (GEO) of fluorescence (indicative of antibody binding) were observed on day 9 ( P = 0.0001 female, and P = 0.0011 male) in H over both L and M ( P < 0.0001 for each). In males, GEO of the Poly I:C group was higher than that of the ITP group ( P = 0.0374), with no difference observed among females over all days. By using a FCM algorithm for population comparisons of fluorescence to investigate binding within H, the T(x) scores indicated higher fluorescence in control and Poly I:C groups over ITP ( P = 0.0001). Unlike chickens, Gallus gallus , which express CD4 primarily on L, kestrels bound the commercial antibody primarily within the gated H subpopulation, suggesting an immunophenotypic difference between taxa, despite a ~60% identity of Falco CD4 amino acid sequences with chicken CD4. The emergent cell subset within the gated H presented dendritic-like cell (DLC) morphological and functional properties, apparently serving as an effector cell. This study adds interpretive context to ecological investigations of infection and of potential immunomodulation by emerging compounds, whereby the early innate responses are mediated by the various cell subsets serving as useful quantitative markers of immunological condition. Data showed that dietary exposure to ITP was immunosuppressive for male and female kestrels over the course of the experiment, reducing DLC frequency compared to the Poly I:C controls. Heterophils and DLC were important in facilitating innate immunological responses.

Environment International↗
USGS2009

Thermal characteristics of amphibian microhabitats in a fire-disturbed landscape

Disturbance has long been a central issue in amphibian conservation, often regarding negative effects of logging or other forest management activities, but some amphibians seem to prefer disturbed habitats. After documenting increased use of recently burned forests by boreal toads (Bufo boreas), we hypothesized that burned habitats provided improved thermal opportunities in terrestrial habitats. We tested this hypothesis by conducting a radio telemetry study of habitat use (reported previously) and by using physical models that simulated the temperature of adult toads. We deployed 108 physical models in and adjacent to a 1-year old burn using a fully-replicated design with three burn severities (unburned, partial, high severity) and four microhabitats (open surface, under vegetation, under log, in burrow). Model temperatures were compared to a range of preferred temperatures in published studies. We found 70% more observations within the preferred temperature range of B. boreas in forests burned with high severity than in unburned areas. Burned forest was warmer than unburned forest across all microhabitats, but the largest relative difference was in burrows, which averaged 3 ??C warmer in high-severity burn areas and remained warmer though the night. More than twice as many observations were within the preferred temperature range in high-severity burrows than in unburned burrows. Areas burned with high severity were still warmer than unburned forest 3 years after the fire. Habitat use of toads during the concurrent radio telemetry study matched that predicted by the physical models. These results suggest there are fitness-linked benefits to toads using burned habitats, such as increased growth, fertility, and possibly disease resistance. However, increased soil temperatures that result from wildfire may be detrimental to other amphibian species that prefer cooler temperatures and stable environments. More broadly, our data illustrate the use of physical models to measure and interpret changes that amphibians may experience from disturbance, and highlight the need for research linking vital rates such as growth and survival to disturbance.

Forest Ecology and Management↗
USGS2013

Rapid chemical evolution of tropospheric volcanic emissions from Redoubt Volcano, Alaska, based on observations of ozone and halogen-containing gases

We report results from an observational and modeling study of reactive chemistry in the tropospheric plume emitted by Redoubt Volcano, Alaska. Our measurements include the first observations of Br and I degassing from an Alaskan volcano, the first study of O 3 evolution in a volcanic plume, as well as the first detection of BrO in the plume of a passively degassing Alaskan volcano. This study also represents the first detailed spatially-resolved comparison of measured and modeled O 3 depletion in a volcanic plume. The composition of the plume was measured on June 20, 2010 using base-treated filter packs (for F, Cl, Br, I, and S) at the crater rim and by an instrumented fixed-wing aircraft on June 21 and August 19, 2010. The aircraft was used to track the chemical evolution of the plume up to ~ 30 km downwind (2 h plume travel time) from the volcano and was equipped to make in situ observations of O 3 , water vapor, CO 2 , SO 2 , and H 2 S during both flights plus remote spectroscopic observations of SO 2 and BrO on the August 19th flight. The airborne data from June 21 reveal rapid chemical O 3 destruction in the plume as well as the strong influence chemical heterogeneity in background air had on plume composition. Spectroscopic retrievals from airborne traverses made under the plume on August 19 show that BrO was present ~ 6 km downwind (20 min plume travel time) and in situ measurements revealed several ppbv of O 3 loss near the center of the plume at a similar location downwind. Simulations with the PlumeChem model reproduce the timing and magnitude of the observed O 3 deficits and suggest that autocatalytic release of reactive bromine and in-plume formation of BrO were primarily responsible for the observed O 3 destruction in the plume. The measurements are therefore in general agreement with recent model studies of reactive halogen formation in volcanic plumes, but also show that field studies must pay close attention to variations in the composition of ambient air entrained into volcanic plumes in order to unambiguously attribute observed O 3 anomalies to specific chemical or dynamic processes. Our results suggest that volcanic eruptions in Alaska are sources of reactive halogen species to the subarctic troposphere.

Alaska↗
USGS2025

What is the lowest latitude of discrete aurorae during superstorms?

From a survey of published accounts of visual sightings of aurorae, a compilation is presented of the lowest identified geomagnetic latitude at which discrete aurorae were seen at local zenith during magnetic storms having intensities with maximum − Dst > 200 nT. The compilation includes data for the superstorms of 2 September 1859, 4 February 1872, and 15 May 1921. A statistical model is developed representing the equatorward boundary of discrete aurorae versus storm intensity. The model indicates that a once-per-century storm would likely induce discrete aurorae at zenith down to a geomagnetic latitude of 34 ° . Insofar as aurorae can be taken as a proxy for electrojet currents, such a storm would expose many nighttime electric-power systems, in the contiguous United States or Europe, to high levels of geomagnetic disturbance. A Carrington-class storm would induce discrete aurorae down to 24 ° . These exposures are much greater than those indicated in recent numerical simulations of extreme magnetic storms. Using the model to infer storm intensity from reports of low-latitude aurorae, a storm on 28 August 1859, likely had maximum − Dst = 673 nT. That this storm occurred just a few days before the Carrington storm of 2 September (maximum − Dst = 964 nT) deserves attention. A storm that occurred on 17 September 1770 is estimated to have had maximum − Dst = 928 nT. The vision of Ezekiel could have been inspired by aurorae from a storm with maximum − Dst = 550 nT.

Space Weather↗
USGS2006

Continental crustal evolution observations

How has the continental crust evolved? What are the primary processes responsible for its composition, structure, and mode of deformation? What role do fluids play in deep crustal processes? In the last dozen years, geophysicists have obtained images of the deep continental crust that can be used to examine these questions and refine geologic models of crustal evolution. In this report we summarize recent progress in geophysical studies of the deep continental crust and highlight some of the more important implications of deep crustal processes.

Eos, Earth and Space Science News↗
USGS2019

Persistence of intense, climate-driven runoff late in Mars history

Mars is dry today, but numerous precipitation-fed paleo-rivers are found across the planet’s surface. These rivers’ existence is a challenge to models of planetary climate evolution. We report results indicating that, for a given catchment area, rivers on Mars were wider than rivers on Earth today. We use the scale (width and wavelength) of Mars paleo-rivers as a proxy for past runoff production. Using multiple methods, we infer that intense runoff production of >(3–20) kg/m 2 per day persisted until <3 billion years (Ga) ago and probably <1 Ga ago, and was globally distributed. Therefore, the intense runoff production inferred from the results of the Mars Science Laboratory rover was not a short-lived or local anomaly. Rather, precipitation-fed runoff production was globally distributed, was intense, and persisted intermittently over >1 Ga. Our improved history of Mars’ river runoff places new constraints on the unknown mechanism that caused wet climates on Mars.

Science Advances↗
USGS2020

Dendritic reidite from the Chesapeake Bay impact horizon, Ocean Drilling Program Site 1073 (offshore northeastern USA): A fingerprint of distal ejecta?

High-pressure minerals provide records of processes not normally preserved in Earth’s crust. Reidite, a quenchable polymorph of zircon, forms at pressures >20 GPa during shock compression. However, there is no broad consensus among empirical, experimental, and theoretical studies on the nature of the polymorphic transformation. Here we decipher a multistage history of reidite growth recorded in a zircon grain in distal impact ejecta (offshore northeastern United States) from the ca. 35 Ma Chesapeake Bay impact event which, remarkably, experienced near-complete conversion (89%) to reidite. The grain displays two distinctive reidite habits: (1) intersecting sets of planar lamellae that are dark in cathodoluminescence (CL); and (2) dendritic epitaxial overgrowths on the lamellae that are luminescent in CL. While the former is similar to that described in literature, the latter has not been previously reported. A two-stage growth model is proposed for reidite formation at >40 GPa in Chesapeake Bay impact ejecta: formation of lamellar reidite by shearing during shock compression, followed by dendrite growth, also at high pressure, via recrystallization. The dendritic reidite is interpreted to nucleate on lamellae and replace damaged zircon adjacent to lamellae, which may be amorphous ZrSiO 4 or possibly an intermediate phase, all before quenching. These results provide new insights on the microstructural evolution of the high-pressure polymorphic transformation over the microseconds-long interval of reidite stability during meteorite impact. Given the formation conditions, dendritic reidite may be a unique indicator of distal ejecta.

Virginia↗
USGS2010

Normal-faulting slip maxima and stress-drop variability: A geological perspective

We present an empirical estimate of maximum slip in continental normal-faulting earthquakes and present evidence that stress drop in intraplate extensional environments is dependent on fault maturity. A survey of reported slip in historical earthquakes globally and in latest Quaternary paleoearthquakes in the Western Cordillera of the United States indicates maximum vertical displacements as large as 6–6.5 m. A difference in the ratio of maximum-to-mean displacements between data sets of prehistoric and historical earthquakes, together with constraints on bias in estimates of mean paleodisplacement, suggest that applying a correction factor of 1.4±0.3 to the largest observed displacement along a paleorupture may provide a reasonable estimate of the maximum displacement. Adjusting the largest paleodisplacements in our regional data set (∼6 m) by a factor of 1.4 yields a possible upper-bound vertical displacement for the Western Cordillera of about 8.4 m, although a smaller correction factor may be more appropriate for the longest ruptures. Because maximum slip is highly localized along strike, if such large displacements occur, they are extremely rare. Static stress drop in surface-rupturing earthquakes in the Western Cordillera, as represented by maximum reported displacement as a fraction of modeled rupture length, appears to be larger on normal faults with low cumulative geologic displacement (<2 km) and larger in regions such as the Rocky Mountains, where immature, low-throw faults are concentrated. This conclusion is consistent with a growing recognition that structural development influences stress drop and indicates that this influence is significant enough to be evident among faults within a single intraplate environment.

Bulletin of the Seismological Society of America↗
USGS2011

Natural phosphorus sources for the Pacific Northwest

Phosphorus is a naturally occurring element found in all rocks; the amount varies by the type of rock. The amount of phosphorus in sediments is expected to be correlated with the amount of phosphorus in the parent rocks. Streambed sediment collected by the National Uranium Resource Evaluation (NURE) Program were used to estimate the variation of phosphorus across the Pacific Northwest. This file provides an estimate of the mean concentration of phosphorus in soils for each incremental catchment of the USGS Pacific Northwest SPARROW model.

Idaho;Oregon;Washington↗
USGS2001

Summary of trends and status analysis for flow, nutrients, and sediments at selected nontidal sites, Chesapeake Bay basin, 1985-99

Water-quality and flow data from 31 sites in nontidal portions of the Chesapeake Bay Basin were analyzed to document annual nutrient and sediment loads and trends for the period 1985 through 1999 as part of an annual reevaluation and reporting for the Chesapeake Bay Program. Annual loads were estimated by use of the U.S. Geological Survey ESTIMATOR model. Trends were estimated using linear regression. Trends were reported for monthly mean flow, monthly load, flow-adjusted concentration, and flow-weighted concentration. Median yields and concentrations were calculated to help facilitate comparisons between basins. The drought of 1999 had pronounced effects on trend results. The trend in flow increased at 4 of the 31 sites, 8 fewer sites than in 1998. Ten less significant trends were estimated for nutrient and sediment loads compared to 1985-98. Trends in flow-weighted and flow-adjusted concentrations varied little by nutrient species and geographic location. Trends were generally downward or not significant for both the nitrogen and phosphorus species throughout the Chesapeake Bay Basin. Trends in flow-adjusted concentration indicated downward trends at most sites for nutrients and about half the sites for sediments, an indication that management actions are reducing nutrient and sediment concentrations.

Open-File Report↗
USGS2002

Progress on geoenvironmental models for selected mineral deposit types

Since the beginning of economic geology as a subdiscipline of the geological sciences, economic geologists have tended to classify mineral deposits on the basis of geological, mineralogical, and geochemical criteria, in efforts to systematize our understanding of mineral deposits as an aid to exploration. These efforts have led to classifications based on commodity, geologic setting (Cox and Singer, 1986), inferred temperatures and pressures of ore formation (Lindgren, 1933), and genetic setting (Park and MacDiarmid, 1975; Jensen and Bateman, 1979). None of these classification schemes is mutually exclusive; instead, there is considerable overlap among all of these classifications. A natural outcome of efforts to classify mineral deposits is the development of &ldquo;mineral deposit models.&rdquo; A mineral deposit model is a systematically arranged body of information that describes some or all of the essential characteristics of a selected group of mineral deposits; it presents a concept within which essential attributes may be distinguished and from which extraneous, coincidental features may be recognized and excluded (Barton, 1993). Barton (1993) noted that the grouping of deposits on the basis of common characteristics forms the basis for a classification, but the specification of the characteristics required for belonging to the group is the basis for a model. Models range from purely descriptive to genetic. A genetic model is superior to a descriptive model because it provides a basis to distinguish essential from extraneous attributes, and it has flexibility to accommodate variability in sources, processes, and local controls. In general, a descriptive model is a necessary prerequisite to a genetic model.

Open-File Report↗
USGS2004

Death Valley regional ground-water flow system, Nevada and California -- hydrogeologic framework and transient ground-water flow model

A numerical three-dimensional (3D) transient ground-water flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the ground-water flow system and previous less extensive ground-water flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect ground-water flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley Regional Ground-Water Flow System (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the ground-water flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural ground-water discharge occurring through evapotranspiration and spring flow; the history of ground-water pumping from 1913 through 1998; ground-water recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional ground-water flow system. Ground-water flow in the Death Valley region is composed of several interconnected, complex ground-water flow systems. Ground-water flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional ground-water flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and evapotranspiration (ET) are the dominant natural ground-water discharge processes. Ground water also is withdrawn for agricultural, commercial, and domestic uses. Ground-water flow in the DVRFS was simulated using MODFLOW-2000, a 3D finite-difference modular ground-water flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 m on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins, also simulated by drains. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model gener

Scientific Investigations Report↗
USGS2013

Application of sediment characteristics and transport conditions to resource management in selected main-stem reaches of the Upper Colorado River, Colorado and Utah, 1965-2007

The Colorado River Basin provides habitat for 14 native fish, including 4 endangered species protected under the Federal Endangered Species Act of 1973. These endangered fish species once thrived in the Colorado River system, but water-resource development, including the building of numerous diversion dams and several large reservoirs, and the introduction of non-native fish, resulted in large reductions in the numbers and range of the four species through loss of habitat and stream function. Understanding how stream conditions and habitat change in response to alterations in streamflow is important for water administrators and wildlife managers and can be determined from an understanding of sediment transport. Characterization of the processes that are controlling sediment transport is an important first step in identifying flow regimes needed for restored channel morphology and the sustained recovery of endangered fishes within these river systems. The U.S. Geological Survey, in cooperation with the Upper Colorado River Endangered Fish Recovery Program, Bureau of Reclamation, U.S. Fish and Wildlife Service, Argonne National Laboratory, Western Area Power Administration, and Wyoming State Engineer’s Office, began a study in 2004 to characterize sediment transport at selected locations on the Colorado, Gunnison, and Green Rivers to begin addressing gaps in existing datasets and conceptual models of the river systems. This report identifies and characterizes the relation between streamflow (magnitude and timing) and sediment transport and presents the findings through discussions of (1) suspended-sediment transport, (2) incipient motion of streambed material, and (3) a case study of sediment-transport conditions for a reach of the Green River identified as a razorback sucker spawning habitat (See report for full abstract).

Colorado;Utah↗