Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Engineering Analysis”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56Linked to original sources

Pressure sensitivity of low permeability sandstones

Detailed core analysis has been carried out on 32 tight sandstones with permeabilities ranging over four orders of magnitude (0.0002 to 4.8 mD at 5000 psi confining pressure). Relationships between gas permeability and net confining pressure were measured for cycles of loading and unloading. For some samples, permeabilities were measured both along and across bedding planes. Large variations in stress sensitivity of permeability were observed from one sample to another. The ratio of permeability at a nominal confining pressure of 500 psi to that at 5000 psi was used to define a stress sensitivity ratio. For a given sample, confining pressure vs permeability followed a linear log-log relationship, the slope of which provided an index of pressure sensitivity. This index, as obtained for first unloading data, was used in testing relationships between stress sensitivity and other measured rock properties. Pressure sensitivity tended to increase with increase in carbonate content and depth, and with decrease in porosity, permeability and sodium feldspar. However, scatter in these relationships increased as permeability decreased. Tests for correlations between pressure sensitivity and various linear combinations of variables are reported. Details of pore structure related to diagenetic changes appears to be of much greater significance to pressure sensitivity than mineral composition.

Journal of Petroleum Science and Engineering↗

Physical and hydrochemical evidence of lake leakage near Jim Woodruff Lock and Dam and of ground-water inflow to Lake Seminole, and an assessment of karst features in and near the lake, southwestern Georgia and northwestern Florida

Hydrogeologic data and water-chemistry analyses indicate that Lake Seminole leaks into the Upper Floridan aquifer near Jim Woodruff Lock and Dam, southwestern Georgia and northwestern Florida, and that ground water enters Lake Seminole along upstream reaches of the lake’s four impoundment arms (Chattahoochee and Flint Rivers, Spring Creek, and Fishpond Drain). Written accounts by U.S. Army Corps of Engineers geologists during dam construction in the late 1940s and early 1950s, and construction-era photographs, document karst-solution features in the limestone that comprise the lake bottom and foundation rock to the dam, and confirm the hydraulic connection of the lake and aquifer. More than 250 karst features having the potential to connect the lake and aquifer were identified from preimpoundment aerial photographs taken during construction. An interactive map containing a photomosaic of 53 photographic negatives was orthorectfied to digital images of 1:24,000-scale topographic maps to aid in identifying karst features that function or have the potential to function as locations of water exchange between Lake Seminole and the Upper Floridan aquifer. Some identified karst features coincide with locations of mapped springs, spring runs, and depressions that are consistent with sinkholes and sinkhole ponds. Hydrographic surveys using a multibeam echosounder (sonar) with sidescan sonar identified sinkholes in the lake bottom along the western lakeshore and in front of the dam. Dye-tracing experiments indicate that lake water enters these sinkholes and is transported through the Upper Floridan aquifer around the west side of the dam at velocities of about 500 feet per hour to locations where water "boils up" on land (at Polk Lake Spring) and in the channel bottom of the Apalachicola River (at the "River Boil"). Water discharging from Polk Lake Spring joins flow from a spring-fed ground-water discharge zone located downstream of the dam; the combined flow disappears into a sinkhole located on the western floodplain of the river and is transmitted through the Upper Floridan aquifer, eventually discharging to the Apalachicola River at the River Boil. Acoustic Doppler current profiling yielded flow estimates from the River Boil in the range from about 140 to 220 cubic feet per second, which represents from about 1 to 3 percent of the average daily flow in the river. Binary mixing-model analysis using naturally occurring isotopes of oxygen and hydrogen (oxygen-18 and deuterium) indicates that discharge from the River Boil consists of a 13-to-1 ratio of lake water to ground water and that other sources of lake leakage and discharge to the boil probably exist. Analyses of major ions, nutrients, radon-222, and stable isotopes of hydrogen and oxygen contained in water samples collected from 29 wells, 7 lake locations, and 5 springs in the Lake Seminole area during 2000 indicate distinct chemical signatures for ground water and surface water. Ground-water samples contained higher concentrations of calcium and magnesium, and higher alkalinity and specific conductance than surface-water samples, which contained relatively high concentrations of total organic carbon and sulfate. Solute and isotopic tracers indicate that, from May to October 2000, springflow exhibited more ground-water qualities (high specific conductance, low dissolved oxygen, and low temperature) than surface water; however, the ratio of ground water to surface water of the springs was difficult to quantify from November to April because of reduced springflow and rapid mixing of springflow and lake water during sampling. The saturation index of calcite in surface-water samples indicates that while surface water is predominately undersaturated with regard to calcite year-round, a higher potential for dissolution of the limestone matrix exists from late fall through early spring than during summer. The relatively short residence time (5–7 hours) and rapid flow velocity (nearly 500 feet per hour) of lake water leaking into the Upper Floridan aquifer and exiting at the River Boil in the Apalachicola River implies that calcite-undersaturated water is in constant contact with the limestone, increasing the potential for limestone dissolution and enlargement of flow pathways by erosion. A relatively low potential exists, however, for limestone dissolution to cause sudden sinkhole collapse followed by catastrophic lake drainage because ground-water levels close to the lake, except near the dam, are nearly the same as lake stage, resulting in low vertical and lateral hydraulic gradients and low flow between the lake and aquifer. An increased potential for lake leakage and sinkhole formation and collapse exists near some in-lake springs during colder months of the year, as density differences and the hydraulic potential between lake water and ground water establish the conditions for calcite-undersaturated lake water to enter nonflowing springs and contact limestone.

Florida, Georgia↗

Maximum entropy derived statistics of sound speed structure in a fine-grained sediment inferred from sparse broadband acoustic measurements on the New England continental shelf

Marginal probability distributions for parameters representing an effective sound-speed structure of a fine-grained sediment are inferred from a data ensemble maximum entropy method that utilizes a sparse spatially distributed set of received pressure time series resulting from multiple explosive sources in a shallow-water ocean environment possessing significant spatial variability of the seabed. A remote sensing seabed acoustics experiment undertaken in March 2017 off the New England Shelf was designed so that multiple independent analyses could infer the statistical properties of the seabed. The current analysis incorporates the measured horizontal variability from interpretations of a subbottom profiling survey of the experimental area. An idealized range- and azimuth-dependent parameterization of the seabed is derived from identification of horizons within the seabed that define multiple sediment layers. A sparse set of explosive charges were deployed on circular tracks with radii of about 2, 4, and 6.5 km with an acoustic array at the center to correlate a set of random measurements to physical acoustic processes that characterize the seabed. The mean values of a surface sound speed ratio and a linear sound speed gradient for the fine-grained sediment layer derived from 12 data samples processed in the 25–275-Hz band provide an estimate of the effective sound-speed structure in a 130-km $^2$ area. The inferred sediment sound speed values are evaluated by predicting measured time series data not used in the statistical inference, and are also compared to historical measurements. Finally, the low-frequency maximum entropy estimate of the sediment sound speed along with physical measurements derived from piston core measurements are utilized to estimate the sediment grain bulk modulus. predictions made by the viscous grain shearing model.

Maine, Vermont, New Hampshire, Massachusetts, Rhod↗

A comparison of synthetic flowpaths derived from light detection and ranging topobathymetric data and National Hydrography Dataset High Resolution Flowlines

Bathymetric and topobathymetric light detection and ranging (lidar) digital elevation models created for the Delaware River were provided to the National Geospatial Program and used to evaluate synthetic flowpath extraction from bathymetric/topobathymetric lidar survey data as a data source for improving the density, distribution, and connectivity of the National Hydrography Dataset High Resolution Flowline Network. As the surface-water component of The National Map, the National Hydrography Dataset maintains the Nation’s drainage network flow information and geometries for surface-water features used in hydrologic, hydraulic, and other science and engineering disciplines. The regional lidar survey for the Delaware River between Hancock, New York, and Trenton, New Jersey, was collected for the U.S. Geological Survey using the Experimental Advanced Airborne Research Lidar sensor system and processed by the Coastal National Elevation Database Applications Program. Using 1 percent of the maximum flow accumulation value for the surveyed Delaware River corridor as the flow accumulation threshold for grid cells at 1-, 5-, and 10-meter resolution created 223 to 283 kilometers of synthetic flowpaths potentially representing the river channel thalweg, which is the deepest point in a riverbed cross-section. There was potential for improving the High Resolution National Hydrography Dataset (HR NHD) Flowline network in places where the Delaware River channel, depicted as an Artificial Path in the HR NHD, is offset from the extracted synthetic river flowpath which sometimes appeared better positioned than the Artificial Path to represent the river thalweg. For the same area, using 0.05 percent of the maximum flow accumulation at the 1-, 5-, and 10-meter resolutions extracted 744 to 1,317 kilometers of synthetic flowpaths, with extracted synthetic flowpaths representing the main river channel and additional synthetic flowpaths representing tributaries or streams adjacent to the main channel. Overlaying these results with the HR NHDFlowline Network indicates that some of the additional synthetic flowpaths are connected to or extend HR NHD stream/river feature types. Some disconnected or isolated synthetic flowpaths not included in stream/river feature types were validated in orthoimagery and U.S. Topo Maps and provide examples of how extracted synthetic flowpaths could be used to delineate new stream/river features. Other additional extracted synthetic flowpaths depict linear features such as canals, tree lines, roads, or linear topographic depressions. For some river reaches where obstructions to flow or where low-relief topographic or bathymetric surfaces alter the flow direction, the software tool used to develop the flow direction grid did not calculate a primary flowpath for the river channel. Based on the results of this analysis, site conditions for the Delaware River corridor did not affect the quality of lidar bathymetric survey data. However, depending on the resolution of the lidar bathymetric digital elevation models (BDEMs), site conditions do have different effects on results for extracted synthetic flowpaths. We found that synthetic flowpaths extracted from 1-meter resolution lidar DEMs had more varied flow directions around in-channel landforms that obstructed flow than synthetic flowpaths extracted from 5- or 10-meter resolution lidar DEMs. As a result the 1-meter resolution DEM created some isolated or discontinuous synthetic flowpath segments where the 5- and 10-meter DEMs developed more continuous flowpaths. In this case the river bed upstream from the in-channel obstruction is shallower than the river bed downstream. Under these conditions the 1-meter resolution DEM provided synthetic flowpaths delineating a potential river thalweg. In this same area, the software solution modified (virtually raised) the river bed in the 5- and 10-meter resolution DEMs and flattened the bathymetric surface to create a continuous downstream flow direction, which caused trellis-patterned synthetic flowpaths to form. Under different site conditions and converse to the above development of synthetic flowpaths at different resolutions, at an abandoned river flood plain (terrace) with low relief that is adjacent to the river channel, the flow direction grid for the 1-meter resolution DEM developed continuous synthetic flowpath corresponding to a HR NHD Flowline network stream/river feature that connected to the main river channel but the larger resolution DEMs created isolated or disconnected synthetic flowpaths. A project to continue an evaluation of benefits of or issues caused by extracting synthetic flowpaths to enhance the HR NHD could include a study to assess the potential for merging surface-water flowpaths extracted from lidar topobathymetry and 3D Elevation Program digital elevation models. The merged DEM approach to synthetic flowpath extraction could extend the HR NHDFlowline network and enhance flow accumulations that might develop better flow direction grids in low-relief areas. Because of the confined lateral extent of the Delaware River, the lidar DEMs were not used to create catchments or watersheds; however, the merged DEM approach could also be tested as a resource for enhancing HR NHD catchments and watersheds. This lidar DEM synthetic flowpath extraction project supports the National Geospatial Program efforts to collect and produce high-quality lidar data to provide 3-dimensional representations of natural feature and aligns with the National Spatial Data Infrastructure to improve utilization of geospatial data. The results also can be useful for understanding strategies that can help maintain quality data in the HR NHD programs. KEYWORDS: bathymetric, digital elevation model, extracted synthetic flowpath, lidar, High Resolution National Hydrography Dataset, topobathymetric

New Jersey↗

Validation of gridded precipitation datasets for flood-typing in select conterminous U.S. basins

Gridded precipitation datasets are required for flood-typing historical annual peak streamflow events in basins across the Conterminous United States. Selected gridded precipitation datasets were validated over the period 1981–2013 through comparisons with gage data from the NOAA Global Historical Climatology Network daily (GHCNd). The ability of each gridded dataset to capture the spatiotemporal characteristics of daily precipitation, including multi-day extremes over six selected regions, was assessed using the Kling-Gupta Efficiency metric and its component statistics. Overall, the Parameter-elevation Regression on Independent Slopes Model and Livneh-unsplit were found to best match the spatiotemporal variability of the GHCNd precipitation data, including extremes. The Analysis of Record for Calibration was found to be the third best-performing dataset in most regions except in the western U.S. The performance of reanalysis datasets evaluated appears to be poor compared to gage-based datasets. The reanalysis datasets might not be able to skillfully capture precipitation amounts at the correct location and time. Gage- and radar-based datasets were found to have relatively small biases (within +/-10% on an annual basis), while reanalysis datasets were found to have larger positive apparent biases, especially in winter and spring in most regions. It is possible that the apparent overestimation of winter and spring precipitation in the reanalysis datasets might reflect snow undercatch at gages especially in the central U.S. An overall deterioration of performance for correlation and/or variability was also observed for the summer season compared to other seasons in the reanalysis datasets. Various precipitation datasets might need to be used for flood-typing during different periods from the late 19th century to present. Datasets from different sources have different biases and errors and might have to be homogenized using downscaling and bias-adjustment methods. Alternatively, precipitation thresholds used in some flood-typing schemes might have to be adjusted as a function of time.

conterminous United States↗

Greenup and evapotranspiration following the Minute 319 pulse flow to Mexico: An analysis using Landsat 8 Normalized Difference Vegetation Index (NDVI) data

In the southwestern U.S., many riparian ecosystems have been altered by dams, water diversions, and other anthropogenic activities. This is particularly true of the Colorado River, where numerous dams and agricultural diversions have affected this water course, especially south of the U.S.–Mexico border. In the spring of 2014, 130 million cubic meters of water was released to the lower Colorado River Delta in Mexico. To understand the impact of this pulse flow release on vegetation in the delta’s riparian corridor, we analyzed a modified form of Landsat 8 Operational Land Imager (OLI) Normalized Difference Vegetation Index (NDVI*) data. We assessed greenup during the growing period and estimated actual evapotranspiration (ET a ) for the period prior to (yr. 2013) and following (i.e., yr. 2014 and 2015) the pulse flow. We found a significant increase in NDVI* from 2013 to 2014 (P < 0.05) and a decrease from 2014 to 2015; however, 2015 levels were still significantly higher than in 2013. ET a was also higher in 2014 vs. 2013, with an estimated 74.5 million cubic meters in 2013 and 88.9 in 2014. The most intense greening occurred in the zone of inundation but also extended into the non-flooded part of the riparian zone, indicating replenishment of groundwater. These findings suggest the peak response by vegetation to the flow lasted about one year, followed by a decrease in NDVI*. As a long term solution to the declining condition of vegetation, additional pulse releases are likely needed for restoration and survival of riparian plant communities in the Colorado River Delta.

Ecological Engineering↗

Comparison of index velocity measurements made with a horizontal acoustic Doppler current profiler

The State of Illinois' annual withdrawal from Lake Michigan is limited by a U.S. Supreme Court decree, and the U.S. Geological Survey (USGS) is responsible for monitoring flows in the Chicago Sanitary and Ship Canal (CSSC) near Lemont, Illinois as a part of the Lake Michigan Diversion Accounting overseen by the U.S. Army Corps of Engineers, Chicago District. Every 5 years, a technical review committee consisting of practicing engineers and academics is convened to review the U.S. Geological Survey's streamgage practices in the CSSC near Lemont, Illinois. The sixth technical review committee raised a number of questions concerning the flows and streamgage practices in the CSSC near Lemont and this report provides answers to many of those questions. In addition, it is the purpose of this report to examine the index velocity meters in use at Lemont and determine whether the acoustic velocity meter (AVM), which is now the primary index velocity meter, can be replaced by the horizontal acoustic Doppler current profiler (H-ADCP), which is currently the backup meter. Application of the AVM and H-ADCP to index velocity measurements in the CSSC near Lemont, Illinois, has produced good ratings to date. The site is well suited to index velocity measurements in spite of the large range of velocities and highly unsteady flows at the site. Flow variability arises from a range of sources: operation of the waterway through control structures, lockage-generated disturbances, commercial and recreational traffic, industrial withdrawals and discharges, natural inflows, seiches, and storm events. The influences of these factors on the index velocity measurements at Lemont is examined in detail in this report. Results of detailed data comparisons and flow analyses show that use of bank-mounted instrumentation such as the AVM and H-ADCP appears to be the best option for index velocity measurement in the CSSC near Lemont. Comparison of the rating curves for the AVM and H-ADCP demonstrates that the H-ADCP is a suitable replacement for the AVM as the primary index velocity meter in the CSSC near Lemont. A key component to Lake Michigan Diversion Accounting is the USGS gaging station on the CSSC near Lemont, Illinois. The importance of this gaging station in monitoring withdrawals from Lake Michigan has made it one of the most highly scrutinized gaging stations in the country. Any changes in streamgaging practices at this gaging station requires detailed analysis to ensure the change will not adversely affect the ability of the USGS to accurately monitor flows. This report provides a detailed analysis of the flow structure and index velocity measurements in the CSSC near Lemont, Illinois, to ensure that decisions regarding the future of this streamgage are made with the best possible understanding of the site and the characteristics of the flow.

Illinois↗

Lateral and vertical channel movement and potential for bed-material movement on the Madison River downstream from Earthquake Lake, Montana

The 1959 Hebgen Lake earthquake caused a massive landslide (Madison Slide) that dammed the Madison River and formed Earthquake Lake. The U.S. Army Corps of Engineers excavated a spillway through the Madison Slide to permit outflow from Earthquake Lake. In June 1970, high streamflows on the Madison River severely eroded the spillway channel and damaged the roadway embankment along U.S. Highway 287 downstream from the Madison Slide. Investigations undertaken following the 1970 flood events concluded that substantial erosion through and downstream from the spillway could be expected for streamflows greater than 3,500 cubic feet per second (ft 3 /s). Accordingly, the owners of Hebgen Dam, upstream from Earthquake Lake, have tried to manage releases from Hebgen Lake to prevent streamflows from exceeding 3,500 ft 3 /s measured at the U.S. Geological Survey (USGS) gaging station 0638800 Madison River at Kirby Ranch, near Cameron, Montana. Management of flow releases from Hebgen Lake to avoid exceeding the threshold streamflow at USGS gaging station 06038800 is difficult, and has been questioned for two reasons. First, no road damage was reported downstream from the Earthquake Lake outlet in 1993, 1996, and 1997 when streamflows exceeded the 3,500-ft 3 /s threshold. Second, the 3,500-ft 3 /s threshold generally precludes releases of higher flows that could be beneficial to the blue-ribbon trout fishery downstream in the Madison River. In response to concerns about minimizing streamflow downstream from Earthquake Lake and the possible armoring of the spillway, the USGS, in cooperation with the Madison River Fisheries Technical Advisory Committee (MADTAC; Bureau of Land Management; Montana Department of Environmental Quality; Montana Fish, Wildlife and Parks; PPL-Montana; U.S. Department of Agriculture Forest Service - Gallatin National Forest; and U.S. Fish and Wildlife Service), conducted a study to determine movement of the Madison River channel downstream from Earthquake Lake and to investigate the potential for bed material movement along the same reach. The purpose of this report is to present information about the lateral and vertical movement of the Madison River from 1970 to 2006 for a 1-mile reach downstream from Earthquake Lake and for Raynolds Pass Bridge, and to provide an analysis of the potential for bed-material movement so that MADTAC can evaluate the applicability of the previously determined threshold streamflow for initiation of damaging erosion. As part of this study channel cross sections originally surveyed by the USGS in 1971 were resurveyed in 2006. Incremental channel-movement distances were determined by comparing the stream centerlines from 14 aerial photographs taken between 1970 and 2006. Depths of channel incision and aggregation were determined by comparing the 2006 and 1971 cross-section and water-surface data. Particle sizes of bed and bank materials were measured in 2006 and 2008 using the pebble-count method and sieve analyses. A one-dimensional hydraulic-flow model (HEC-RAS) was used to calculate mean boundary-shear stresses for various streamflows; these calculated boundary-shear stresses were compared to calculated critical-shear stresses for the bed materials to determine the potential for bed-material movement. A comparison of lateral channel movement distances with annual peak streamflows shows that streamflows higher than the 3,500-ft 3 /s threshold were followed by lateral channel movement except from 1991 to 1992 and possibly from 1996 to 1997. However, it was not possible to discern whether the channel moved gradually or suddenly, or in response to one peak flow, to several peak flows, or to sustained flows. The channel moved between 2002 and 2005 even when streamflows were less than the threshold streamflow of 3,500 ft 3 /s. Comparisons of cross sections and aerial photographs show that the channel has moved laterally and incised and aggraded to varying degrees. The channel has developed meander bends and has incised as much as 5&ndash;12 feet (ft) through the upstream part of the Madison Slide (cross sections 1400&ndash;800). Near cross section 800, the stream has eroded into the steep right bank between the stream and the road where fill was mechanically placed after 1970. Channel movement also was noted downstream from the Madison Slide. Near Raynolds Pass Bridge, about 3 miles (mi) downstream from Earthquake Lake, elevations across the channel have changed by -1.4 ft to +1.9 ft, but these changes were local in nature and could represent a few rocks or depressions in the bed. Overall, it does not appear that the materials eroded from the Madison Slide are causing aggradation in the subreach near the Raynolds Pass Bridge. Comparisons of critical shear stresses to mean boundary-shear stresses indicate that the D50 particle sizes (median size) along the right side of the bed between cross sections 400 and 500 and along the right side of the bed between cross sections 1300 and 1400 could move at the threshold streamflow. In contrast, most of the D84 particle sizes at those two locations probably will not move at the threshold streamflow. This lack of movement for the larger particles at the threshold streamflow could lead to further armoring of the bed as the D50 and smaller-sized particles are removed from the bed and transported downstream. The Shields parameter values from 0.04 to 0.08 that were used to calculate critical shear stresses could be conservative for a high-gradient stream such as the Madison. A higher, less conservative, Shields parameter would result in higher critical shear stresses, meaning that higher streamflows would be required to move material than those reported herein. In addition, because materials in the channel thalweg are exposed to higher boundary-shear stresses than the materials along the sides of the channel, larger, more erosion-resistant materials likely exist in the deeper parts of the channel where high-flow depths and velocities prevented sediment sampling. Movement of these materials might require higher critical shear stresses than estimated in this report. Characterization of sediment sizes in the center of the stream and observation of bed-material movement for a range of streamflows could provide information to help refine the Shields parameter and critical-shear stress estimates for bed materials in the Madison River downstream from Earthquake Lake. Furthermore, resurveying cross sections and water-surface elevations more frequently (either annually or after high streamflows) could better define the relation between streamflow and lateral and vertical channel movement.

Montana↗

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗

Concentrations, distribution, and sources of selenium from irrigated lands in western United States

Concentrations, distribution, and sources of selenium from irrigated lands were studied between 1986 and 1990 at 20 reconnaissance project areas in 17 western states under the Department of Interior's National Irrigation Water Quality Program. Samples of water, bottom sediment, whole‐body fish, and bird livers for analysis of selenium concentrations were collected before, during, and after irrigation season from streams, canals, lakes, and ground water in each project area. Selenium concentrations in water ranged from less than the detection limit of 1 μg/L in 42% of the 586 samples collected to 4,800 μg/L from a well in the Pine River Area, in southern Colorado. Selenium concentrations in 223 samples of bottom sediment ranged from less than the detection limit of 0.1 μg/g to 85 μg/g in a sample from the Middle Green River Basin, in Utah. Selenium concentrations in whole‐body fish (all species) ranged from 0.1 μg to 50 μg per gram dry weight, with the maximum concentration observed in a carp from the Gunnison River Basin, in western Colorado. Selenium concentrations in bird livers (all species) ranged from less than 0.32 μg to 170 μg per gram dry weight with the maximum concentration observed in the liver of an avocet from the Kendrick Reclamation Project, in Wyoming. Cretaceous‐age marine shales probably are the original sources for selenium in 16 of the 20 areas studied. In several areas, evaporative concentration of applied irrigation water and drainage of this water to wetlands, canals, streams, and lakes appears to be responsible for elevated concentrations of selenium.

Colorado, Utah, Wyoming↗

The influence of frequency and duration of seismic ground motion on the size of triggered landslides—A regional view

Observation, theory, and intuition all suggest that larger earthquakes should trigger larger landslides. Many factors could contribute to this, including depth-dependent shear strength or non-linearity of ground motion in soils and rock, but we hypothesize that the key characteristics of large earthquakes causing this phenomenon are (in addition to magnitude) the frequency and duration of the strong ground motion. Because of the paucity of site-specific data for detailed analysis, we take a regional approach to this question by analyzing strong-motion records and earthquake-induced landslide (EQIL) inventories from six well-documented earthquakes. Ground motion is characterized using earthquake magnitude and the median durations and frequencies (mean periods) of subsets of strong-motion records relevant to landslide triggering. EQIL inventories are characterized using the median landslide area of the entire inventory as well as the median areas of the largest 1% of the landslides and the largest 10 landslides. We then compare ground-motion characteristics with landslide size statistics to determine possible correlations. Comparisons of all earthquake- and landslide-size statistics show strong positive correlations between landslide size and (1) magnitude, (2) ground-motion duration, and (3) mean period. Although all the ground-motion measures yield highly correlated regressions, mean period appears to be the best overall predictor of landslide size. Landslide modeling using Newmark's sliding-block method also shows that longer mean periods and durations and larger magnitudes correlate strongly with increases in modeled displacements. These results support our hypothesis that increasing period and duration of seismic ground motion are the physical factors driving increased landslide sizes for larger earthquakes. Additional studies including data from a much larger set of earthquakes is needed to confirm the results of this initial study.

Engineering Geology↗

Development, Testing, and Sensitivity and Uncertainty Analyses of a Transport and Reaction Simulation Engine (TaRSE) for Spatially Distributed Modeling of Phosphorus in South Florida Peat Marsh Wetlands

Alterations to the predevelopment delivery of water and nutrients into the Everglades of southern Florida have been occurring for nearly a century. Major regional drainage projects, large-scale agricultural development, and changes to the hydrology of the Kissimmee River-Lake Okeechobee watershed have resulted in substantial phosphorus transport increases by surface waters. Excess phosphorus has accumulated in the soils of northern Everglades marshes to levels that have impaired the natural resources of the region. Regulations now limit the amount of phosphorous that enters the Everglades through an extensive network of water-control structures. This study involved the development and application of water-quality modeling components that may be applied to existing hydrologic models of southern Florida to evaluate the effects of different management scenarios. The result of this work is a spatially distributed water-quality model for phosphorus transport and cycling in wetlands. The model solves the advection-dispersion equation on an unstructured triangular mesh and incorporates a wide range of user-selectable mechanisms for phosphorus uptake and release parameters. In general, the phosphorus model contains transfers between stores; examples of stores that can be included are soil, water column (solutes), pore water, macrophytes, suspended solids (plankton), and biofilm. Examples of transfers are growth, senescence, settling, diffusion, and so forth, described with first order, second order, and Monod types of transformations. Local water depths and velocities are determined from an existing two-dimensional, overland-flow hydrologic model. The South Florida Water Management District Regional Simulation Model was used in this study. The model is applied to three case studies: intact cores of wetland soils with water, outdoor mesocosoms, and a large constructed wetland; namely, Cell 4 of Stormwater Treatment Area 1 West (STA-1W Cell 4). Different levels of complexity in the phosphorus cycling mechanisms were simulated in these case studies using different combinations of phosphorus reaction equations. Changes in water column phosphorus concentrations observed under the controlled conditions of laboratory incubations, and mesocosm studies were reproduced with model simulations. Short-term phosphorus flux rates and changes in phosphorus storages were within the range of values reported in the literature, whereas unknown rate constants were used to calibrate the model output. In STA-1W Cell 4, the dominant mechanism for phosphorus flow and transport is overland flow. Over many life cycles of the biological components, however, soils accrue and become enriched in phosphorus. Inflow total phosphorus concentrations and flow rates for the period between 1995 and 2000 were used to simulate Cell 4 phosphorus removal, outflow concentrations, and soil phosphorus enrichment over time. This full-scale application of the model successfully incorporated parameter values derived from the literature and short-term experiments, and reproduced the observed long-term outflow phosphorus concentrations and increased soil phosphorus storage within the system. A global sensitivity and uncertainty analysis of the model was performed using modern techniques such as a qualitative screening tool (Morris method) and the quantitative, variance-based, Fourier Amplitude Sensitivity Test (FAST) method. These techniques allowed an in-depth exploration of the effect of model complexity and flow velocity on model outputs. Three increasingly complex levels of possible application to southern Florida were studied corresponding to a simple soil pore-water and surface-water system (level 1), the addition of plankton (level 2), and of macrophytes (level 3). In the analysis for each complexity level, three surface-water velocities were considered that each correspond to residence times for the selected area (1-kilometer long) of 2, 10, and 20

Scientific Investigations Report↗

Geohydrology and water quality of stratified-drift aquifers in the lower Merrimack and coastal river basins, southeastern New Hampshire

Communities in the lower Merrimack River basin and coastal river basins of southeastern New Hampshire are experiencing increased demands for water because of a rapid increase in population. The population in 1987 was 225,495 and is expected to increase by 30 percent during the next decade. As of 1987, five towns used the stratified-drift aquifers for municipal supply and withdrew an estimated 6 million gallons per day. Four towns used the bedrock aquifer for municipal supply and withdrew an average of 1 .6 million gallons per day. Stratified-drift deposits cover 78 of the 327 square miles of the study area. These deposits are generally less than 10 square miles in areal extent, and their saturated thickness ranges front less than 20 feet to as much as 100 feet. Transmissivity exceeds 4,000 square feet per day in several locations. Stratified-drift aquifers in the eastern part are predominantly small ice-contact deposits surrounded by marine sediments or till of low hydraulic conductivity. Stratified-drift aquifers in the western part consist of ice-contact and proglacial deposits that are large in areal extent and are commonly in contact with surface-water bodies. Five stratified-drift aquifers, in the towns of Derry, Windham, Kingston, North Hampton, and Greenland, have the greatest potential to supply additional amounts of water. Potential yields and contributing areas of hypothetical supply wells were estimated for an aquifer in Windham near Cobbetts Pond and for an aquifer in Kingston along the Powwow River by use of a method analogous to superposition in conjunction with a numerical ground-waterflow model. The potential yield is estimated to be 0 .6 million gallons per day for the Windham-Cobbetts Pond aquifer and 4 .0 million gallons per day for the Kingston-Powwow River aquifer. Contributing recharge area for supply wells is estimated to be 1.6 square miles in the Windham-Cobbetts Pond aquifer and 4.9 square miles in the Kingston-Powwow River aquifer. Analyses of water samples from 30 wells indicate that the water quality in the basins studied is generally suitable for drinking and other domestic purposes. Concentrations of iron and manganese exceeded the U.S. Environmental Protection Agency's (USEPA) and the New Hampshire Water Supply Engineering Bureau's secondary maximum contaminant levels for drinking water in 20 samples. With one exception, concentrations of volatile organic compounds at all wells sampled met New Hampshire Water Supply and Engineering Bureau's drinking-water standards. At one well, trichloroethylene was detected at a concentration of 5.7 micrograms per liter. Ground-water contamination has been detected at several hazardous-waste sites in the study area. Currently, 5 sites are on the USEPA's National Priority List of superfund sites, 10 sites are Resource Conservation and Recovery Act of 1976 sites, and 1 site is a Department of Defense hazardous-waste site of stratigraphic layers is a product of a material's density and the velocity at which sound travels through that material. The reflected signals return to the hydrophones at the water surface and are then filtered, amplified, and displayed graphically on the chart recorder to allow interpretation of aquifer stratigraphy and bedrock depths. Lithologic data from nearby wells and test holes were used as control points to check the interpretation of the reflection profiles. Test drilling was done at 66 locations (pls. 1-3) to determine sediment grain size, stratigraphy, depth to water table, depth to bedrock, and ground water quality. A 6-inch-diameter, hollow-stem auger was used for test drilling. Split-spoon samples of subsurface materials collected at specific depths were used to evaluate the grain-size characteristics and identify the stratigraphic sequence of materials comprising the aquifers. Thirty-eight test holes cased with a 2-inch-diameter polyvinyl-chloride (PVC) pipe and slotted screens were used to make ground-water-level measurements and collect ground-water-quality samples. Surface-water-discharge measurements were made at 16 sites during low flow when the surface water is primarily ground-water discharge. These low-flow measurements indicate quantities of ground water potentially available from aquifers. Hydraulic conductivities of aquifer materials were estimated from grain-size-distribution data from 61 samples of stratified drift. Transmissivity was estimated from well logs by assigning hydraulic conductivity to specific well-log intervals, multiplying by the saturated thickness of the interval, and summing the results. Additional transmissivity values were obtained from an analysis of specific capacity and aquifer-test data. Long-term aquifer yields and contributing areas to hypothetical supply wells were estimated by application of a method that is analogous to super position and incorporates a ground-water-flow model developed by McDonald and Harbaugh (1988). This method was applied to two aquifers judged to have the best potential for providing additional ground-water supplies. Samples of ground water from 26 test wells and 4 municipal wells were collected in March and August 1987 for analysis of common inorganic, organic, and volatile organic constituents. Methods for collecting and analyzing the samples are described by Fishman and Freidman (1989). The water-quality results from the well samples were used to characterize background water quality in the stratified-drift aquifers.

New Hampshire↗

Mitigating flood risks in urban estuaries: Tidal dynamics, shoreline hardening, nature-based solutions, and floodgates in San Francisco Bay

Hydrodynamic models are valuable tools for understanding the primary factors influencing daily and peak water levels and for guiding discussions on potential adaptation strategies for managing flood risk in coastal areas. This analysis uses the Delft3D San Francisco Bay-Delta Community Model to simulate water levels and incorporates the effects of a number of adaptation measures in the urban San Francisco Bay estuary, California. In particular, we examine the influence of shoreline hardening, nature-based solutions, and subregional floodgates on regional water levels. The result shows that under present conditions, tidal amplification is responsible for generating a wide distribution of extreme water levels across San Francisco Bay. Tidal amplification is found to decrease under sea level rise, thereby producing a relative damping effect on extremes. A comparison of different shoreline scenarios demonstrates that hard frontal shorelines result in higher tidal amplification, whereas restored (soft) shorelines lower amplification. The current shoreline configuration has both hard and soft characteristics and results in an intermediate tidal response. In some areas, wetland restoration reduces extreme water levels by as much as 20 cm, whereas hard-shoreline addition elevates them by as much as 10 cm for 1.5 m of sea level rise. Furthermore, local floodgates can significantly reduce high water levels without major adverse effects elsewhere in San Francisco Bay. These findings point toward the justification for a range of adaptive measures across political boundaries, weighing hard and soft options in addressing the mounting danger of sea level rise.

California↗

Earthquake scenario and probabilistic ground-shaking hazard maps for the Albuquerque-Belen-Santa Fe, New Mexico, corridor

New Mexico's population is concentrated along the corridor that extends from Belen in the south to Espa&ntilde;ola in the north and includes Albuquerque and Santa Fe. The Rio Grande rift, which encompasses the corridor, is a major tectonically, volcanically, and seismically active continental rift in the western U.S. Although only one large earthquake (moment magnitude (M) &ge; 6) has possibly occurred in the New Mexico portion of the rift since 1849, paleoseismic data indicate that prehistoric surface-faulting earthquakes of M 6.5 and greater have occurred on aver- age every 400 yrs on many faults throughout the Rio Grande rift. We have developed a series of nine scenario and probabilistic hazard maps that portray the ground shaking that could occur in the Albuquerque-Belen-Santa Fe corridor from future earthquakes in New Mexico. These maps, at a scale of 1:500,000, display color-contoured ground-motion values in terms of the parameters of peak horizontal acceleration and horizontal spectral accelerations at 0.2 and 1.0 second (sec) periods. The maps depict surficial ground shaking and incorporate the site-response effects at locations underlain by unconsolidated sediments. The scenario maps are for a M 7.0 earthquake rupturing the Sandia-Rincon faults, which are adjacent to and dip west beneath Albuquerque. The probabilistic maps are for the two annual exceedance probabilities of building code relevance, 10% and 2% exceedance probabilities in 50 yrs (corresponding to return periods of 500 and 2,500 yrs, respectively). We included 57 Quaternary faults, all located within the Rio Grande rift, in the probabilistic seismic hazard analysis. These faults were characterized in terms of their geometry, rupture behavior (including possible segmentation), maximum expected earthquake magnitude, recurrence model, probability of activity, and slip rate. Preferred maximum magnitude values for these faults ranged from M 6.1 to 7.4 and preferred slip rates from 0.01 to approximately 0.12 mm/yr. Regional source zones and Gaussian smoothing of the historical seismicity were also included in the probabilistic hazard analysis to account for the hazard from background earthquakes (M &le; 6.5). A numerical ground-motion modeling approach and empirical attenuation relation- ships appropriate for extensional tectonic regimes were used to compute the scenario earthquake and probabilistic ground motions on rock. Amplification factors were then used to modify the rock motions and hence to incorporate site response into the hazard maps. These factors were based on three generalized geologic site-response categories (hard rock, soft rock, and firm/stiff soil) and were adopted from similar California-based categories because insufficient subsurface geologic and geotechnical data are available for the map area. The resulting hazard maps indicate that from both scenario and probabilistic perspectives, the ground-shaking hazard in the Albuquerque&ndash;Belen&ndash;Santa Fe corridor from future earthquakes could be severe, damaging, and potentially disastrous. In the event of a M 7.0 earthquake occurring on the Sandia&ndash;Rincon faults, ground shaking as characterized by peak ground acceleration could reach 0.7 g in much of the eastern half of the Albuquerque metropolitan area. (1 g = 980 cm/sec, the rate of gravitational acceleration.) These high ground motions will be attributable to the city&rsquo;s location directly over the Sandia&ndash;Rincon faults and the amplifying effect of the unconsolidated sediments within the Albuquerque Basin. These levels of ground shaking will probably result in severe damage to traditional adobe construction and even to modern buildings. Long- period ground motions (> 1.0 sec), which are significant to long and tall structures (e.g., tall buildings, long bridges, and highway overpasses), will also be high (> 1.0 g). Injuries and loss of life will be likely. For the 500- and 2,500-yr return period maps, the highest peak accelerations are predicted to be at the damaging levels of 0.3 g and 0.6 g, respectively. All maps show dramatically the frequency-dependent amplification of unconsolidated sediments in the basins along the Rio Grande valley (e.g., Albuquerque Basin). The pattern of amplification and deamplification is clearly a function of the distribution of unconsolidated sediments. These maps are not intended to be a substitute for site-specific studies for engineering design nor to replace standard maps commonly referenced in building codes. Rather, we hope that these maps will be used as a guide by government agencies; the engineering, urban planning, emergency preparedness, and response communities; and the general public as part of an overall program to reduce earthquake risk and losses in New Mexico.

New Mexico↗

Fishway Entrance Palisade

This technical report summarizes the work that was conducted by the University of Massachusetts Amherst and the United States Geological Survey (USGS), along with other project partners, on the Fishway Entrance Palisade (EP), a projected funded through the Department of Energy’s (DOE) funding opportunity titled ‘Innovative Solutions for Fish Passage at Hydropower Dams’ (DE‐FOA‐0001662). The period of performance ranged from September 1, 2018 through September 30, 2021. The EP is a novel fish passage engineering technology designed to provide more favorable entry conditions for fish and to reduce costs relative to conventional fishway auxiliary water systems (AWS). The EP project has four primary components. First, the Northeast United States Auxiliary Water Systems Database was created (Northeast Fishway Auxiliary Water Systems Database Section). The database, developed with material provided by the U.S. Fish and Wildlife Service, contains information on fishway type (e.g., lift, Denil, pool and weir) and Auxiliary Water System (AWS) details (e.g., water conveyance method, diffuser type) for 60 hydroelectric sites in the region. Findings indicate that nearly 4 out of every 10 fishway in the region is a fish lift and approximately 1 out of every 4 is a Denil ladder. The remainder are a mix of vertical slot fishways, pool and weirs, and Ice Harbor fishways. Furthermore, over half of all AWS systems use floor diffusers to discharge the auxiliary (or attraction) water into the entrance of a fishway, whereas only 14% use wall diffusers. Second, limited experiments on a conventional AWS with live, actively migrating fish were conducted at the USGS Easter Ecological Science Center (EESC) S.O. Conte Research Laboratory (Conventional Auxiliary Water System Experiments Section). This study determined how water velocity through a wall diffuser, without turning vanes or timber baffles to distribute the flow, affects the behavior and passage of adult American shad, a conservative surrogate species for migratory fish on the East Coast. Two gross diffuser velocity treatments were examined, 0.5 ft/s and 1.0 ft/s. These wall diffuser velocities represented current (0.5 ft/s) and past (1.0 ft/s) design criteria guidelines set forth by the USFWS North Atlantic-Appalachian Region (Rojas 2020; USFWS 2019). Six trials with a total of 151 American Shad were conducted in June of 2019 for the two treatments. No differences in American shad passage efficiency were discovered between the two treatments, while approximately 3 in every 4 attempts were successful at passing the diffuser. While these results may appear to indicate that the generally accepted gross wall diffuser velocity criteria for American shad of 0.5 ft/s could be safely increased to 1.0 ft/s, further analysis is warranted. Furthermore, it is unknown how other migratory and resident fish species that traverse these structures would be impacted by such a change. Studying the wall diffuser hydraulics led to an important AWS observation. Without turning vanes or timber baffles in this study, doubling the diffuser area was insufficient at producing the type of flow field change one may expect by halving the gross diffuser velocity. Instead, the flow fields throughout each treatments study area were similar, which led to similar results in shad performance. This not only highlights the importance of installing flow guidance devices like turning vanes, but also to the importance of properly maintaining them, which can be costly. Third, more expansive experiments on the novel EP were conducted in the spring of 2019 and 2021 (Fishway Entrance Palisade Experiments). The goal of this study was to determine how adult American shad responded to a variety of conditions at a full-scale EP. A total of six treatments were examined by changing the average auxiliary channel velocity between 1.0 and 5.0 ft/s in intervals of 1.0 ft/s and by inserting/removing an entrance gate at the opening of the fishway. Thirty trials with a total of 1,273 shad were conducted over the two years. In all treatments, at least ~7 out of every 10 fish successfully passed the EP diffuser and swam into the entrance channel within the 3.5-hour long trial, highlighting the general effectiveness of the novel AWS technology. In both study years, lower velocities through the EP diffuser led to increased shad performance, though performance peaked for the 2 ft/s velocity treatment. This treatment condition represents an approximate six-fold increase in gross diffuser velocity relative to conventional auxiliary water systems, which in turn presents opportunities for cost savings (e.g., reduction in diffuser size). Shad performance, in general, was worse in 2019 than in 2021, potentially due to the different run timing when our trials were conducted (2019 trials occurred near the end of the migration season, unlike in 2021). Treatments in 2019 had approximately a 20% reduction in entrance efficiency by the trial end, including a 16.7% drop for the 3 ft/s velocity treatment in 2019 relative to 2021 (the only carryover treatment between years). Lastly, adding an entrance gate caused a significant delay to entry. The time to 25% entry raised ~20 minutes from the near instantaneous 25% entry that was reported for the other treatments conducted in the same year (2021). Though by the end of the 3.5-hour trial, the overall entrance efficiency nearly matched those of the other 2021 treatments. The fourth and final component of the EP project was an economic analysis that focused on the cost of attraction and environmental flows (Modeling Power Generation Losses Due to Environmental and Fish Passage Attraction Flows at a Run-Of-River Hydroelectric Operation in the Northeast). The study assessed the economic impact of meeting environmental flow requirements at a representative hydroelectric facility and fish lift in the Northeast. An initial finding of the study was that there is a paucity of published data on the costs of meeting attraction and environmental flows. This is due, in part, to the proprietary nature of this data. To explore the costs associated with these flows, three types of environmental flows were assessed: upstream fishway attraction flows, downstream fishway attraction flows, and habitat maintenance minimum flows. A physics-based model was developed and calibrated with three years of hourly generation and flow data as inputs. Gage flow inputs were adjusted and used to calculate power generated. To address hydrologic variability, the model was executed to simulate 30 years of historical flows. Results indicate that both interannual and seasonal climatic factors impact the costs of meeting environmental flow requirements. Generation potential is most strongly curtailed during dry years in terms of maximizing the capacity factor (the percent of time a plant generates at capacity). Dry years, and especially dry summers, have the most significant costs associated with mitigation flows. Of the three types of flows, habitat flows are most costly in terms of power production, followed by upstream attraction flows. Downstream attraction flows are least costly. This finding is the likely result of differences in both flow rates and duration of the seasonal requirement for each flow. Overall, environmental flows represented a 2-12% loss in annual generation, but losses during a dry summer can reach over 20%.

Final Technical Report↗

Field determination of optimal dates for the discrimination of invasive wetland plant species using derivative spectral analysis

Mapping invasive plant species in aquatic and terrestrial ecosystems helps to understand the causes of their progression, manage some of their negative consequences, and control them. In recent years, a variety of new remote-sensing techniques, like Derivative Spectral Analysis (DSA) of hyperspectral data, have been developed to facilitate this mapping. A number of questions related to these techniques remain to be addressed. This article attempts to answer one of these questions: Is the application of DSA optimal at certain times of the year? Field radiometric data gathered weekly during the summer of 1999 at selected field sites in upstate New York, populated with purple loosestrife ( Lythrum salicaria L.), common reed ( Phragmites australis (Cav.)) and cattail ( Typha L.) are analyzed using DSA to differentiate among plant community types. First, second and higher-order derivatives of the reflectance spectra of nine field plots, varying in plant composition, are calculated and analyzed in detail to identify spectral ranges in which one or more community types have distinguishing features. On the basis of the occurrence and extent of these spectral ranges, experimental observations suggest that a satisfactory differentiation among community types was feasible on 30 August, when plants experienced characteristic phenological changes (transition from flowers to seed heads). Generally, dates in August appear optimal from the point of view of species differentiability and could be selected for image acquisitions. This observation, as well as the methodology adopted in this article, should provide a firm basis for the acquisition of hyperspectral imagery and for mapping the targeted species over a broad range of spatial scales.

Photogrammetric Engineering and Remote Sensing↗

Methods for estimating selected low-flow frequency statistics and mean annual flow for ungaged locations on streams in North Georgia

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources, Environmental Protection Division, developed regional regression equations for estimating selected low-flow frequency and mean annual flow statistics for ungaged streams in north Georgia that are not substantially affected by regulation, diversions, or urbanization. Selected low-flow frequency statistics and basin characteristics for 56 streamgage locations within north Georgia and 75 miles beyond the State’s borders in Alabama, Tennessee, North Carolina, and South Carolina were combined to form the final dataset used in the regional regression analysis. Because some of the streamgages in the study recorded zero flow, the final regression equations were developed using weighted left-censored regression analysis to analyze the flow data in an unbiased manner, with weights based on the number of years of record. The set of equations includes the annual minimum 1- and 7-day average streamflow with the 10-year recurrence interval (referred to as 1Q10 and 7Q10), monthly 7Q10, and mean annual flow. The final regional regression equations are functions of drainage area, mean annual precipitation, and relief ratio for the selected low-flow frequency statistics and drainage area and mean annual precipitation for mean annual flow. The average standard error of estimate was 13.7 percent for the mean annual flow regression equation and ranged from 26.1 to 91.6 percent for the selected low-flow frequency equations. The equations, which are based on data from streams with little to no flow alterations, can be used to provide estimates of the natural flows for selected ungaged stream locations in the area of Georgia north of the Fall Line. The regression equations are not to be used to estimate flows for streams that have been altered by the effects of major dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, or wastewater discharges. The regression equations should be used only for ungaged sites with drainage areas between 1.67 and 576 square miles, mean annual precipitation between 47.6 and 81.6 inches, and relief ratios between 0.146 and 0.607; these are the ranges of the explanatory variables used to develop the equations. An attempt was made to develop regional regression equations for the area of Georgia south of the Fall Line by using the same approach used during this study for north Georgia; however, the equations resulted with high average standard errors of estimates and poorly predicted flows below 0.5 cubic foot per second, which may be attributed to the karst topography common in that area. The final regression equations developed from this study are planned to be incorporated into the U.S. Geological Survey StreamStats program. StreamStats is a Web-based geographic information system that provides users with access to an assortment of analytical tools useful for water-resources planning and management, and for engineering design applications, such as the design of bridges. The StreamStats program provides streamflow statistics and basin characteristics for U.S. Geological Survey streamgage locations and ungaged sites of interest. StreamStats also can compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in Georgia.

Georgia↗