Geology ReportsSearch

SEARCH · Geology Reports

Results for “Lateral Lines”

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.

269 records · Page 15Linked to original sources

Saudi Arabian seismic deep-refraction profiles; final project report

In February 1978 a seismic deep-refraction profile was recorded by the U.S. Geological Survey along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat-Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used: five on land, with most charges placed below the water table in drill holes, and one at sea, with charges placed on the sea floor and detonated from a ship. Slightly more than 61 metric tons of explosives were used in 19 discrete firings. Seismic energy was recorded by 100 newly-developed portable seismic stations deployed in approximately 200 km-long arrays for each firing. Each station consisted of a standard 2-Hz vertical component geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. In this final report, we fully document the field and data-processing procedures and present the final seismogram data set as both a digital magnetic tape and as record sections for each shot point. Record sections include a normalized set of seismograms, reduced at 6 km/s, and a true-amplitude set, reduced at 8 km/s, which have been adjusted for amplifier gain, individual shot size, and distance from the shot point. Appendices give recorder station and shot information, digital data set descriptions, computer program listings, arrival times used in the interpretation, and a bibliography of reports published as a result of this project. We used two-dimensional ray-tracing techniques in the data analysis, and our interpretation is based primarily on horizontally layered models. The Arabian Shield is composed, to first-order, of two layers, each about 20 km thick, with average velocities of 6.3 km/s and 7.0 km/s, respectively. At the western shield margin the crust thins to less than 20 km total thickness, beyond which the Red Sea shelf and coastal plain are interpreted to be underlain by oceanic crust. A major crustal lateral velocity inhomogeneity northeast of Sabhah in the Shammar Tectonic Province is interpreted as the suture zone of two crustal blocks of different composition. Several high-velocity anomalies in the upper crust correlate with mapped gneissic dome structures. Two intra-crustal reflectors at13 km depth are interpreted as the tops of mafic intrusives. The Mohorovicic discontinuity beneath the shield varies from 43 km depth in the northeast with 8.2 km/s mantle velocity to 38 km depth in the southwest with 8.0 km/s mantle velocity. Two velocity discontinuities are identified in the upper mantle, at 59 and 70 km depth. We suggest further work, including refined analyses of the data employing filtering and synthetic seismogram techniques, as well as consideration of attenuation properties. Extension of the seismic refraction profile to the Arabian Gulf and some short profiles perpendicular to the existing profile would be fruitful areas for future field work.

Open-File Report

Geology and ground-water resources of the islands of Lanai and Kahoolawe, Hawaii

Lanai lies 59 miles southeast of Honolulu, Oahu, has an area of 141 square miles, and is 3,370 feet high. (See fig. 1 and pl. 1.) Lanai City is the only town of importance. The island produces pineapples and cattle. The surface above about 1,200 feet is generally covered with lateritic soil, which reaches a maximum depth of about 50 feet. Below this level the island is partly devoid of vegetation and is strewn with boulders, the result of having been once submerged by the ocean to this depth. Traces of various emerged and submerged shore lines are described, the highest fossiliferous marine deposits being 1,070 feet above sea level. Lanai is an eroded extinct basaltic volcano built during one period of activity. No secondary eruptions occurred as on most of the other islands. It has three rift zones and a summit caldera. The summit plateau has resulted from collapse along the northwest rift zone. Elsewhere there is much evidence of faulting. About 100 faults and 275 dikes were recorded, but they are so close together in places that it was not possible to show them all on the map. The climate is semitropical, the mean annual temperature of Lanai City, altitude 1,620 feet, being 68° F. Because Lanai lies to the lee of Maui Island it is dry. The mean annual rainfall ranges from 38 inches on the summit to less than 10 inches on the coast. The windward (northeast) side is carved by streams into deep canyons. Maunalei Gulch has the only perennial stream, and it does not reach the sea. Ground water, the lifeblood of Lanai is scarce. Lanai City obtains some of its water supply by a tunnel from gravel in Maumalei Gulch. This water apparently rises from the dike complex in this gulch. The rest of the supply comes from a recently constructed shaft tapping the dike complex not far downstream. The total quantity of high-level ground water discharged by springs and tunnels ranges from about 600,000 gallons a day in wet weather to about 250,000 gallons a day in dry weather. The basal water, although potable, is fairly high in salt. Several sites are recommended for developing and conserving ground water. Kahoolawe Island is 11 miles long, 6 miles wide, 1,491 feet high, covers 45 square miles, and lies 94 miles southeast of Honolulu and 6 3 / 4 miles southwest of Maui. It is a shield-shaped extinct volcano composed chiefly of thin flows of primitive basalt poured in rapid succession from three rift zones and a vent at their intersection. At one stage the volcano was indented with a caldera about 3 miles across which was later completely filled. A graben led southwestward from it. The rocks are divided into Late Tertiary (?) or early Pleistocene(?) pre-caldera basalts, caldera-filling basalts and basaltic andesites, post-caldera basalts and andesites, and Recent post erosional basalts. A few thin vitric tuff beds and cinder cones were found. Marine erosion has cut cliffs as high as 800 feet along the east and south shores and exposed a cross section of the caldera. Only shallow ephemeral gulches exist. The entire summit has been eroded to a hard-pan surface by the wind as a result of the vegetation being destroyed by livestock. The island is semi-arid and well water is needed for stook: "The stook is now supplied entirely from storage of rain and flood waters. During droughts water is hauled by boat from the island of Maui. All the wells dug so far yield water that is too brackish for stock except at the fairly inaccessible south side of Kanapou Bay. The resistivity survey indicates a water table 1.5 feet or less above sea level for 2.25 miles inland. A few sites for wells are recommended in the dike complex where small supplies of water suitable for stock might be found. Petrographic studies by Gordon A. Macdonald indicate that the pre-caldera and caldera-filling lavas are largely normal olivine basalt of the type which forms the bulk of all Hawaiian volcanoes thus far investigated. It represents the undifferentiated magma of the Hawaiian petrographic province. Toward the close of the caldera-filling epoch the vent became less active, and magmatic differentiation produced basaltic andesites, which are interbedded with normal basalts. The post-caldera lavas are largely basaltic andesites and andesites. The much younger lavas, erupted after a period of extensive erosion, are olivine basalts similar in composition to the pre-caldera flows. The mineralogy of the Kahoolawe rocks is described in detail.

Hawaii

Ground water in the Crow Creek-Sand Lake area, Brown and Marshall Counties, South Dakota

The Crow Creek Sand Lake area is in the northeastern part of South Dakota and consists of flat to gently rolling land. It is roughly 25 miles square and is bounded on the west by the James River, on the north by the North Dakota State line, and on the south and east by arbitrary land lines. Natural drainage, which is very poor, is augmented by a now very dilapidated system of man made drains; undrained or poorly drained potholes and sloughs are common. Most of the report area is in the drainage basin of Crow Creek, a tributary of the James River. About 60 square miles along the west border of the area is drained directly by the James River; and about 36 square miles in the northeastern corner of the area is drained by the Wild Rice River, a tributary of the Red River of the North. Granite and other related rocks of Precambrian age underlie the entire area below a depth of about 1,100 to 1,200 feet. Under part of the area these rocks are overlain by the Sioux quartzite, also of Precambrian age. Next above is the Dakota sandstone of Cretaceous age, which is about 250 feet thick and is the aquifer tapped by more than 600 deep wells in the area. The Dakota sandstone in turn is overlain by younger formations of Cretaceous age, which include the Graneros shale, Greenhorn limestone, Carlile shale, Niobrara formation, and Pierre shale. These younger formations have an aggregate thickness of about 700 to 950 feet and generally do not yield water. Mantling, the consolidated stratified rocks are unconsolidated deposits of Quaternary age that are absent locally but in places are as much as 200 feet thick. These deposits consist principally of glacial outwash, glacial till, and lake sediments and are the source of water for all shallow wells in the area. Precipitation is the principal source of recharge to the ground-water reservoir tapped by the shallow wells. The water table ranges in depth from less than 1 foot to a little more than 30 feet below the land surface; in July 1951 the depth to water was less than 5 feet in 15 percent of the area and 5 to 10 feet in an additional 25 percent of the area. Because the hydraulic gradient in most of the area is less than 10 feet per mile and because much of the water-bearing material is only slightly permeable, lateral movement of the ground water in the unconsolidated deposits is exceedingly slow. Transpiration by vegetation and evaporation account for most of the ground water discharged in the area; relatively little ground water is discharged by streams, wells, and springs, or as underflow out of the area. Much of the low-lying land is waterlogged. However, by improving and cleaning both the natural and artifical drains and by lowering the water table by pumping, waterlogging can be materially decreased and its recurrence prevented.

North Dakota

Geophysical studies of the Crump Geyser known geothermal resource area, Oregon, in 1975

The U.S. Geological Survey (USGS) conducted geophysical studies in support of the resource appraisal of the Crump Geyser Known Geothermal Resource Area (KGRA). This area was designated as a KGRA by the USGS, and this designation became effective on December 24, 1970. The land classification standards for a KGRA were established by the Geothermal Steam Act of 1970 (Public Law 91-581). Federal lands so classified required competitive leasing for the development of geothermal resources. The author presented an administrative report of USGS geophysical studies entitled 'Geophysical background of the Crump Geyser area, Oregon, KGRA' to a USGS resource committee on June 17, 1975. This report, which essentially was a description of geophysical data and a preliminary interpretation without discussion of resource appraisal, is in Appendix 1. Reduction of sheets or plates in the original administrative report to page-size figures, which are listed and appended to the back of the text in Appendix 1, did not seem to significantly degrade legibility. Bold print in the text indicates where minor changes were made. A colored page-size index and tectonic map, which also show regional geology not shown in figure 2, was substituted for original figure 1. Detailed descriptions for the geologic units referenced in the text and shown on figures 1 and 2 were separately defined by Walker and Repenning (1965) and presumably were discussed in other reports to the committee. Heavy dashed lines on figures 1 and 2 indicate the approximate KGRA boundary. One of the principal results of the geophysical studies was to obtain a gravity map (Appendix 1, fig. 10; Plouff, and Conradi, 1975, pl. 9), which reflects the fault-bounded steepness of the west edge of sediments and locates the maximum thickness of valley sediments at about 10 kilometers south of Crump Geyser. Based on the indicated regional-gravity profile and density-contrast assumptions for the two-dimensional profile, the maximum sediment thickness was estimated at 820 meters. A three-dimensional gravity model would have yielded a greater thickness. Audiomagnotelluric measurements were not made as far south as the location of the gravity low, as determined in the field, due to a lack of communication at that time. A boat was borrowed to collect gravity measurements along the edge of Crump Lake, but the attempt was curtailed by harsh, snowy weather on May 21, 1975, which shortly followed days of hot temperature. Most of the geophysical data and illustrations in Appendix 1 have been published (Gregory and Martinez, 1975; Plouff, 1975; and Plouff and Conradi, 1975), and Donald Plouff (1986) discussed a gravity interpretation of Warner Valley at the Fall 1986 American Geophysical Union meeting in San Francisco. Further interpretation of possible subsurface geologic sources of geophysical anomalies was not discussed in Appendix 1. For example, how were apparent resistivity lows (Appendix 1, figs. 3-6) centered near Crump Geyser affected by a well and other manmade electrically conductive or magnetic objects? What is the geologic significance of the 15-milligal eastward decrease across Warner Valley? The explanation that the two-dimensional gravity model (Appendix 1, fig. 14) was based on an inverse iterative method suggested by Bott (1960) was not included. Inasmuch as there was no local subsurface rock density distribution information to further constrain the gravity model, the three-dimensional methodology suggested by Plouff (1976) was not attempted. Inasmuch as the associated publication by Plouff (1975), which released the gravity data, is difficult to obtain and not in digital format, that report is reproduced in Appendix 2. Two figures of the publication are appended to the back of the text. A later formula for the theoretical value of gravity for the given latitudes at sea level (International Association of Geodesy, 1971) should be used to re-compute gravity anomalies. To merge t

Open-File Report

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

Documentation of the Surface-Water Routing (SWR1) Process for modeling surface-water flow with the U.S. Geological Survey Modular Ground-Water Model (MODFLOW-2005)

A flexible Surface-Water Routing (SWR1) Process that solves the continuity equation for one-dimensional and two-dimensional surface-water flow routing has been developed for the U.S. Geological Survey three-dimensional groundwater model, MODFLOW-2005. Simple level- and tilted-pool reservoir routing and a diffusive-wave approximation of the Saint-Venant equations have been implemented. Both methods can be implemented in the same model and the solution method can be simplified to represent constant-stage elements that are functionally equivalent to the standard MODFLOW River or Drain Package boundary conditions. A generic approach has been used to represent surface-water features (reaches) and allows implementation of a variety of geometric forms. One-dimensional geometric forms include rectangular, trapezoidal, and irregular cross section reaches to simulate one-dimensional surface-water features, such as canals and streams. Two-dimensional geometric forms include reaches defined using specified stage-volume-area-perimeter (SVAP) tables and reaches covering entire finite-difference grid cells to simulate two-dimensional surface-water features, such as wetlands and lakes. Specified SVAP tables can be used to represent reaches that are smaller than the finite-difference grid cell (for example, isolated lakes), or reaches that cannot be represented accurately using the defined top of the model. Specified lateral flows (which can represent point and distributed flows) and stage-dependent rainfall and evaporation can be applied to each reach. The SWR1 Process can be used with the MODFLOW Unsaturated Zone Flow (UZF1) Package to permit dynamic simulation of runoff from the land surface to specified reaches. Surface-water/groundwater interactions in the SWR1 Process are mathematically defined to be a function of the difference between simulated stages and groundwater levels, and the specific form of the reach conductance equation used in each reach. Conductance can be specified directly or calculated as a function of the simulated wetted perimeter and defined reach bed hydraulic properties, or as a weighted combination of both reach bed hydraulic properties and horizontal hydraulic conductivity. Each reach can be explicitly coupled to a single specific groundwater-model layer or coupled to multiple groundwater-model layers based on the reach geometry and groundwater-model layer elevations in the row and column containing the reach. Surface-water flow between reservoirs is simulated using control structures. Surface-water flow between reaches, simulated by the diffusive-wave approximation, can also be simulated using control structures. A variety of control structures have been included in the SWR1 Process and include (1) excess-volume structures, (2) uncontrolled-discharge structures, (3) pumps, (4) defined stage-discharge relations, (5) culverts, (6) fixed- or movable-crest weirs, and (7) fixed or operable gated spillways. Multiple control structures can be implemented in individual reaches and are treated as composite flow structures. Solution of the continuity equation at the reach-group scale (a single reach or a user-defined collection of individual reaches) is achieved using exact Newton methods with direct solution methods or exact and inexact Newton methods with Krylov sub-space methods. Newton methods have been used in the SWR1 Process because of their ability to solve nonlinear problems. Multiple SWR1 time steps can be simulated for each MODFLOW time step, and a simple adaptive time-step algorithm, based on user-specified rainfall, stage, flow, or convergence constraints, has been implemented to better resolve surface-water response. A simple linear- or sigmoid-depth scaling approach also has been implemented to account for increased bed roughness at small surface-water depths and to increase numerical stability. A line-search algorithm also has been included to improve the quality of the Newton-step upgrade vector, if possible. The SWR1 Process has been benchmarked against one- and two-dimensional numerical solutions from existing one- and two-dimensional numerical codes that solve the dynamic-wave approximation of the Saint-Venant equations. Two-dimensional solutions test the ability of the SWR1 Process to simulate the response of a surface-water system to (1) steady flow conditions for an inclined surface (solution of Manning's equation), and (2) transient inflow and rainfall for an inclined surface. The one-dimensional solution tests the ability of the SWR1 Process to simulate a looped network with multiple upstream inflows and several control structures. The SWR1 Process also has been compared to a level-pool reservoir solution. A synthetic test problem was developed to evaluate a number of different SWR1 solution options and simulate surface-water/groundwater interaction. The solution approach used in the SWR1 Process may not be applicable for all surface-water/groundwater problems. The SWR1 Process is best suited for modeling long-term changes (days to years) in surface-water and groundwater flow. Use of the SWR1 Process is not recommended for modeling the transient exchange of water between streams and aquifers when local and convective acceleration and other secondary effects (for example, wind and Coriolis forces) are substantial. Dam break evaluations and two-dimensional evaluations of spatially extensive domains are examples where acceleration terms and secondary effects would be significant, respectively.

Techniques and Methods

Trapping and migration of methane associated with the gas hydrate stability zone at the Blake Ridge Diapir: New insights from seismic data

The Blake Ridge Diapir is the southernmost of a line of salt diapirs along the Carolina trough. Diapirs cause faulting of the superjacent sediments, creating pathways for migration of fluids and gas to the seafloor. We analyzed reflection seismic data from the Blake Ridge Diapir, which is located in a region with known abundant gas hydrate occurrence. A striking feature in these data is a significant shallowing of the base of gas hydrate stability (BGHS) over the center of the diapir. The seafloor is warped up by about 100 m above the diapir, from about 2300 m to about 2200 m. The BGHS, as indicated by a bottom simulating reflection (BSR), is about 4.5 s off the flanks of the diapir, rising to about 4.15 s at the center. Above the diapir, a fault system appears to rise vertically from the BGHS to about 0.05 s below the seafloor (40-50 m); it then diverges into several steeply dipping faults that breach the seafloor and cover an area ~700 m in diameter. Other secondary faults diverge from the main fault or emerge directly from the BGHS near the crest of the diapir. Gas and other fluids may migrate upward through the faults. We performed complex trace analysis to compare the reflection strength and instantaneous frequency along individual reflections. A low-frequency anomaly over the center of the diapir indicates high seismic attenuation. This is interpreted to be caused by migration of fluids (probably methane) along fault pathways. The migration of gas (i.e. probably mainly methane) through the gas hydrate stability zone is not yet understood. We speculate that pore fluids in the faults may be too warm and too salty to form gas hydrate, even at depths where gas hydrate is stable away from the diapir. Alternatively, gas hydrates may seal the fault walls such that water supply is too low to transform all the gas into gas hydrates. The shallowing of the BSR may reflect increased heatflow above the diapir either caused by the high thermal conductivity of the underlying salt or by advective heat transport along with fluids. High pore water salinity shifts the gas hydrate stability to lower temperatures and may also play a significant role in BSR shallowing. We, therefore, investigated the possible effect of pore water salinity on shallowing of the BSR. We found that BSR shallowing may theoretically be entirely caused by increased salinity over the diapir, although geologically this would not be reasonable. This observation demonstrates the potential importance of pore water salinity for lateral variations of BSR depths, in particular, above salt structures: (C) 2000 Elsevier Science B.V.

Blake Ridge Diapir

Gold placers of the historical Fortymile River region, Alaska

The Fortymile River region in east-central Alaska has a long and colorful history as the site of the first major gold discovery in interior Alaska. Placer gold has been mined in the region nearly every year since its original discovery in 1886. Total gold production is approximately 500,000 troy ounces. Although many of the rich deposits have been mined, there still exist areas that contain gold. Areas of mined and unmined gold-bearing creek and terrace gravels are outlined on the accompanying geologic map. The early history of the Fortymile area centered on the small frontier settlement of Fortymile City located at the junction of the Fortymile and Yukon Rivers in Canadian territory. This was the supply and jumping-off point for prospectors who worked their way into Alaska up the Fortymile River and found gold on many of its tributaries. Hand mining, both underground and surface, using sluice boxes and (or) rockers were the earliest methods; later, hydraulicking, dredging, and draglining methods were used. More recently, bulldozers and elevated trammels have been used, as well as very portable floating suction dredges. The rich mining lore of the area is closely associated with events of the nearby world-famous Klondike District. Bedrock and placer geology and mining history of individual gold-rich creeks are herein updated. The Fortymile area, which is part of the Yukon-Tanana Upland, contains quartzite, schist, gneiss, amphibolite, marble, serpentinite, and granite overlain by basalt, sandstone, conglomerate, shale, tuff, and coal; overlying these rocks are several deposits of varying ages consisting of gold-bearing gravel and colluvium. The close spatial association of creeks containing placer gold and the gneiss, schist, amphibolite, and marble unit strongly suggests this metamorphic unit is the gold source. High terrace gravels record a time from the late Tertiary to early Pleistocene when the ancestral Fortymile River and its major tributaries, the North and South Forks, had floodplains roughly 1 to 2 miles (2-3 kilometers) wide and gradients of about 4 feet per mile (0.75 meters per kilometer). Base-level lowering during the post-early Pleistocene caused the rivers to cut into their floodplains and to develop the youthful characteristics they have today such as V-shaped canyons, narrow floodplains, and gradients of at least twice those of the old river. Colluvium marginal to creek deposits in steep-sided valleys is often gold bearing. Much of the unconsolidated gravel within the major drainages of the Fortymile River, South Fork, North Fork, and Mosquito Fork is colluvium. Heavy-mineral-concentrate samples from the gold-producing creeks and high terrace gravels contain varying amounts of magnetite (20 to 80 percent) and ilmenite (10 to 30 percent), and samples from creeks draining areas principally composed of metamorphic rocks contain abundant garnet (10 to 30 percent). Gold fineness ranges from 620 to 927, but it is difficult to attach any geologic significance to the fineness data. Most placer gold in the Fortymile River area has been recovered at, or near, the gravelbedrock contact. The lowermost 3.3 feet (1 meter) of gravel and the uppermost 1.6 feet (0.5 meter) of bedrock may contain as much as 80 to 90 percent of the gold that is ultimately recovered. Gold nuggets are rare and most of the gold recovered is in the form of flattened fragments less than .2 inches (5 millimeters) in greatest dimension. However, large gold nuggets have been found on Wade Creek; examples are ones of 25,33,56, and 70 ounces. Occasionally, large nuggets may still be found in the tailing piles along the creek. The Fortymile River and its tributaries the South Fork, Walker Fork, and Mosquito Fork, all of which at one time were the sites of bucket-line dredge operations, now are almost exclusively mined using floating suction dredges. Unmined gold-bearing gravel is present in the floodplain of the Walker Fork valley below Cherry Creek and in low (about 100 to 130 feet or 30 to 40 meters) terraces along the north side of Walker Fork and east side of Cherry Creek. Considering the locations of where most gold has been found in the South Fork valley both by the older bucket dredges and the modern suction dredges, it seems likely that the tributary drainages of Lost Chicken, Napoleon, Franklin, and Buckskin Creeks have supplied the bulk of the gold to the South Fork valley. A quarter acre (0.10 hectare), 130-foot-thick ( 40 meters) section of the high terrace gravels on the north side of Napoleon Creek was mined for placer gold and yielded values estimated to be $8.50 per cubic yard (or $6.50 per cubic meter) at $350 per troy ounce. The unmined high terrace gravels on the south side of Buckskin Creek contain gold; however, this gravel is only 3 to 6.5 feet (1 to 2 meters) thick. The search for a lode gold source in the Fortymile River region may be in vain, because substantially more gold than has been recovered from the placers can be derived by the gradual erosion of large volumes of source rocks that contain background mean gold amounts. Using Leon's mass balance equation, 5,167 metric tons of gold may exist in the placers of the Fortymile River region, less than 1 percent of the recovered amount of 15.6 tons. The largest gold resource remaining in the Fortymile River region is probably in the high terrace gravels exposed along many of the creeks and rivers. Until there is exploratory drilling or a comprehensive sampling program, the amount of gold in these gravels will remain unknown. Environmental constraints imposed by Federal and State agencies have slowed, but not stopped, placer mining in the Fortymile River area, and a significant gold price rise would result in more mining.

Alaska

Geophysical analysis of the Salmon Peak Formation near Amistad Reservoir Dam, Val Verde County, Texas, and Coahuila, Mexico, March 2006, to aid in piezometer placement

Since 1992, numerous sinkholes have developed northwest of the Amistad Reservoir dam on the Rio Grande. Increases in the discharge of springs south of the dam, on the western side of the Rio Grande, in Coahuila, Mexico, have been documented. In 1995 the Mexico Section of the International Boundary and Water Commission (IBWC) completed a study of the western embankment (Coahuila, Mexico) of the dam that included surface geophysics, borehole geophysics, and installation of piezometers to learn more about subsurface conditions. As part of a 5-year safety inspection in 2005, technical advisors recommended that one line of similarly constructed piezometers be installed on the eastern embankment (Val Verde County, Texas) of the dam for comparison of water levels (potentiometric head) on both the western and eastern embankments of Amistad Reservoir dam. To provide technical assistance for the horizontal and vertical placement of piezometers on the eastern embankment of Amistad Reservoir dam, the U.S. Geological Survey, in cooperation with the U.S. Section of the IBWC, conducted a study along both the western and eastern embankments of Amistad Reservoir dam. The study involved an integrated approach using surface and borehole geophysical methods. In the western embankment investigation, geological and geophysical characteristics that indicate relatively large water-yielding properties of the Salmon Peak Formation were identified. The direct-current (DC) resistivity method was selected as the surface geophysical reconnaissance technique to correlate relatively large water-yielding properties of the Salmon Peak Formation, identified from analysis of borehole geophysical logs, with variations in subsurface resistivity. The dipole-dipole array and the reciprocal Schlumberger array were selected as the most applicable DC resistivity arrays. Two resistivity units were identified in both the dipole-dipole array data and the reciprocal Schlumberger array data along DC resistivity profiles on both embankments. Resistivity unit 1 generally is of relatively low resistivity, ranging from 45 to 150 ohm-meters compared with resistivity unit 2, which ranges from 120 to 345 ohm-meters (depending on the DC array type). The presence of mapped sinkholes in the reservoir north of the western embankment study area and the zone of increased water content (as indicated by zones of low neutron log count rates in nearby piezometers) leads to the conclusion that resistivity unit 1 is a preferential flow path where surface water from Amistad Reservoir is forced into the ground-water system (because of increased head from the reservoir). In the eastern embankment investigation, trends in the spatial distribution of sinkholes and the occurrence of weathered zones were identified from geologic descriptions of cores. The correlation of surface geophysical DC resistivity, historical lithologic data, and general trend of documented sinkholes along the eastern end of the eastern embankment profile were used to justify further exploration (drilling of piezometers) in the eastern expression of resistivity unit 1. The spatial location of the piezometers and the screened intervals were selected to best match the locations of the screened intervals of the western embankment piezometers. Six piezometers were installed on the eastern embankment and logged using borehole geophysical techniques. Surface DC resistivity sections superimposed on the resistivity logs for two piezometers indicate three discernible resistivity units that correlate with resistivity units 2, 1, and 2, respectively, identified in the western embankment study area. Resistivity units 1 and 2 in the DC resistivity profiles generally correspond with low and high resistivity zones, respectively, on the normal and lateral resistivity logs collected in the nearby piezometers at the time of installation.

Scientific Investigations Report

Effects of runoff controls on the quantity and quality of urban runoff at two locations in Austin, Texas

Rapid urban development in the Austin metropolitan area, Texas, is causing concern about increasing peak discharges from storm runoff and the degradation of the quality of water in receiving streams, lakes, and aquifers. In an attempt to decrease peak discharges and improve water quality, runoff controls are being required in some watersheds. This report summarizes the precipitation, streamflow, and water-quality data collected from September 1982 to September 1984 upstream and downstream from runoff controls at two locations, and presents the effects of these runoff controls on streamflow and the quality of runoff water. The two controls are a detention and filtering pond near Barton Creek Square Shopping Center, a large shopping center southwest of downtown Austin, and a grass-swale control in the Al ta Vista Planned Unit Development, a multipie-family housing area. At Barton Creek Square Shopping Center, rainfall for the storms analyzed ranged from 0.14 to 2.88 inches. The rainfall rate for the September 7, 1983, storm exceeded the 100-year return period for the 5- and 10-minute duration and was equal to the 50-year return period for the 15-minute duration. Peak discharge at the inflow station to the detention pond was closely related to the maximum rainfall during a 5-minute period and occurred about 10 minutes later. The maximum inflow at this station was 185 cubic feet per second and appeared to be the limit of the storm sewer system. For small- and moderate-sized storms, the runoff is contained in the detention pond and passes through a filter system. Runoff from large storms overflows into the drop outlet. For storms contained in the pond, peak discharges at the outflow station generally were less than 3.1 cubic feet per second. As time passed, the outflow peak discharges tended to decrease as a result of reduced permeability of the filter. Cleaning the filter appeared to increase the peak flows but did not restore them to the previous level. The runoff-rainfall ratio averaged 0.85 at the inflow station and 0.36 at the outflow station. A water budget shows unexplained losses to average 20 percent. At the Barton Creek Square Shopping Center, discharge-weighted densities of fecal-col iform and fecal-streptococci bacteria and discharge-weighted concentrations of biochemical oxygen demand, chemical oxygen demand, total organic carbon, suspended solids, total ammonia plus organic nitrogen, and total phosphorus generally were larger in the inflow than in the outflow. Discharge-weighted concentrations of dissolved lead, dissolved iron, and dissolved zinc generally were small in both the inflow and outflow; however, the larger discharge-weighted concentrations of these constituents generally were found in the inflow. Discharge-weighted concentrations of volatile dissolved solids were smaller in the inflow than in the outflow for 10 of the 22 storms analyzed. Discharge-weighted concentrations of total nitrite plus nitrate nitrogen and dissolved solids generally were much smaller in the inflow than in the outflow. It is likely that organic and ammonia nitrogen trapped in the pond from previous storms and in the inflow water as it flows through the pond is being oxidized to nitrite and nitrate nitrogen. Similarly, dissolved solids retained in the filter or on the bed of the pond from previous storms are being leached to the outflow. Measured peak concentrations or densities of most constituents in the inflow were significantly larger than those in the outflow for most constituents. An exception was noted for concentrations of total nitrite plus nitrate which were larger in the outflow than the inflow as indicated by discrete sample analysis for six storms. Loads of most constituents and total numbers of bacteria were significantly larger in the inflow than in the outflow. The total numbers of bacteria were reduced by approximately 80 percent. Average removal efficiencies for suspended solids, biochemical oxygen demand, total phosphorus, total organic carbon, chemical oxygen demand, and dissolved zinc ranged between 60 and 80 percent. The average loads of dissolved solids were approximately 13 percent larger in the outflow than the inflow. Average loads of total nitrite plus nitrate nitrogen were approximately 110 percent larger in the outflow than in the inflow. The increase in loads of these constituents is due to material being leached from the bed of the pond or from the filter system. At Al ta Vista, rainfall for the storms analyzed ranged from 0.25 to 2.00 inches. The maximum rainfall intensity was 0.30 inch for a 5-minute interval. The runoff-rainfall ratio averaged 0.42 and appeared to be evenly distributed about the mean ratio line. The peak discharge at the inflow station to the grass-covered swale area was 0.93 cubic foot per second. Inaccuracies of discharge at the outflow station and variations in the ungaged drainage area with the size of the storm prevented a hydrologic analysis of the basin above this station. Discharge-weighted concentrations of total phosphorus were larger in the outflow than in the inflow for each of the 19 storms analyzed. Discharge-weighted concentrations of dissolved solids, volatile dissolved solids, biochemical oxygen demand, chemical oxygen demand, and total organic carbon were larger in the outflow than in the inflow for at least 12 of the 19 storms analyzed. Discharge-weighted densities of fecal streptococci were decreased between the inflow and outflow, with discharge-weighted densities of fecal streptococci being less in the outflow for 15 of the 19 storms analyzed. Because of the relatively small variations in concentrations and densities of constituents between the inflow and outflow sites, and because of the errors in discharge at the outflow gage, it is not feasible to determine the effect of the grass-covered swales on discharge-weighted concentrations and densities of water-quality constituents. Discrete concentrations or densities of most constituents were not decreased. Peak concentrations of dissolved solids in the outflow exceeded peak concentrations in the inflow for all five of the storms analyzed with discrete samples. Peak concentrations of suspended solids, total ammonia plus organic nitrogen, total nitrite plus nitrate nitrogen, total nitrogen, and dissolved iron were larger in the outflow than in the inflow for four of the five storms analyzed. Load-removal efficiencies of water-quality constituents could not be determined because of inaccuracies in measuring discharge at the outflow site,

Texas

Seismic refraction profile, Kingdom of Saudi Arabia: field operations, instrumentation, and initial results

In February 1978 a seismic deep-refraction profile was recorded by the USGS along a 1000-km line across the Arabian Shield in western Saudi Arabia. The line begins in Paleozoic and Mesozoic cover rocks near Riyadh on the Arabian Platform, leads southwesterly across three major Precambrian tectonic provinces, traverses Cenozoic rocks of the coastal plain near Jizan (Tihamat Asir), and terminates at the outer edge of the Farasan Bank in the southern Red Sea. More than 500 surveyed recording sites were occupied, including 19 in the Farasan Islands. Six shot points were used--five on land, with charges placed mostly below water table in drill holes, and one at sea, with charges placed on the sea floor and fired from a ship. The total charge consumed was slightly in excess of 61 metric tons in 21 discrete firings. Seismic energy was recorded by means of a set of 100 newly developed portable seismic stations. Each station consists of a standard 2-Hz vertical geophone coupled to a self-contained analog recording instrument equipped with a magnetic-tape cassette. The stations were deployed in groups of 20 by five observer teams, each generally consisting of two scientist-technicians and a surveyor-guide. On the day prior to deployment, the instruments were calibrated and programmed for automatic operation by means of a specially designed device called a hand-held tester. At each of ten pre-selected recording time windows on a designated firing day, the instruments were programmed to turn on, stabilize, record internal calibration signals, record the seismic signals at three levels of amplification, and then deactivate. After the final window in the firing sequence, all instruments were retrieved and their data tapes removed for processing. A specially designed, field tape- dubbing system was utilized at shot point camps to organize and edit data recorded on the cassette tapes. The main functions of this system are to concatenate all data from each shot on any given day onto a single shot tape, and to provide hard copy for monitoring recorder performance so that any problems can be corrected prior to the next deployment. Composite digital record sections were produced from the dubbed tapes for each shot point by a portable processing and plotting system. The heart of this system is a DEC PDP 11VO3 computer, which controls a cassette playback unit identical to those used in the recorders and dubbers, a set of discriminators, a time-code translator, a digitizer, and a digital plotter. The system was used to maintain various informational data sets and to produce tabulations and listings of various sorts during the field operations, in addition to its main task of producing digital record sections. Two master clocks, both set to time signals broadcast by the British Broadcasting Corporation, provided absolute time for the recording operations. One was located on the ship and the other was stationed at a base camp on the mainland. The land-based master clock was used to set three additional master clocks located at the other active shot points a few days in advance of each firing, and these clocks were then used to set the internal clocks in the portable seismic stations via the hand-held tester. A master clock signal was also linked to the firing system at each shot point for determination of the absolute shot instant. It is possible to construct a generalized crustal model from examination of the six shot point composite record sections obtained in the field. Such a model rests upon a number of simplifying assumptions and will almost certainly be modified at a later stage of interpretation. The main assumptions are that the crust consists of two homogeneous isotropic layers having no velocity inversion,, that the Mohorovicic discontinuity is sharp, and that effects of surface inhomogeneities and elevation changes can be ignored. The main characteristics of the tentative model are the following: (1) The thickness of th

Open-File Report

Multi-resource analysis: A proof of concept study of natural resource tradeoffs in the Piceance Basin, Colorado, using the net resources assessment (NetRA) decision support tool

Executive Summary The U.S. Geological Survey (USGS) is developing a multi-resource analysis (MRA) line of products to inform land-use decision makers. Specifically, MRA products will integrate scientific information, include considerations for natural resource interrelations, and quantify the effects of resource management decisions in biophysical, economic, and societal terms. As part of the establishment of the MRA, the USGS, in collaboration with the University of New Mexico, has developed the Net Resources Assessment (NetRA) decision support tool. As a proof of concept analysis, the NetRA was applied to the Piceance basin in Colorado in a hypothetical example to illustrate how resource managers could use the NetRA to consider tradeoffs of natural resources among alternative development plans and land cover patterns within a geographic region. The NetRA is a policy-relevant approach to assess the availability of multiple natural resources. It is an analytical toolset that may be used to examine the spatiotemporal relations between development of energy and mineral resources and delivery of biological natural resources. The NetRA operates at multiple map scales and contains a set of integrated, compatible submodels with specific data requirements for natural resource stocks, engineering economics, biophysical, and ecological data for ecosystem services stocks, market prices, regulations, and nonmarket values. The NetRA includes an explicit process to consider the interdependence between development and conservation, which is a crucial consideration in land-management and land-use decisions. The NetRA is used to estimate an expected net resource value (NRV). The NRV is the expected, present value, economic benefit from the extraction of a resource (for example, natural gas) minus the total cost of production, which is the aggregation of the development, production, and social costs. Social costs include private costs plus any external costs. There can be external social benefits associated with natural gas production, such as increased demand for locally produced goods and increased employment in the local area through backward and forward linkages of natural gas production. The NRV is used to compare development outcomes (scenarios) from a range of exploration and development plans for cumulative energy production. The Piceance basin application of the NetRA uses the NRV to assess the tradeoff between continuous natural gas extraction and the effects to the local populations of Odocoileus hemionus (mule deer) and aquatic species and to consumptive water uses for an area the size and resolution of a USGS energy resource assessment unit. In the proof of concept simulation, the 2.9-square-mile-area of USGS oil and gas assessment unit 50200263 (Piceance basin continuous gas unit of the Mesaverde Total Petroleum System) was gridded into 588 cells. From this area, seven clusters with potential for development and three that cannot be developed were identified; the three clusters that cannot be developed were identified as wilderness study areas, areas of critical environmental concern, and national forests. On the basis of these criteria, there are 118 cells unsuitable for development in the oil and gas assessment unit: 84 are in national forests, 23 are areas of critical environmental concern, and 11 are wilderness study areas. The remaining cells in the oil and gas assessment unit can be developed on both private and public lands. Two scenarios were considered that are distinguished as plan 1 and plan 2. Plan 1 keeps the amount of land disturbance unchanged and limits the number of development locations to 140 grid cells for the production period, which constrains the amount of the energy resources available for development; the plan requires the usage of the Bureau of Land Management (BLM) unsuitability criteria. Plan 2 also limits the number of development locations to 140 grid cells for the production period but provides a constant volume of energy production by increasing the density of well pads within the cells. The effects of plan 2 to the NRV when there are five wells per pad and five pads per square mile happen mostly in the first 5 years of development, even though the effects on the population of mule deer continue in later years. This outcome is the result of the upfront development and investment costs and the initial effect to the ecosystem services.

Colorado

Hydrogeology, simulation of regional ground-water flow, and saltwater intrusion, Potomac-Raritan-Magothy Aquifer System, Northern Coastal Plain of New Jersey

The Potomac-Raritan-Magothy aquifer system in Middlesex and Monmouth Counties in the northern Coastal Plain of New Jersey consists primarily of unconsolidated Cretaceous sediments, which are divided into the upper and middle aquifers and confining units. These units, which strike northeastsouthwest along the Fall Line, dip and thicken to the southeast. The upper aquifer consists primarily of the Old Bridge Sand Member of the Magothy Formation, which is composed of coarse-grained sands, localized thin clay beds, and younger surficial sands and gravels in and near the outcrop. Transmissivity ranges from 1,760 to 19,400 ft 2 /d (feet squared per day) and tends to be higher in updip areas. Estimated withdrawals from the upper aquifer in the northern Coastal Plain were approximately 42 Mgal/d (million gallons per day) in 1986. Cones of depression whose centers range from 36 to 42 ft (feet) below sea level have developed as a result of these withdrawals. The upper aquifer is confined throughout most of the northern New Jersey Coastal Plain by clays and silts of the Cretaceous Woodbury Clay and Merchantville Formation and younger sediments of the Magothy Formation. This confining unit generally is greater than 200 ft thick. The simulated vertical hydraulic conductivity for the confining unit ranges from 8.4 x 10 -5 to 5.6 x 10 -3 feet per day; interpreted vertical hydraulic conductivities generally are lower except in southwestern Middlesex County, where the vertical hydraulic conductivities of the confining unit are higher. The middle aquifer consists primarily of the Farrington Sand Member of the Cretaceous Raritan Formation and surficial Holocene and Miocene sands and gravels in its outcrop area. It also can include the uppermost sands of the Cretaceous Potomac Group in parts of Monmouth County. The middle aquifer is composed of fine to coarse sand that contains some lignite and pyrite, and, locally, some clay beds. It pinches out in the northern part of Sayreville Township, near Raritan River. The transmissivity of the aquifer ranges from 2,140 to 13,800 ft 2/ d and tends to decrease in the northern part of the northern Coastal Plain of New Jersey where the aquifer thins. A poorly permeable confining unit composed mostly of clays and silts of the Woodbridge Clay Member of the Raritan Formation overlies the aquifer in most of this area. The confining unit generally is greater than 100 ft thick, although it thins and is sandy in the southwestern part of Middlesex County, where a good hydraulic connection exists between the middle and upper aquifers. Estimated withdrawals from the middle aquifer in the northern Coastal Plain were about 22 Mgal/d in 1986. These withdrawals have caused cones of depression whose centers range from 77 to 93 ft below sea level. A finite-difference, quasi-three-dimensional ground-water flow model was developed to simulate ground-water flow in the aquifer system. The confined and unconfined areas of the upper and middle aquifers were modeled as separate layers. The model was calibrated primarily by adjusting vertical hydraulic conductivity in the confining units and horizontal hydraulic conductivity in the aquifers, then matching simulated and measured groundwater levels for the period 1896-1986 and simulated and interpreted potentiometric surfaces under predevelopment conditions and in 1984. For the predevelopment period, the total flow into and out of the upper and middle aquifers is 35 and 21 Mgal/d, respectively. Recharge to the aquifer system is from direct recharge in the unconfined areas and from vertical leakage through overlying confining units. The main recharge areas are the topographically high areas in southwestern Middlesex County for both aquifers, in the eastern Sayreville area for the upper aquifer, and north of the Raritan River for the middle aquifer. Most ground water discharges to low-lying regional surface-water drains (streams), which flow into the South River. For 1984 transient conditions, the total ground-water flow into and out of the upper and middle aquifers is 61 and 34 Mgal/d, respectively. The largest amount of recharge is from direct recharge in the unconfined areas, but some recharge also is derived from vertical leakage through the Merchantville-Woodbury confining unit, captured ground-water discharge to streams, and induced inflow at artificial-recharge facilities. Regional flow is from recharge areas toward major cones of depression. Sensitivity analysis showed that the model was useful for representing flow in the system, especially in the confined-aquifer areas. Model representation of lateral and vertical boundary conditions was judged acceptable. Simulation results were less sensitive to changes in aquifer properties in the unconfined areas of the aquifers and to changes in storage in the confining units. Sensitivity analysis and calibration of hydraulic parameters and conditions showed that the distribution of hydraulic head was sensitive to changes in horizontal hydraulic conductivity in the aquifers, vertical hydraulic conductivity in the confining units, magnitudes of ground-water withdrawals, and initial hydraulic head in aquifer outcrop areas. Two scenarios were simulated to determine the effects of ground-water withdrawals from 1986 through 2019. For the scenario in which ground-water withdrawals increase to about 69 Mgal/d in the upper aquifer and 37 Mgal/d in the middle aquifer, centers of cones of depression are as deep as 100 ft below sea level in the upper aquifer and 170 ft below sea level in the middle aquifer. For this scenario, most of the additional water comes from captured surface-water discharge, induced cross-formational flow from overlying aquifers, and increases in induced flow from artificial-recharge areas. Induced flow from Raritan Bay also increases. For the scenario in which ground water withdrawals are reduced to 42.5 Mgal/d in the upper aquifer and 15 Mgal/d in the middle aquifer, water levels recover to above sea level nearly everywhere. In each aquifer, ground-water discharge to streams increases and induced flow through the confining units and from the overlying sediments decreases, and discharge of ground water to Raritan Bay in the upper aquifer exceeds the induced recharge from Raritan Bay. Reversal of ground-water gradients has caused saltwater intrusion in the two aquifers. Chloride concentrations in water from the upper aquifer in Keyport and Union Beach Boroughs were as high as 2,100 mg/L (milligrams per liter) in 1986. The intrusion has not increased significantly since well fields in the area were closed in the late 1970's. Elevated chloride concentrations also were measured in Keanesburg Borough in 1986. In both of these areas, saltwater has entered the upper aquifer from the Bay because of movement of the freshwater-saltwater interface in response to increasing ground-water withdrawals. Chloride concentrations in well-water samples from the middle aquifer were as high as 6,000 mg/L in Sayreville Borough in 1987; concentrations in samples from drive-point wells from the same aquifer near the Washington Canal, the main source of saltwater, were as high as 7,100 mg/L. The migration of the saltwater front at about 470 feet per year to the southeast is influenced mainly by a thinning of the middle aquifer, which constrains flow, and by the locations of regional cones of depression caused by groundwater withdrawals.

New Jersey

Hydrogeology of the unsaturated zone, North Ramp area of the Exploratory Studies Facility, Yucca Mountain, Nevada

Yucca Mountain, in southern Nevada, is being investigated by the U.S. Department of Energy as a potential site for a repository for high-level radioactive waste. This report documents the results of surface-based geologic, pneumatic, hydrologic, and geochemical studies conducted during 1992 to 1996 by the U.S. Geological Survey in the vicinity of the North Ramp of the Exploratory Studies Facility (ESF) that are pertinent to understanding multiphase fluid flow within the deep unsaturated zone. Detailed stratigraphic and structural characteristics of the study area provided the hydrogeologic framework for these investigations. Multiple lines of evidence indicate that gas flow and liquid flow within the welded tuffs of the unsaturated zone occur primarily through fractures. Fracture densities are highest in the Tiva Canyon welded (TCw) and Topopah Spring welded (TSw) hydrogeologic units. Although fracture density is much lower in the intervening nonwelded and bedded tuffs of the Paintbrush nonwelded hydrogeologic unit (PTn), pneumatic and aqueous-phase isotopic evidence indicates that substantial secondary permeability is present locally in the PTn, especially in the vicinity of faults. Borehole air-injection tests indicate that bulk air-permeability ranges from 3.5x10 -14 to 5.4x10 -11 square meters for the welded tuffs and from 1.2x10 -13 to 3.0x10 -12 square meters for the non welded and bedded tuffs of the PTn. Analyses of in-situ pneumatic-pressure data from monitored boreholes produced estimates of bulk permeability that were comparable to those determined from the air-injection tests. In many cases, both sets of estimates are two to three orders of magnitude larger than estimates based on laboratory analyses of unfractured core samples. The in-situ pneumatic-pressure records also indicate that the unsaturated-zone pneumatic system consists of four subsystems that coincide with the four major hydrogeologic units of the unsaturated zone at Yucca Mountain. In descending order, these hydrogeologic units are the Tiva Canyon welded (TCw), Paintbrush nonwelded (PTn), Topopah Spring welded (TSw ), and Calico Hills nonwelded (CHn). Deep percolation takes place as episodic pulses of inflow that propagate rapidly to depth and apparently bypass most of the rock matrix. Field-scale and core-scale water potentials throughout much of the PTn and TSw are very high, generally greater than -0.3 megapascals, and are nearly depth invariant. Thus, the imbibition capacity of the densely welded tuffs, at least near fractures, is very small because of low matrix permeabilities and low water-potential gradients across the fracture-matrix interface. The combination of high fracture permeability, high water potentials, high matrix saturations, and low matrix permeabilities results in a percolation environment that favors deep fracture flow. The episodic pulses of inflow are evidenced in the sporadic but nevertheless commonplace occurrence of water with concentrations of radioactive isotopes indicative of origins postdating the atmospheric testing of nuclear weapons. High concentrations of tritium have been detected at many horizons within the PTn and in the top of the TSw. Much lower concentrations of tritium, indicating the mixing of a bomb-pulse component with older water, have been detected in the deeper sections of the TSw and in the CHn. Evidence for fracture flow also is apparent in the widespread occurrence of perched water with chemical and isotopic signatures that indicate a fracture-flow origin for at least some of this water. In the North Ramp area, perched water has been detected at the base of the Topopah Spring Tuff or in the top of the underlying non welded to partially welded tuffs of the Calico Hills Formation in every dry-drilled borehole of sufficient depth to penetrate the Topopah Spring Tuff-Calico Hills Formation contact. The concentrations of the major ions of the perched water are similar to that of TSw pore water at borehole UZ-14, CHn pore water, and saturated-zone water at boreholes NRG-7 a and SD-9. The absolute chloride concentration of the perched water, however, is much lower than the chloride concentration of pore water from either the PTn or the TSw. The chemical and isotopic compositions of perched water indicate that this water was derived primarily from fracture flow, with little or no contribution from water in the matrix of the overlying rock. Carbon-14 ages of perched water range from 3,000 to 7,000 years. Strontium-87 isotope ratios indicate dissolution of surficial pedogenic calcite and calcite fracture fillings, which supports a fracture-flow origin for perched water. Moreover, carbon-13 and deuterium isotope values indicate rapid infiltration into fractures with little or no prior evaporation. Evidence for deep fracture flow into the Calico Hills Formation at UZ-14 is indicated by carbon-14 values that are from 65 and 95 percent modem carbon, equivalent to apparent ages of about 3,500 to 500 years. Some of these ages are younger than age estimates for perched water in the overlying Topopah Spring Tuff and are much younger than any that could be derived from a matrix-flow model. Evidence is lacking for extensive lateral flow within the PTn or for interception and diversion of this flow downward along structural pathways (faults), two key features of the original conceptual model for unsaturated flow at Yucca Mountain. Where data are available to infer lateral flow in the PTn, it is not certain that fracture flow could not have produced the same results. Pneumatic data, derived primarily from analysis of the interference effects from excavation of the North Ramp tunnel, indicate that faults within the Topopah Spring Tuff are open over substantial distances and are very permeable. Tunnel-boring-induced pneumatic disturbances have been propagated along these faults over distances that exceed 500 meters. These disturbances also have been detected in the pneumatic-pressure record of the overlying PTn in the vicinity of these faults. In spite of the apparent high permeability of faults, the existing data have neither confirmed nor refuted the hypothetical role of these faults in intercepting lateral flow from within or from above the PTn and diverting this flow downward into the deeper subsurface. On the basis of measured temperature gradients within the TSw, deep percolation appears to be greatest beneath active channels of major drainages, diminishing toward the margins and hillslopes bordering these channels. Numerical simulations indicate that this downward percolation is accompanied by lateral spreading as the percolation front moves downward through the PTn and across the contact between the PTn and underlying TSw. Temperature data from a well-documented site in Pagany Wash indicate the presence of a significant heat-flow deficit between the PTn and underlying TSw that most likely is due to nonconductive heat-flow processes with substantial capacity to extract heat. Percolation fluxes on the order of 10 to 20 millimeters per year beneath the Pagany Wash channel and on the order of 5 millimeters per year or less beneath the hillslopes bordering this drainage accounted for the apparent heat-flow deficit. Analyses of borehole temperature gradients in Drill Hole Wash indicate similar percolation fluxes and flux distributions within that drainage. An analysis of residence times estimated from uncorrected carbon-14 activities of perched-water samples and estimates for the volume of the structurally controlled reservoir, however, showed that the perched-water reservoir intersected by borehole UZ-14 under Drill Hole Wash could be sustained by percolation fluxes through the TSw of as little as 0.001 to 0.29 millimeter per year. The significance and implications of these findings with respect to waste isolation are discussed in the appendix of this report.

Nevada

A field guide to the central, creeping section of the San Andreas fault and the San Andreas Fault Observatory at Depth

This field trip is along the central section of the San Andreas fault and consists of eight stops that illustrate surface evidence of faulting, in general, and features associated with active fault creep, in particular. Fault creep is slippage along a fault that occurs either in association with small-magnitude earthquakes or without any associated large-magnitude earthquakes. Another aspect of the trip is to highlight where there are multiple fault traces along this section of the San Andreas fault zone in order to gain a better understanding of plate-boundary processes. The first stop is along the Calaveras fault, part of the San Andreas fault system, at a location where evidence of active fault creep is abundant and readily accessible. The stops that follow are along the San Andreas fault and at convenient locations to present and discuss rock types juxtaposed across the fault that have been transported tens to hundreds of kilometers by right-lateral motion along the San Andreas fault. Stops 6 and 7 are examples of recent studies of different aspects of the fault: drilling into the fault at the depth of repeating magnitude (M) 2 earthquakes with the San Andreas Fault Observatory at Depth (SAFOD) and the geological, geophysical, and seismological study of M 6 earthquakes near the town of Parkfield. Along with the eight official stops on this field trip are 12 “rolling stops”—sites of geologic interest that add to the understanding of features and processes in the creeping section of the fault. Many of the rolling stops are located where stopping is difficult to dangerous; some of these sites are not appropriate for large vehicles (buses) or groups; some sites are not appropriate for people at all. We include photographs of or from many of these sites to add to the reader's experience without adding too many stops or hazards to the trip. An extensive set of literature is available for those interested in the San Andreas fault or in the creeping section, in particular. For more scientifically oriented overviews of the fault, see Wallace ( 1990 ) and Irwin ( 1990 ); for a more generalized overview with abundant, colorful illustrations, see Collier ( 1999 ). Although the presence of small sections of the San Andreas fault was known before the great 1906 San Francisco earthquake, it was only after that event and subsequent geologic investigations reported in Lawson ( 1908 ) that showed the fault as a long structure, extending all the way from east of Los Angeles into northern California. Prentice ( 1999 ) described the importance of the 1908 “Lawson report” and how it pivotally influenced the understanding of the San Andreas. Hill ( 1981 ) presented a wonderful introduction to the evolution of thought on the San Andreas. Geologic maps and maps of the most recently active fault trace in the creeping section, or large parts of it, include those by Brown ( 1970 ), Dibblee ( 1971 , 1980 ), and Wagner et al. ( 2002 ); detailed geologic maps are discussed at various stops in this guide. Various aspects of the creeping section of the San Andreas fault have been the focus of many geologic field trips in the past few decades. Guidebooks for some of those trips include those by Gribi ( 1963a , 1963b ), Brabb et al. ( 1966 ), Rogers ( 1969 ), Bucknam and Haller ( 1989 ), Harden et al. ( 2001 ), and Stoffer ( 2005 ). The creeping section of the San Andreas fault zone lies between areas that experienced large-displacement surface breakage during great earthquakes in 1857 and 1906 (Fig. 1 inset). Burford and Harsh ( 1980 ) divided the creeping section into three segments: (1) a northwest section where the creep rate increases to the southeast in step-like increments, (2) a central section where the creep rate is relatively constant at a maximum value of ∼30 mm/yr (∼1.2 in/yr), and (3) a southeast section where the creep rate decreases to the southeast (Fig. 2 ). The rate of slip along the creeping section of the fault zone has been measured using creepmeters, alignment arrays, and laser distance-measuring devices. The aperture of measurements over which these measurements are made ranges from 10 m (∼33 ft) (creepmeters) to 100 m (∼330 ft) (alignment arrays) to kilometers and tens of kilometers (laser measuring devices). Creepmeter and alignment-array measurements are here termed “near-fault” measurements; laser measurements over distances of 1–2 km (∼0.6–1.2 mi) are termed “intermediate-scale” measurements; laser measurements over tens of kilometers (miles) are termed “broadscale” measurements. Comparisons among near-fault, intermediate-scale, and broadscale measurements and geologic maps show that the northwest part of the creeping section of the fault is composed of two narrow zones of active deformation, one along the San Andreas fault and one along the Calaveras-Paicines fault, whereas the central and southeast sections are both composed of a single relatively narrow zone of deformation. The southeast section is transitional to a locked zone southeast of Cholame; a locked fault is one that slips only in association with a moderate to large earthquake. Throughout the creeping section of the San Andreas fault zone, broadscale measurements generally indicate more deformation than near-fault and intermediate-scale measurements, which are in reasonably close agreement except at Monarch Peak (Mustang Ridge), near the center of the creeping section and our Stop 5 ((Figs. 1 ) and 2 ). Features that we see on this trip include offset street curbs, closed depressions (sag ponds), fault scarps (steep slopes formed by movement along a fault), a split and displaced tree, offset fence lines, fresh fractures, and offset road lines (Fig. 3 is a sketch showing some of the landforms that represent deformation by an active fault). We also see evidence of long-term maturity of the San Andreas fault, as indicated by fault features and displaced rock types (Fig. 4 ). Finally, we will visit sites of ongoing research into the processes associated with earthquakes and their effects. Discussions include drilling into the San Andreas fault at the SAFOD drill site and the 2004 Parkfield earthquake and its effects and implications.

California

Geophysics- and geochemistry-based assessment of the geochemical characteristics and groundwater-flow system of the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system in Doña Ana County, New Mexico, and El Paso County, Texas, 2010–12

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.” Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system. The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin. The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits. Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the south end of the midbasin uplift. Groundwater from the Aden Hills and the East and West Potrillo Mountains flows to the south end of the midbasin uplift and then continues east towards the Paso del Norte. Throughout most of the Mesilla Valley, the vertical hydraulic gradient was downward because the water-level altitude in the Rio Grande alluvium was higher than it was in the Santa Fe Group, but in some areas (typically in the middle and southern parts of the Mesilla Valley), the vertical hydraulic gradient was substantially reduced or even reversed to an upward hydraulic gradient. The geochemistry data indicate that there was a complex system of multiple geochemical endmembers and mixing between these endmembers with recharge to the Rio Grande alluvium and Santa Fe Group composed mostly of seepage from the Rio Grande, inflows from deeper or neighboring water systems, and mountain-front recharge. Five distinct geochemical groups were identified in the Mesilla Basin study area: (1) ancestral Rio Grande (pre-Pleistocene) geochemical group, (2) modern Rio Grande (Pleistocene to present) geochemical group, (3) mountain-front geochemical group, (4) deep groundwater upwelling geochemical group, and (5) unknown freshwater geochemical group. The ancestral Rio Grande groundwater was water that recharged into the system as seepage losses from the ancestral Rio Grande; this groundwater generally flows from north to south-southeast towards the Paso del Norte. Groundwater on the west side of the midbasin uplift generally flows south until it reaches the southern part of the study area; from the southern part of the study area, the groundwater flows east towards the Paso del Norte. Groundwater on the east side of the midbasin uplift flows south-southeast towards the Paso del Norte where it mixes with groundwater from the modern Rio Grande, uplifted areas in the west, and the deep saline source. The water type of the modern Rio Grande geochemical group ranged from calcium-sulfate water type in the northern part of the study area to sodium-chloride-sulfate water type in the southern part of the study area; from north to south there was a substantial increase in specific conductance, strontium-87/strontium-86 ratio, potassium, and the trace metals of iron and lithium, changing the water chemistry such that it became similar to the water chemistry of the deep groundwater upwelling geochemical group. From age-dating results, water in the modern Rio Grande geochemical group was recharged to the Rio Grande alluvium within the past 10 years. The mountain-front geochemical group was generally old water (apparent age was greater than 10,000 carbon-14 years before present) that was somewhat mineralized and has relatively high concentrations of fluoride and silica, which might indicate longer exposure to volcanic and siliciclastic rocks or aluminosilicate minerals. There were five different locations of recharge determined from the groundwater geochemistry within the mountain-front geochemical group, all having a slightly different geochemical signature: (1) the Rough and Ready Hills, Robledo Mountains, and the Sleeping Lady Hills, (2) the Doña Ana Mountains, (3) the Aden Hills and West Potrillo Mountains, (4) the East Potrillo Mountains, and (5) the Sierra Juárez in Mexico. The groundwater from the Rough and Ready Hills, Robledo Mountains, the Sleeping Lady Hills, and the Doña Ana Mountains generally flows toward the Rio Grande and eventually mixes together and with the modern Rio Grande groundwater. The groundwater originating from the Aden Hills and East and West Potrillo Mountains generally flows east to southeast at a slow rate and eventually mixes and continues east, where it mixes with groundwater from the ancestral Rio Grande geochemical group and with the groundwater from the Sierra Juárez. The groundwater from the Sierra Juárez flows north and then east towards the Paso del Norte where it mixes with groundwater from the uplifted areas in the west, ancestral and modern Rio Grande groundwater, and the upwelling groundwater from a deep saline source. The deep groundwater upwelling geochemical group had the highest concentrations of bicarbonate, potassium, silica, aluminum, iron, and lithium within the study area, indicating that it had been in contact with carbonate and siliciclastic rocks for a much longer period of time and at higher temperatures compared to the other geochemical groups, and was most likely ancient marine groundwater originating from the Paleozoic and Cretaceous carbonate rocks which was upwelling into the Mesilla Basin aquifer system in the southeastern part of the study area through the extensive fault systems. Direct-current resistivity and time-domain electromagnetic soundings support the interpretation of ancient marine groundwater upwelling into the Mesilla Basin aquifer system, as do the analytical results from wells, and the helicopter frequency domain electromagnetic data collected along the Rio Grande. The hydrogen-2/hydrogen-1 ratio and oxygen-18/oxygen-16 ratio isotopic results for samples in the unknown freshwater geochemical group did not plot on the Rio Grande evaporation line, indicating this group did not have a Rio Grande signature (that is, there was no isotopic evidence of a component of Rio Grande water) and it also had the lowest mineralized content of any geochemical group in the study area.

New Mexico, Texas

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida