Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Natural Resource Modeling”

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

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

At least 1,675 records · Page 93Linked to original sources

Mineral potential for sediment-hosted copper deposits in the Islamic Republic of Mauritania (phase V, deliverable 75)

The presence of Neoproterozoic through Cambrian, continental, siliciclastic sedimentary rocks interbedded with dolomitic carbonates, shales, and glacial tillites similar to the Katanga Supergroup host rocks of the Central African Copperbelt and other sediment-hosted copper-bearing Proterozoic sequences worldwide, is first order criteria for consideration of the Neoproterozoic units of the Taoudeni Basin in Mauritania as prospective for sediment-hosted copper deposits. Review of the National Mineral Occurrences Database (Marsh and Anderson, 2015) and previous literature suggest that only a handful of small sediment-hosted copper occurrences have been found to date in Mauritania and that the resource potential for this deposit type is low. In the northern Taoudeni Basin, the most important occurrence is at Taradent. This occurrence consists of three mineralized horizons in the lower Neoproterozoic Char Group in three outcrop areas separated by alluvium over a strike length of 12 kilometers (km). The most extensively mineralized horizon consists of malachite and disseminated copper sulfides, and is concentrated at the base of a dolomitic interval, consistent with a reduced faciestype sediment-hosted copper deposit model. Additional and poorly described copper occurrences in the Taoudeni Basin margin sedimentary rocks in northeastern Mauritania, such as Chegga Guettatira and Sidi Bara, may be sediment-hosted copper occurrences and extend the potential throughout this portion of the Basin. Additional potential for sediment-hosted copper deposits is indicated in reports by the British Geological Survey (BGS; Gunn and others, 2004) and BRGM (Salpeteur, 2005) in the Neoproterozoic-Cambrian sedimentary rocks (Teniagouri Group) of the central and southern Taoudeni Basin at a number of occurrences in the vicinity of Fouges, Weringuel, and Gig Zig. These occurrences exhibit characteristics consistent with a reduced-facies-type sediment-hosted copper deposit model and based on correlation of similar host rocks extend the potential throughout the entire Taoudeni Basin margin. Very poorly described copper occurrences in the northeastern and southeastern part of the Basin at Kreb en Naga, Kreb en Naga SW, Berbere, Dhar Nema, Dhar Nema 2, Saile, and Nejam-Medroume all are described as occurring in sedimentary rocks and are hosted by or are in close proximity to voluminous outcrops of Jurassic microgabbros (also referred to as dolerite or diabase). Uncertainty exists as to whether they are sediment-hosted copper occurrences. They are likely more closely analogous to copper occurrences associated with mafic igneous activity. The determination of whether reported copper occurrences and associated sedimentary rocks at Kreb en Naga and Kreb en Naga SW in the northeastern Taoudeni Basin and the Dhar Nema occurrences in the southeastern Taoudeni are permissive of sediment-hosted copper mineralization was inconclusive due to the absence of copper showings. While sedimentary host rocks at these locations must be broadly considered permissive, the absence of other favorable geologic features in these areas indicates that the potential for mineral deposits of this type is low. The widespread presence of mafic sills and dikes of Jurassic age throughout these areas and descriptions of associated minor copper mineralization suggests that many of the copper occurrences in these areas are not of the sediment-hosted copper type. Highly speculative potential exists in the Coastal Basin for Revett-type copper deposits in unexposed Jurassic-Cretaceous sandstones that may contain structural traps with gas accumulations. Eocene continental red bed sandstones of the Gorgol Formation, where overlain by the shallow marine sedimentary rocks of the Rinndiao Formation, represent an additional low potential rock sequence within the southeastern portion of the Coastal Basin that is permissive of sediment-hosted copper deposits. Field examinations of outcrop areas in the vicinity of many of the known occurrences were performed in late 2007 to provide additional criteria for evaluating the potential for this type of mineralization in Mauritania. In general, the geographic coordinates of occurrences are imprecise and many of the occurrences could not be located. Examination of outcrops of the Char Group in the Taradent area confirms the permissive nature of these rocks for reduced-facies-type sediment-hosted copper deposits. Inspection of weak copper mineralization hosted by black shales at several occurrences in the Gadel-Maylime area confirms that these occurrences are also of the reduced-faciestype and that the black and green shales of the Bouly and Ould Yenje Formations are permissive. Reported copper analyses up to 0.75 percent copper at Taradent and a single analysis of mineralized shale in the Gadel-Maylime area containing 0.55 percent copper indicate that low grade ore material exists at these locations. However indications that such copper grades are continuous over significant intervals and present in large tonnages are lacking. Although mineral occurrence data and descriptive geological information are adequate to delineate areas favorable for sediment-hosted copper deposits, this review indicates that potential for this type of deposit in Mauritania is low.

Open-File Report↗

Residence times and nitrate transport in ground water discharging to streams in the Chesapeake Bay Watershed

One of the major water-quality problems in the Chesapeake Bay is an overabundance of nutrients from the streams and rivers that discharge to the Bay. Some of these nutrients are from nonpoint sources such as atmospheric deposition, agricultural manure and fertilizer, and septic systems. The effects of efforts to control nonpoint sources, however, can be difficult to quantify because of the lag time between changes at the land surface and the response in the base-flow (ground water) component of streams. To help resource managers understand the lag time between implementation of management practices and subsequent response in the nutrient concentrations in the base-flow component of streamflow, a study of ground-water discharge, residence time, and nitrate transport in springs throughout the Chesapeake Bay Watershed and in four smaller watersheds in selected hydrogeomorphic regions (HGMRs) was conducted. The four watersheds were in the Coastal Plain Uplands, Piedmont crystalline, Valley and Ridge carbonate, and Valley and Ridge siliciclastic HGMRs. A study of springs to estimate an apparent age of the ground water was based on analyses for concentrations of chlorofluorocarbons in water samples collected from 48 springs in the Chesapeake Bay Watershed. Results of the analysis indicate that median age for all the samples was 10 years, with the 25th percentile having an age of 7 years and the 75th percentile having an age of 13 years. Although the number of samples collected in each HGMR was limited, there did not appear to be distinct differences in the ages between the HGMRs. The ranges were similar between the major HGMRs above the Fall Line (modern to about 50 years), with only two HGMRs of small geographic extent (Piedmont carbonate and Mesozoic Lowland) having ranges of modern to about 10 years. The median values of all the HGMRs ranged from 7 to 11 years. Not enough samples were collected in the Coastal Plain for comparison. Spring samples showed slightly younger water under wet conditions than under dry conditions. The apparent age of water from wells, springs, and other ground-water discharge points in the four targeted watersheds was modern to 60 years, which was similar to the apparent ages from the spring study. In the Pocomoke River Watershed in the Coastal Plain Uplands HGMR, the apparent age of ground-water samples ranged from 0 to 60 years; the ages in the vicinity of the streams ranged from 0 to 23 years. The apparent ages of ground water in the Polecat Creek Watershed in the Piedmont crystalline HGMR ranged from 2 to 30 years. The apparent ages of water from wells in the Muddy Creek Watershed in the Valley and Ridge carbonate HGMR ranged from 10 to 20 years (except for a single sample that was 45 years). The ages in the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR ranged from 0 to 50 years. The distribution in apparent age of water from wells in the targeted watersheds, however, generally is older than that for water from the springs. The median age of water from wells in the Muddy Creek Watershed, for example, was 15 years, compared to 11 years for the water from the springs in that watershed, and less than 10 years for water from all springs in the spring study. The similarity in the ranges in apparent age of water from the wells and from the springs shows that the samples from the targeted watersheds and springs have bracketed the range of apparent ages that would be expected in the shallow ground-water-flow systems throughout the Chesapeake Bay Watershed. The apparent age of water from individual wells does not necessarily represent the entire distribution of ages of the discharging ground water, and it is this distribution of ages that affects the response of nutrient concentrations in stream base flow. Nutrient-reduction scenarios were modeled for two watersheds for which the distribution of apparent ground-water ages was available, the East Mahantango Creek Watershed in the Valley and Ridge siliciclastic HGMR and the Locust Grove Watershed in the Coastal Plain Uplands HGMR. A nutrient-reduction scenario was created for East Mahantango Creek, where the average residence time was determined to be approximately 10 years on the basis of the output of particle tracking from a ground-water-flow model. This scenario showed decreases of nearly 50 percent in base-flow concentrations of nitrate in streams within the first year after the reduction in nitrogen input; smaller reductions in nitrate concentration occurred in each subsequent year. A second scenario for that same watershed, in which the same 10-year average residence time was assumed and an exponential model was used for analysis, showed that a 50-percent reduction in base-flow concentrations of nitrate could take up to 5 years. For the Locust Grove Watershed, in which an average residence time of 32 years was assumed, simulation with the exponential model showed that it may take more than 20 years to achieve a 50-percent reduction in base-flow concentra-tions of nitrate. Although it was not possible to construct such scenarios for all watersheds, these examples show the range of possible responses to changes in nutrient inputs in two very different types of watersheds. Findings from this study include information on factors that affect ground-water age, spatial distribution of ages, and nitrogen transport. In the East Mahantango Creek Watershed and the Polecat Creek Watershed, the residence time varied spatially depending on the position of the flow path, and temporally depending on the recharge conditions. Generally, ground water in areas near the stream had short residence times and the water in upland areas had longer residence times. Water traveling through deep layers had longer residence times than water traveling through shallow layers, and residence times were faster under high recharge conditions than low recharge conditions. Ground water in the Pocomoke Watershed exhibits a similar pattern: younger water discharges to small order streams in headwater basins and older water discharges to larger streams near the basin outlet. Factors affecting nitrogen transport in ground water include spatial and temporal variation in input sources, ground-water age, and aquifer processes that lead to denitrification. Spatial and temporal variations in nitrogen sources affect all the watersheds. Tributaries with higher inputs of nitrogen have higher concentrations in stream base flow. Areas where nitrogen application rates have increased over time show an age-nitrate relation in ground-water samples. The age-nitrate relation can be affected by denitrification, which occurs in Pocomoke and East Mahantango Creeks but is not evident in Polecat and Muddy Creeks. In East Mahantango Creek, the level of denitrification is significant in water with residence times greater than 20 years, but because this is a small component of overall ground-water discharge to a stream, it may not remove a significant quantity of nitrogen from the system. Denitrification in Pocomoke Creek is significant and appears to affect mostly older water discharging to streams. Therefore, if most of the nitrogen entering these two streams is associated with the discharge of younger ground water, denitrification may not greatly affect the overall nitrogen delivery to these streams. Other findings of this study show that nitrate in ground water discharging along preferential flow paths may not be affected by natural processes, such as denitrification or uptake by riparian vegetation. Seeps to swales and ditches beneath the north uplands at Polecat Creek indicate a shallow water table and discharge of young ground water whereas the absence of such seeps on the south side indicates a deep water table and a lack of young ground water. Similarly, discharge at the base of the slope and to the valley wetland south of the creek but not north of the creek indicates a different role for the riparian forest on the two sides of the creek. In many of the systems where water discharges at the base of slopes to wetlands, ditches have been dug to drain the valley. Such drainage circumvents possible removal of nitrate by riparian vegetation. Because ground-water residence times do not appear directly related to the HGMRs, the targeting of management practices will achieve the most rapid response in water quality if directed at 1) watersheds with large agricultural sources of nitrate, 2) areas with the shortest ground-water-flow paths and 3) areas not affected by significant denitrification. The fastest response in stream base-flow concentrations of nitrogen to implementation of management practices would be to implement practices in those areas with the highest loads rather than attempt to target practices on the basis of HGMR stratification. Overall findings of the study indicate that 1) ground-water contributions to nitrogen in streamflow are significant, 2) some response to management practices should be evident in base-flow concentrations of nitrogen and loads within 1 to 5 years in watersheds with the shortest average residence times, but response time may be closer to 20 years in watersheds with longer average ground-water residence times, 3) the majority of the response in ground-water discharge to any changes in management practices will be distributed over a 10-year time period even in the watersheds with the fastest response times, and 4) given that half the streamflow is from ground-water discharge and the other half is runoff or soil water, about 90 percent of total water being discharged to a stream will be less than about a decade old; therefore, full implementation of nutrient reductions may result in improved streamwater quality in about a decade. In the more-likely scenario of gradual source reduction, the reduction in concentrations of nitrate in streams and aquifers would take longer than the examples shown here.

Chesapeake Bay watershed↗

Geophysical studies of the West Antarctic rift system

The West Antarctic rift system extends over a 3000 × 750 km, largely ice covered area from the Ross Sea to the base of the Antarctic Peninsula, comparable in area to the Basin and Range and the East African rift system. A spectacular rift shoulder scarp along which peaks reach 4–5 km maximum elevation marks one flank and extends from northern Victoria Land-Queen Maud Mountains to the Ellsworth-Whitmore-Horlick Mountains. The rift shoulder has maximum present physiographic relief of 5 km in the Ross Embayment and 7 km in the Ellsworth Mountains-Byrd Subglacial Basin area. The Transantarctic Mountains part of the rift shoulder (and probably the entire shoulder) has been interpreted as rising since about 60 Ma, at episodic rates of ∼1 km/m.y., most recently since mid-Pliocene time, rather than continuously at the mean rate of 100 m/m.y. The rift system is characterized by bimodal alkaline volcanic rocks ranging from at least Oligocene to the present. These are exposed asymmetrically along the rift flanks and at the south end of the Antarctic Peninsula. The trend of the Jurassic tholeiites (Ferrar dolerites, Kirkpatric basalts) marking the Jurassic Transantarctic rift is coincident with exposures of the late Cenozoic volcanic rocks along the section of the Transantarctic Mountains from northern Victoria Land to the Horlick Mountains. The Cenozoic rift shoulder diverges here from the Jurassic tholeiite trend, and the tholeiites are exposed continuously (including the Dufek intrusion) along the lower- elevation (1–2 km) section of Transantarctic Mountains to the Weddell Sea. Widely spaced aeromagnetic profiles in West Antarctica indicate the absence of Cenozoic volcanic rocks in the ice covered part of the Whitmore-Ellsworth-Mountain block and suggest their widespread occurrence beneath the western part of the ice sheet overlying the Byrd Subglacial Basin. A German Federal Institute for Geosciences and Natural Resources (BGR)-U.S. Geological Survey (USGS) aeromagnetic survey over the Ross Sea continental shelf indicates rift fabric and suggests numerous submarine volcanoes along discrete NNW trending zones. A Bouguer anomaly range of approximately 200 (+50 to −150) mGal having 4–7 mGal/km gradients where measured in places marks the rift shoulder from northern Victoria Land possibly to the Ellsworth Mountains (where data are too sparse to determine maximum amplitude and gradient). The steepest gravity gradients across the rift shoulder require high density (mafic or ultramafic?) rock within the crust as well as at least 12 km of thinner crust beneath the West Antarctic rift system in contrast to East Antarctica. Sparse land seismic data reported along the rift shoulder, where velocities are greater than 7 km/s, and marine data indicating velocities above 7 km/s beneath the Ross Sea continental shelf support this interpretation. The maximum Bouguer gravity range in the Pensacola Mountains area of the Transantarctic Mountains is only about 130 mGal with a maximum 2 mGal/km gradient, which can be explained solely by 8 km of crustal thickening. Large offset seismic profiles over the Ross Sea shelf collected by the German Antarctic North Victoria Land Expedition V (GANOVEX V) combined with earlier USGS and other results indicate 17–21 km thickness for the crust beneath the Ross Sea shelf which we interpret as evidence of extended rifted continental crust. A regional positive Bouguer anomaly (0 to +50 mGal), the width of the rift, extends from the Ross Sea continental shelf throughout the Ross Embayment and Byrd Subglacial Basin area of the West Antarctic rift system and indicates that the Moho is approximately 20 km deep tied to the seismic results (probably coincident with the top of the asthenosphere) rather than the 30 km reported in earlier interpretations. The interpretation of horst and graben structures in the Ross Sea, made from marine seismic reflection data, probably can be extended throughout the rift (i.e., the Ross Ice shelf and the Byrd Subglacial Basin areas). The near absence of earthquakes in the West Antarctic rift system probably results from a combination of primarily sparse seismograph coverage and, secondarily, suppression of earthquakes by the ice sheet (e.g., Johnston, 1987) and very high seismicity shortly after deglaciation in the Ross Embayment followed by abnormally low seismicity at present (e.g., Muir Wood, 1989). The evidence of high temperatures at shallow depth beneath the Ross Sea continental shelf and adjacent Transantarctic Mountains is supportive of thermal uplift of the mountains associated with lateral heat conduction from the rift and can possibly also explain the volcanism, rifting, and high elevation of the entire rift shoulder to the Ellsworth-Horlick-Whitmore Mountains. We infer that the Gondwana breakup and the West Antarctic rift are part of a continuously propagating rift that started in the Jurassic when Africa separated from East Antarctica (including the failed Jurassic Transantarctic rift). Rifting proceeded clockwise around East Antarctica to the separation of New Zealand and the Campbell Plateau about 85–95 Ma and has continued (with a spreading center jump) to its present location in the Ross Embayment and West Antarctica. The Cenozoic activity of the West Antarctic rift system appears to be continuous in time with rifting in the same area that began only in the late Mesozoic. Although the mechanism for rifting is not completely explained, we suggest a combination of the flexural rigidity model (Stem and ten Brink, 1989) proposed for the Ross Embayment and the thermal plume or hot spot concepts. The propagating rift may have been “captured” by the thermal plume.

Tectonics↗

Hydrologic analysis of the Mojave River, California, using a mathematical model

The channel of the Mojave Rive'r in California is normally dry and is highly permeable over much of its length, and large quantities of water from natural floodflows in the channel infiltrate through the channel bed to the underlying ground-water body. From 1930 to 1972 only 18 floods at The Forks produced flow at Barstow, 55 miles (88 kilometres) downstream from The Forks. Peak discharges at Barstow from these floods ranged from 180 cubic feet per second (5 cubic metres per second) in 1967 to 64,300 cubic feet per second (1,820 cubic metres per second) in 1938. Total stream infiltration, primarily as ground-water recharge, ranged from 3,600 acre-feet (4.40 cubic hectometres) in 1935 to 50,400 acre-feet (62.1 cubic hectometres) in 1969 between The Forks and Victorville and from about 7,000 acre-feet (8.61 cubic hectometres) in 1935 to 128,000 acre-feet (158 cubic hectometres) in 1969 between Victorville and Barstow. The Mojave Water Agency is considering the use of the channel of the river to convey water imported from northern California through Silverwood Reservoir (5 miles or 8 kilometres upstream from The Forks) downstream to the Barstow area. The imported water would be used to replenish aquifers underlying the Barstow area. A mathematical model was developed that simulates the advance of discharge down the initially dry channel of the Mojave River, and the model was used to evaluate the potential of the channel to move imported water downstream to Barstow. Results of simulation by modeling indicate that the channel of the Mojave River can be used to efficiently convey imported water to Barstow only when the absorption capacity of the channel has been reduced by an antecedent flood. The volume of imported water that can reach Barstow depends on the volume and duration of the antecedent flood, on the volume of imported water released trom Silverwood Reservoir, and on the rate at which imported water is released. A release of 20,000 acre-feet (24.6 cubic hectometres) of imported water may produce at Barstow a maximum volume of imported water of 2,500, 8,000, 11,000, or 15,000 acre-feet (3.08, 9.86, 13.6, or 18.4 cubic hectometres) for a release rate of 500, 750, 1,000, or 2,000 cubic feet per second (14.2, 21.2, 28.3, or 56.8 cubic metres per second). For planning purposes in evaluating some of the hydraulic effects of recharge on the aquifer, a simulation of the aquifer near the Barstow area using an electrical analog model showed that a combination of no pumping and a yearly recharge of 5,000 acre-feet (6.17 cubic hectometres) for 10 years could raise ground-water levels at least 10 feet (3 metres) over an area of about 10 square miles (25 square kilometres). To obtain the water-level changes due to the combined effects of pumping and recharge, the above water-level changes should be superimposed on the separate effects of pumping.

California↗

Regional evaluation of evapotranspiration in the Everglades

Nine sites in the Florida Everglades were selected and instrumented for collection of data necessary for evapotranspiration-determination using the Bowen-ratio energy-budget method. The sites were selected to represent the sawgrass or cattail marshes, wet prairie, and open-water areas that constitute most of the natural Everglades system. At each site, measurements necessary for evapotranspiration (ET) calculation and modeling were automatically made and stored on-site at 15- or 30-minute intervals. Data collected included air temperature and humidity at two heights, wind speed and direction, incoming solar radiation, net solar radiation, water level and temperature, soil moisture content, soil temperature, soil heat flux, and rainfall. Data summarized in this report were collected from January 1996 through December 1997, and the development of site-specific and regional models of ET for this period is described. Latent heat flux is the energy flux density equivalent of the ET rate. Modified Priestley-Taylor models of latent heat flux as a function of selected independent variables were developed at each site. These models were used to fill in periods of missing latent heat flux measurement, and to develop regional models of the entire Everglades region. The regional models may be used to estimate ET in wet prairie, sawgrass or cattail marsh, and open-water portions of the natural Everglades system. The models are not applicable to forested areas or to the brackish areas adjacent to Florida Bay. Two types of regional models were developed. One type of model uses measurements of available energy at a site, together with incoming solar energy and water depth, to estimate hourly ET. This available-energy model requires site data for net radiation, water heat storage, and soil heat flux, as well as data for incoming solar radiation and water depth. The other type of model requires only incoming solar energy, air temperature, and water depth data to provide estimates of hourly ET. The second model thus uses data that are more readily available than the data required for the available-energy model. Computed ET mean annual totals for all nine sites for the 1996-97 period ranged from 42.4 inches per year at a site where the water level is below land surface for several months each year to 57.4 inches per year at an open-water site with no emergent vegetation. Although the density of photosynthetically-active plant leaves has been shown to relate directly to ET in some studies, it does not appear to relate directly to ET in the Everglades, based on comparison of annual ET data with leaf-area index, defined as the Normalized Difference Vegetation Index (NDVI), data from satellite imagery. NDVI and ET appear to be inversely related in the Everglades. The greatest ET rates occurred at open-water sites where the NDVI data indicated the lowest leaf-area index. Among the remaining vegetated sites, there is no clear relation between ET and NDVI, though the highest ET rate corresponded to the lowest NDVI and one of the lowest ET rates corresponded to the highest NDVI value. The variation in ET follows a seasonal pattern, with lowest monthly ET totals occurring in December through February, and highest ET occurring in May through August. The greatest range in monthly ET among all nine sites for the 2-year period occurred at site 3: from 1.81 inches in December 1997 to 6.84 inches in July 1996. A study to compare the Bowen-ratio/energy balance method of ET measurement with the eddy-correlation method was done at one site from June 22, 1998, through September 28, 1998. This comparison indicated that both methods gave comparable values of the Bowen ratio, but there was a considerable difference in available energy measured by the two methods. The mean of all 30-minute measured turbulent heat fluxes from the eddy-correlation apparatus for June 22 through September 29, 1998, was 137.4 watts per square meter, and the mean of the corresponding measured energy was 163.6 watts per square meter, or about 20 percent greater. The disagreement in mean energy fluxes measured by the two methods is problematical and is not fully understood. Although the difference seems to be related to friction velocity, and is practically non-existent at values of friction velocity greater than 0.3 meter per second, the "correctness" of either method cannot be determined with the data available.

Florida↗

Relations between hydrology, water quality, and taste-and-odor causing organisms and compounds in Lake Houston, Texas, April 2006-September 2008

Lake Houston is a surface-water-supply reservoir and an important recreational resource for the city of Houston, Texas. Growing concerns over water quality in Lake Houston prompted a detailed assessment of water quality in the reservoir. The assessment focused on water-quality constituents that affect the aesthetic quality of drinking water. The hydrologic and water-quality conditions influencing the occurrence of taste-and-odor causing organisms and compounds in Lake Houston were assessed using discrete and continuously monitored water-quality data collected during April 2006– September 2008. The hydrology of Lake Houston is characterized by rapidly changing conditions. During inflow events, water residence time can change by orders of magnitude within a matter of hours. Likewise, the reservoir can stratify and destratify over a period of several hours, even during non-summer and at relatively short water residence times, given extended periods with warm temperatures and little wind. The rapidly changing hydrology likely influences all other aspects of water quality in Lake Houston, including the occurrence of taste-and-odor causing organisms and compounds. Water quality in Lake Houston varied with respect to season and water residence time but typically was indicative of turbid, eutrophic to hypereutrophic conditions. In general, turbidity and nutrient concentrations were largest during non-summer (October–May) and when water residence times were relatively short (less than 100 days), which reflects the influence of inflow events on water-quality conditions. Large inflow events can cause substantial changes in water-quality conditions over relatively short periods of time (hours). The taste-and-odor causing organisms cyanobacteria and actinomycetes bacteria were always present in Lake Houston. Cyanobacterial biovolume was largest during summer (June– September) and when water residence time was greater than 100 days. Annual maxima in cyanobacterial biovolume occurred during July-September of each year, when temperatures were larger than 27 degrees Celsius and water residence times were longer than 400 days. In contrast, actinomycetes bacteria were most abundant during non-summer and when water residence times were less than 100 days, reflecting the close association between these organisms and transport of suspended sediments. Geosmin and 2-methylisoborneol are the taste-and-odor causing compounds most commonly produced by cyanobacteria and actinomycetes bacteria. Geosmin was detected more frequently (62 percent of samples) than 2-methylisoborneol (29 percent of samples) in Lake Houston. Geosmin exceeded the human detection threshold (10 nanograms per liter) only once during the study period and 2-methylisoborneol exceeded the human detection threshold twice. Manganese is a naturally occurring trace element that can occasionally cause taste-andodor problems in drinking water. Manganese concentrations exceeded the human detection threshold (about 50 micrograms per liter) in about 50 percent of samples collected near the surface and 84 percent of samples collected near the bottom. The cyanotoxin microcystin was detected relatively infrequently (16 percent of samples) and at small concentrations (less than or equal to 0.2 micrograms per liter). The abundance of the taste-and-odor causing organisms cyanobacteria and actinomycetes bacteria in Lake Houston was coupled with inflow events and subsequent changes in water-quality conditions. Cyanobacterial biovolume (biomass) in Lake Houston was largest during warm periods with little inflow and relatively small turbidity values. In contrast, actinomycetes bacteria were most abundant following inflow events when turbidity was relatively large. Severe taste-and-odor problems were not observed during the study period, precluding quantification of the hydrologic and water-quality conditions associated with large concentrations of taste-and-odor causing compounds and development of predictive models. Reservoir inflow (water residence time) and turbidity, variables related to the abundance of potential taste-andodor causing organisms, are currently (2011) continuously measured in Lake Houston, and predictive models could be developed in the future when the hydrologic and water-quality conditions associated with taste-and-odor problems have been better quantified. Seasonal and water residence time influences on water-quality conditions altered relations between hydrologic and water-quality conditions and taste-and-odor causing organisms and compounds. Future data collection and development of predictive models need to account for the variability associated with season and water residence time.

Texas↗

Analysis of steady-state flow and advective transport in the eastern Snake River Plain aquifer system, Idaho

Quantitative estimates of ground-water flow directions and traveltimes for advective flow were developed for the regional aquifer system of the eastern Snake River Plain, Idaho. The work included: (1) descriptions of compartments in the aquifer that function as intermediate and regional flow systems, (2) descriptions of pathlines for flow originating at or near the water table, and (3) quantitative estimates of traveltimes for advective transport originating at or near the water table. A particle-tracking postprocessing program was used to compute pathlines on the basis of output from an existing three-dimensional steady-state flow model. The flow model uses 1980 conditions to approximate average annual conditions for 1950-80. The advective transport model required additional information about the nature of flow across model boundaries, aquifer thickness, and porosity. Porosity of two types of basalt strata has been reported for more than 1,500 individual cores from test holes, wells, and outcrops near the south side of the Idaho National Engineering Laboratory. The central 80 percent of samples had porosities of 0.08 to 0.25, the central 50 percent of samples, O. 11 to 0.21. Calibration of the model involved choosing a value for porosity that yielded the best solution. Two radiologic contaminants, iodine-129 and tritium, both introduced to the flow system about 40 years ago, are relatively conservative tracers. Iodine- 129 was considered to be more useful because of a lower analytical detection limit, longer half-life, and longer flow path. The calibration value for porosity was 0.21. Most flow in the aquifer is contained within a regional-scale compartment and follows paths that discharge to the Snake River downstream from Milner Dam. Two intermediate-scale compartments exist along the southeast side of the aquifer and near Mud Lake.One intermediate-scale compartment along the southeast side of the aquifer discharges to the Snake River near American Fails Reservoir and covers an area of nearly 1,000 square miles. This compartment, which receives recharge from an area of intensive surface-water irrigation, is apparently fairly stable. The other intermediate-scale compartment near Mud Lake covers an area of 300 square miles. The stability and size of this compartment are uncertain, but are assumed to be in a state of change. Traveltimes for advective flow from the water table to discharge points in the regional compartment ranged from 12 to 350 years for 80 percent of the particles; in the intermediate-scale flow compartment near American Falls Reservoir, from 7 to 60 years for 80 percent of the particles; and in the intermediate-scale compartment near Mud Lake, from 25 to 100 years for 80 percent of the particles. Traveltimes are sensitive to porosity and assumptions regarding the importance of the strength of internal sinks, which represent ground-water pumpage. A decrease in porosity results in shorter traveltimes but not a uniform decrease in traveltime, because the porosity and thickness is different in each model layer. Most flow was horizontal and occurred in the top 500 feet of the aquifer. An important limitation of the model is the assumption of steady-state flow. The most recent trend in the flow system has been a decrease in recharge since 1987 because of an extended drought and changes in land use. A decrease in flow through the system will result in longer traveltimes than those predicted for a greater flow. Because the interpretation of the model was limited to flow on a larger scale, and did not consider individual wells or well fields, the interpretations were not seriously limited by the discretization of well discharge. The interpretations made from this model also were limited by the discretization of the major discharge areas. Near discharge areas, pathlines might not be representative at the resolution of the grid. Most improvement in the estimates of ground-waterflow directions and travelt

Water-Resources Investigations Report↗

Simulated effects of projected pumping on the availability of freshwater in the Evangeline Aquifer in an area southwest of Corpus Christi, Texas

This study is an investigation of the continued availability of freshwater in the Evangeline aquifer along the Texas Gulf Coast and the potential for degradation of the water quality by salinewater intrusion. Recharge to the aquifer occurs by the infiltration of precipitation in the outcrop area and by cross-formational flow from deeper aquifers. The predevelopment recharge rate is about 6 to 8 cubic feet per second. The predevelopment flow is toward the coast. The flow is semiconfined in the outcrop area and confined underneath the Chicot aquifer in the eastern two-thirds of the study area. Discharge, under natural conditions, is upward into the Chicot aquifer and to the Nueces River or Gulf of Mexico. Intensive pumping by irrigators, industries, and municipalities over the last 80 years has created a cone of depression as deep as 219 feet below sea level under the city of Kingsville in Kleberg County. The total rate of pumpage in 1982 was 29.6 cubic feet per second. A mathematical model of the flow and water quality in the Evangel ine aquifer was developed using available data to simulate the historical effect of pumping on the potentiometric surface and water quality, and to simulate the effect of projected pumping on the potentiometric surface and water quality to the year 2020. The water quality in the aquifer is only marginally suitable for drinking water. The chloride concentration before development in the 1930's and 1940's, ranged from 9 to 1,971 milligrams per liter. The mean chloride concentration was 353 (standard deviation 262) milligrams per liter. The potential sources of water-quality degradation on a regional scale are: Salinewater intrusion from under the Gulf of Mexico; movement of poor quality water within outlying sections of the aquifer; and downward leakage from the overlying Chicot aquifer. Leakage from the Chicot is the most likely to cause serious regional water-quality degradation. Other local potential sources of contamination are: Leaky well casings, oil-field brine disposal, water movement along faults, and in-situ uranium mining. These sources might create some local water-quality degradation. The results of the historical period simulation indicate, as do current field data, that little or no significant deterioration has occurred in the water quality of the Evangeline aquifer. The simulations and the sensitivity tests of the aquifer properties, conditions, and assumptions indicate that vertical conductivity of the Chicot aquifer is the most sensitive and least well known part of the system. The storage coefficient of the Evangeline aquifer and the aggregate thickness of high-conductivity sand layers within the aquifer as well as the vertical distribution of these layers are also important properties that are not well known. Two simulations of the projected pumping a low estimate, as much as 46.2 cubic feet per second during 2011-20; and a high estimate, as much as 60.0 cubic feet per second during the same period indicate that no further regional water-quality deterioration is likely to occur. Many important properties and conditions are estimated from poor or insufficient field data, and possible ranges of these properties and conditions are tested. In spite of the errors and data deficiencies, the results are based on the best estimates currently available. The reliability of the conclusions rests on the adequacy of the data and the demonstrated sensitivity of the model results to errors in estimates of these properties.

Texas↗

Hydrogeology of the area near the J4 test cell, Arnold Air Force Base, Tennessee

The U.S. Air Force operates a major aerospace systems testing facility at Arnold Engineering Development Center (AEDC) in Coffee County, Tennessee. Dewatering operations at one of the test facilities, the J4 test cell, has affected the local ground-water hydrology. The J4 test cell is approximately 100 feet in diameter, extends approximately 250 feet below land surface, and penetrates several aquifers. Ground water is pumped continuously from around the test cell to keep the cell structurally intact. Because of the test cell's depth, dewatering has depressed water levels in the aquifers surrounding the site. The depressions that have developed exhibit anisotropy that is controlled by zones of high permeability in the aquifers. Additionally, contaminants - predominately volatile organic compounds - are present in the ground-water discharge from the test cell and in ground water at several other Installation Restoration Program (IRP) sites within the AEDC facility. The dewatering activities at J4 are drawing these contaminants from the nearby sites. The effects of dewatering at the J4 test cell were investigated by studying the lithologic and hydraulic characteristics of the aquifers, investigating the anisotropy and zones of secondary permeability using geophysical techniques, mapping the potentiometric surfaces of the underlying aquifers, and developing a conceptual model of the ground-water-flow system local to the test cell. Contour maps of the potentiometric surfaces in the shallow, Manchester, and Fort Payne aquifers (collectively, part of the Highland Rim aquifer system) show anisotropic water-level depressions centered on the J4 test cell. This anisotropy is the result of features of high permeability such as chert-gravel zones in the regolith and fractures, joints, and bedding planes in the bedrock. The presence of these features of high permeability in the Manchester aquifer results in complex flow patterns in the Highland Rim aquifers near the J4 test cell. The occurrence, distribution, and orientation of these features has a great effect on ground-water flow to the J4 test cell. The depression caused by dewatering extends out horizontally through the aquifers along the most permeable pathways. Since the aquifers above the Chattanooga Shale are not separated by distinct confining units, areas in adjacent aquifers above and below these zones of high permeability in the Manchester aquifer are also dewatered. Conditions in all Highland Rim aquifers approximate steady-state equilibrium because ground-water withdrawal at the test cell has been continuous since the late 1960's. The average ground-water discharge from the dewatering system at the J4 test cell was 105 gallons per minute, for 1992-95. The ground-water capture areas in each aquifer extend into all or parts of landfill #2 and leaching pit #2 (IRP site 1), the main testing area (IRP site 7), and the old fire training area (IRP site 10). IRP sites 8 and 12 are outside the ground-water capture areas. Of the 35 sampled wells in the J4 area, 10 produced water samples containing chlorinated organic compounds such as 1,2-dichloroethane (1,2-DCA), 1,1-dichloroethylene (1,1-DCE), and trichloroethylene (TCE) in concentrations which exceeded the Tennessee Department of Environment and Conservation Maximum Contaminant Levels (MCL's) for public water-supply systems. The highest concentrations were detected in samples from well AEDC-274 with 45 micrograms per liter (mg/L) 1,2-DCA, 320 mg/L 1,1-DCE, and 1,200 mg/L TCE. These compounds are synthetic and do not occur naturally in the environment. A sample of the ground-water discharge from the J4 test cell also contained concentrations of these compounds that exceed MCL's. Chlorinated organic compounds, including 1,2-DCA; 1,1-DCE; and TCE also have been detected at IRP sites 1, 7, 8, nd 10. The six dewatering wells surrounding the J4 test cell penetrate the Chattanooga Shale and are open to the Highland Rim aquifer system, there

Water-Resources Investigations Report↗

Simulation of Ground-Water Flow and Effects of Ground-Water Irrigation on Base Flow in the Elkhorn and Loup River Basins, Nebraska

Irrigated agriculture is vital to the livelihood of communities in the Elkhorn and Loup River Basins in Nebraska, and ground water is used to irrigate most of the cropland. Concerns about the sustainability of ground-water and surface-water resources have prompted State and regional agencies to evaluate the cumulative effects of ground-water irrigation in this area. To facilitate understanding of the effects of ground-water irrigation, a numerical computer model was developed to simulate ground-water flow and assess the effects of ground-water irrigation (including ground-water withdrawals, hereinafter referred to as pumpage, and enhanced recharge) on stream base flow. The study area covers approximately 30,800 square miles, and includes the Elkhorn River Basin upstream from Norfolk, Nebraska, and the Loup River Basin upstream from Columbus, Nebraska. The water-table aquifer consists of Quaternary-age sands and gravels and Tertiary-age silts, sands, and gravels. The simulation was constructed using one layer with 2-mile by 2-mile cell size. Simulations were constructed to represent the ground-water system before 1940 and from 1940 through 2005, and to simulate hypothetical conditions from 2006 through 2045 or 2055. The first simulation represents steady-state conditions of the system before anthropogenic effects, and then simulates the effects of early surface-water development activities and recharge of water leaking from canals during 1895 to 1940. The first simulation ends at 1940 because before that time, very little pumpage for irrigation occurred, but after that time it became increasingly commonplace. The pre-1940 simulation was calibrated against measured water levels and estimated long-term base flow, and the 1940 through 2005 simulation was calibrated against measured water-level changes and estimated long-term base flow. The calibrated 1940 through 2005 simulation was used as the basis for analyzing hypothetical scenarios to evaluate the effects of ground-water irrigation on stream base flow for 1940 through 2005 and for 2006 through 2045. Simulated base flows were compared for scenarios that alternately did or did not include a representation of the effects of ground-water irrigation. The difference between simulated base flows for the two scenarios represents the predicted effects of ground-water irrigation on base flow. Comparison of base flows between simulations with ground-water irrigation and no ground-water irrigation indicated that ground-water irrigation has cumulatively reduced streamflows from 1940 through 2005 by 888,000 acre-feet in the Elkhorn River Basin and by 2,273,000 acre-feet in the Loup River Basin. Generally, predicted cumulative effects of ground-water irrigation on base flow were 5 to 10 times larger from 2006 through 2045 than from 1940 through 2005, and were 7,678,000 acre-feet for the Elkhorn River Basin and 14,784,000 acre-feet for the Loup River Basin. The calibrated simulation also was used to estimate base-flow depletion as a percentage of pumping volumes for a 50-year future time period, because base-flow depletion percentages are used to guide the placement of management boundaries in Nebraska. Mapped results of the base-flow depletion analysis conducted for most of the interior of the study area indicated that pumpage of one additional theoretical well simulated for a future 50-year period generally would result in more than 80 percent depletion when it was located close to the stream, except in areas where depletion was partly offset by reduced ground-water discharge to evapotranspiration in wetland areas. In many areas, depletion for the 50-year future period composed greater than 10 percent of the pumped water volume for theoretical wells placed less than 7 or 8 miles from the stream, though considerable variations existed because of the heterogeneity of the natural system represented in the simulation. For a few streams, predicted future simulated base flows dec

Scientific Investigations Report↗

Aquatic assemblages and their relation to temperature variables of least-disturbed streams in the Salmon River basin, central Idaho, 2001

In the late 1990s, Idaho’s established stream temperature criteria for the protection of coldwater biota and salmonid spawning were considered inadequate because the criteria did not agree with observed biological conditions in many instances and did not allow for variability in environmental condition or species diversity across a broad area such as the entire State of Idaho. In 2001, benthic invertebrate and fish assemblages in 34 least-disturbed streams in the Salmon River Basin, central Idaho, were evaluated in relation to stream temperature and other environmental variables. The Salmon River Basin retains watersheds that are minimally affected by human activities. These “natural” stream conditions provide a basis for deriving attainable stream temperatures that can be used to set new, and revise existing, water-quality criteria for stream habitats affected by human activities. During July through September 2001, data were collected to document the thermal regime of least-disturbed streams, characterize the distribution of aquatic biota at streams representing a gradient of temperature, and describe the relations between environmental variables and benthic invertebrate and fish assemblages. Nine stream temperature metrics were compared with the U.S. Environmental Protection Agency’s criterion of 10 ° C (degrees Celsius) for bull trout spawning and juvenile rearing. The maximum weekly-maximum temperature at all 33 sites where temperature data were available exceeded this criterion. The maximum daily-maximum temperature (MDMT) at 30 of the sites exceeded the Idaho Department of Environmental Quality (IDEQ) criterion of 13.0 ° C for the protection of salmonid spawning; and the maximum daily-average temperature at all 33 sites exceeded the 9.0 ° C criterion for the protection of salmonid spawning. The thermal regime at most sites did not exceed the IDEQ criteria for the protection of coldwater biota. Nine environmental variables—water-surface gradient, discharge, wetted channel width, width:depth ratio, aspect, total seasonal thermal input, open canopy, riparian canopy density, and elevation were selected for correlation with the nine stream temperature metrics. Elevation showed the strongest inverse correlation with the stream temperature metrics. Two hundred and one benthic invertebrate taxa from the 34 sampling sites were identified. The most abundant taxa were Baetis tricaudatus , Oligochaeta, Tvetenia bavarica gr., Acari, Rhithrogena , Cinygmula , Heterlimnius , Micropsectra , Eukiefferiella devonica gr., Drunella doddsi , and Cricotopus . Of the 201 benthic invertebrate taxa collected during this study, 57 taxa (present at a minimum of 5 sampling sites) were significantly correlated with one or more of the stream temperature metrics. Among the invertebrate taxa collected, 32 were identified as coldwater taxa. Of the coldwater taxa collected, Zapada oregonensis gr. showed the strongest inverse correlation with the stream temperature metrics and was collected at sites where maximum weekly-maximum temperature (based on date of sample and 6 days prior) ranged from 11.3 ° to 18.5 °C. Ten species of fish in the families Salmonidae, Cottidae, and Cyprinidae were collected. Two species (bull trout and chinook salmon) listed as threatened under the U.S. Fish and Wildlife Service Endangered Species Act were collected. Among all species, bull trout showed the strongest inverse correlation between relative fish abundance and stream temperature. Bull trout and juvenile bull trout densities were inversely correlated with stream temperature. The probability of occurrence of juvenile bull trout was significantly correlated with MDMT on the basis of results from a logistic regression model developed during this study. However, this model differed from a similar model developed by the U.S. Forest Service on the basis of regional data collected in the Pacific Northwest. The regression model based on data collected during this study showed higher probabilities of occurrence of bull trout at colder stream temperatures (10 ° to 15 ° C) and lower probabilities of occurrence at warmer stream temperatures (16 ° to 21 ° C) than did the model based on regional data. The model comparisons suggest that regional or local differences need to be considered when deriving stream temperature criteria.

Idaho↗

Simulated effects of ground-water augmentation on the hydrology of Round and Halfmoon Lakes in northwestern Hillsborough County, Florida

Pumpage from the Upper Floridan aquifer in northwest Hillsborough County near Tampa, Florida, has induced downward leakage from the overlying surficial aquifer and lowered the water table in many areas. Leakage is highest where the confining layer separating the aquifers is breached, which is common beneath many of the lakes in the study area. Leakage of water to the Upper Floridan aquifer has lowered the water level in many lakes and drained many wetlands. Ground water from the Upper Floridan aquifer has been added (augmented) to some lakes in an effort to maintain lake levels, but the resulting lake-water chemistry and lake leakage patterns are substantially different from those of natural lakes. Changes in lake-water chemistry can cause changes in lake flora, fauna, and lake sediment composition, and large volumes of lake leakage are suspected to enhance the formation of sinkholes near the shoreline of augmented lakes. The leakage rate of lake water through the surficial aquifer to the Upper Floridan aquifer was estimated in this study using ground-water-flow models developed for an augmented lake (Round Lake) and non-augmented lake (Halfmoon Lake). Flow models developed with MODFLOW were calibrated through nonlinear regression with UCODE to measured water levels and monthly net ground-water-flow rates from the lakes estimated from lake-water budgets. Monthly estimates of ground-water recharge were computed using an unsaturated flow model (LEACHM) that simulated daily changes in storage of water in the soil profile, thus estimating recharge as drainage to the water table. Aquifer properties in the Round Lake model were estimated through transient-state simulations using two sets of monthly recharge rates computed during July 1996 to February 1999, which spanned both average conditions (July 1996 through October 1997), and an El Niño event (November 1997 through September 1998) when the recharge rate doubled. Aquifer properties in the Halfmoon Lake model were estimated through steady-state simulations of average conditions in July 1996. Simulated hydrographs computed by the Round and Halfmoon Lake models closely matched measured water-level fluctuations, except during El Niño, when the Halfmoon Lake model was unable to accurately reproduce water levels. Possibly, potential recharge during El Niño was diverted through ground-water-flow outlets that were not represented in the Halfmoon Lake model, or a large part of the rainfall was diverted into runoff before it could become recharge. Solute transport simulations with MT3D indicate that leakage of lake water extended 250 to 400 feet into the surficial aquifer around Round Lake, and from 75 to 150 feet around Halfmoon Lake before flowing to the underlying Upper Floridan aquifer. These results are in agreement with concentrations of stable isotopes of oxygen-18 (d18O) and deuterium (dD) in the surficial aquifer. Schedules of monthly augmentation rates to maintain constant stages in Round and Halfmoon Lakes were computed using an equation that accounted for changes in the Upper Floridan aquifer head and the deviation from the mean recharge rate. Resulting lake stages were nearly constant during the first half of the study, but increased above target lake stages during El Niño; modifying the computation of augmentation rates to account for the higher recharge rate during El Niño resulted in lake stages that were closer to the target lake stage. Substantially more lake leakage flows to the Upper Floridan aquifer from Round Lake than from Halfmoon Lake, because the estimated vertical hydraulic conductivities of lake and confining layer sediments and breaches in the confining layer beneath Round Lake are much greater. Augmentation rates required to maintain the low guidance stages in Round Lake (53 feet) and Halfmoon Lake (42 feet) under average Upper Floridan aquifer heads are estimated as 33,850 cubic feet per day and 1,330 to 10,000 cubic feet per day, respectively. These rates equate to 26 inches per month of water applied to the entire surface of Round Lake and 0.34 to 2.5 inches per month of water applied to Halfmoon Lake.

Florida↗

Hydrogeologic setting and preliminary estimates of hydrologic components for Bull Run Lake and the Bull Run Lake drainage basin, Multnomah and Clackamas counties, Oregon

The hydrogeologic setting was described and preliminary estimates of hydrologic components prepared for the Bull Run Lake and for the Bull Run Lake drainage basin, in the Cascade Range of northwestern Oregon. The 0.73-square-mile lake and the 3.44-square-mile drainage basin lie within the Bull Run Watershed, the principal water supply for the Portland, Oregon, metropolitan area. During periods of high demand or low inflows to the watershed, the City of Portland, Bureau of Water Works, releases water from Bull Run Lake to augment the supply. Bull Run Lake is impounded by a natural dam formed by a landslide. Outflow of ground water from the lake through the landslide emerges as springflow at the toe of the landslide and forms the headwaters of the Bull Run River. The approximately 4,300-Mgal (million gallons) discharge of the Bull Run River measured below the springs during the 1993 water year is composed of (1) outflow of ground water from Bull Run Lake through the landslide (approximately 60 percent), (2) ground water originating from the contributing drainage area between the lake and the springs (approximately 34 percent), (3) streamflow from Bull Run Lake (approximately 5 percent), and (4) surface runoff (streamflow and overland flow) from the contributing drainage area between the lake and the springs (approximately 1 percent). Estimated ranges for inflows to the Bull Run Lake drainage basin during the 1993 water year were about 3,400 to 9,200 Mgal from precipitation from rain and snow, and about 0 to 3,300 Mgal from fog drip. Estimated ranges for outflows from the lake basin, listed from largest to smallest, were about 1,800 to 3,400 Mgal for ground-water outflow through the landslide; about 600 to 1,800 Mgal for evapotranspiration from the land surface; about 170 to 410 Mgal for lake evaporation; and about 0 to 400 Mgal for streamflow from the lake. Ground- water outflow through the consolidated rocks could not be evaluated owing to the lack of data. The lake storage increased by a range of from about 1,700 to 1,900 Mgal. Changes in ground-water storage and soil-moisture storage could not be evaluated as a result of insufficient data. Estimated inflows to Bull Run Lake from precipitation on the lake surface during the 1993 water year ranged from about 600 to 1,600 Mgal. Inflows from ground water and surface runoff could not be evaluated owing to the lack of data. Estimated ranges for outflows from the lake were about 1,800 to 3,400 Mgal from ground-wateroutflow through the landslide, about 170 to 410 Mgal from lake evaporation, and about 0 to 400 Mgal from streamflow. Outflow of ground water through the consolidated rocks could not be evaluated owing to the lack of data. Lake storage increased by a range of from about 1,700 to 1,900 Mgal. Suggestions for further study include (1) evaluation of the surface-runoff component of inflow to the lake; (2) use of a cross-sectional ground-water flow model to estimate ground-water inflow, outflow, and storage; (3) additional data collection to reduce the uncertainties of the hydrologic components that have large relative uncertainties; and (4) determination of long-term trends for a wide range of climatic and hydrologic conditions.

Oregon↗

Geochemical conditions and the occurrence of selected trace elements in groundwater basins used for public drinking-water supply, Desert and Basin and Range hydrogeologic provinces, 2006-11: California GAMA Priority Basin Project

The geochemical conditions, occurrence of selected trace elements, and processes controlling the occurrence of selected trace elements in groundwater were investigated in groundwater basins of the Desert and Basin and Range (DBR) hydrogeologic provinces in southeastern California as part of the Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA PBP is designed to provide an assessment of the quality of untreated (raw) groundwater in the aquifer systems that are used for public drinking-water supply. The GAMA PBP is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The DBR hydrogeologic provinces consist of 141 defined groundwater basins separated by mountain ranges, faults, and other features. This report presents analyses of data collected from nine study areas within the DBR hydrogeologic provinces: Antelope Valley, Borrego Valley, the Central Desert area, Coachella Valley, Colorado River, Indian Wells Valley, Low-Use Basins of the Mojave and Sonoran Deserts, the Mojave, and Owens Valley. Collectively, these nine study areas are referred to as the DBR study unit. The study unit covers approximately 7,000 square miles and includes the 63 groundwater basins in the DBR hydrogeologic provinces in which groundwater is used for public drinking-water supply. The vast majority of the 223 wells sampled for this study were long-screened production wells used primarily for public supply. Uncorrected carbon-14 ( 14 C) groundwater ages for samples collected in the DBR study unit ranged from less than (<) 100 to 33,700 years before present (BP). Sixty-six percent of sample ages were greater than (>) 100 years BP, and 40 percent were >3,800 years BP. Samples collected from wells located adjacent to mountain-front recharge areas or major surface-water features generally had younger groundwater ages than did samples collected from wells located away from mountain fronts or towards the distal ends of basin groundwater flow paths. Most groundwater sampled in the DBR study unit had alkaline pH: 89 percent of sample pH values ranged from 7.1 to 9.8, with 37 percent greater than or equal to (&ge;) 7.9. Groundwater age was significantly correlated (positively) with pH, likely because silicate weathering is a primary control on groundwater pH and is a slow process. The oxidation-reduction (redox) condition of the groundwater sampled in the DBR study unit was predominantly oxic (71 percent), except in the Colorado River study area where organic-rich fluvial aquifers provide the electron donors necessary to support iron-reducing (anoxic-Fe) redox processes. The cation type of 78 percent of the samples was either sodium- or mixed-type, and the anion type of 83 percent of the samples was either bicarbonate- or mixed-type. Sodium-type groundwaters generally were older and more alkaline than calcium-type groundwaters, consistent with the change in water chemistry expected from cation exchange between groundwater and aquifer sediments over long periods of time. Because of the correlation with young groundwater, calcium-type groundwater was predominantly from wells located adjacent to mountain-front recharge areas. Arsenic (As), boron (B), fluoride (F), molybdenum (Mo), strontium (Sr), uranium (U), and vanadium (V) were selected for assessment in this study because they occurred at concentrations greater than California Department of Public Health or U.S. Environmental Protection Agency regulatory or non-regulatory drinking-water-quality benchmarks in more than 2 percent of the 223 samples collected in the DBR study unit. As and F were detected most commonly (18 and 13 percent, respectively) at concentrations above associated water-quality benchmarks and Sr and V least frequently (both at 3 percent). Given that 14 C groundwater ages are predominantly >100 years BP, land use in the study unit is primarily undeveloped, and chemicals derived from anthropogenic sources, such as volatile organic compounds, were infrequently detected, high concentrations of these trace elements in groundwater were most likely the result of natural factors and not anthropogenic factors. As, F, Mo, and V concentrations showed significant positive correlations to groundwater age and to pH. This relation is partly due to the sources of trace elements likely being the weathering of primary minerals, such as silicate minerals, which is a slow process that takes place over hundreds to thousands of years. This relation also reflects the positive correlation between groundwater age and pH. Geochemical modeling predicted that the dominant species of As, Mo, and V in solution were oxyanions (HAsO 4 2&ndash; , MoO 4 2&ndash; , and H 2 VO 4&ndash; ), which are likely to be mobile in alkaline groundwater because mineral surfaces composing aquifer matrices have a predominantly negative surface charge under alkaline conditions. F also exists predominantly as a negatively charged ion (F &ndash; ). At pH values >7.5, saturation indices generated by the geochemical modeling program PHREEQC indicated that F solubility may be somewhat limited by the precipitation of the mineral fluorapatite [Ca 5 (PO 4 ) 3 F]. Speciation modeling of As in anoxic-Fe groundwater (iron-reducing conditions) showed that samples were supersaturated with orpiment (As 2 S 3 ), indicating that mineral precipitation may be responsible for low As concentrations observed in reducing groundwater. In contrast, U concentrations showed significant negative correlations to groundwater age and to pH. Higher U concentrations generally occurred in samples for which geochemical modeling indicated that the uncharged ternary complex Ca 2 UO 2 (CO 3 ) 3 was the dominant aqueous U species. This uncharged complex is not attracted to the charged surfaces of minerals and thus increases U solubility. Formation of Ca 2 UO 2 (CO 3 ) 3 was greater in younger groundwaters because calcium and uranium concentrations generally were lower in older groundwaters, likely due to cation-exchange processes and precipitation of the mineral calcite as groundwater pH increased. Co-precipitation of U with the calcite (CaCO 3 ) may remove U from the aqueous phase. Saturation indices indicated that the anoxic-Fe groundwaters from the Colorado River study area were supersaturated with the mineral uraninite (UO 2 ), suggesting that UO 2 precipitation may be responsible for the low concentrations of U observed in these samples. Concentrations of strontium, which exists primarily in a cationic form (Sr 2+ ), were not significantly correlated with either groundwater age or pH. Strontium concentrations showed a strong positive correlation with total dissolved solids (TDS). Dissolved constituents, such as Sr, that interact with mineral surfaces through outer-sphere complexation become increasingly soluble with increasing TDS concentrations of groundwater. Boron concentrations also showed a significant positive correlation with TDS, indicating the B may interact to a large degree with mineral surfaces through outer-sphere complexation.

California↗

Groundwater movement and interaction with surface water near the confluence of the Platte and Elkhorn rivers, Nebraska, 2016–18

The State of Nebraska requires a sustainable balance between long-term water supplies and uses of groundwater and surface water and requires Natural Resources Districts to include the effect of groundwater use on surface-water systems as part of their respective integrated management plans. Recent droughts in Nebraska (2000–6; 2012–13) have amplified concerns about the long-term sustainability of groundwater and surface-water resources in the state, and concerns about the effect of groundwater irrigation on both streamflow and the water supplies needed to meet wildlife, recreational, and municipal needs. The lower Platte River provides nearly 100 percent of drinking-water supplies to Lincoln, Nebraska, 40 to 60 percent of drinking-water supplies to Omaha, Nebr., and critical aquatic and riparian habitat for threatened and endangered species. The Lower Platte River Basin-wide Management Plan has been jointly developed by the Nebraska Department of Natural Resources and seven Natural Resources Districts to address some of these concerns by managing groundwater and surface-water resources conjunctively. To sustain flows in the lower Platte River that are needed for municipal water supplies, water managers have proposed projects aimed at temporary storage of surface water in upstream parts of the basin to mitigate periods of low flow in the lower Platte River. To increase scientific understanding and provide support for any potential future streamflow augmentation projects, the Papio-Missouri River Natural Resources District, the Lower Platte North Natural Resources District, and the Nebraska Department of Natural Resources, in cooperation with the U.S. Geological Survey, initiated this study to examine groundwater/surface-water interaction along the lower Platte and Elkhorn Rivers upstream from their confluence. The study design described herein focused on understanding seasonal characteristics of groundwater movement and interaction with surface water during periods of high groundwater demand (June through August) and low groundwater demand (all other months). Understanding how groundwater movement and interaction with surface water are affected by streamflow conditions and local groundwater demand is critical to the development of any streamflow augmentation project intended to sustain streamflow and mitigate periods of low flow in the lower Platte River. The characteristics of groundwater movement and interaction with surface water are affected by hydrologic and local climatic conditions. For the study area, 2016–18 conditions can be broadly characterized as above normal precipitation. The flows measured at the Elkhorn River at Waterloo, Nebr., streamflow-gaging station (U.S. Geological Survey station 06800500) were above the long-term median, and the streamflow of the Platte River near Leshara, Nebr., streamflow-gaging station (06796500) remained normal or slightly above normal for the duration of this study. Continuous streamflow and water-level data were interpreted to examine differences in groundwater movement and interaction with surface water between the Platte and Elkhorn Rivers during high and low groundwater demand periods. Although the streamflow for the Platte and Elkhorn Rivers and their tributaries was less during the high groundwater demand period, the hydraulic gradient along a transect of recorder wells was identical (0.0012 foot per foot) during the high and low groundwater demand synoptic water-level and streamflow surveys. The hydraulic gradient between the Platte and Elkhorn Rivers generally remained between 0.0011 and 0.0012 foot per foot. It can be inferred that the hydraulic gradient, which is the only temporally variable factor in Darcy’s Law, is consistent throughout the study period and that groundwater flow does not vary appreciably along this transect. The northern part of the study area (north of the transect of recorder wells) has consistent groundwater and tributary flow from Big Slough, Rawhide Creek (Old Channel), and Rawhide Creek for low and high groundwater demand periods. In the southern part of the study area (south of the transect of recorder wells), tributary flow is more variable and dependent on local groundwater demand and flow conditions of the Platte River. Small decreases (less than 2 feet) in the groundwater levels, such as those measured during the high groundwater demand period, can have substantial changes in the streamflow in an unnamed tributary to the Elkhorn River. The streamflow measured during the high groundwater demand synoptic water-level and streamflow survey had decreased by nearly a factor of 20 when compared to the low groundwater demand period. The volume of groundwater discharge received by the Elkhorn River was estimated by examining the changes in streamflow between measurement locations. Streamflow measurements indicate that the groundwater discharge received by the Elkhorn River in the southern part of the study area was seasonably variable, making it difficult if not impossible to estimate an annual value. In the Elkhorn River, between the Elkhorn River at Waterloo, Nebr., streamflow-gaging station and the Q Street Bridge, streamflow measurements collected during the low groundwater demand period indicated a gain of 80 cubic feet per second, which is comparable to the gain estimated using aerial thermal infrared imagery and water temperature data. Streamflow measurements collected during the high groundwater demand period indicate a loss of 80 cubic feet per second across this same reach. In assessing water supply conditions in the lower Platte River system, the term “loss” in reference to streamflow in the Elkhorn River should be used with caution. Most likely, flow from the Elkhorn River which is “lost” to the groundwater system will later discharge to surface water closer to the confluence of the Platte and Elkhorn Rivers as underflow. A calibrated groundwater flow model of the study area likely is required to predict the fate of this water and to quantify groundwater discharge during varying hydrologic conditions along this reach. Aerial thermal infrared imagery indicated that much of the groundwater discharge in the southern part of the study area is focused across a 3-mile reach where the Elkhorn River turns southwest, perpendicular to the regional groundwater flow direction. Points of focused groundwater discharge were not detected with aerial thermal infrared imagery, indicating that groundwater discharge is diffuse rather than concentrated at focused points. Temperature-based streambed flux estimates indicated that strong regional groundwater gradients are not driving groundwater discharge and hyporheic flow is the dominant groundwater/surface-water exchange process.

Nebraska↗

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2&mdash;two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗

Connectivity of Mojave Desert tortoise populations—Management implications for maintaining a viable recovery network

Executive Summary The historic distribution of Mojave desert tortoises ( Gopherus agassizii ) was relatively continuous across the range, and the importance of tortoise habitat outside of designated tortoise conservation areas (TCAs) to recovery has long been recognized for its contributions to supporting gene flow between TCAs and to minimizing impacts and edge effects within TCAs. However, connectivity of Mojave desert tortoise populations has become a concern because of recent and proposed development of large tracts of desert tortoise habitat that cross, fragment, and surround designated conservation areas. This paper summarizes the underlying concepts and importance of connectivity for Mojave desert tortoise populations by reviewing current information on connectivity and providing information to managers for maintaining or enhancing desert tortoise population connectivity as they consider future proposals for development and management actions. Maintaining an ecological network for the Mojave desert tortoise, with a system of core habitats (TCAs) connected by linkages, is necessary to support demographically viable populations and long-term gene flow within and between TCAs. There are four points for wildlife and land-management agencies to consider when making decisions that could affect connectivity of Mojave desert tortoise populations (for example, in updating actions in resource management plans or amendments that could help maintain or restore functional connectivity in light of the latest information): Management of all desert tortoise habitat for persistence and connectivity . Desert tortoise populations continue to decline within most TCAs, and it is unlikely that trends are better in populations outside protected areas. Fragmentation exacerbates negative population trends by breaking large continuous populations into smaller isolated populations. Connectivity within large populations can enhance resilience to localized disturbances due to rescue by neighboring individuals. In contrast, smaller fragmented populations are resistant to rescue by their isolation and thus could suffer irreversible declines to extirpation from a variety of threats and stochastic events. Enhanced threat reduction to reverse declines within TCAs and to maintain occupied habitat in the surrounding matrix would help reduce the variability in population growth rates and improve the resilience of protected populations even while implementing efforts to improve connectivity. Each TCA has unique strengths and weaknesses regarding its ability to support minimum sustainable populations based on areal extent and its ability to support population increases based on landscape connection with adjacent populations. Considering how proposed projects (inside or outside of TCAs) affect connectivity and the ability of TCAs to support at least 5,000 adult tortoises (the numerical goal for each TCA) could help managers to maintain the resilience of TCAs to population declines. The same project, in an alternative location, could have very different impacts on local and regional populations. For example, within the habitat matrix surrounding TCAs, narrowly delineated corridors may not allow for natural population dynamics if they do not accommodate overlapping home ranges along most of their widths so that tortoises reside, grow, find mates, and produce offspring that can replace older tortoises. In addition, most habitat outside TCAs may receive more surface disturbance than habitat within TCAs. Therefore, managing the entire remaining matrix of desert tortoise habitat for permeability may be better than delineating fixed corridors. These concepts apply, especially given uncertainty about long-term condition of habitat, within and outside of TCAs under a changing climate. Ultimately, questions such as “ What are the critical linkages that need to be protected ?” could be better framed as “ How can we manage the remaining habitat matrix in ways that sustain ecological processes and habitat suitability for special status species ?” Land-management decisions made in the context of the latter question may be more conducive to maintenance of a functional ecological network. Limitations on landscape-level disturbance across habitat managed for the desert tortoise Clearly delineating habitat linkages and differentiating them from non-delineated areas by the uses that are permitted or prohibited within them by specific management guidelines can help achieve functional connectivity. Such guidelines would be most effective if they considered and accounted for all surface disturbances (for example, temporary disturbances such as fiberoptic lines or off-highway vehicle routes, right-of-ways, utility-scale solar development, urbanization) to the extent possible. A weighted framework that varies with the permanence or severity of the disturbance, and can be additive to quantify cumulative effects, could be useful (Xiong, 2020). For example, minor roads can alter tortoise movements independently of other features (Peaden and others, 2017; Hromada and others, 2020), but if the isolated dirt road is accompanied by a powerline that encourages raven predation (Xiong, 2020), then the two features together may be additive. Ignoring minor or temporary disturbance on the landscape could result in a cumulatively large impact that is not explicitly acknowledged (Goble, 2009); therefore, understanding and quantifying all surface disturbance on a given landscape is prudent. In California, the Bureau of Land Management established 0.1–1.0 percent caps on new surface-disturbance for TCAs and mapped linkages that address the issues described in number 1 of this list. Nevada, Utah, and Arizona currently do not have surface-disturbance limits. Limits comparable to those in the Desert Renewable Energy Conservation Plan (DRECP) would be 0.5 percent within TCAs and 1 percent within the linkages modeled by Averill-Murray and others (2013). Limits in some areas of California within the Desert Renewable Energy Conservation Plan, such as Ivanpah Valley, are more restrictive, at 0.1 percent. Continuity across the state line in Nevada could be achieved with comparable limits in the adjacent portion of Ivanpah Valley, as well as the Greater Trout Canyon Translocation Area and the Stump Springs Regional Augmentation Site. These more restrictive limits would help protect remaining habitat in the major interstate connectivity pathway through Ivanpah Valley and focal areas of population augmentation that provide additional population connectivity along the western flank of the Spring Mountains. In a recent study that analyzed 13 years of desert tortoise monitoring data, nearly all desert tortoise observations were at sites in which 5 percent or less of the surrounding landscape within 1 kilometer was disturbed (Carter and others, 2020a). To help maintain tortoise habitability and permeability across all other non-conservation-designated tortoise habitat, all surface disturbance could be limited to less than 5-percent development per square kilometer because the 5-percent threshold for development is the point at which tortoise occupation drops precipitously (Carter and others, 2020a). However, although individual desert tortoises were observed at development levels up to 5 percent, we do not know the fitness or reproductive characteristics of these individuals. This level of development also may not allow for long-term persistence of healthy populations that are of adequate size needed for demographic or functional connectivity; therefore, a conservative interpretation suggests that, ideally, development could be lower. Lower development levels would be particularly useful in areas within the upper 5th percentile of connectivity values modeled by Gray and others (2019). Reducing ancillary threats in places where connectivity is restricted to narrow strips of habitat, for example, narrow mountain passes or vegetated strips between solar development, could enhance the functionality of these vulnerable linkages. In such areas, maintaining multiple, redundant linkages could further enhance overall connectivity. Minimization of mortality from roads and maximization of passage under roads . Roads pose a significant threat to the long-term persistence of local tortoise populations, and roads of high traffic volume lead to severe population declines, which ultimately fragments populations farther away from the roads. Three points (a.–c.) pertain to reducing direct mortality of tortoises on the many paved roads that cross desert tortoise habitat and to maintaining a minimal level of permeability across these roads: Tortoise-exclusion fencing tied into culverts, underpasses, overpasses, or other passages below roads in desert tortoise habitat, would limit vehicular mortality of tortoises and provide opportunities for movement across the roads. Installation of shade structures on the habitat side of fences installed in areas with narrow population-depletion zones would limit overheating of tortoises that may pace the fence. Passages below highways could be maintained or retrofitted to ensure safe tortoise access, for example, by filling eroded drop-offs or modifying erosion-control features such as rip-rap to make them safer and more passable for tortoises. Wildlife management agencies could work with transportation departments to develop construction standards that are consistent with hydrologic/erosion management goals, while also incorporating a design and materials consistent with tortoise survival and passage and make the standards widely available. The process would be most effective if the status of passages was regularly monitored and built into management plans. Healthy tortoise populations along fenced highways could be supported by ensuring that land inside tortoise-exclusion fences is not so degraded that it leads to degradation of tortoise habitat outside the exclusion areas. For example, severe invasive plant infestations inside a highway exclusion could cause an increase of invasive plants outside the exclusion area and degrade habitat; therefore, invasive plants inside road rights of way could be mown or treated with herbicide to limit their spread into adjacent tortoise habitat and minimize the risk of these plants carrying wildfires into adjacent habitat. Adaptation of management based on new information . Future research will continue to build upon and refine models related to desert tortoise population connectivity and develop new ones. New models could consider landscape levels of development and be constructed such that they share common foundations to support future synthesis efforts. If model development was undertaken in partnership with entities that are responsible for management of desert tortoise habitat, it would facilitate incorporation of current and future modeling results into their land management decisions. There are specific topics that may be clarified with further evaluation: The effects of climate change on desert tortoise habitat, distribution, and population connectivity; The effects of large-scale fires, especially within repeatedly burned habitat, on desert tortoise distribution and population connectivity; The ability of solar energy facilities or similar developments to support tortoise movement and presence by leaving washes intact; leaving native vegetation intact whenever possible, or if not possible, mowing the site, allowing vegetation to re-sprout, and managing weeds; and allowing tortoises to occupy the sites; and The design and frequency of underpasses necessary to maintain functional demographic and genetic connectivity across linear features, like highways.

Arizona, California, Nevada↗

Analytical methods, numerical modeling, and monitoring strategies for evaluating the effects of ground-water withdrawals on unconfined aquifers in the New Jersey Coastal Plain

Analytical and numerical solutions of ground-water withdrawals in the unconfined part of the Kirkwood-Cohansey aquifer system of the Coastal Plain of New Jersey were evaluated for their usefulness in predicting the area of influence of a pumped well and in determining hydraulic characteristics of an aquifer. Additionally, simulations of ground-water withdrawal using a finite-difference model provided information on the ways in which prudent well-location strategies can disperse the local effects of withdrawal over a larger part of an aquifer system. The design of a monitoring network that is sensitive to the ground-water hydraulics of streams and wetlands of the Coastal Plain of New Jersey also was considered for its utility in providing hydrologic data necessary to establish the baseline hydrologic conditions near wetlands and streams and in signaling when ground-water levels are being adversely affected by withdrawals elsewhere in the system. The application of methods based on the Theis analytical solution to ground-water flow in unconfined aquifers can lead to erroneous estimates of the size of the area of influence generated by ground-water withdrawals. Analysis oftime-drawdown data from an unconfined aquifer system are best evaluated by means of the Neuman solution, which accounts for the effects of gravity drainage; however, the pumped well must be far enough from streams so that ground water is not drawn from nearby streams. Time-drawdown data from a test well in Winslow County, N.J., were analyzed by means of the Neuman solution. Results indicate that the aquifer has a relatively high vertical to horizontal anisotropy of 1:198, and a specific yield of 0.028, an indication that the area of influence of a pumped well at the test site would be relatively large. Results from a finite-difference ground-water-flow model of the northeastern part of the Mullica River Basin near Chesilhurst, N.J., show that the area influenced by a long-term withdrawal is best estimated from a steady state ground-water-flow analysis that includes the effects of average areal recharge. Withdrawal simulations indicate an order-of-magnitude difference between the size of the area of influence generated from a 3-day (72 hour) withdrawal and the size of the area produced under steady-state conditions. An aquifer characterized by a low specific yield will cause the area of influence to extend farther away from the pumped well. The contributing area of flow to the pumped well includes areas on the water table that would, under natural conditions, be incorporated into the contributing areas of flow to streams. Ground water that is drawn to a pumped well is diverted from nearby streams; the withdrawal decreases the size of the contributing areas of flow to streams by an amount equal to the contributing area of flow to the well. Withdrawals made from a well close to a stream divert ground water that would, under natural conditions, flow to the stream. The diverted ground water causes the area of influence of the well to be smaller than it would if the well were far from the stream. Water-table declines caused by withdrawals near streams are, to some degree, mitigated by ground-water diversion from streams. However, the withdrawals can significantly reduce ground-water seepage to nearby streams, especially along stream reaches and wetlands close to the well. Alternatively, these effects can be dispersed over a large part of the aquifer if wells are located on surface-water divides. Measurements of seasonal water-level fluctuations in the Mullica River Basin indicate that the greatest fluctuations in water levels are found in upland areas, where the average fluctuation is 3.4 feet. Fluctuations in hydraulic head in the wetland areas averages 1.3 feet. The bimodal average of ranges in water levels show that upland areas are more sensitive to recharge than lowland areas. The pattern of yearly mean water levels fluctuates irregularly about a long-term mean value. Abnormally low or high yearly average values that are brought on by periods of drought or excess recharge are short lived; over time, hydrologic conditions shift back to average levels under natural conditions. Wetland areas in the New Jersey Coastal Plain are characterized by ground-water seepage into wide, shallow depressions. Periods of inundation are longest in the deepest part of the depression, whereas inundation of areas near the fringes of wetlands due to ground-water seepage is only seasonal. The seepage face in the fringe areas expand and contract in response to seasonal variation in water-table elevation and in response to precipitation. Values of the aquifer storage coefficient and transmissivity can, in some cases, be determined by use of hydraulic head or streamflow recession analysis as an alternative to aquifer testing. The recession curves developed from hydro graphs of Middle Branch and McDonalds Branch in the New Jersey Coastal Plain indicate that the aquifer near McDonalds Branch has about 2.6 times the storage capacity of the aquifer adjacent to Middle Branch; this finding is consistent with the relatively small ranges of water-level changes measured in McDonalds Branch compared to those measured in Middle Branch.

New Jersey↗