Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Names”

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 1,765 records · Page 98Linked to original sources

Hydrogeologic characterization of the Yucaipa groundwater subbasin

Introduction Water management in the Santa Ana River watershed in San Bernardino and Riverside Counties in southern California (fig. A1) is complex with various water purveyors navigating geographic, geologic, hydrologic, and political challenges to provide a reliable water supply to stakeholders. As the population has increased throughout southern California, so has the demand for water. The Yucaipa groundwater subbasin (hereafter referred to as “Yucaipa subbasin”), one of nine groundwater subbasins in what the California Department of Water Resources (DWR) refers to as the Upper Santa Ana Valley groundwater basin (California Department of Water Resources, 2016; fig. A1; the DWR naming convention is used within this report), is no exception; steady population growth since the 1940s and changes in water use have forced local water purveyors to regularly adapt their water infrastructure. Water demands within the Yucaipa subbasin have historically been supplied by groundwater, but water imported via the California State Water Project has augmented the total water supply through direct use and through anthropogenic recharge at the Wilson Creek and Oak Glen Creek spreading basins since 2002. Overall demand for groundwater continues to rise, and local water managers are concerned that despite the influx of imported water, groundwater levels may decline to a point where producing water will be uneconomical, severely limiting the ability of local agencies to meet water-supply demand. To better understand the hydrogeology and water resources in the Yucaipa subbasin, the U.S. Geological Survey (USGS) initiated a study in cooperation with the San Bernardino Valley Municipal Water District (SBVMWD) to characterize and model the hydrologic system of the Yucaipa subbasin and the surrounding Yucaipa Valley watershed (YVW; fig. A2). To gain this comprehensive understanding, a three-dimensional (3D) hydrogeologic framework model (HFM; Cromwell and Matti, 2022) was constructed to quantify the structure and extent of hydrogeologic units in the YVW; the hydrologic system was conceptualized and quantified (described in chapter A); and the Yucaipa Integrated Hydrological Model (YIHM; described in chapter B ) was developed to simulate the integrated surface-water and aquifer systems, including natural and anthropogenic recharge and discharge throughout the study area during 1947–2014.

Scientific Investigations Report↗

Geology and hydrogeology of the Yucaipa groundwater subbasin, San Bernardino and Riverside Counties, California

The Yucaipa groundwater subbasin (referred to in this report as the Yucaipa subbasin) is located about 75 miles (mi) east of of Los Angeles and about 12 mi southeast of the City of San Bernardino. In the Yucaipa subbasin, as in much of southern California, limited annual rainfall and large water demands can strain existing water supplies; therefore, understanding local surface water and groundwater conditions is essential for managing these resources. To better understand the hydrogeology and water resources in the Yucaipa subbasin, especially groundwater, the San Bernardino Valley Municipal Water District and the U.S. Geological Survey initiated a cooperative study to evaluate the hydrogeologic system of the Yucaipa subbasin and the encompassing Yucaipa Valley watershed. Previous studies of the area provided information on general geologic and hydrologic conditions, but this study provides the first comprehensive definition of the hydrogeology of the subsurface throughout the entire subbasin. The Yucaipa subbasin is located between the northwest trending San Andreas fault zone and San Jacinto fault. Several northeast-trending dip-slip faults dissect the Yucaipa subbasin, providing the mechanism for structural relief within the sediment-filled subbasin and between the subbasin and surrounding mountains and highlands. Several of these dip-slip faults have been previously identified as potential barriers to groundwater flow. This report provides a synthesis of previous studies and a discussion of the geologic interpretations that were used as the foundation for hydrogeologic classification of the Yucaipa subbasin. Notably, this report (1) adopts the recently named and classified sedimentary deposits of Live Oak Canyon geologic formation and extends the mapped distribution of the formation into the Yucaipa subbasin, and (2) adopts the interpretation that activity along the Banning fault predates the deposition of most basin-fill sedimentary materials in the Yucaipa subbasin. Four hydrogeologic units were classified in the Yucaipa subbasin: (1) crystalline basement, (2) consolidated sedimentary materials, (3) unconsolidated sediment, and (4) surficial materials. The crystalline basement unit forms the bottom boundary of the aquifer system, and the three other units comprise the basin-fill aquifer system. The four hydrogeologic units vary in extent, thickness, and structural relief across the subbasin, with the unconsolidated sediment unit serving as the primary aquifer unit. A three-dimensional hydrogeologic framework model was developed for the Yucaipa subbasin and surrounding area to characterize the thickness, extent, and hydrogeologic variability of the aquifer system. Geologic maps, borehole geophysical logs, drillers’ lithology logs, and depth-to-basement gravity data were used to map and interpolate the subsurface extent and structure of the hydrogeologic units within the subbasin. Faults and structures of geologic and (or) hydrogeologic importance were included in the model for future evaluation of their potential effects on groundwater flow. The resulting hydrogeologic framework is consistent with existing geologic concepts and the tectonic and structural history of the Yucaipa subbasin and surrounding area. The framework is also suitable for use in basin-scale hydrogeologic investigations.

California↗

Revisions to the Virginia Coastal Plain hydrogeologic framework southwest of the James River

New drilling information reveals that altitudes of some hydrogeologic units of the Virginia Coastal Plain aquifer system differ by as much as 50 feet (ft) from those previously known, namely the Aquia and Potomac aquifers, the Potomac confining zone, and the Nanjemoy-Marlboro and Saint Marys confining units. In addition, the lateral margins of some hydrogeologic units are located as much as several miles from previously estimated locations. The largest revisions to unit margins were for the Aquia aquifer and the Nanjemoy-Marlboro and Saint Marys confining units. Interpretation of new geophysical logs, sediment core, and cuttings as well as revised interpretations to existing data indicate channels and embayments are also preserved on eroded top surfaces of the shallowest hydrogeologic units including the Yorktown confining zone, Yorktown-Eastover aquifer, Saint Marys confining unit, Potomac confining zone, and Potomac aquifer. Enhanced details on the configuration of part of the aquifer system southwest of the James River are provided by sediment cores and cuttings as well as geophysical logs from 36 recently drilled boreholes. These, along with reinterpretation of data from 93 preexisting boreholes, form the basis for revised top-surface altitudes and margins of hydrogeologic units beneath parts of Prince George, Surry, Sussex, Isle of Wight, and Southampton Counties and the cities of Franklin and Suffolk. Groundwater withdrawals in the Virginia Coastal Plain cause widespread water-level declines, create the potential for saltwater intrusion, and contribute to regionwide land subsidence. A description of the aquifer system, termed a hydrogeologic framework, was developed by the U.S. Geological Survey in 2006 and provides information needed to base withdrawal-permitting decisions by the Virginia Department of Environmental Quality. This revision of part of the hydrogeologic framework southwest of the James River is based on interpretations of both new and previously analyzed borehole data. The revision is strictly confined to the study area extent and hydrogeologic units not found within the study area were not revised and are not discussed in this report. The newly determined hydrogeologic-unit altitudes and margins have implications for groundwater-withdrawal permitting. New interpretations have found that the Yorktown Eastover aquifer is absent in the southwestern part of the City of Suffolk, owing to what is most likely an isolated area of sediment-texture facies change. Most notably, the top-surface altitudes of the Aquia and Potomac aquifers have been lowered by as much as 50 ft from previous interpretations. This means that wells previously believed to be screened in the top of the Potomac aquifer could, based on these new interpretations, be screened in the bottom of the Aquia aquifer. These changes to aquifers in which wells are screened means that there is potentially more room in the groundwater withdrawal permitting for the Potomac aquifer, the largest and most productive aquifer in Virginia, and overpumping occurring in the Aquia aquifer.

Virginia↗

Virginia Bridge Scour Pilot Study—Hydrological Tools

Hydrologic and geophysical components interact to produce streambed scour. This study investigates methods for improving the utility of estimates of hydrologic flow in streams and rivers used when evaluating potential pier scour over the design-life of highway bridges in Virginia. Recent studies of streambed composition identify potential bridge design cost savings when attributes of cohesive soil and weathered rock unique to certain streambeds are considered within the bridge planning design. To achieve potential cost savings, however, attributes and effects of scour forces caused by water movement across the streambed surface must be accurately described and estimated. This study explores the potential for improving estimates of the hydrologic component, namely hydrologic flow, afforded by empirically based deterministic, probabilistic, and statistical modeling of flows using streamgage data from 10 selected sites in Virginia. Methods are described and tools are provided that may assist with estimating hydrological components of flow duration and potential cumulative stream power for bridge designs in specific settings, and calculation of comprehensive projections of anticipated individual bridge pier scour rates. Examples of hydrologic properties needed to determine the rates of streambed scour are described for sites spanning a range of basin sizes and locations in Virginia. Deterministic, probabilistic, and statistical modeling methods are demonstrated for estimating hydrological components of streambed scour over a bridge design lifespan. Eight tools provide examples of streamflow analysis using daily and instantaneous streamflow data collected at 10 study sites in Virginia. Tool 1 provides a generalized system dynamics model of streamflow and sediment motion that may be used to estimate hydrologic flow over time. Tool 2 illustrates at-a-station hydraulic geometry using methods pioneered by Leopold and others. Tool 3 provides a system dynamics model developed to test the use of Monte-Carlo sampling of instantaneous streamflow measurements to augment and increase precision of site-specific period-of-record daily-flow values useful for driving stream-power and streambed scour estimates. Tool 4 integrates deterministic modeling, maximum likelihood logistic regression, and Monte-Carlo sampling to identify probable hydrologic flows. Tool 5 provides instantaneous flow hydrologic envelope profiles, using measured instantaneous flow data integrated with measured daily-flow value data. Tool 6 provides precise estimates of hydrologic flow over entire data time-series suitable for driving scour simulation models. Tool 7 provides a threshold of flow and probability of time-under-load interactive calculator that allows selection of a desired bridge design lifespan, ranging from 1 to 250 years, and identification of a flow interval of interest. Tool 8 provides a flow-random sampling interactive tool, developed to facilitate easy access to large datasets of randomly sampled flow data measurements from unique locations for purposes of computing and testing future models of bridge pier scour.

Virginia↗

Modeling the water-quality effects to the Klamath River from recirculation in drains and canals, Oregon and California, 2006–15

The potential recirculation of Klamath Strait Drain (hereafter called by its local name, “Klamath Straits Drain”) water into Ady Canal to reduce the drain discharge of high nutrient loads into the Klamath River was assessed by the U.S. Geological Survey for the Bureau of Reclamation. To study the feasibility of recirculation, this investigation evaluated three recirculation scenarios over a 10-year period from 2006 to 2015, as a series of 1-year model simulations. A combination of two existing hydrodynamic, water-temperature, and water-quality models (CE-QUAL-W2) were used, including (1) the Link-Keno reach of the Klamath River, using Klamath Straits Drain as a tributary and for calendar years 2006–11, and (2) the same Link-Keno model used for calendar years 2012–15 in combination with an independent Klamath Straits Drain model from 2012 to 2015. Model simulations using the water-quality models were configured for the base case conditions and three different sets of recirculation scenarios: the maximum year-round recirculation without limits (scenario 1), limited year-round recirculation fixed by the current pipe flow configuration (scenario 2), and limited seasonal recirculation (May–September) also fixed by the current pipe flow configuration (scenario 3). In the base case, estimates of annual average daily total nitrogen loads and daily total phosphorus loads exported to the Klamath River from the Klamath Straits Drain were as much as 3,060 and 457 pounds per day (lbs/day), respectively. Currently (2023), the Total Maximum Daily Loads allocations for the Klamath Straits Drain are 21 and 268 lbs/day for total phosphorus and total nitrogen, respectively, so these maximum estimates exceed the current Total Maximum Daily Loads by greater than an order of magnitude. With scenario 1, load reductions occurred year-round for all constituents evaluated (total nitrogen, total phosphorus, 5-day biochemical oxygen demand [BOD5], 5-day carbonaceous biochemical oxygen demand) for the Klamath Straits Drain discharging to the Klamath River. Scenario 2 also had large reductions in total nitrogen, total phosphorus, and BOD5 loads. Substantial reductions did occur for scenario 3 but were constrained to only the active recirculation period from May through September. Despite the restricted period, the average reductions in the annual average daily load for total phosphorus and total nitrogen were 32.1 percent and 26.5 percent, respectively. The Ady Canal diverts high nutrient loads from the Klamath River, so the loading tradeoffs to the Klamath River between no recirculation and the recirculation scenarios were calculated. On an annual basis, the overall net balance between the Klamath Straits Drain and Ady Canal resulted in more total nitrogen and total phosphorus load reductions to the Klamath River for the three recirculation scenarios than the base case, for most years. In contrast, the net balance for BOD5 loads was higher to the Klamath River for the three recirculation scenarios than the base case, for most years. With the recirculation scenarios, the optimal recirculation periods to benefit Ady Canal, Klamath River, and Klamath Straits Drain did not always coincide. Recirculation would be most effective at reducing loads toward the Klamath Straits Drain Total Maximum Daily Load allocations in the spring (March–May) of each year. However, recirculation during these months would also increase salinity in the Ady Canal. In summer, recirculation would reduce Klamath Straits Drain loads toward the Total Maximum Daily Load allocations, though recirculation could decrease Klamath River water quality mostly because of decreased withdrawals of Klamath River water by the Ady Canal. Scenario 3 avoided recirculation into Ady Canal in the early spring months when salinity concerns would be the highest, while still decreasing nutrient loads exported from the Klamath Straits Drain to the Klamath River in the summer months.

California, Oregon↗

Evaluation of the characteristics, discharge, and water quality of selected springs at Fort Irwin National Training Center, San Bernardino County, California

Eight springs and seeps at Fort Irwin National Training Center were described and categorized by their general characteristics, discharge, geophysical properties, and water quality between 2015 and 2017. The data collected establish a modern (2017) baseline of hydrologic conditions at the springs. Two types of springs were identified: (1) precipitation-fed upland springs (Cave, Desert King, Devouge, No Name, and Panther Springs) and (2) groundwater discharge-fed basin springs (Garlic, Bitter, and Jack Springs). Comparison of electrical resistivity tomography data collected at groundwater basin springs from 2015 to 2017 indicated that spring discharge and connection to the underlying groundwater system is highly focused, although the springs themselves appear diffuse and are spread out over a large area. Spring discharge was consistently less than reported by Thompson (1929), except at Garlic Spring where discharges and vegetation have increased in recent years. Multiple discrete flume and seepage meter measurements taken between October 2015 and April 2016 indicated that discharge changed predictably on diurnal and seasonal timescales in response to evapotranspiration. These preliminary results and the lush vegetation noted at some of the springs, particularly at Bitter, Garlic, and Jack Springs, indicated plant evapotranspiration accounts for a substantial part of the discharge from these springs. The quality of water ranges from fresh in precipitation-fed upland springs (Cave, Desert King, Devouge, and Panther Springs) to slightly saline (Garlic and Jack Springs) and moderately saline (Bitter Spring) in groundwater-fed discharge springs. Nitrate concentrations from water at most of the springs were less than 3 milligrams per liter, except for samples from Devouge and Desert King Springs and one sample from Jack Spring. An analysis of delta nitrogen-15 in nitrate (δ 15 N-NO 3 ) and delta oxygen-18 in nitrate (δ 18 O-NO 3 ) indicates high nitrate concentrations in excess of the U.S. Environmental Protection Agency maximum contaminant level at Jack Spring and Desert King Spring resulting from the dissolution of nitrate-bearing caliche deposits; nitrate concentrations at Devouge Spring are a result of algal growth within the spring, and the source of nitrate concentrations in Garlic Spring are consistent with a treated wastewater origin from Langford Valley-Irwin subbasin upgradient. The source of water in upland springs, indicated by values of delta oxygen-18 (δ 18 O) and delta deuterium (δD) are consistent with recharge from winter precipitation. In groundwater basin springs, values of δ 18 O and δD are consistent with groundwater sampled from nearby wells. Summer monsoonal precipitation appears to contribute little water to spring flow. Most springs contain low levels of tritium and appear to be primarily older (pre-1950s) groundwater. Groundwater basin springs with detectable tritium may result from occasional streamflow in nearby washes. These springs could be susceptible to decreases in flow during extended dry periods when the localized recharge may be reduced due to the loss of focused recharge through nearby washes. Groundwater samples from Garlic and Bitter Springs contained arsenic concentrations above the U.S. Environmental Protection Agency maximum contaminant level. Groundwater samples from all springs, except Cave, Desert King, and Devouge Springs, exceeded the State of California maximum contaminant level for fluoride. Garlic Spring was the only sampled spring that contained vanadium concentrations that exceeded the State of California notification level. Only a single water sample from Jack Spring contained uranium at a concentration that exceeded the U.S. Environmental Protection Agency maximum contaminant level. Many other constituents of concern were analyzed, including those from anthropogenic sources that may be a result of military activities. Most of these constituents were not detected above their respective reporting levels in spring water; only 15 were detected in spring waters. Diesel and gasoline degradants, many of which also occur naturally, were the most commonly detected compounds. Several other organic compounds, primarily solvents or their degradants, were detected in groundwater basin springs. These constituents, in order of decreasing detection frequency, were carbon disulfide; perchlorate; mercury; acetone; methylnaphthalene; toluene; methyl ethyl ketone; cyanide; and styrene; 4-iso-propyl-toluene; isopropylbenzene; methyl salicylate; and phenol. Except for Garlic Spring, which is affected by discharges of treated wastewater, the quality of water from most springs appears to be relatively unaffected by activities at the Fort Irwin National Training Center.

California↗

Hydrogeologic framework and extent of saltwater intrusion in Kings, Queens, and Nassau Counties, Long Island, New York

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gamma-ray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco-Magothy aquifer in southwestern Nassau County appears to be increasing.

New York↗

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon↗

Evaluating hydrologic data products for scientific and management applications related to potential future streamflow conditions in the Upper Mississippi and Illinois Rivers

The hydrology of the Upper Mississippi and Illinois Rivers is a fundamental driver of ecosystem patterns and processes across a large portion of the United States. Quantitative hydrologic data for the main stems of these rivers underlie numerous scientific investigations, statistical models, and decision-making processes for local, State, and Federal agencies involved in the Upper Mississippi River Restoration program. Although historical hydrologic data exist, data representing potential future conditions of the Upper Mississippi and Illinois Rivers lack the resolution necessary to anticipate biotic and abiotic responses to altered hydrology and to determine resilient management actions. A source of future hydrologic scenarios is the readily available LOCA–VIC–mizuRoute hydrologic data products (named for the chain of models the data are produced from—localized constructed analogs, Variable Infiltration Capacity macroscale hydrological model, and the mizuRoute hydrologic routing model—that we shorten further to LVM in this report) that include simulated discharges for historic and future timeframes. The objective of this study is to assess the reliability of the hydrologic data products for their use in Upper Mississippi River Restoration program applications. Key study questions are (1) do the hydrologic data products reproduce characteristics of hydrology necessary to support ecological modeling and restoration decision-making applications within the Upper Mississippi River Restoration program? and (2) are there geographic differences in the reliability of the hydrologic data products? Seven characteristics of river hydrology were selected related to flow magnitude, seasonality, and regime for evaluation. The seven characteristics were calculated using observed and historical simulated hydrologic data at 19 U.S. Geological Survey streamgages throughout the basins of the Upper Mississippi and Illinois Rivers; two streamgages are located on the main stem of the Mississippi River and two streamgages are located on the main stem of the Illinois River. Statistical comparisons between observed and historical simulated characteristics indicated that the hydrologic data products did not reliably represent historical hydrologic conditions in the basin or main stem. The hydrologic data products we evaluated could not reliably capture the overall hydrologic regime or flow magnitudes; the latter is evidenced by substantial underestimates of discharge at most streamgages. Seasonal hydrologic characteristics were captured more reliably than flow magnitude, but overall correspondence was low for most streamgages. A weak latitudinal pattern in seasonal characteristics indicated the hydrologic data products poorly represent streamflow timing in snow-affected regions of the basin. Discrepancies in magnitude, seasonality, and regime indicate the potential for multiple sources of error. Because poor correspondence was present across all 19 streamgages, it was not possible to identify specific drivers of poor performance (that is, drainage area or geography). The modeling chain should be evaluated for biases associated with meteorologic forcing data, as well as hydrologic model formulation and calibration. We conclude that the hydrologic data products we evaluated appear unsuitable for applications tied to habitat and ecosystem restoration and management in the Upper Mississippi and Illinois Rivers. Plans to develop a future hydrology dataset for the Upper Mississippi River Restoration program would benefit from ongoing work to improve global climate model output downscaling methods, to improve hydrologic models, to make use of innovations in machine-learning approaches for projecting hydrology, and other efforts. The framework developed herein to evaluate hydrometeorological outputs generated using global climate models for a specific water resources application is a transferrable approach that could be applied to other data products and river systems.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗

Geology and uranium deposits of the Caribou area, Boulder County, Colorado

Pitchblende was discovered in the Caribou mine, Boulder County. Colo., in 1948 by Consolidated Caribou Silver Mines, incorporated during the reopening of the 1040-level. Since 1948, a small quantity of pitchblende has been produced from the Radium vein during exploratory work that was done under contract with the Atomic Energy Commission. The Caribou mine, primarily a silver mine, develops five veins (Caribou. No Name, Poorman, Sherman and Silver Dollar) that were formerly worked as separate mines. Three other veins, the Radium, Elmer, and Nelson, are exposed in the Caribou mine only at and below the 920-level. The general geology of the Caribou area, known since 1869 as a silver district, has been discussed by Bastin and Hill (1917) and more recently by Lovering and Goddard (1950). In 1937 a detailed study of the petrography of the Caribou stock was made by Smith (1938). In 1948 and 1949, R.U. King (1950) of the U. S., Geological Survey examined the Caribou mine to evaluate the uranium deposit. In 1950, D.M. Sheridan and E.N. Hinrichs of the U. S. Geological Survey made a radioactivity reconnaissance of the mine dumps at Caribou and vicinity. This reconnaissance included examination of most of the accessible workings of the Caribou mine. Other reports on the area are listed in the bibliography. As a result of King’s preliminary examination and the later exploration for uranium on the property, the U.S. Geological Survey in 1951 began a detailed investigation of the Caribou area on behalf of the Division of Raw Materials of the Atomic Energy Commission, the principal purpose of which was to evaluate the uranium occurrences in the area. The study was begun by E. P. Kaiser and W. S., Cavender, who established a triangulation net, using a base line in Caribou Park for horizontal control and a U.S. Geological Survey bench mark on the top of Caribou Hill for vertical control. Upon the assignment of Kaiser to another project in August 1951 the investigation was continued until November 1951, by F.B. Moore and W. S. Cavender. During the study of the district the writers prepared a geologic map at a scale of 1 inch to 100 feet of an area of about one square mile. The mapped area includes most of the mines in the Caribou district and embraces the eastern half of Caribou Hill and adjacent areas to the east and north. In addition, all accessible underground workings were mapped at a scale of 1 inch to 40 feet. The writers wish to acknowledge the cooperation of the staff of the Consolidated Caribou Silver Mines incorporated, who made the Caribou mine available for examination at all times and who furnished maps and suggestions that were of great assistance. Mr. A.E. Blakesley, owner of the Comstock mine, was also most cooperative in making possible the examination of his mine. Thanks are due Dr., E.E. Wahlstrom of the University of Colorado and to the Boulder Daily Camera for the use of their files containing information on the Caribou mine. To S. Lovering of the U.S. Geological Survey made many valuable suggestions on the identification of alteration products in thin section.

Colorado↗

Montroseite, a new vanadium oxide from the Colorado plateaus

Montroseite, a new vanadium mineral named from Montrose County, Colorado, has been found in four mines in western Colorado and in two mines in eastern Utah. It is black, opaque, submetallic, and occurs in microscopic bladed crystals of the orthorhombic dipyramidal class. The axial ratio is a:b:c = 0.509:1:0.310, the common forms are b {010}, m {110}, p {121}, and a large vicinal form is approximately {0, 10, 1}. The observed specific gravity is 4.0 and the calculated specific gravity is 4.15. The composition is essentially VO(OH), with some iron commonly substituted for vanadium. Partial oxidation to VO 2 has taken place. Chemical analyses and X-ray diffraction data are given. Single crystal study indicated that the space group symmetry is Pbnm(D 2h 16 ).

Colorado;Utah↗

Selected annotated bibliography of the geology of uraniferous and radioactive native bituminous substances, exclusive of coals, in the United States

Native bituminous substances are divided into two groups, 1) bitumens and, 2) pyrobitumens. Bitumens are composed principally of hydrocarbons substantially free from oxygenated bodies, are fusible, and are soluble in carbon disulfide. Native bitumens occur in liquid and solid forms. The native liquid bitumens include all petroleums or crude oils. Native solid bitumens include native waxes such as ozocerite, asphalts or petroleum tars, and asphaltites such as gilsonite and grahamite. Pyrobitumens are composed principally of hydrocarbons which may contain oxygenated bodies. They are infusible and are insoluble, or nearly insoluble, in carbon disulfide. Native pyrobitumens are divided into an oxygen-containing group including peats, lignites, and coals, and an essentially oxygen-free, asphaltic group including such substances as wurtzilite, albertite, impsonite, and ingramite. Thucholites, which are carbonaceous substances that may contain uranium, thorium, and rare earths, commonly are considered to be pyrobitumens. Their compositions are variable and may fall into either the oxygen-containing or oxygen-free group. All varieties of native bituminous substances may be associated with mineral matter. The nomenclature of bitumens and pyrobitumens is used very loosely in the literature. This circumstance arises from the difficulty in recognizing many of these substances by visual examination, and because many of them can be identified accurately only by chemical methods. Inasmuch as some of the chemical procedures are time-consuming and satisfactory analytical methods have not been devised for all these substances, geologists generally have not obtained precise identifications but rather have used names that appeared most appropriate to the circumstances. It is expected that future research will show many substances called "asphaltite," "thucholite," etc., to be incorrectly identified. The nomenclature used by the authors of the various references of this bibliography is followed without deviation or further discussion. The stratigraphic nomenclature also is that used by the authors. In this bibliography emphasis is placed on reports dealing with the uranium contents and radioactivity of native bituminous substances rather than on mineralogical and chemical studies of these substances. The distribution of the substances described in the references is shown on the accompanying map. The indicated presence of these substances does not infer that they contain sufficient radioactive elements to constitute ores.

Trace Elements Investigations↗

US Topo Product Standard

This document defines a U.S. Geological Survey (USGS) digital topographic map. This map product series, named “US Topo,” is modeled on the now historical USGS 7.5-minute (1:24,000 scale) topographic map series produced and printed by the USGS from 1947 to 2006. US Topo maps have the same extent, scale, and general layout as the historical topographic maps. US Topo maps incorporate an orthorectified image (hereinafter referred to as “orthoimage”) and shaded relief image along with a selection of data that were included in the historical 7.5-minute topographic maps. Between June and September of 2017, the USGS transitioned the format of US Topo maps to be published, by using a geospatial extension, in an International Organization for Standardization (ISO) 32000-compliant Adobe® portable document format (PDF) that is called a “geospatial PDF.” Previously, US Topo maps were published, by using geospatial extensions patented by TerraGo® Technologies, in PDF in a format called a “GeoPDF®.” The geospatial PDF design allows a user to zoom in and out in a georeferenced environment, turn layers on and off, view or print any combination of layers, and print any portion of the map at the published scale. US Topo maps are intended to serve conventional map users by providing geographic information system (GIS) information in symbolized form in the customary topographic map layout. The maps are not intended for advanced GIS analysis applications. These products are built on standard coordinate systems and include full U.S. National Grid (USNG) lines, making US Topo maps particularly useful for emergency first-response operations. These maps are also used by traditional topographic map users, such as resource managers, planners, and recreational users who continue to have a need for the symbolized feature data contained in the 7.5-minute quadrangle maps. Full-size style sheet templates in PDF defining the placement of map elements, marginalia, and font sizes and styles accompany this standard. US Topo maps published as geospatial PDFs are fashioned to conform to these style sheets so that a user can print out a map at the 1:24,000, 1:25,000, or 1:20,000 scale using the dimensions of the traditional standard 7.5-minute quadrangle. Symbology and type specifications for feature content and detailed requirements for geospatial content will be published separately. This document is an update of the US Topo Product Standard published in 2011 (Cooley and others, 2011). It is applicable to all US Topo maps. Updates in this version include the introduction of an ISO 32000-compliant geospatial PDF as a new file format for published maps; new style sheet templates for 1:24,000-scale maps (conterminous United States and Hawaii), 1:25,000-scale maps (Alaska), and 1:20,000-scale maps (Puerto Rico and U.S. Virgin Islands); an updated US Topo Map Symbol attachment; minor updates to text, including changes to the features and layers included in the US Topo product and the sheet size of the US Topo maps; updated figures demonstrating the US Topo product; an updated metadata file containing map-specific information.

Techniques and Methods↗

Kendall-Theil Robust Line (KTRLine--version 1.0)-A Visual Basic Program for Calculating and Graphing Robust Nonparametric Estimates of Linear-Regression Coefficients Between Two Continuous Variables

The Kendall-Theil Robust Line software (KTRLine-version 1.0) is a Visual Basic program that may be used with the Microsoft Windows operating system to calculate parameters for robust, nonparametric estimates of linear-regression coefficients between two continuous variables. The KTRLine software was developed by the U.S. Geological Survey, in cooperation with the Federal Highway Administration, for use in stochastic data modeling with local, regional, and national hydrologic data sets to develop planning-level estimates of potential effects of highway runoff on the quality of receiving waters. The Kendall-Theil robust line was selected because this robust nonparametric method is resistant to the effects of outliers and nonnormality in residuals that commonly characterize hydrologic data sets. The slope of the line is calculated as the median of all possible pairwise slopes between points. The intercept is calculated so that the line will run through the median of input data. A single-line model or a multisegment model may be specified. The program was developed to provide regression equations with an error component for stochastic data generation because nonparametric multisegment regression tools are not available with the software that is commonly used to develop regression models. The Kendall-Theil robust line is a median line and, therefore, may underestimate total mass, volume, or loads unless the error component or a bias correction factor is incorporated into the estimate. Regression statistics such as the median error, the median absolute deviation, the prediction error sum of squares, the root mean square error, the confidence interval for the slope, and the bias correction factor for median estimates are calculated by use of nonparametric methods. These statistics, however, may be used to formulate estimates of mass, volume, or total loads. The program is used to read a two- or three-column tab-delimited input file with variable names in the first row and data in subsequent rows. The user may choose the columns that contain the independent (X) and dependent (Y) variable. A third column, if present, may contain metadata such as the sample-collection location and date. The program screens the input files and plots the data. The KTRLine software is a graphical tool that facilitates development of regression models by use of graphs of the regression line with data, the regression residuals (with X or Y), and percentile plots of the cumulative frequency of the X variable, Y variable, and the regression residuals. The user may individually transform the independent and dependent variables to reduce heteroscedasticity and to linearize data. The program plots the data and the regression line. The program also prints model specifications and regression statistics to the screen. The user may save and print the regression results. The program can accept data sets that contain up to about 15,000 XY data points, but because the program must sort the array of all pairwise slopes, the program may be perceptibly slow with data sets that contain more than about 1,000 points.

Techniques and Methods↗

TracerLPM (Version 1): An Excel® workbook for interpreting groundwater age distributions from environmental tracer data

TracerLPM is an interactive Excel® (2007 or later) workbook program for evaluating groundwater age distributions from environmental tracer data by using lumped parameter models (LPMs). Lumped parameter models are mathematical models of transport based on simplified aquifer geometry and flow configurations that account for effects of hydrodynamic dispersion or mixing within the aquifer, well bore, or discharge area. Five primary LPMs are included in the workbook: piston-flow model (PFM), exponential mixing model (EMM), exponential piston-flow model (EPM), partial exponential model (PEM), and dispersion model (DM). Binary mixing models (BMM) can be created by combining primary LPMs in various combinations. Travel time through the unsaturated zone can be included as an additional parameter. TracerLPM also allows users to enter age distributions determined from other methods, such as particle tracking results from numerical groundwater-flow models or from other LPMs not included in this program. Tracers of both young groundwater (anthropogenic atmospheric gases and isotopic substances indicating post-1940s recharge) and much older groundwater (carbon-14 and helium-4) can be interpreted simultaneously so that estimates of the groundwater age distribution for samples with a wide range of ages can be constrained. TracerLPM is organized to permit a comprehensive interpretive approach consisting of hydrogeologic conceptualization, visual examination of data and models, and best-fit parameter estimation. Groundwater age distributions can be evaluated by comparing measured and modeled tracer concentrations in two ways: (1) multiple tracers analyzed simultaneously can be evaluated against each other for concordance with modeled concentrations (tracer-tracer application) or (2) tracer time-series data can be evaluated for concordance with modeled trends (tracer-time application). Groundwater-age estimates can also be obtained for samples with a single tracer measurement at one point in time; however, prior knowledge of an appropriate LPM is required because the mean age is often non-unique. LPM output concentrations depend on model parameters and sample date. All of the LPMs have a parameter for mean age. The EPM, PEM, and DM have an additional parameter that characterizes the degree of age mixing in the sample. BMMs have a parameter for the fraction of the first component in the mixture. An LPM, together with its parameter values, provides a description of the age distribution or the fractional contribution of water for every age of recharge contained within a sample. For the PFM, the age distribution is a unit pulse at one distinct age. For the other LPMs, the age distribution can be much broader and span decades, centuries, millennia, or more. For a sample with a mixture of groundwater ages, the reported interpretation of tracer data includes the LPM name, the mean age, and the values of any other independent model parameters. TracerLPM also can be used for simulating the responses of wells, springs, streams, or other groundwater discharge receptors to nonpoint-source contaminants that are introduced in recharge, such as nitrate. This is done by combining an LPM or user-defined age distribution with information on contaminant loading at the water table. Information on historic contaminant loading can be used to help evaluate a model's ability to match real world conditions and understand observed contaminant trends, while information on future contaminant loading scenarios can be used to forecast potential contaminant trends.

Techniques and Methods↗

Advective transport observations with MODPATH-OBS--documentation of the MODPATH observation process

The MODPATH-OBS computer program described in this report is designed to calculate simulated equivalents for observations related to advective groundwater transport that can be represented in a quantitative way by using simulated particle-tracking data. The simulated equivalents supported by MODPATH-OBS are (1) distance from a source location at a defined time, or proximity to an observed location; (2) time of travel from an initial location to defined locations, areas, or volumes of the simulated system; (3) concentrations used to simulate groundwater age; and (4) percentages of water derived from contributing source areas. Although particle tracking only simulates the advective component of conservative transport, effects of non-conservative processes such as retardation can be approximated through manipulation of the effective-porosity value used to calculate velocity based on the properties of selected conservative tracers. This program can also account for simple decay or production, but it cannot account for diffusion. Dispersion can be represented through direct simulation of subsurface heterogeneity and the use of many particles. MODPATH-OBS acts as a postprocessor to MODPATH, so that the sequence of model runs generally required is MODFLOW, MODPATH, and MODPATH-OBS. The version of MODFLOW and MODPATH that support the version of MODPATH-OBS presented in this report are MODFLOW-2005 or MODFLOW-LGR, and MODPATH-LGR. MODFLOW-LGR is derived from MODFLOW-2005, MODPATH 5, and MODPATH 6 and supports local grid refinement. MODPATH-LGR is derived from MODPATH 5. It supports the forward and backward tracking of particles through locally refined grids and provides the output needed for MODPATH_OBS. For a single grid and no observations, MODPATH-LGR results are equivalent to MODPATH 5. MODPATH-LGR and MODPATH-OBS simulations can use nearly all of the capabilities of MODFLOW-2005 and MODFLOW-LGR; for example, simulations may be steady-state, transient, or a combination. Though the program name MODPATH-OBS specifically refers to observations, the program also can be used to calculate model prediction of observations. MODPATH-OBS is primarily intended for use with separate programs that conduct sensitivity analysis, data needs assessment, parameter estimation, and uncertainty analysis, such as UCODE_2005, and PEST. In many circumstances, refined grids in selected parts of a model are important to simulated hydraulics, detailed inflows and outflows, or other system characteristics. MODFLOW-LGR and MODPATH-LGR support accurate local grid refinement in which both mass (flows) and energy (head) are conserved across the local grid boundary. MODPATH-OBS is designed to take advantage of these capabilities. For example, particles tracked between a pumping well and a nearby stream, which are simulated poorly if a river and well are located in a single large grid cell, can be simulated with improved accuracy using a locally refined grid in MODFLOW-LGR, MODPATH-LGR, and MODPATH-OBS. The locally-refined-grid approach can provide more accurate simulated equivalents to observed transport between the well and the river. The documentation presented here includes a brief discussion of previous work, description of the methods, and detailed descriptions of the required input files and how the output files are typically used.

Techniques and Methods↗

Description of input and examples for PHREEQC version 3: A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations

PHREEQC version 3 is a computer program written in the C and C++ programming languages that is designed to perform a wide variety of aqueous geochemical calculations. PHREEQC implements several types of aqueous models: two ion-association aqueous models (the Lawrence Livermore National Laboratory model and WATEQ4F), a Pitzer specific-ion-interaction aqueous model, and the SIT (Specific ion Interaction Theory) aqueous model. Using any of these aqueous models, PHREEQC has capabilities for (1) speciation and saturation-index calculations; (2) batch-reaction and one-dimensional (1D) transport calculations with reversible and irreversible reactions, which include aqueous, mineral, gas, solid-solution, surface-complexation, and ion-exchange equilibria, and specified mole transfers of reactants, kinetically controlled reactions, mixing of solutions, and pressure and temperature changes; and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for differences in composition between waters within specified compositional uncertainty limits. Many new modeling features were added to PHREEQC version 3 relative to version 2. The Pitzer aqueous model ( pitzer.dat database, with keyword PITZER ) can be used for high-salinity waters that are beyond the range of application for the Debye-Hückel theory. The Peng-Robinson equation of state has been implemented for calculating the solubility of gases at high pressure. Specific volumes of aqueous species are calculated as a function of the dielectric properties of water and the ionic strength of the solution, which allows calculation of pressure effects on chemical reactions and the density of a solution. The specific conductance and the density of a solution are calculated and printed in the output file. In addition to Runge-Kutta integration, a stiff ordinary differential equation solver (CVODE) has been included for kinetic calculations with multiple rates that occur at widely different time scales. Surface complexation can be calculated with the CD-MUSIC (Charge Distribution MUltiSIte Complexation) triple-layer model in addition to the diffuse-layer model. The composition of the electrical double layer of a surface can be estimated by using the Donnan approach, which is more robust and faster than the alternative Borkovec-Westall integration. Multicomponent diffusion, diffusion in the electrostatic double layer on a surface, and transport of colloids with simultaneous surface complexation have been added to the transport module. A series of keyword data blocks has been added for isotope calculations— ISOTOPES, CALCULATE_VALUES, ISOTOPE_ALPHAS, ISOTOPE_RATIOS, and NAMED_EXPRESSIONS . Solution isotopic data can be input in conventional units (for example, permil, percent modern carbon, or tritium units) and the numbers are converted to moles of isotope by PHREEQC. The isotopes are treated as individual components (they must be defined as individual master species) so that each isotope has its own set of aqueous species, gases, and solids. The isotope-related keywords allow calculating equilibrium fractionation of isotopes among the species and phases of a system. The calculated isotopic compositions are printed in easily readable conventional units. New keywords and options facilitate the setup of input files and the interpretation of the results. Keyword data blocks can be copied (keyword COPY ) and deleted (keyword DELETE ). Keyword data items can be altered by using the keyword data blocks with the _MODIFY extension and a simulation can be run with all reactants of a given index number (keyword RUN_CELLS ). The definition of the complete chemical state of all reactants of PHREEQC can be saved in a file in a raw data format ( DUMP and _RAW keywords). The file can be read as part of another input file with the INCLUDE$ keyword. These keywords facilitate the use of IPhreeqc, which is a module implementing all PHREEQC version 3 capabilities; the module is designed to be used in other programs that need to implement geochemical calculations; for example, transport codes. Charting capabilities have been added to some versions of PHREEQC. Charting capabilities have been added to Windows distributions of PHREEQC version 3. (Charting on Linux requires installation of Wine.) The keyword data block USER_GRAPH allows selection of data for plotting and manipulation of chart appearance. Almost any results from geochemical simulations (for example, concentrations, activities, or saturation indices) can be retrieved by using Basic language functions and specified as data for plotting in USER_GRAPH . Results of transport simulations can be plotted against distance or time. Data can be added to a chart from tab-separated-values files. All input for PHREEQC version 3 is defined in keyword data blocks, each of which may have a series of identifiers for specific types of data. This report provides a complete description of each keyword data block and its associated identifiers. Input files for 22 examples that demonstrate most of the capabilities of PHREEQC version 3 are described and the results of the example simulations are presented and discussed.

Techniques and Methods↗