Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geologic Notes”

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 577 records · Page 32Linked to original sources
USGS2007

Revisiting the cape cod bacteria injection experiment using a stochastic modeling approach

Bromide and resting-cell bacteria tracer tests conducted in a sandy aquifer at the U.S. Geological Survey Cape Cod site in 1987 were reinterpreted using a three-dimensional stochastic approach. Bacteria transport was coupled to colloid filtration theory through functional dependence of local-scale colloid transport parameters upon hydraulic conductivity and seepage velocity in a stochastic advection - dispersion/attachment - detachment model. Geostatistical information on the hydraulic conductivity (K) field that was unavailable at the time of the original test was utilized as input. Using geostatistical parameters, a groundwater flow and particle-tracking model of conservative solute transport was calibrated to the bromide-tracer breakthrough data. An optimization routine was employed over 100 realizations to adjust the mean and variance ofthe natural-logarithm of hydraulic conductivity (InK) field to achieve best fit of a simulated, average bromide breakthrough curve. A stochastic particle-tracking model for the bacteria was run without adjustments to the local-scale colloid transport parameters. Good predictions of mean bacteria breakthrough were achieved using several approaches for modeling components of the system. Simulations incorporating the recent Tufenkji and Elimelech (Environ. Sci. Technol. 2004, 38, 529-536) correlation equation for estimating single collector efficiency were compared to those using the older Rajagopalan and Tien (AIChE J. 1976, 22, 523-533) model. Both appeared to work equally well at predicting mean bacteria breakthrough using a constant mean bacteria diameter for this set of field conditions. Simulations using a distribution of bacterial cell diameters available from original field notes yielded a slight improvement in the model and data agreement compared to simulations using an average bacterial diameter. The stochastic approach based on estimates of local-scale parameters for the bacteria-transport process reasonably captured the mean bacteria transport behavior and calculated an envelope of uncertainty that bracketed the observations in most simulation cases.

Massachusetts↗
USGS1990

Studies of angiospermous woods in Australian brown coal by nuclear magnetic resonance and analytical pyrolysis: new insight into early coalification

Many Tertiary coals contain abundant fossilized remains of angiosperms, which commonly dominated the ancient peat-swamp environments; modern analogs of such swamps can be found in tropical and subtropical regions of the world. Comparisons of angiospermous wood from Australian brown coal with similar wood buried in modern peat swamps of Indonesia have provided some new insights into coalification reactions. These comparisons were made by using solid-state 13 C nuclear magnetic resonance (NMR) techniques and pyrolsis-gas chromatography-mass spectrometry (py-gc-ms). These two modern techniques are especially suited for detailed structural evaluation of the complex macromolecules in coal. The earliest transformation (peatification) of organic matter in angiospermous wood is the degradation and removal of cellulosic components and the concomitant selective preservation of lignin-derived components. The angiospermous lignin that becomes enriched in wood as a result of cellulose degradation also is modified by coalification reactions; this modification, however, does not involve degradation and removal of the lignin. Rather, the early coalification process transforms the lignin phenols (guiacyl and syringyl) to eventually yield the aromatic structures typically found in brown coal. One such transformation, which is determined from NMR data, involves the cleavage of aryl-ether bonds that link guaiacyl and syringyl units in lignin, and this transformation leads to the formation of free lignin phenols. Another transformation, which is also determined from the NMR data, involves the loss of methoxyl groups, probably via demethylation, to produce catechol-like structures. Coincident with ether-cleavage and demethylation, the aromatic rings derived from lignin phenols become more carbon-substituted and cross linked, as determined by dipolar-dephasing NMR studies. This cross linking is probably responsible for preventing the lignin phenols, which are freed from the lignin macromolecule by ether cleavage, from being removed from the coal by dissolution. Pyrolysis data suggest that the syringyl units are altered more readily than are guaiacyl units, and this difference in resistance leads to an enrichment of the guaiacyl units in fossil angiospermous woods. Many of the coalification reactions noted above occur to some degree in all angiospermous fossil wood examined; however, some significant differences are observed in the degree of coalification of the fossil wood samples from the same burial depth in the brown coal. These differences indicate that the depth and duration of burial are probably not entirely responsible for the variations in degree of coalification. Different rates of degradation in peat may have contributed to the variations in the apparent degree of coalification; some woods may have been altered more rapidly at the peat stage than others. Although preliminary, this systematic study of botanically related wood in peat and coal results in a more detailed differentiation of coalification reactions than have previous investigations. The combined use of solid-state 13 C NMR and py-gc-ms has facilitated this detailed insight into the coalification of angiospermous wood.

International Journal of Coal Geology↗
USGS1955

Radioactivity and uranium content of the Sharon Springs member of the Pierre shale and associated rocks in western Kansas and eastern Colorado

As a part of the Geological Survey's program of investigating uranium-bearing carbonaceous rocks on behalf of the Division of Raw Materials of the U.S. Atomic Energy Commission, a reconnaissance of the Sharon Springs member of the Pierre shale in western Kansas and eastern Colorado was conducted during 1954. The Sharon Springs member of the Pierre shale and its lateral equivalents ranges from 155 to about 500 feet in thickness and generally contains about 0.001 percent uranium, but some beds contain larger amounts. A 6-foot thick shale bed in Cheyenne County, Colo., contains about 0.006 percent uranium, a 4 1/2-foot thick sequence of beds in Crowley County, Colo., is estimated to contain between 0.004 and 0.005 percent uranium, and a 3 1/2-foot thick sequence of beds in Kiowa County, Colo., contains about 0.004 percent uranium. At several outcrop localities, sequences of beds as much as 9 1/2 feet thick contain about 0.003 percent uranium. Data from wells indicate that the 4 1/2-foot thick sequence of beds in Crowley County, Colo., may have a lateral extent of at least 5 1/2 miles. A gamma-ray log of a well in Yuma County, Colo., indicates the presence of a sequence of beds 66 feet thick which contains 0.005 to 0.010 percent equivalent uranium. No definite pattern of areal distribution of radioactivity and uranium content in the Sharon Springs is indicated by available data. Lateral variation in uranium content of individual beds was not noted in outcrops, which seldom extend more than 150 feet, but subsurface data from gamma-ray logs of wells indicate that both the maximum radioactivity and the thickness of radioactive beds are variable within distances of a few miles. Vertical variation in radioactivity and uranium content of the more radioactive beds is usually abrupt, but in the rocks as a whole the range of uranium content is so small that large variations in content are absent. In most of the gamma-ray logs examined there is only part of the sequence of rocks comprising the Pierre shale and Niobrara formation that exhibits radioactivity in excess of the average radioactivity of the two formations. Comparison of features of gamma-ray logs of wells in north-eastern Colorado suggests that the most radioactive part referred to above is a laterally correlatable sequence of beds. The stratigraphic position of the radioactive unit relative to the Pierre shale-Niobrara formation contact in oil industry scout reports, as identified from electric logs and wells, is variable within short distances. This may indicate that some of the Pierre-Niobrara contacts picked from electric logs may not correspond to the boundary that would be selected by examination of the rocks themselves, or may indicate that there is a facies relationship between teh lowermost part of the Pierre shale and the uppermost part of the Niobrara formation.

Colorado;Kansas↗
USGS1986

Preliminary map showing the thickness of glacial deposits in Ohio

This map was compiled as part of a U.S. Geological Survey project to portray the character and thickness of glacial deposits east of the Rocky Mountains. It is a preliminary map, compiled at a scale of 1:250,000 from sources of differing type and quality (fig. I, table 1). County drift-thickness maps produced by the Ohio Geological Survey are the most reliable source of data used. These county maps, either published or in manuscript form, were available for 36 of the 72 counties that contain glacial or glacially related deposits more than 50 ft thick. An earlier, reconnaissance drift-thickness map of northwestern Ohio (Ohio Department of Natural Resources, 1960) provided data for parts of 12 additional counties. Published maps of bedrock topography were utilized for 9 counties and for parts of 6 counties. The computation of drift thickness in these counties is less reliable then the drift-thickness maps published by the Ohio Geological Survey because the bedrock topography maps, most of which are at least 20 years old, were compiled from a more limited set of data. In the remaining counties, only well-log data were available. Drift-thickness data in these counties are the least reliable; the determination of thickness trends was commonly guided by the configuration of preglacial drainage channels (as drawn by Stout and others, 1943) and by topography. Throughout the glaciated areas of Ohio, and especially near the glacial margin, a network of drainage channels is buried by glacial drift. Some of these channels arc of a preglacial age while others were formed during interglacial periods; some were cut into the bedrock by rivers— either preglacial rivers or meltwater streams flowing along the glacial margin—while others, such as the lower Cuyahoga River valley, may have been overdeepened by the erosive power of glacial ice. The thickest deposits of drift in the state occur in the lower Cuyahoga River valley near Lake Erie, and above the buried valley in west-central Ohio that has been historically considered a part of the Teays Valley. The thickest drift reliably measured in Ohio to date is 602 ft in a drill hole near Cleveland. In contrast to the extreme variations in drift thickness encountered in the vicinity of buried channels, drift on the upland arcus is generally thinner and the variations in thickness are much less pronounced. Worthy of note, however, are three large areas where the drift sheet is relatively thick. In northwestern Ohio, a large volume of drift was deposited along the flanks of the Erie ice lobe (fig. 2) near the interlobate position with the Saginaw lobe to the northwest; drift thickness there exceeds 200 ft. Thick drift was also deposited in a roughly cast-west band across the Miami lobe. The mechanism that produced this band of thick drift is not obvious, but it may have been influenced in part by bedrock topography. Bedrock control of drift thickness is more clearly indicated to the cast of Columbus, along the eastern flank of the Scioto lobe, where ice slow was resisted by rocks of the Allegheny plateau. The edge of the plateau, or the Allegheny escarpment, is obscured by glacial deposits but its likely position (Fenneman, 1938; Stout and others, 1913; Dove, 1960; and Root and others, 1961) is shown on the map. Southward from the ice margin's reentrant position in southern Richland County, ice flowing eastward from the Scioto lobe encountered the topographically higher plateau, which constrained the ice and caused drift to accumulate in significant thicknesses just to the west of the escarpment.

Ohio↗
USGS2012

Use of ASTER and MODIS thermal infrared data to quantify heat flow and hydrothermal change at Yellowstone National Park

The overarching aim of this study was to use satellite thermal infrared (TIR) remote sensing to monitor geothermal activity within the Yellowstone geothermal area to meet the missions of both the U.S. Geological Survey and the Yellowstone National Park Geology Program. Specific goals were to: 1) address the challenges of monitoring the surface thermal characteristics of the > 10,000 spatially and temporally dynamic thermal features in the Park (including hot springs, pools, geysers, fumaroles, and mud pots) that are spread out over ~ 5000 km 2 , by using satellite TIR remote sensing tools (e.g., ASTER and MODIS), 2) to estimate the radiant geothermal heat flux (GHF) for Yellowstone's thermal areas, and 3) to identify normal, background thermal changes so that significant, abnormal changes can be recognized, should they ever occur (e.g., changes related to tectonic, hydrothermal, impending volcanic processes, or human activities, such as nearby geothermal development). ASTER TIR data (90-m pixels) were used to estimate the radiant GHF from all of Yellowstone's thermal features and update maps of thermal areas. MODIS TIR data (1-km pixels) were used to record background thermal radiance variations from March 2000 through December 2010 and establish thermal change detection limits. A lower limit for the radiant GHF estimated from ASTER TIR temperature data was established at ~ 2.0 GW, which is ~ 30–45% of the heat flux estimated through geochemical thermometry. Also, about 5 km 2 of thermal areas was added to the geodatabase of mapped thermal areas. A decade-long time-series of MODIS TIR radiance data was dominated by seasonal cycles. A background subtraction technique was used in an attempt to isolate variations due to geothermal changes. Several statistically significant perturbations were noted in the time-series from Norris Geyser Basin, however many of these did not correspond to documented thermal disturbances. This study provides concrete examples of the strengths and limitations of current satellite TIR monitoring of geothermal areas, highlighting some specific areas that can be improved. This work provides a framework for future satellite-based thermal monitoring at Yellowstone and other volcanic and geothermal systems.

Idaho, Montana, Oregon, Wyoming↗
USGS2005

Simulation of flow and sediment mobility using a multidimensional flow model for the White Sturgeon critical-habitat reach, Kootenai River near Bonners Ferry, Idaho

In 1994, the Kootenai River white sturgeon (Acipenser transmontanus) was listed as an Endangered Species as a direct result of two related observations. First, biologists observed that the white sturgeon population in the Kootenai River was declining. Second, they observed a decline in recruitment of juvenile sturgeon beginning in the 1950s with an almost total absence of recruitment since 1974, following the closure of Libby Dam in 1972. This second observation was attributed to changes in spawning and (or) rearing habitat resulting from alterations in the physical habitat, including flow regime, sediment-transport regime, and bed morphology of the river. The Kootenai River White Sturgeon Recovery Team was established to find and implement ways to improve spawning and rearing habitat used by white sturgeon. They identified the need to develop and apply a multidimensional flow model to certain reaches of the river to quantify physical habitat in a spatially distributed manner. The U.S. Geological Survey has addressed these needs by developing, calibrating, and validating a multidimensional flow model used to simulate streamflow and sediment mobility in the white sturgeon critical-habitat reach of the Kootenai River. This report describes the model and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to biological or other habitat data. This study was conducted in cooperation with the Kootenai Tribe of Idaho along a 23-kilometer reach of the Kootenai River, including the white sturgeon spawning reach near Bonners Ferry, Idaho that is about 108 to 131 kilometers below Libby Dam. U.S. Geological Survey's MultiDimensional Surface-Water Modeling System was used to construct a flow model for the critical-habitat reach of the Kootenai River white sturgeon, between river kilometers 228.4 and 245.9. Given streamflow, bed roughness, and downstream water-surface elevation, the model computes the velocity field, water-surface elevations, and boundary shear stress throughout the modeled reach. The 17.5 kilometer model reach was subdivided into two segments on the basis of predominant grain size: a straight reach with a sand, gravel, and cobble substrate located between the upstream model boundary at river kilometer 245.9 and the upstream end of Ambush Rock at river kilometer 244.6, and a meandering reach with a predominately sand substrate located between upstream end of Ambush Rock and the downstream model boundary at river kilometer 228.4. Model cell size in the x and y (horizontal) dimensions is 5 meters by 5 meters along the computational grid centerline with 15 nodes in the z (vertical) dimension. The model was calibrated to historical streamflows evenly distributed between 141.6 and 2,548.9 cubic meters per second. The model was validated by comparing simulated velocities with velocities measured at 15 cross sections during steady streamflow. These 15 cross sections were each measured multiple (7-13) times to obtain velocities suitable for comparison to the model results. Comparison of modeled and measured velocities suggests that the model does a good job of reproducing flow patterns in the river, although some discrepancies were noted. The model was used to simulate water-surface elevation, depth, velocity, bed shear stress, and sediment mobility for Kootenai River streamflows of 170, 566, 1,130, 1,700, and 2,270 cubic meters per second (6,000, 20,000, 40,000, 60,000, and 80,000 cubic feet per second). The three lowest streamflow simulations represent a range of typical river conditions before and since the construction of Libby Dam, and the highest streamflow simulation (2,270 cubic meters per second) is approximately equal to the annual median peak streamflow prior to emplacement of Libby Dam in 1972. Streamflow greater than 566 cubic meters per second were incrementally increased by 570 cubic meters per second. For each

Scientific Investigations Report↗
USGS1985

Isotopic studies of the late Archean plutonic rocks of the Wind River Range, Wyoming

Isotopic studies of the Rb-Sr and U-Th-Pb systems in whole-rock samples and the U-Pb systematics for zircons document the existence of two late Arehean intrusive events in the Wind River Range. All of the systems examined indicate an age of ∼2,630 ± 20 m.y. for the Louis Lake batholith. Apparent ages for the Bears Ears pluton range from 2,504 ± 40 m.y. to 2,575 ± 50 m.y. The scatter in apparent ages for the Bears Ears pluton does not appear to be primarily the result of disturbance by postintrusive events, but it may be explained by an isotopically inhomogenous magma at the time of intrusion. Data for a few samples indicate that the Wind River Range was affected locally by a postmagmatic hydrothermal event that was approximately Tertiary in age. This event lowered δ 18 O values and disturbed parent-daughter relationships in most of the isotopic systems investigated, but it was recent enough that there is no demonstrable effect in the Pb-Pb system. The Bears Ears pluton has some chemical and petrologic features that are similar to those reported for the granites in the Granite Mountains to the east. These granites are spatially associated with low-temperature uranium deposits of Tertiary age and have been shown to have lost large amounts of uranium during the early to middle Tertiary. U-Pb systematics indicate, however, that the low to moderate uranium contents and highly variable Th/U values noted for the Bears Ears pluton are best interpreted as being primary features. If uranium was lost after magma generation, the loss most likely occurred at the time of intrusion. Such a loss could account for uraniferous Precambrian pegmatites southwest of the main part of the Range. The two intrusive units apparently were derived from different protoliths that were formed during early to middle Archean. Initial isotopic ratios and petrochemistry for the Louis Lake batholith are consistent with an early Archean trondhjemitic to tonalitic source. The protolith for the Bears Ears pluton must have been more evolved and somewhat younger. Inconsistencies as to the degree of evolution of this protolith, as inferred from isotopic and trace-element data, suggest that the protolith may have been subjected to high-grade meta-morphism that caused loss of Rb and U prior to generation of the magma.

Wyoming↗
USGS2002

Topographic map of Mars M 25M RKN

NOTES ON BASE This map is based on data from the Mars Orbiter Laser Altimeter (MOLA) (Smith and others, 2001), an instrument on NASA's Mars Global Surveyor (MGS) spacecraft (Albee and others, 2001). The image used for the base of this map represents more than 600 million measurements gathered between 1999 and 2001, adjusted for consistency (Neumann and others, 2001, 2002) and converted to planetary radii. These have been converted to elevations above the areoid as determined from a martian gravity field solution GMM2 (Lemoine and others, 2001), truncated to degree and order 50, and oriented according to current standards (see below). The average accuracy of each point is originally ~100 meters in horizontal position, and ~1 meter in radius (Neumann, 2001). However, the total elevation uncertainty is at least +3 m due to the global error in the areoid (+1.8 meters according to Lemoine and others, 2001) and regional uncertainties in its shape (communication from Neumann, 2002). The measurements were converted into a digital elevation model (DEM) (communication from Neumann, 2002; Neumann and others, 2001; Smith and others, 2001) using Generic Mapping Tools software (Wessel and Smith, 1998), with a resolution of 0.015625 degrees per pixel or 64 pixels per degree. In projection, the pixels are 926.17 meters in size at the equator. PROJECTION The Mercator projection is used between latitudes +57°, with a central meridian at 0° and latitude equal to the nominal scale at 0°. The Polar Stereographic projection is used for the polar regions north of the +55° parallel and south of the –55° parallel with a central meridian set for both at 0°. The adopted equatorial radius is 3,396.19 km (Duxbury and others 2002; Seidelmann and others 2002). COORDINATE SYSTEM The MOLA data were initially referenced to an internally consistent inertial coordinate system, derived from tracking of the MGS spacecraft. By adopting appropriate values for the orientation of Mars as defined by the International Astronomical Union (IAU) and the International Association of Geodesy (IAG) (Seidelmann and others, 2002), these inertial coordinates were converted into the planet-fixed coordinates (longitude and latitude) used on this map. These values include the orientation of the north pole of Mars (including the effects of precession), the rotation rate of Mars, and a value for W0 of 176.630°, where W0 is the angle along the equator to the east, between the 0° meridian and the equator's intersection with the celestial equator at the standard epoch J2000.0 (Seidelmann and others, 2002). This value of W0 was chosen (Duxbury and others, 2002) in order to place the 0° meridian through the center of the small (~500 m) crater Airy-0, located in the crater Airy (de Vaucouleurs and others, 1973; Seidelmann and others, 2002). Longitude increases to the east, and latitude is planetocentric as allowed by IAU/IAG standards (Seidelmann and others, 2002) and in accordance with current NASA and USGS standards (Duxbury and others, 2002). A secondary grid (printed in red) has been added to the map as a reference to the west longitude/planetographic latitude system that is also allowed by IAU/IAG standards (Seidelmann and others, 2002) and has also been used for Mars. The figure adopted to compute this secondary grid is an oblate spheroid with an equatorial radius of 3,396.19 km and a polar radius of 3,376.2 km (Duxbury and others, 2002; Seidelmann and others, 2002). MAPPING TECHNIQUES To create the topographic base image, the original DEM produced by the MOLA team in Simple Cylindrical projection with a resolution of 64 pixels per degree was projected into the Mercator and Polar Stereographic pieces. A shaded relief was generated from each DEM with a sun angle of 30° from horizontal and a sun azimuth of 270°, as measured clockwise from north, and a vertical exaggeration of 100%. Illumination is from the west, which follows a long-standing USGS tradition for planetary maps. This allows for continuity in the shading between maps and quadrangles, and most closely resembles lighting conditions found on imagery. The DEM values were then mapped to a smooth global color look-up table. Note that the chosen color scheme simply represents elevation changes and is not intended to imply anything about surface characteristics (e.g., past or current presence of water or ice). These two files were then merged and scaled to 1:25 million for the Mercator portion and 1:15,196,708 for the two Polar Stereographic portions, with a resolution of 300 dots per inch. The projections have a common scale of 1:13,923,113 at +56° latitude. NOMENCLATURE Names on this sheet are approved by the IAU and have been applied for features clearly visible at the scale of this map. For a complete list of the IAU-approved nomenclature for Mars, see the Gazetteer of Planetary Nomenclature at http://planetarynames.wr.usgs.gov/. Font color was chosen for readability. Names followed by an asterisk are provisionally approved. *M 25M RKN: Abbreviation for Mars, 1:25,000,000 series, shaded relief (R), with color (K) and nomenclature (N) (Greeley and Batson, 1990)

Open-File Report↗
USGS2025

Grammar to graph—An approach for semantic transformation of annotations to triples

Data annotation is the process of labeling data to show the outcome that a related data model should predict. In this study, annotation data were transformed into semantic graph triples, mainly for use with the Resource Description Framework (RDF), a type of entity-relationship-attribute data model for graph databases. The transformation of annotation data to semantic graph triples provides complex linguistic meaning with data handling advantages such as reduced data storage needs, improved logical specification of relations between objects, and reusable classes and properties that support logic and inference. A grammar-based framework in graph form supports user questions and queries. The words defining approximately 334 topographic feature types compiled by the U.S. Geological Survey were tokenized as units of analysis and grouped by part of speech. Their dependency relations were identified for this study using natural language processing libraries. Dependency concepts are used as structured semantic relations among part-of-speech classes. Tokens, units equivalent to words, form instances of classes and were quantified within a tabular output format using PostgreSQL data storage software. Table data were logically aligned as triples following a mapping file and stored with an ontology file using Ontop virtual triplestore software. A grammar ontology schema for the data was synchronized to match queries whose results validated the graph’s structure. The text analysis produced 8 part-of-speech classes of content words for object representations and 4 classes of function words for operational applications. Dependency relations formed 27 ontology properties for topographic subgraph structures. Token occurrences shaped overall ontology salience and formed a lexicon of syntactic terms for subgraph objects and properties. The schema ontology of class and property population shapes formed the lexicon of English terms. SPARQL Protocol and RDF Query Language (SPARQL) was used with the lexicon to conform data to RDF guidelines. This study confirms the hypothesis that although linguistic logic varies from description logic, its approximation applies to ontology design. Property and query use case patterns extracted from the analysis support queries concerning complex topographic relations and patterns normally embedded within text definitions. The method used in this study could be applied to text forms in other domains, such as survey notes.

Scientific Investigations Report↗
USGS2016

Geology and geologic history of the Moscow-Pullman basin, Idaho and Washington, from late Grande Ronde to late Saddle Mountains time

The Moscow-Pullman basin, located on the eastern margin of the Columbia River flood basalt province, consists of a subsurface mosaic of interlayered Miocene sediments and lava flows of the Imnaha, Grande Ronde, Wanapum, and Saddle Mountains Basalts of the Columbia River Basalt Group. This sequence is ~1800 ft (550 m) thick in the east around Moscow, Idaho, and exceeds 2300 ft (700 m) in the west at Pullman, Washington. Most flows entered from the west into a topographic low, partially surrounded by steep mountainous terrain. These flows caused a rapid rise in base level and deposition of immature sediments. This field guide focuses on the upper Grande Ronde Basalt, Wanapum Basalt, and sediments of the Latah Formation. Late Grande Ronde flows terminated midway into the basin to begin the formation of a topographic high that now separates a thick sediment wedge of the Vantage Member to the east of the high from a thin layer to the west. Disrupted by lava flows, streams were pushed from a west-flowing direction to a north-northwest orientation and drained the basin through a gap between steptoes toward Palouse, Washington. Emplacement of the Roza flow of the Wanapum Basalt against the western side of the topographic high was instrumental in this process, plugging west-flowing drainages and increasing deposition of Vantage sediments east of the high. The overlying basalt of Lolo covered both the Roza flow and Vantage sediments, blocking all drainages, and was in turn covered by sediments interlayered with local Saddle Mountains Basalt flows. Reestablishment of west-flowing drainages has been slow. The uppermost Grande Ronde, the Vantage, and the Wanapum contain what is known as the upper aquifer. The water supply is controlled, in part, by thickness, composition, and distribution of the Vantage sediments. A buried channel of the Vantage likely connects the upper aquifer to Palouse, Washington, outside the basin. This field guide locates outcrops; relates them to stratigraphic well data; outlines paleogeographic basin evolution from late Grande Ronde to the present time; and notes structures, basin margin differences, and features that influence upper aquifer water supply.

Idaho, Nevada, Oregon, Washington↗
USGS2002

CMGTooL user's manual

During the past several years, the sediment transport group in the Coastal and Marine Geology Program (CMGP) of the U. S. Geological Survey has made major revisions to its methodology of processing, analyzing, and maintaining the variety of oceanographic time-series data. First, CMGP completed the transition of the its oceanographic time-series database to a self-documenting NetCDF (Rew et al., 1997) data format. Second, CMGP’s oceanographic data variety and complexity have been greatly expanded from traditional 2-dimensional, single-point time-series measurements (e.g., Electro-magnetic current meters, transmissometers) to more advanced 3-dimensional and profiling time-series measurements due to many new acquisitions of modern instruments such as Acoustic Doppler Current Profiler (RDI, 1996), Acoustic Doppler Velocitimeter, Pulse-Coherence Acoustic Doppler Profiler (SonTek, 2001), Acoustic Bacscatter Sensor (Aquatec, 1001001001001001001). In order to accommodate the NetCDF format of data from the new instruments, a software package of processing, analyzing, and visualizing time-series oceanographic data was developed. It is named CMGTooL. The CMGTooL package contains two basic components: a user-friendly GUI for NetCDF file analysis, processing and manipulation; and a data analyzing program library. Most of the routines in the library are stand-alone programs suitable for batch processing. CMGTooL is written in MATLAB computing language (The Mathworks, 1997), therefore users must have MATLAB installed on their computer in order to use this software package. In addition, MATLAB’s Signal Processing Toolbox is also required by some CMGTooL’s routines. Like most MATLAB programs, all CMGTooL codes are compatible with different computing platforms including PC, MAC, and UNIX machines (Note: CMGTooL has been tested on different platforms that run MATLAB 5.2 (Release 10) or lower versions. Some of the commands related to MAC may not be compatible with later releases of MATLAB). The GUI and some of the library routines call low-level NetCDF file I/O, variable and attribute functions. These NetCDF exclusive functions are supported by a MATLAB toolbox named NetCDF, created by Dr. Charles Denham . This toolbox has to be installed in order to use the CMGTooL GUI. The CMGTooL GUI calls several routines that were initially developed by others. The authors would like to acknowledge the following scientists for their ideas and codes: Dr. Rich Signell (USGS), Dr. Chris Sherwood (USGS), and Dr. Bob Beardsley (WHOI). Many special terms that carry special meanings in either MATLAB or the NetCDF Toolbox are used in this manual. Users are encouraged to read the documents of MATLAB and NetCDF for references.

Open-File Report↗
USGS1989

Heavy-mineral reconnaissance off the coast of the Apalachicola river delta, northwest Florida: A summary and new interpretations

Two-hundred-and-fifty sediment samples were collected for heavy-mineral and textural analysis along the northwest Florida coastline from approximately 24 km offshore of Apalachee Bay to the same distance offshore of Pensacola Bay. The heavy-mineral suite characterizing sediments within this region consists of opaque minerals, kyanite, staurolite, tourmaline, zircon and rutile. Minor constituents of this suite include epidote, sphene, amphibole, sillimanite, garnet and leucoxene. The average heavy-mineral concentration within these sediments is approximately 0.12 wt.%. Specifically, the 2 to 3∅ grain-size interval contains an average of 0.51 wt.%, whereas the 3 to 4∅ interval contains an average of 4.39 wt.% heavy minerals. Note that the 3 to 4∅ interval typically represents only 4% of the sample volume. There is a general westward increase in heavy-mineral concentrations throughout the study area. Superimposed on this regional trend, areas of maximum heavy-mineral concentration occur within sediments offshore of St. George and Santa Rosa islands. The primary source of sediments in the region is the crystalline rocks of the southern Appalachians. Granulometric analyses of these sediments reveal a westward increase in values of sample mean grain size, and decrease in standard deviation, and percent fines. It is postulated from these data, in addition to the interpretation of sample grain-size distributions, skewness, and kurtosis, that these inner continental shelf sediments are primarily fluvial in origin. These sediments have been transported to the shelf by the Apalachicola and surrounding major rivers during Pleistocene low sea-level stands. Data also indicate evidence of reworking by coastal or marine offshore wave processes.

Marine Geology↗
USGS1981

Quality of surface water at selected sites in the Suwannee River basin, Florida, April 1978 to December 1979

This report presents the results of analyses of water-quality samples collected from 13 surface-water sites in the Suwannee River basin in Florida from April 1978 to December 1979. The U.S. Geological Survey, in a cooperative investigative program with the Suwannee River Authority, has sampled a network of nine water-quality monitoring stations since 1968 to detect short- and long-term changes in surface-water quality in the Suwannee River. Two additional stations, Camp Branch near Genoa and Rocky Creek near Houston, were established in October 1978. In addition to this 11-station network, Deep Creek and Robinson Creek near Suwannee Valley, which originate in the Osceola National Forest and have been supported since April 1976 through a cooperative agreement with the U.S. Bureau of Land Management, have been included in this report. The analyses of samples collected routinely included: nutrients, total organic carbon, and 5-day biological oxygen demand, bimonthly; trace metals, annually. The array of constituents sampled was expanded in October 1978 at three of the original nine stations to provide quality-of-water information for streams draining an industrial area: Rocky Creek near Belmont, Hunter Creek near Belmont, and Swift Creek at Facil. Data collected at these three sites now include: major chemical constituents, six times per year; radium-226, two times per year; and trace metals, one time per year. These constituents are determined in addition to nutrients, total organic carbon, and biological oxygen demand which continue to be collected six times per year. All results of analyses of the water-quality samples collected from April 1978 to December 1979 remained within, or near, previously measured ranges and water-quality fluctuations were similar to those noted from data collected since 1970.

Florida↗
USGS2022

Sources and characteristics of dissolved organic carbon in the McKenzie River, Oregon, related to the formation of disinfection by-products in treated drinking water

Executive Summary This study characterized the concentration and quality of dissolved organic carbon (DOC) in the McKenzie River, a relatively undeveloped watershed in western Oregon, and its link to forming disinfection by-products (DBPs) in treated drinking water. The study aimed to identify the primary source(s) of DOC in source water for the Eugene Water & Electric Board’s (EWEB) conventional treatment plant on the McKenzie River near river mile 11, upstream of Hayden Bridge. The two classes of regulated compounds examined—trihalomethanes (THMs) and haloacetic acids (HAAs)—form when organic carbon in raw source water reacts with chlorine and (or) bromine during water treatment. The objectives of the study were to: characterize the amount and quality of DOC in the McKenzie River and select tributaries during storms; identify the most common types of carbon using UV-vis spectroscopy and other methods; evaluate optical properties for predicting DBP precursors in surface water; and identify land cover classes or vegetation types that may be important sources of organic carbon and DBP precursors in EWEB’s source water. Eleven storms were sampled synoptically in upstream-to-downstream fashion to provide a “snapshot” of water quality conditions at four sites on the McKenzie River from Frissell Bridge (6 miles downstream from Trail Bridge Reservoir) to the EWEB water treatment plant at Hayden Bridge and nine contributing tributaries. Storms included late summer and early autumn “first flush” events and late autumn, winter, and spring storms spanning a range in streamflows from 3,000 to 26,000 cubic feet per second as measured in the main stem McKenzie River at the EWEB water intake. Water samples were analyzed for DOC concentrations and optical properties (fluorescence and ultraviolet absorbance [UVA]) across a range of wavelengths to characterize the quantity and quality of dissolved organic matter (DOM) in the McKenzie River at the drinking water intake and upstream locations. Paired sets of source and finished water samples were collected at the EWEB treatment plant to identify DOC quality parameters in raw source water that might predict DBP concentrations in finished drinking water. DOC concentrations were relatively low in the McKenzie River (0.4–3 milligrams per liter [mg/L]; average 1.5 mg/L) but much higher in the tributaries. The highest DOC concentrations occurred during “first flush” storms in October 2012 and September 2013; the highest value (16 mg/L) was measured at the 52nd Street stormwater outfall. The average DOC concentration in the lower basin-tributaries was 3.8 mg/L; three middle basin tributaries—Quartz, Gate, and Haagen Creeks, which drain private forestland with less coniferous forest compared with other higher elevation tributaries— had slightly lower average DOC concentrations (2.8 mg/L). These middle-basin watersheds may be important sources of DOC and DBP precursors to the McKenzie River, even more so than the lower basin tributaries, depending on their flows (and loads). This is particularly true after the September 2020 Holiday Farm fire, which burned much of this area. DOC concentrations increased 68 percent in the McKenzie River between the uppermost reference site at Frissell Bridge and Vida; this includes drainage from Quartz Creek, Blue River Lake and Cougar Reservoir, which all contributed DOC to the main stem. In contrast, the lowermost tributaries draining most of the agricultural and urban land did not have a large effect on DOC in the McKenzie River despite their higher DOC concentrations because of their presumed relatively low streamflows and, consequently, DOC loads. Apart from the continuous flow monitors in the McKenzie River and some tributaries (Blue River and South Fork McKenzie River, and streamflow at Hayden Bridge and Vida, Camp Creek and some other locations), streamflow was not assessed during sample collection for this study. This lack of streamflow data precludes a detailed analysis of loads, which is discussed in the future studies section. All DBP concentrations in finished drinking water were less than EPA maximum contaminant levels (MCLs) of 0.080 mg/L for the four trihalomethanes (THM4) and 0.060 mg/L for five haloacetic acids (HAA5). During the 11 storm sampling events the maximum summed concentrations were about 0.040 mg/L for both THM4 and HAA5. Compliance monitoring samples, collected separately by EWEB, yielded some higher concentrations—0.046 mg/L THM4 and 0.047 HAA5—during the December 2012 storm. The corresponding benchmark quotient (BQ) values, which indicate how close a measured DBP concentration is to the MCL, were 0.58 and 0.78, respectively, for THM4 and HAA5. Compared with a similar 2007–08 McKenzie River study that did not target storm events, concentrations of THM4 and HAA5 in finished water were 68 percent and 33 percent higher, respectively, during the current study. Due to the high dilution rates in the McKenzie River main stem, many of the individual fluorescence excitation-emission measurements were low (<0.1 Raman units) and approached analytical detection limits. Parallel factor analysis (PARAFAC) resulted in a five-component model (C1–C5) that represents five unique organic fluorophores. Components C1, C2, and C3 represent DOM associated with soil-derived, humic-like, more degraded organic matter. In contrast, components C4 and C5 represent “fresher” DOM, derived from terrestrial and aquatic plants, including algae and cyanobacteria that are common in the McKenzie River and its tributaries and reservoirs. The fluorescence data and PARAFAC modeling suggest that most of the DOC in the McKenzie River originated from terrestrial sources (primarily components C1 and C2). The largest increases in DOC in the main stem occurred in the reach upstream of Vida, from inflows by Quartz Creek, Blue River, South Fork McKenzie River, and other tributaries. Concentrations of DBPs in EWEB’s finished drinking water were positively correlated with DOC concentrations in raw source water (THM4, p <0.05; HAA5, p <0.01) for paired samples collected 12−24 hours apart. DOC concentrations were significantly positively correlated ( p <0.001) with laboratory-based fluorescent dissolved organic matter (fDOM) measurements, suggesting fDOM as a useful parameter for monitoring and predicting DOC concentration in surface water and DBP concentrations in finished water. Of all the PARAFAC components in surface water, C5 had the highest correlations with DBPs in finished water (rho = 0.77–0.84, p <0.01), followed by components C1 and C2 (rho = 0.75 and 0.71, respectively, p <0.01). This C5 carbon is associated with recently produced DOM, possibly from decomposed terrestrial and aquatic vegetation. Model loadings of these three components were considerably higher in the sampled tributaries relative to the main stem McKenzie River, with most of the observed increases in the main stem apparent at Vida. This points to Quartz Creek or other tributaries in the reach between Frissell Bridge and the sampling site near Vida (South Fork McKenzie and Blue Rivers) as potentially key contributors of DOM source material that leads to the production of DBPs in treated drinking water. A limited load analysis showed that the reservoirs contributed 8–37 percent of the instantaneous DOC loads observed at Vida at the time of sampling, which suggests other sources such as Quartz Creek and other streams in the reach between Frissell Bridge and Vida are more important. Random forest analyses identified PARAFAC components C1 and C5 and fluorescence peaks A, C, M, T and N as the best predictors for HAA5 concentrations in finished drinking water, explaining 62.5 percent of the variation. The best predictors for THM4 were C1, C4 + C5, and peaks T, A, and N, which explained 33 percent of the variation. Several land cover and vegetation classes were correlated with DOC concentration and other optical measurements. The percentage of evergreen forest in each of the subwatersheds sampled was negatively correlated ( p <0.001) with DOC concentration and many optical indicators of DOM quantity: UVA 254 , fDOM, and all of the fluorescence peaks. In contrast, mixed (deciduous) forest was positively correlated ( p <0.001) with DOC, fDOM, UVA 254 , and several fluorescence peaks, demonstrating the importance of deciduous leaf fall in generating DOC and DBP precursors. The high level of human activities in the middle and lower portion of the basin—including timber harvesting and road construction on private forestland, agricultural, rural, industrial, and urban development—have resulted in the greatest loss in native coniferous and mixed deciduous forests in the basin. DOC loading from these tributaries and reservoir releases, which contain DOC from terrestrial and aquatic productivity, both enrich the McKenzie River. Concentrations of DOC increased an average of 71 percent (range 30–120 percent) in the McKenzie River between Frissell Bridge, the upstream reference site, and Vida. PARAFAC components C1, C2, and C5—which were correlated with DBPs in finished water—increased, on average, 109–136 percent (range 20–250 percent) in this same Frissell-to-Vida reach. These increases occur from input of tributaries in the middle basin such as Quartz Creek and others, as noted above. Future monitoring, field, and lab studies can improve our understanding of seasonal and spatial sources of organic carbon contributing DBP precursors to the McKenzie River and allow detection of long-term trends resulting from the recent Holiday Farm Fire, which burned 173,393 acres of forestland, including riparian areas along the main stem, and numerous structures, homes, and outbuildings in September 2020. Future studies could examine DOC fluxes and flushing of carbon from the watershed, investigate the role of precipitation amount and intensity in mobilizing carbon and sediment, and evaluate impacts to aquatic communities and human health as part of a post-fire assessment. Other areas ripe for study include evaluating the impacts of potential temperature increases on carbon sequestration and decomposition in the burned and unburned forests and identifying practices that foster sequestration of carbon in forest soils. The use of fluorescence sensors such as fDOM to monitor the concentration and composition of raw water supplies may be improved for detection of specific DBP precursors, to provide continuous and real-time information to treatment plant operators. Future studies that monitor DOM amount and quality, and DBP Formation Potential (FP), particularly during storm events, paired with streamflow measurements, as suggested above, could help identify areas that contribute high DOC loads and thus help managers identify the key areas to focus restoration activities. Other studies could examine treatment options for currently regulated DBPs and potentially unregulated compounds, including advanced biological treatments for their removal. This study was a collaboration between the U.S. Geological Survey (USGS) and EWEB in Eugene, Oregon, with additional funding provided from USGS Cooperative Matching Funds Program.

Oregon↗
USGS1975

Notes on the origin of Colluma Crater, Bolivia

Colluma Crater, on the semiarid Altiplano (high plain) of western Bolivia, is an oval structure having overall dimensions of 6.7 km by 6.0 km. The structure has two almost concentric cuestaform rims (the inner rim is 3.6 km by 3.1 km) composed of poorly consolidated clastic sediments that dip outward. The center of the crater is about 80 m below the maximum height of the rims and about at the altitude of the surrounding plain. Because of the double rim, the centripetal drainage, and the absence of volcanic rocks, this structure is considered a collapsed dome. We believe it was probably formed by the doming of lower Quaternary (?) sediments by a subjacent igneous intrusion, partial retreat of the magma, collapse of the central part of the dome, erosional etching of the two rims, and partial filling of the center by detritus from the walls. Evidence for origin by impact (nickel-iron materials, shock structures, ejecta, and so forth) is lacking or was unrecognized, but this mode of origin is not rejected at this time. Geophysical surveys are recommended to determine whether the structure continues in depth and if an igneous, plug is below the crater.

Colluma Crater↗
USGS1971

Hydrogeologic data for the Upper Connecticut River Basin, Connecticut

This report contains geologic, ground-water, and quality-of-water data collected and compiled for a water resources investigation of the upper Connecticut River basin, Connecticut by the U.S. Geological Survey in financial cooperation with the Connecticut Department of Environmental Protection. These data, together with surface-water data, were collected at sites shown on plate A and are, with exceptions discussed below, presented herein. The interpretation of these data will be published separately in Connecticut Water Resources Bulletin No. 24. As defined for this investigation, the upper Connecticut River basin is a 508-square-mile area in north-central Connecticut (see front cover and plate A) drained by the Connecticut River and its tributaries including the Farmington River downstream from the community of Tariffville. The basin includes all or parts of the following 30 towns: Avon, Bloomfield, Bolton, Cromwell, East Granby, East Hampton, East Hartford, East Windsor, Ellington, Enfield, Farmington, Glastonbury, Hartford, Manchester, New Brltain, Newington, Plainville, Portland, Rocky Hill, Simsbury, Somers, South Windsor, Stafford, Suffield, Tolland, Vernon, West Hartford, Wethersfield, Windsor, and Windsor Locks. The data compiled for this investigation consist of those collected specifically for the study during the period July 1967 - December 1969 and those collected during previous investigations. Many records of wells and test holes in thus report were published in Connecticut Water Resources Bulletin No. 4. This report is out of print, therefore these previously published records, together with newer ones, are published in this report to provide a readily available source of ground-water information. Previously published well numbers have been retained and are those published in this report. However, test holes previously published as wells have been assigned test-hole numbers; the former well number Is noted in table 3.

Connecticut↗
USGS2012

Lidar base specification

In late 2009, a $14.3 million allocation from the American Recovery and Reinvestment Act (ARRA) for new light detection and ranging (lidar) elevation data acquisition prompted the U.S. Geological Survey (USGS) National Geospatial Program (NGP) to develop a common minimum specification for all lidar data acquired for The National Map. Released as a working draft in 2010 and formally published in 2012, the USGS–NGP Lidar Base Specification (LBS) was quickly embraced by numerous States, counties, and foreign countries as the foundation for their own lidar specifications. Prompted by a growing appreciation for the wide applicability and inherent value of lidar, a consortium of Federal agencies commissioned the National Enhanced Elevation Assessment (NEEA) study in 2010 to quantify the costs and benefits of a national lidar program. Published in 2012, the NEEA report documented a substantial return on such an investment, defined five quality levels (QL) for elevation data, and recommended an 8-year collection cycle of QL2 lidar data as the optimum balance of benefit and affordability. In response to the study, the USGS–NGP established the 3D Elevation Program (3DEP) in 2013 as the interagency vehicle through which the NEEA recommendations could be realized. Lidar is a quickly evolving technology and much has changed in the industry since the previous version of the Lidar Base Specification (LBS) was published. Lidar data have improved in accuracy and spatial resolution, the American Society for Photogrammetry and Remote Sensing has revised the geospatial accuracy standards, industry standard file formats have been expanded, additional applications for lidar have become accepted, and the need for interoperable data across collections has been realized. This revision to the LBS addresses some of those changes and provides continued guidance towards a nationally consistent lidar dataset. Note Future versions of the Lidar base specification will be released online at https://www.usgs.gov/3DEP/lidarspec .

Techniques and Methods↗
USGS2013

Determination of flow losses in the Cape Fear River between B. Everett Jordan Lake and Lillington, North Carolina, 2008-2010

During 2008-2010, the U.S. Geological Survey conducted a hydrologic investigation in cooperation with the Triangle J Council of Governments Cape Fear River Flow Study Committee and the North Carolina Division of Water Resources to collect hydrologic data in the Cape Fear River between B. Everett Jordan Lake and Lillington in central North Carolina to help determine if suspected flow losses occur in the reach. Flow loss analyses were completed by summing the daily flow releases at Jordan Lake Dam with the daily discharges at Deep River at Moncure and Buckhorn Creek near Corinth, then subtracting these values from the daily discharges at Cape Fear River at Lillington. Examination of long-term records revealed that during 10,227 days of the 1983-2010 water years, 408 days (4.0 percent) had flow loss when conditions were relatively steady with respect to the previous day's records. The flow loss that occurred on these 40 days ranged from 0.49 to 2,150 cubic feet per second with a median flow loss of 37.2 cubic feet per second. The months with the highest number of days with flow losses were June (16. percent), September (16.9 percent), and October (19.4 percent). A series of synoptic discharge measurements made on six separate days in 2009 provided "snapshots" of overall flow conditions along the study reach. The largest water diversion is just downstream from the confluence of the Haw and Deep Rivers, and discharges substantially decrease in the main stem downstream from the intake point. Downstream from Buckhorn Dam, minimal gain or loss between the dam and Raven Rock State Park was noted. Analyses of discharge measurements and ratings for two streamgages-one at Deep River at Moncure and the other at Cape Fear River at Lillington-were completed to address the accuracy of the relation between stage and discharge at these sites. The ratings analyses did not indicate a particular time during the 1982-2011 water years in which a consistent bias occurred in the computations of discharge records that would indicate false flow losses. A total of 34 measured discharges at a streamgage on the Haw River below B. Everett Jordan Lake near Moncure were compared with the reported hourly flow releases from Jordan Lake Dam. Because 28 of 34 measurements were within plus or minus 10 percent of the hourly flow releases reported by the U.S Army Corps of Engineers, use of the current discharge computation tables for reporting Jordan Lake Dam flow releases is generally supported. A stage gage was operated on the Cape Fear River at Buckhorn Dam near Corinth to collect continuous stage-only records. Throughout the study period, flow over the dam was observed along its length, and flow loss within the study reach is not attributed to river-level fluctuations at the dam. Water-use information and (or) data were obtained for five industrial facilities, a regional power utility, two municipalities, one small hydropower facility on the Deep River, and one quarry operation also adjacent to the Deep River. The largest water users are the regional power producer, a small hydropower operation, and the two municipalities. The total water-use diversions for these facilities range from almost 25.5 to 38.5 cubic feet per second (39.5 to 59.5 million gallons per day) during the winter and summer periods, respectively. This range is equivalent to 69 to 104 percent of the 37 cubic feet per second median flow loss. The Lockville hydropower station is on the Deep River about 1 mile downstream from the streamgage near Moncure. Run-of-river operations at the facility do not appear to affect flow losses in the study reach. The largest water user in the study area is a regional power producer at a coal-fired power-generation plant located immediately adjacent to the Cape Fear River just downstream from the confluence of the Haw an Deep Rivers. Comparisons of daily water withdrawals, sup-plied by the regional power producer, and discharge records at a streamgage on the diversion canal indicated many days when consumption exceeded the producer's estimates for the cooling towers. Uncertainty surrounding reasonable estimates of consumption remained in effect at the end of the study. Data concerning evaporative losses were compiled using two approaches-an analysis of available pan-evaporation data from a National Weather Service cooperative observer station in Chapel Hill, North Carolina; and a compilation of reference open-water evaporation computed by the State Climate Office of North Carolina. The potential flow loss by evaporation from the main stem and the Deep River was estimated to be in the range of 4 to 14 cubic feet per second during May through October, equivalent to 10 to 38 percent of the 37 cubic feet per second median flow loss. Daily water-use diversions and evaporation losses were compared to flow-loss occurrences during the period April 2008 through September 2010. In comparing the surface-water, water-use, and evaporation data compiled for 2008-2010, it is evident that documented water diversions combined with flow losses by open-water evaporation can exceed the net flow gain in the study area and result in flow losses from the reach. Analysis of data from a streamgage downstream from the regional power plant on the diversion canal adjacent to the Cape Fear River provided insight into the occurrence of an apparent flow loss at the streamgage at Lillington. Assessment of the daily discharges and subsequent hydrographs for the canal streamgage indicated at least 24 instances during the study when the flows suddenly changed by magnitudes of 100 to more that 200 cubic feet per second, resulting in a noted time-lag effect on the downstream discharges at the Lillington streamgage, beginning 8 to 16 hours after the sudden flow change. A fiber-optic distributed temperature-sensing survey was conducted on the Cape Fear River at the Raven Rock State Park reach August 12-14, 2009, to determine if the presence of diabase dikes were preferentially directing groundwater discharge. No temperature anomalies of colder water were measured during the survey, which indicated that at the time of the survey that particular reach of the Cape Fear River was a "no-flow" or losing stream. An aerial thermal-infrared survey was conducted on the Haw and Cape Fear Rivers on February 27, 2010, from Jordan Lake Dam to Lillington to qualitatively delineate areas of groundwater discharge on the basis of the contrast between warm groundwater discharge and cold surface-water temperatures. Dis-charge generally was noted as diffuse seepage, but in a few cases springs were detected as inflow at a discrete point of discharge. Two reaches of the Cape Fear River (regional power plant and Bradley Road reaches) were selected for groundwater monitoring with a transect of piezometers installed within the flood plain. Groundwater-level altitudes at these reaches were analyzed for 1 water year (October 1, 2009, to September 30, 2010). Data collected as part of this study represent only a brief period of time and may not represent all conditions and all years; however, the data indicate that, during the dry summer months, the Cape Fear River within the study area is losing an undetermined quantity of water through seepage. Analyses completed during this investigation indicate a study reach with complex flow patterns affected by numerous concurrent factors resulting in flow losses. The causes of flow loss could not be solely attributed to any one factor. Among the factors considered, the occurrences of water diversions and evaporative losses were determined to be sufficient on some days (particularly during the base-flow period) to exceed the net gain in flows between the upstream and downstream ends of the study area. Losses by diversions and evaporation can exceed the median flow loss of 3 cubic feet per second, which indicates that flow loss from the study reach is real. Groundwater data collected during 2009-2010 indicate the possibility of localized flow loss during the summer, particularly in the impounded reach above Buckhorn Dam. However, no indication of unusual patterns was noted that would cause substantial flow loss by groundwater and surface-water interaction at the river bottom.

North Carolina↗