Search USGSSearch

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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 829 records · Page 46Linked to original sources

Development of a regionally consistent geospatial dataset of agricultural lands in the Upper Colorado River Basin, 2007-10

Irrigation in arid environments can alter the natural rate at which salts are dissolved and transported to streams. Irrigated agricultural lands are the major anthropogenic source of dissolved solids in the Upper Colorado River Basin (UCRB). Understanding the location, spatial distribution, and irrigation status of agricultural lands and the method used to deliver water to agricultural lands are important to help improve the understanding of agriculturally derived dissolved-solids loading to surface water in the UCRB. Irrigation status is the presence or absence of irrigation on an agricultural field during the selected growing season or seasons. Irrigation method is the system used to irrigate a field. Irrigation method can broadly be grouped into sprinkler or flood methods, although other techniques such as drip irrigation are used in the UCRB. Flood irrigation generally causes greater dissolved-solids loading to streams than sprinkler irrigation. Agricultural lands in the UCRB mapped by state agencies at varying spatial and temporal resolutions were assembled and edited to represent conditions in the UCRB between 2007 and 2010. Edits were based on examination of 1-meter resolution aerial imagery collected between 2009 and 2011. Remote sensing classification techniques were used to classify irrigation status for the June to September growing seasons between 2007 and 2010. The final dataset contains polygons representing approximately 1,759,900 acres of agricultural lands in the UCRB. Approximately 66 percent of the mapped agricultural lands were likely irrigated during the study period.

Arizona, Colorado, New Mexico, Utah, Wyoming

Titanium mineral resources in heavy-mineral sands in the Atlantic coastal plain of the southeastern United States

This study examined titanium distribution in the Atlantic Coastal Plain of the southeastern United States; the titanium is found in heavy-mineral sands that include the minerals ilmenite (Fe 2+ TiO 3 ), rutile (TiO 2 ), or leucoxene (an alteration product of ilmenite). Deposits of heavy-mineral sands in ancient and modern coastal plains are a significant feedstock source for the titanium dioxide pigments industry. Currently, two heavy-mineral sands mining and processing operations are active in the southeast United States producing concentrates of ilmenite-leucoxene, rutile, and zircon. The results of this study indicate the potential for similar deposits in many areas of the Atlantic Coastal Plain. This study used the titanium analyses of 3,457 stream sediment samples that were analyzed as part of the U.S. Geological Survey’s National Geochemical Survey program. This data set was analyzed by an integrated spatial modeling technique known as Bayesian hierarchical modeling to map the regional-scale, spatial distribution of titanium concentrations. In particular, clusters of anomalous concentrations of titanium occur: (1) along the Fall Zone, from Virginia to Alabama, where metamorphic and igneous rocks of the Piedmont region contact younger sediments of the Coastal Plain; (2) a paleovalley near the South Carolina and North Carolina border; (3) the upper and middle Atlantic Coastal Plain of North Carolina; (4) the majority of the Atlantic Coastal Plain of Virginia; and (5) barrier islands and stretches of the modern shoreline from South Carolina to northeast Florida. The areas mapped by this study could help mining companies delimit areas for exploration.

Scientific Investigations Report

Characterization of peak streamflows and flooding in select areas of Pennsylvania from the remnants of Hurricane Ida, September 1–2, 2021

Pennsylvania experienced heavy rainfall on September 1 and 2, 2021, as the remnants of Hurricane Ida swept over parts of the State. Much of eastern and south-central Pennsylvania received 5 to 10 inches of rain, and most of the rainfall fell within little more than 6 hours. Southeastern Pennsylvania experienced widespread, substantial flooding, and the city of Philadelphia and surrounding areas were particularly affected by the flooding. U.S. Geological Survey (USGS) streamgages registered peak streamflows of record at 19 locations, and 52 locations experienced top 5 peak streamflows for the period of record and an annual exceedance probability estimate of at least 10 percent. During this September 2021 flood event, USGS personnel made over 60 streamflow measurements at streamgages in Pennsylvania using direct and indirect methods. Many of those streamflow measurements were made to verify or improve the accuracy, extent, or development of new stage-streamflow relations at streamgages operated by the USGS. After the floodwaters receded, USGS personnel identified and documented a total of 338 high-water marks in Pennsylvania, noting such things as their general description, location, height above land surface, and quality. Many of these high-water marks were used to create five flood-documentation maps for selected communities in southeastern Pennsylvania that experienced substantial flooding because of the remnants of Hurricane Ida. Digital datasets of the inundated areas, mapped boundaries, and water depth are available (Stuckey and Conlon, 2023).

Pennsylvania

A model to locate potential areas for lake sturgeon spawning habitat construction in the St. Clair–Detroit River System

In response to a need for objective scientific information that could be used to help remediate loss of fish spawning habitat in the St. Clair River and Detroit River International Areas of Concern, this paper summarizes a large-scale geographic mapping investigation. Our study integrates data on two variables that many riverine fishes respond to in selecting where to spawn in these waters (water flow velocity and water depth) with available maps of the St. Clair–Detroit River System (SC–DRS). Our objectives were to locate and map these two physical components of fish habitat in the St. Clair and Detroit rivers and Lake St. Clair using a geographic information system (GIS) and to identify where, theoretically, fish spawning habitat could be remediated in these rivers. The target fish species to which this model applies is lake sturgeon (Acipenser fulvescens), but spawning reefs constructed for lake sturgeon in this system have been used for spawning by 17 species of fish. Our analysis revealed areas in each river that possessed suitable water velocity and depth for fish spawning and therefore could theoretically be remediated by the addition of rock-rubble substrate like that used at two previously remediated sites in the Detroit River at Belle Isle and Fighting Island. Results of our analysis revealed that only 3% of the total area of the SC–DRS possesses the necessary combination of water depth and high flow velocity to be indicated by the model as potential spawning habitat for lake sturgeon.

Michigan, Ontario

Occurrence, Distribution, Sources, and Trends of Elevated Chloride Concentrations in the Mississippi River Valley Alluvial Aquifer in Southeastern Arkansas

Water-quality data from approximately 2,500 sites were used to investigate the distribution of chloride concentrations in the Mississippi River Valley alluvial aquifer in southeastern Arkansas. The large volume and areal distribution of the data used for the investigation proved useful in delineating areas of elevated (greater than 100 milligrams per liter) chloride concentrations, assessing potential sources of saline water, and evaluating trends in chloride distribution and concentration over time. Irrigation water containing elevated chloride concentrations is associated with negative effects to rice and soybeans, two of the major crops in Arkansas, and a groundwater chloride concentration of 100 milligrams per liter is recommended as the upper limit for use on rice. As such, accurately delineating areas with high salinity ground water, defining potential sources of chloride, and documenting trends over time is important in assisting the agricultural community in water management. The distribution and range of chloride concentrations in the study area revealed distinct areas of elevated chloride concentrations. Area I includes an elongated, generally northwest-southeast trending band of moderately elevated chloride concentrations in the northern part of the study area. This band of elevated chloride concentrations is approximately 40 miles in length and varies from approximately 2 to 9 miles in width, with a maximum chloride concentration of 360 milligrams per liter. Area II is a narrow, north-south trending band of elevated chloride concentrations in the southern part of the study area, with a maximum chloride concentration of 1,639 milligrams per liter. A zone of chloride concentrations exceeding 200 milligrams per liter is approximately 25 miles in length and 5 to 6 miles in width. In Area I, low chloride concentrations in samples from wells completed in the alluvial aquifer next to the Arkansas River and in samples from the upper Claiborne aquifer, which underlies the alluvial aquifer, indicate that leakage from the river and upward flow of saline water in underlying aquifers are not likely sources for the saline water in the alluvial aquifer in Area I. A good comparison was noted for chloride concentrations in Area I and surface geomorphology. In the majority of cases, elevated chloride concentrations occurred in backswamp deposits, with low concentrations (less than 50 milligrams per liter) in areas of active or abandoned channel deposits. The fine-grained, clay-rich deposits associated with backswamp areas likely restrict recharge, induce increased ratios between evapotranspiration and recharge, and experience minimal flushing of salts concentrated during evapotranspiration. In Area II, chloride isoconcentration maps of the underlying upper Claiborne aquifer, in addition to samples from wells completed in the middle and lower Claiborne aquifers, showed a similar chloride distribution to that of the alluvial aquifer with decreasing chloride concentrations to the east of the zone of elevated chloride concentrations, which suggests a deeper source of saline water that affects Tertiary and Quaternary aquifer systems. Mixing curves developed from bromide/chloride ratios in water samples from the alluvial aquifer, Tertiary aquifers, and samples of brine water from the Jurrasic Smackover Formation additionally discounted upward flow of saline water from underlying Tertiary formations as a potential mechanism for salinity in the alluvial aquifer in Area II. A review of information on oil exploration wells in Chicot County revealed that most of these wells were drilled from 1960 to 1980, after the elevated chloride concentrations were detected in the early 1950s. The elongated nature of the zone of elevated chloride concentrations in Area II suggests a line source or linear conduit connection with the source. Maps of a fractured limestone in the Smackover Formation in Arkansas, Mississippi, and Louisiana for purpose

Scientific Investigations Report

Wyoming Landscape Conservation Initiative—A case study in partnership development

The Wyoming Landscape Conservation Initiative (WLCI) is a successful example of collaboration between science and natural resource management at the landscape scale. In southwestern Wyoming, expanding energy and mineral development, urban growth, and other changes in land use over recent decades, combined with landscape-scale drivers such as climate change and invasive species, have presented compelling challenges to resource managers and a diverse group of Federal, State, industry, and non-governmental organizations, as well as citizen stakeholders. To address these challenges, the WLCI was established as a collaborative forum and interagency partnership to develop and implement science-based conservation actions. About a decade after being established, this report documents the establishment and history of the WLCI, focusing on the path to success of the initiative and providing insights and details that may be useful in developing similar partnerships in other locations. Not merely retrospective, the elements of the WLCI that are presented herein are still in play, still evolving, and still contributing to the resolution of compelling conservation challenges in the Western United States. The U.S. Geological Survey has developed many successful longstanding partnerships, of which the WLCI is one example. “As the Nation’s largest water, earth, and biological science and civilian mapping agency, the U.S. Geological Survey collects, monitors, analyzes, and provides scientific understanding about natural resource conditions, issues, and problems. The diversity of our scientific expertise enables us to carry out large-scale, multi-disciplinary investigations and provide impartial scientific information to resource managers, planners, and other customers” (U.S. Geological Survey, 2016).

Wyoming

Development of regional skews for selected flood durations for the Central Valley Region, California, based on data through water year 2008

Flood-frequency information is important in the Central Valley region of California because of the high risk of catastrophic flooding. Most traditional flood-frequency studies focus on peak flows, but for the assessment of the adequacy of reservoirs, levees, other flood control structures, sustained flood flow (flood duration) frequency data are needed. This study focuses on rainfall or rain-on-snow floods, rather than the annual maximum, because rain events produce the largest floods in the region. A key to estimating flood-duration frequency is determining the regional skew for such data. Of the 50 sites used in this study to determine regional skew, 28 sites were considered to have little to no significant regulated flows, and for the 22 sites considered significantly regulated, unregulated daily flow data were synthesized by using reservoir storage changes and diversion records. The unregulated, annual maximum rainfall flood flows for selected durations (1-day, 3-day, 7-day, 15-day, and 30-day) for all 50 sites were furnished by the U.S. Army Corps of Engineers. Station skew was determined by using the expected moments algorithm program for fitting the Pearson Type 3 flood-frequency distribution to the logarithms of annual flood-duration data. Bayesian generalized least squares regression procedures used in earlier studies were modified to address problems caused by large cross correlations among concurrent rainfall floods in California and to address the extensive censoring of low outliers at some sites, by using the new expected moments algorithm for fitting the LP3 distribution to rainfall flood-duration data. To properly account for these problems and to develop suitable regional-skew regression models and regression diagnostics, a combination of ordinary least squares, weighted least squares, and Bayesian generalized least squares regressions were adopted. This new methodology determined that a nonlinear model relating regional skew to mean basin elevation was the best model for each flood duration. The regional-skew values ranged from -0.74 for a flood duration of 1-day and a mean basin elevation less than 2,500 feet to values near 0 for a flood duration of 7-days and a mean basin elevation greater than 4,500 feet. This relation between skew and elevation reflects the interaction of snow and rain, which increases with increased elevation. The regional skews are more accurate, and the mean squared errors are less than in the Interagency Advisory Committee on Water Data's National skew map of Bulletin 17B.

California

Mean annual, seasonal, and monthly precipitation and runoff in Arkansas, 1951-2011

This report describes long-term annual, seasonal, and monthly means for precipitation and runoff in Arkansas for the period from 1951 through 2011. Precipitation means were estimated using data from the Parameter-elevation Regressions on Independent Slopes Model database; while total runoff, groundwater runoff, and surface runoff means were estimated using data from 123 active and inactive U.S. Geological Survey continuous-record streamflow-gaging stations located in Arkansas and surrounding States. Annual precipitation in Arkansas for the period from 1951 through 2011 had a mean of 49.8 inches. Of the six physiographic sections in Arkansas, the Ouachita Mountains had the largest mean annual precipitation at 53.0 inches, while the Springfield-Salem plateaus had the smallest mean annual precipitation at 45.5 inches. The mean annual total runoff for Arkansas was 17.8 inches. The Ouachita Mountains had the largest mean annual total runoff at 20.4 inches, while the Springfield-Salem plateaus had the smallest mean annual total runoff at 15.0 inches. Runoff is diminished during the dry season, which is attributed to increased losses from evapotranspiration, consumptive uses including irrigation, and increased withdrawals for public and private water supplies. The decline in runoff during the dry season is observed across the State in all physiographic sections. Spatial results for precipitation and runoff are presented in a series of maps that are available for download from the publication Web page in georeferenced raster formats.

Arkansas

Two-dimensional hydraulic analyses of Joachim Creek, De Soto, Missouri

A two-dimensional hydraulic model; water-surface profiles; and digital maps of water-surface elevation, velocities, and water depths were developed for a 6.7-mile reach of Joachim Creek within and near the city of De Soto, Missouri. Water-surface profiles were generated for the 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probability (10-, 25-, 50-, 100-, and 500-year recurrence interval) flows. Digital maps of water-surface elevation, water depth, and velocity were generated for the 1- and 0.2-percent annual exceedance probability flows. Water-surface elevations and inundation extents of generated profiles and maps were substantially lower than similar products produced for the 2019 flood-insurance study that included the study reach. The differences in water-surface elevations can be attributed to differences in input streamflows and hydraulic simulation techniques. The water-surface elevations generated for the 1- and 0.2-percent annual exceedance probability flows were used to assess the vulnerability and inundation depths of 231 selected structures within the city of De Soto. Results indicate that 157 to 177 of the 231 structures were affected at the 1-percent annual exceedance probability flow, depending on the adjacent grade elevation used for reference. Between 185 and 198 structures were affected at the 0.2-percent annual exceedance probability flow, depending on grade elevation. Inundation depths at the affected structures were 0.02 to 9.28 feet (ft), depending on the flow and adjacent grade reference. Flood elevations were computed for Joachim Creek using a two-dimensional, finite-volume numerical modeling application for river hydraulics. The hydraulic model was calibrated using high-water marks from the April 18, 2013, flood and the maximum measured streamflow at the U.S. Geological Survey streamgage Joachim Creek at De Soto, Mo. (station 07019500), on September 8, 2018. The calibrated model was then used to compute the hydraulic conditions associated with the 10-, 4-, 2-, 1-, and 0.2-percent annual exceedance probability flows. The simulated water-surface elevations and digital elevation model (derived from light detection and ranging data having a 0.60-ft vertical accuracy and a 1.97-ft horizontal resolution) were used to generate products including water-surface profiles and maps of inundated area, water depth, and velocities using model postprocessing software.

Missouri

Bathymetric and underwater video survey of Lower Granite Reservoir and vicinity, Washington and Idaho, 2009-10

The U.S. Geological Survey conducted a bathymetric survey of the Lower Granite Reservoir, Washington, using a multibeam echosounder, and an underwater video mapping survey during autumn 2009 and winter 2010. The surveys were conducted as part of the U.S. Army Corps of Engineer's study on sediment deposition and control in the reservoir. The multibeam echosounder survey was performed in 1-mile increments between river mile (RM) 130 and 142 on the Snake River, and between RM 0 and 2 on the Clearwater River. The result of the survey is a digital elevation dataset in ASCII coordinate positioning data (easting, northing, and elevation) useful in rendering a 3×3-foot point grid showing bed elevation and reservoir geomorphology. The underwater video mapping survey was conducted from RM 107.73 to 141.78 on the Snake River and RM 0 to 1.66 on the Clearwater River, along 61 U.S. Army Corps of Engineers established cross sections, and dredge material deposit transects. More than 900 videos and 90 bank photographs were used to characterize the sediment facies and ground-truth the multibeam echosounder data. Combined, the surveys were used to create a surficial sediment facies map that displays type of substrate, level of embeddedness, and presence of silt.

Washington;Idaho

Flood-inundation maps and updated components for a flood-warning system or the City of Marietta, Ohio and selected communities along the Lower Muskingum River and Ohio River

Digital flood-inundation maps for lower reaches of the Muskingum River and a reach of the Ohio River in southeast Ohio were created by the U.S. Geological Survey (USGS), in cooperation with the Muskingum Watershed Conservancy District and the City of Marietta, Ohio. To complete the inundation maps, Ohio River and lower Muskingum River bathymetry was updated and two streamgages, one on the Ohio River upstream of Marietta near Sardis, Ohio, and one on the Muskingum River in Beverly, Ohio, were added as basic components of the flood-warning system. An updated hydraulic model component also led to the new flood-inundation maps. The maps, which can be accessed through the USGS Flood Inundation Mapping Science Web site at http://water.usgs.gov/osw/flood_inundation/ depict estimates of the areal extent of flooding corresponding to water levels (stages) at one or more of the following USGS streamgages: Muskingum River at McConnelsville, Ohio (03150000); Muskingum River at Beverly, Ohio (03150500); and Ohio River at Marietta, Ohio (03150700). The maps can be used in conjunction with National Weather Service flood-forecast data to show areas of estimated flood inundation associated with forecasted flood-peak stages. Flood profiles for selected reaches were prepared by calibrating steady-state step-backwater models to selected streamgage rating curves. The step-backwater models were used to determine water-surface-elevation profiles for up to 12 flood stages at a streamgage with corresponding stream-flows ranging from approximately the 10- to 0.2-percent chance annual-exceedance probabilities for each of the 3 streamgages that correspond to the flood-inundation maps. Additional hydraulic modeling was used to account for the effects of backwater from the Ohio River on water levels in the Muskingum River. The computed longitudinal profiles of flood levels were used with a Geographic Information System digital elevation model (derived from light detection and ranging) to delineate flood-inundation areas. Digital maps showing flood-inundation areas overlain on digital orthophotographs were prepared for the selected floods.

Ohio

Methods for estimating regional skewness of annual peak flows in parts of the Great Lakes and Ohio River Basins, based on data through water year 2013

Bulletin 17C (B17C) recommends fitting the log-Pearson Type III (LP−III) distribution to a series of annual peak flows at a streamgage by using the method of moments. The third moment, the skewness coefficient (or skew), is important because the magnitudes of annual exceedance probability (AEP) flows estimated by using the LP−III distribution are affected by the skew; interest is focused on the right-hand tail of the distribution, which represents the larger annual peak flows that correspond to small AEPs. For streamgages having modest record lengths, the skew is sensitive to extreme events like large floods, which cause a sample to be highly asymmetrical or “skewed.” For this reason, B17C recommends using a weighted-average skew computed from the station skew for a given streamgage and a regional skew. This report generates an estimate of regional skew for a study area encompassing most of the Great Lakes Basin (hydrologic unit 04) and part of the Ohio River Basin (hydrologic unit 05). A total of 551 candidate streamgages that were unaffected by extensive regulation, diversion, urbanization, or channelization were considered for use in the skew analysis; after screening for redundancy and pseudo record length greater than 36 years, 368 streamgages were selected for use in the study. Flood frequencies for candidate streamgages were analyzed by employing the Expected Moments Algorithm, which extends the method of moments so that it can accommodate interval, censored, and historic/paleo flow data, as well as the Multiple Grubbs-Beck test to identify potentially influential low floods in the data series. Bayesian weighted least squares/Bayesian generalized least squares regression was used to develop a regional skew model for the study area that would incorporate possible variables (basin characteristics) to explain the variation in skew in the study area. Twelve basin characteristics were considered as possible explanatory variables; however, none produced a pseudo coefficient of determination greater than 5 percent; as a result, these characteristics did not help to explain the variation in skew in the study area. Therefore, a constant model having a regional skew coefficient of 0.086 and an average variance of prediction ( AVP new ) (which corresponds to the mean square error [MSE]) of 0.13 at a new streamgage was selected. The AVP new corresponds to an effective record length of 54 years, a marked improvement over the Bulletin 17B national skew map, whose reported MSE of 0.302 indicated a corresponding effective record length of only 17 years.

Great Lakes and Ohio River Basins

Methods for estimating regional skewness of annual peak flows in parts of eastern New York and Pennsylvania, based on data through water year 2013

Bulletin 17C (B17C) recommends fitting the log-Pearson Type III (LP−III) distribution to a series of annual peak flows at a streamgage by using the method of moments. The third moment, the skewness coefficient (or skew), is important because the magnitudes of annual exceedance probability (AEP) flows estimated by using the LP–III distribution are affected by the skew; interest is focused on the right-hand tail of the distribution, which represents the larger annual peak flows that correspond to small AEPs. For streamgages having modest record lengths, the skew is sensitive to extreme events like large floods, which cause a sample to be highly asymmetrical or “skewed.” For this reason, B17C recommends using a weighted-average skew computed from the skew of the annual peak flows for a given streamgage and a regional skew. This report presents an estimate of regional skew for a study area encompassing parts of eastern New York and Pennsylvania. A total of 232 candidate U.S. Geological Survey streamgages that were unaffected by extensive regulation, diversion, urbanization, or channelization were considered for use in the skew analysis; after screening for redundancy and pseudo record length ( P RL ) of at least 36 years, 183 streamgages were selected for use in the study. Flood frequencies for candidate streamgages were analyzed by employing the expected moments algorithm, which extends the method of moments so that it can accommodate interval, censored, and historical/paleo flow data, as well as the multiple Grubbs-Beck test to identify potentially influential low floods in the data series. Bayesian weighted least squares/Bayesian generalized least squares regression was used to develop a regional skew model for the study area that would incorporate possible variables (basin characteristics) to explain the variation in skew in the study area. Ten basin characteristics were considered as possible explanatory variables; however, none produced a pseudo coefficient of determination greater than 1 percent; as a result, these characteristics did not help to explain the variation in skew in the study area. Therefore, a constant model that had a regional skew coefficient of 0.32 and an average variance of prediction at a new streamgage ( AVP new , which corresponds to the mean square error [MSE] of 0.11) was selected. The AVP new corresponds to an effective record length of 68 years, a marked improvement over the Bulletin 17B national skew map, whose reported MSE of 0.302 indicated a corresponding effective record length of only 17 years.

New York, Pennsylvania

Processing, Analysis, and General Evaluation of Well-Driller Logs for Estimating Hydrogeologic Parameters of the Glacial Sediments in a Ground-Water Flow Model of the Lake Michigan Basin

In 2005, the U.S. Geological Survey began a pilot study for the National Assessment of Water Availability and Use Program to assess the availability of water and water use in the Great Lakes Basin. Part of the study involves constructing a ground-water flow model for the Lake Michigan part of the Basin. Most ground-water flow occurs in the glacial sediments above the bedrock formations; therefore, adequate representation by the model of the horizontal and vertical hydraulic conductivity of the glacial sediments is important to the accuracy of model simulations. This work processed and analyzed well records to provide the hydrogeologic parameters of horizontal and vertical hydraulic conductivity and ground-water levels for the model layers used to simulated ground-water flow in the glacial sediments. The methods used to convert (1) lithology descriptions into assumed values of horizontal and vertical hydraulic conductivity for entire model layers, (2) aquifer-test data into point values of horizontal hydraulic conductivity, and (3) static water levels into water-level calibration data are presented. A large data set of about 458,000 well driller well logs for monitoring, observation, and water wells was available from three statewide electronic data bases to characterize hydrogeologic parameters. More than 1.8 million records of lithology from the well logs were used to create a lithologic-based representation of horizontal and vertical hydraulic conductivity of the glacial sediments. Specific-capacity data from about 292,000 well logs were converted into horizontal hydraulic conductivity values to determine specific values of horizontal hydraulic conductivity and its aerial variation. About 396,000 well logs contained data on ground-water levels that were assembled into a water-level calibration data set. A lithology-based distribution of hydraulic conductivity was created by use of a computer program to convert well-log lithology descriptions into aquifer or nonaquifer categories and to calculate equivalent horizontal and vertical hydraulic conductivities (K and KZ, respectively) for each of the glacial layers of the model. The K was based on an assumed value of 100 ft/d (feet per day) for aquifer materials and 1 ft/d for nonaquifer materials, whereas the equivalent KZ was based on an assumed value of 10 ft/d for aquifer materials and 0.001 ft/d for nonaquifer materials. These values were assumed for convenience to determine a relative contrast between aquifer and nonaquifer materials. The point values of K and KZ from wells that penetrate at least 50 percent of a model layer were interpolated into a grid of values. The K distribution was based on an inverse distance weighting equation that used an exponent of 2. The KZ distribution used inverse distance weighting with an exponent of 4 to represent the abrupt change in KZ that commonly occurs between aquifer and nonaquifer materials. The values of equivalent hydraulic conductivity for aquifer sediments needed to be adjusted to actual values in the study area for the ground-water flow modeling. The specific-capacity data (discharge, drawdown, and time data) from the well logs were input to a modified version of the Theis equation to calculate specific capacity based horizontal hydraulic conductivity values (KSC). The KSC values were used as a guide for adjusting the assumed value of 100 ft/d for aquifer deposits to actual values used in the model. Water levels from well logs were processed to improve reliability of water levels for comparison to simulated water levels in a model layer during model calibration. Water levels were interpolated by kriging to determine a composite water-level surface. The difference between the kriged surface and individual water levels was used to identify outlier water levels. Examination of the well-log lithology data in map form revealed that the data were not only useful for model input, but also were useful for understanding th

Scientific Investigations Report

Modified method for estimating petroleum source-rock potential using wireline logs, with application to the Kingak Shale, Alaska North Slope

In 2012, the U.S. Geological Survey completed an assessment of undiscovered, technically recoverable oil and gas resources in three source rocks of the Alaska North Slope, including the lower part of the Jurassic to Lower Cretaceous Kingak Shale. In order to identify organic shale potential in the absence of a robust geochemical dataset from the lower Kingak Shale, we introduce two quantitative parameters, $\Delta DT_\bar{x}$ and $\Delta DT_z$, estimated from wireline logs from exploration wells and based in part on the commonly used delta-log resistivity ($\Delta \text{ }log\text{ }R$) technique. Calculation of $\Delta DT_\bar{x}$ and $\Delta DT_z$ is intended to produce objective parameters that may be proportional to the quality and volume, respectively, of potential source rocks penetrated by a well and can be used as mapping parameters to convey the spatial distribution of source-rock potential. Both the $\Delta DT_\bar{x}$ and $\Delta DT_z$ mapping parameters show increased source-rock potential from north to south across the North Slope, with the largest values at the toe of clinoforms in the lower Kingak Shale. Because thermal maturity is not considered in the calculation of $\Delta DT_\bar{x}$ or $\Delta DT_z$, total organic carbon values for individual wells cannot be calculated on the basis of $\Delta DT_\bar{x}$ or $\Delta DT_z$ alone. Therefore, the $\Delta DT_\bar{x}$ and $\Delta DT_z$ mapping parameters should be viewed as first-step reconnaissance tools for identifying source-rock potential.

Alaska

Precipitation-driven flood-inundation mapping of the Little Blue River at Grandview, Missouri

The U.S. Geological Survey (USGS), in cooperation with the City of Grandview, Missouri, assessed flooding of the Little Blue River at Grandview resulting from varying precipitation magnitudes and durations and expected land-cover changes. The precipitation scenarios were used to develop a library of flood-inundation maps that included a 3.5-mile reach of the Little Blue River and tributaries within and adjacent to the city. A hydrologic model of the upper Little Blue River Basin and a hydraulic model of a selected study reach of the Little Blue River and tributaries were constructed to assess streamflow magnitudes associated with simulated precipitation amounts and the resulting flood-inundation conditions. The U.S. Army Corps of Engineers Hydrologic Engineering Center-Hydrologic Modeling System (HEC–HMS; version 4.4.1) was used to simulate the amount of streamflow produced from a range of rain events. The Hydrologic Engineering Center-River Analysis System (HEC–RAS; version 5.0.7) was then used to construct a steady-state hydraulic model to map resulting areas of flood inundation. Both models were calibrated to the May 28, 2020, high-flow event that produced a peak streamflow approximating a 10-percent annual exceedance probability (10-year flood-frequency recurrence interval) at the Little Blue River at Grandview streamgage (USGS station 06893750). The calibrated HEC–HMS model was used to simulate streamflows from design rainfall events of 1- to 8-hour durations and ranging from a 100- to 0.2-percent annual exceedance probability. Flood-inundation maps were produced for USGS streamflow stages of 17.0 feet (ft), or near bankfull, to 23.0 ft, or a stage exceeding the 0.2-percent annual exceedance interval flood, using the HEC–RAS model. The consequence of each precipitation duration-frequency value was represented by a 1-ft increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding stage at the reference USGS streamgage. Four scenarios were developed with the HEC–HMS hydrologic model: (1) current (2016) land cover, normal antecedent soil-moisture conditions; (2) current land cover, wet antecedent soil-moisture conditions; (3) future land cover, normal antecedent soil-moisture conditions; and (4) future land cover, wet antecedent soil-moisture conditions. The future land-cover condition was estimated based on anticipated development in the basin. All precipitation scenarios were input into each of the four land-cover antecedent moisture conditions and then assigned to a resulting flood-inundation map based on the generated peak flow and corresponding stage at the reference streamgage.

Missouri

Precipitation-driven flood-inundation mapping of Muddy Creek at Harrisonville, Missouri

The U.S. Geological Survey, in cooperation with the city of Harrisonville, Missouri, assessed flooding of Muddy Creek resulting from varying precipitation magnitudes and durations, antecedent runoff conditions, and channel modifications (cleaned culverts and added detention storage). The precipitation scenarios were used to develop a library of flood-inundation maps that included a 3.8-mile reach of Muddy Creek and tributaries within and adjacent to the city. Hydrologic and hydraulic models of the upper Muddy Creek study basin were used to assess streamflow magnitudes associated with simulated precipitation amounts and the resulting flood-inundation conditions. The U.S. Army Corps of Engineers Hydrologic Engineering Center-Hydrologic Modeling System (HEC–HMS; version 4.4.1) was used to simulate the amount of streamflow produced from a range of precipitation events. The Hydrologic Engineering Center-River Analysis System (HEC–RAS; version 5.0.7) was then used to route streamflows and map resulting areas of flood inundation. The hydrologic and hydraulic models were calibrated to the September 28, 2019; May 27, 2021; and June 25, 2021, runoff events representing a range of antecedent runoff conditions and hydrologic responses. The calibrated HEC–HMS model was used to simulate streamflows from design rainfall events of 30-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability. Flood-inundation maps were produced for reference stages of 1.0 foot (ft), or near bankfull, to 4.0 ft, or a stage exceeding the 0.1-percent annual exceedance probability interval precipitation, using the HEC–RAS model. The results of each precipitation duration-frequency value were represented by a 0.5-ft increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding stage at the Muddy Creek reference location. Seven scenarios were developed with the HEC–HMS hydrologic model with resulting streamflows routed in a HEC–RAS hydraulic model, and these scenarios varied by antecedent runoff condition and potential channel modifications. The same precipitation scenarios were used in each of the seven antecedent runoff and channel conditions, and the simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding stage at the Muddy Creek reference location.

Missouri

Influence of geologic setting on ground-water availability in the Lawrenceville area, Gwinnett County, Georgia

Obtaining large quantities of ground water needed for municipal and industrial supply in the Piedmont and Blue Ridge physiographic provinces can be challenging because of the complex geology and the typically low primary permeability of igneous and metamorphic rocks. Areas of enhanced secondary permeability in the bedrock do occur, however, and 'high-yield' wells are not uncommon, particularly where careful site-selection techniques are used prior to test drilling. The U.S. Geological Survey - in cooperation with the City of Lawrenceville, Georgia - conducted this study from 2000 to 2002 to learn more about how different geologic settings influence the availability of ground water in igneous and metamorphic bedrock with the expectation that this knowledge could be used to help identify additional water resources in the area. In compositionally layered-rock settings, wells derive water almost exclusively from lithologically and structurally controlled water-bearing zones formed parallel to foliation and compositional layering. These high-permeability, water-bearing zones - termed foliation-parallel parting systems -combined with high-angle joint systems, are the primary control for the high-yield wells drilled in the Lawrenceville area; yields range from 100 to several hundred gallons per minute (gal/min). Near Lawrenceville, areas with high ground-water yield are present in sequences of amphibolite, biotite gneiss, and button schist where the structural attitude of the rocks is gently dipping, in areas characterized by abundant jointing, and in topographic settings with a continuous source of recharge along these structures. In massive-rock settings, wells derive water mostly from joint systems, although foliation-parallel parting systems also may be important. Wells deriving water primarily from steeply-dipping joint systems typically have low yields ranging from 1 to 5 gal/min. Joint systems in massive-rock settings can be identified and characterized by using many of the methods described in this report. Geologic mapping was the primary method used to determine the distribution, variability, and relative concentrations (intensity) of joint systems. In the subsurface, joints were characterized by taking orientation measurements in the open boreholes of wells using acoustic and/or optical televiewers. In this investigation, the only practical approach found for locating areas of high ground-water potential was first through detailed geologic mapping followed by test drilling, borehole geophysical logging, and aquifer testing. Geologic methods help characterize both large- and small-scale structures and other lithologic and stratigraphic features that influence development of increased secondary permeability in the bedrock. The rock types, discontinuities, depth of weathering, topographic position, and recharge potential - which were the principal factors assessed through detailed geologic mapping - must be evaluated carefully, in relation to one another, to assess the ground-water potential in a given area.

Georgia