USGS Science⌕ Search

SEARCH · USGS Science

Results for “Near Surface Geophysics”

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 127 records · Page 7Linked to original sources
USGS2001

Use of advanced borehole geophysical techniques to delineate fractured-rock ground-water flow and fractures along water-tunnel facilities in northern Queens County, New York

Advanced borehole geophysical methods were used to assess the geohydrology of crystalline bedrock along the course of a new water tunnel for New York City. The logging methods include natural gamma, spontaneous potential, single-point resistance, mechanical and acoustic caliper, focused electromagnetic induction, electromagnetic resistivity, magnetic susceptibility, borehole-fluid temperature and conductance, differential temperature, heat-pulse flowmeter, acoustic televiewer, borehole deviation, optical televiewer, and borehole radar. Integrated interpretation of the geophysical logs from an 825-foot borehole (1) provided information on the extent, orientation, and structure (foliation and fractures) within the entire borehole, including intensely fractured intervals from which core recovery may be poor; (2) delineated transmissive fracture zones intersected by the borehole and provided estimates of their transmissivity and hydraulic head; and (3) enabled mapping of the location and orientation of structures at distances as much as 100 ft from the borehole. Analyses of the borehole-wall image and the geophysical logs from the borehole on Crescent Street, in northern Queens County, are presented here to illustrate the application of the methods. The borehole penetrates gneiss and other crystalline bedrock that has predominantly southeastward dipping foliation and nearly horizontal and southeastward-dipping fractures. The heat-pulse flowmeter logs obtained under pumping and nonpumping conditions, together with the other geophysical logs, indicate five transmissive fracture zones. More than 90 percent of the open-hole transmissivity is associated with a fracture zone 272 feet BLS (below land surface). A transmissive zone at 787 feet BLS that consists of nearly parallel fractures lies within the projected tunnel path; here the hydraulic head is 12 to 15 feet lower than that of transmissive zones above the 315-foot depth. The 60-megahertz directional borehole radar logs indicate the location and orientation of two closely spaced radar reflectors that would intersect the projection of the borehole below its drilled depth. Subsequent excavation of the tunnel past the borehole allowed comparison of the log analysis with conditions observed in the tunnel. The tunnel was found to intersect gneiss with southeastward dipping foliation; many nearly horizontal fractures; and a southeastward dipping fracture zone whose location, character, and orientation was consistent with that of the mapped radar reflectors. The fracture zone produced inflow to the tunnel at a rate of 50 to 100 gallons per minute. All conditions indicated by the logging methods were consistent with those observed within the tunnel.

New York↗
USGS2011

Expansion rate and geometry of floating vegetation mats on the margins of thermokarst lakes, northern Seward Peninsula, Alaska, USA

Investigations on the northern Seward Peninsula in Alaska identified zones of recent (<50 years) permafrost collapse that led to the formation of floating vegetation mats along thermokarst lake margins. The occurrence of floating vegetation mat features indicates rapid degradation of near‐surface permafrost and lake expansion. This paper reports on the recent expansion of these collapse features and their geometry is determined using geophysical and remote sensing measurements. The vegetation mats were observed to have an average thickness of 0.57 m and petrophysical modeling indicated that gas content of 1.5–5% enabled floatation above the lake surface. Furthermore, geophysical investigation provides evidence that the mats form by thaw and subsidence of the underlying permafrost rather than terrestrialization. The temperature of the water below a vegetation mat was observed to remain above freezing late in the winter. Analysis of satellite and aerial imagery indicates that these features have expanded at maximum rates of 1–2 m yr ‐1 over a 56 year period. Including the spatial coverage of floating ‘thermokarst mats’ increases estimates of lake area by as much as 4% in some lakes.

Alaska↗
USGS2014

Characterization of subsurface stratigraphy along the lower American River floodplain using electrical resistivity, Sacramento, California, 2011

In July 2011, the U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, completed a geophysical survey using electrical resistivity along an approximately 6-mile reach of the lower American River in Sacramento, California, to map near-surface lithological variations. This survey is a part of a manifold and comprehensive study of river-flow dynamics and geologic boundary-property knowledge necessary to estimate scour potential and levee erosion risk. Data were acquired on the left (south or west) bank between river mile 5 and 10.7 as well as a short section on the right bank from river mile 5.4 to 6. Thirteen direct-current resistivity profiles and approximately 8.3 miles of capacitively coupled resisistivity data were acquired along accessible areas of the floodplain between the levee and river bank. Capacitively coupled resistivity was used as a reconnaissance tool, because it allowed for greater spatial coverage of data but with lower resolution and depth of investigation than the DC resistivity method. The study area contains Pleistocene-age alluvial deposits, dominated by gravels, sands, silts, and clays, that vary in both lateral extent and depth. Several generations of lithologic logs were used to help interpret resistivity variations observed in the resistivity models.

California↗
USGS2011

Sources and physical processes responsible for OH/H2O in the lunarsoil as revealed by the Moon Mineralogy Mapper (M3)

Analysis of two absorption features near 3 μ m in the lunar reflectance spectrum, observed by the orbiting M 3 spectrometer and interpreted as being due to OH and H 2 O, is presented, and the results are used to discuss the processes producing these molecules. This analysis focuses on the dependence of the absorptions on lunar physical properties, including composition, illumination, latitude, and temperature. Solar wind proton-induced hydroxylation is proposed as the creation process, and its products could be a source for other reported types of hydrogen-rich material and water. The irregular and damaged fine-grained lunar soil seems especially adapted for trapping solar wind protons and forming OH owing to abundant dangling oxygen bonds. The M 3 data reveal that the strengths of the two absorptions are correlated and widespread, and both are correlated with lunar composition but in different ways. Feldspathic material seems richer in OH. These results seem to rule out water from the lunar interior and cometary infall as major sources. There appear to be correlations of apparent band strengths with time of day and lighting conditions. However, thermal emission from the Moon reduces the apparent strengths of the M 3 absorptions, and its removal is not yet completely successful. Further, many of the lunar physical properties are themselves intercorrelated, and so separating these dependencies on the absorptions is difficult, due to the incomplete M 3 data set. This process should also operate on other airless silicate surfaces, such as Mercury and Vesta, which will be visited by the Dawn spacecraft in mid-2011.

Journal of Geophysical Research E: Planets↗
USGS1949

Suggested correlation in Simpson seeps region

During the 1949 drilling season eight core tests were completed in the Simpson Seeps area up to October 25th„ This program, which is still in progress, was undertaken to provide an explanation for the surface seeps near Caps Simpson. The shallowest of the eight holes reached 800 feet; the deepest was drilled to a depth of 1460 feet. Geophysical data on the shallow horizons indicate that the Simpson Seeps Core Tests shown on the section have been drilled in an area of nearly north-south strike. The five core tests that most nearly approximate a section normal to the strike were chosen for the accompanying cross Section. An east-west strike, accompanied by strong anomalous dips of fifteen to twenty degrees to the north is indicated by geophysical data about five miles north of the cross section.

Alaska↗
USGS1993

Geophysical investigations of concealed faults near Yucca Mountain, Southwest Nevada

Detailed gravity and ground magnetic data collected along surveyed traverses across Midway Valley, on the eastern flank of Yucca Mountain, Nevada reveal that these methods can be used to delineate concealed faults. These studies are part of an effort to evaluate faulting in the vicinity of the proposed surface facilities for a potential nuclear waste repository at Yucca Mountain. The largest gravity and magnetic anomaly in the vicinity of Midway Valley is associated with the Paintbrush fault on the west flank of Alice Ridge. Geophysical data infer a vertical offset of about 200 m (650 ft). Another prominent gravity and magnetic anomaly is associated with the Bow Ridge fault in the western part of Midway Valley.

Conference Paper↗
USGS2015

11.12 – Tools and techniques: gravitational method

The gravitational method is used to investigate density variations within the subsurface at depths of several meters to tens of meters, as in depth-to-bedrock investigations, or at depths of several kilometers, as in sedimentary basin thickness investigations. This chapter covers fundamental relations, densities of Earth materials, instruments, field procedures, data reduction, filtering, forward modeling, inversion, and field examples. The focus is on near-surface investigations as distinct from the solid Earth studies found elsewhere in this treatise. The gravitational method is often used in conjunction with other geophysical methods, such as the magnetic method or the seismic method, which target similar physical properties at similar depths.

Book chapter↗
USGS2015

Geologic control on the evolution of the inner shelf morphology offshore of the Mississippi barrier islands, northern Gulf of Mexico, USA

Between 2008 and 2013, high-resolution geophysical surveys were conducted around the Mississippi barrier islands and offshore. The sonar surveys included swath and single-beam bathymetry, sidescan, and chirp subbottom data collection. The geophysical data were groundtruthed using vibracore sediment collection. The results provide insight into the evolution of the inner shelf and the relationship between the near surface geologic framework and the morphology of the coastal zone. This study focuses on the buried Pleistocene fluvial deposits and late Holocene shore-oblique sand ridges offshore of Petit Bois Island and Petit Bois Pass. Prior to this study, the physical characteristics, evolution, and interrelationship of the ridges between both the shelf geology and the adjacent barrier island platform had not been evaluated. Numerous studies elsewhere along the coastal margin attribute shoal origin and sand-ridge evolution to hydrodynamic processes in shallow water (<20 m). Here we characterize the correlation between the geologic framework and surface morphology and demonstrate that the underlying stratigraphy must also be considered when developing an evolutionary conceptual model. It is important to understand this near surface, nearshore dynamic in order to understand how the stratigraphy influences the long-term response of the coastal zone to sea-level rise. The study also contributes to a growing body of work characterizing shore-oblique sand ridges which, along with the related geology, are recognized as increasingly important components to a nearshore framework whose origins and evolution must be understood and inventoried to effectively manage the coastal zone.

Continental Shelf Research↗
USGS2006

Quantitative estimation of minimum offset for multichannel surface-wave survey with actively exciting source

Multichannel analysis of surface waves is a developing method widely used in shallow subsurface investigations. The field procedures and related parameters are very important for successful applications. Among these parameters, the source-receiver offset range is seldom discussed in theory and normally determined by empirical or semi-quantitative methods in current practice. This paper discusses the problem from a theoretical perspective. A formula for quantitatively evaluating a layered homogenous elastic model was developed. The analytical results based on simple models and experimental data demonstrate that the formula is correct for surface wave surveys for near-surface applications. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗
USGS2001

Resolution analysis of finite fault source inversion using one- and three-dimensional Green's functions 2. Combining seismic and geodetic data

Using numerical tests for a prescribed heterogeneous earthquake slip distribution, we examine the importance of accurate Green's functions (GF) for finite fault source inversions which rely on coseismic GPS displacements and leveling line uplift alone and in combination with near-source strong ground motions. The static displacements, while sensitive to the three-dimensional (3-D) structure, are less so than seismic waveforms and thus are an important contribution, particularly when used in conjunction with waveform inversions. For numerical tests of an earthquake source and data distribution modeled after the 1994 Northridge earthquake, a joint geodetic and seismic inversion allows for reasonable recovery of the heterogeneous slip distribution on the fault. In contrast, inaccurate 3-D GFs or multiple 1-D GFs allow only partial recovery of the slip distribution given strong motion data alone. Likewise, using just the GPS and leveling line data requires significant smoothing for inversion stability, and hence, only a blurred vision of the prescribed slip is recovered. Although the half-space approximation for computing the surface static deformation field is no longer justifiable based on the high level of accuracy for current GPS data acquisition and the computed differences between 3-D and half-space surface displacements, a layered 1-D approximation to 3-D Earth structure provides adequate representation of the surface displacement field. However, even with the half-space approximation, geodetic data can provide additional slip resolution in the joint seismic and geodetic inversion provided a priori fault location and geometry are correct. Nevertheless, the sensitivity of the static displacements to the Earth structure begs caution for interpretation of surface displacements, particularly those recorded at monuments located in or near basin environments.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Hydrothermal processes at Gusev Crater: An evaluation of Paso Robles class soils

The Mars Exploration Rover Spirit analyzed multiple occurrences of sulfur-rich, light-toned soils along its traverse within Gusev Crater. These hydrated deposits are not readily apparent in images of undisturbed soil but are present at shallow depths and were exposed by the actions of the rover wheels. Referred to as 'Paso Robles' class soils, they are dominated by ferric iron sulfates, silica, and Mg-sulfates. Ca-sulfates, Ca-phosphates, and other minor phases are also indicated in certain specific samples. The chemical compositions are highly variable over both centimeter-scale distances and between the widely separated exposures, but they clearly reflect the elemental signatures of nearby rocks. The quantity of typical basaltic soil mixed into the light-toned materials prior to excavation by the rover wheels is minimal, suggesting negligible reworking of the deposits after their initial formation. The mineralogy, geochemistry, variability, association with local compositions, and geologic setting of the deposits suggest that Paso Robles class soils likely formed as hydrothermal and famarolic condensates derived from magma degassing and/ or oxidative alteration of crustal iron sulfide deposits. Their occurrence as unconsolidated, near-surface soils permits, though does not require, an age that is significantly younger than that of the surrounding rocks. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2021

Airborne geophysical imaging of weak zones on Iliamna Volcano, Alaska: Implications for slope stability

Water‐saturated, hydrothermally altered rocks reduce the strength of volcanic edifices and increase the potential for sector collapses and far‐traveled mass flows of unconsolidated debris. Iliamna Volcano is an andesitic stratovolcano located on the western side of the Cook Inlet, ∼225 km southwest of Anchorage and is a source of repeated avalanches. The widespread snow and ice cover on Iliamna Volcano make surface alteration difficult to identify. However, intense hydrothermal alteration significantly reduces both the electrical resistivity and magnetization of volcanic rock and can therefore be identified with airborne geophysical measurements. We use airborne electromagnetic and magnetic data to map snow and ice thickness and identify underlying alteration zones at Iliamna Volcano, Alaska. Resistivities were calculated to an average depth of >300 m, and a 3‐D susceptibility model extends from the surface to the base of the volcano, about 3,000 m below the summit. Geophysical models image low resistivity (<30 ohm‐m) and low susceptibilities near the summit of Iliamna and below its older vent complex, with the low susceptibilities indicating alteration up to ∼800 m in thickness. Thin conductors (∼50–100 m thick) on the edifice slopes coincide with recorded locations of repeated debris avalanches over the past ∼60 years and are attributed to saturated zones at high elevation. Three‐dimensional slope stability models based upon the geophysically constrained alteration distribution suggest the edifice of Iliamna is unstable and could lead to collapse scars ∼400 m deep near the current and former vent complexes.

Alaska↗
USGS2012

Upper crustal structure of the Mount Hood, Oregon, region as revealed by time term analysis

Seismic refraction data with a dense areal distribution were collected to study the seismic structure of Mount Hood and the surrounding region. This area is typical of Cascade volcanoes and is geologically quite complex. The prime goals of this project were to search for velocity variations in the upper crustal rocks and to determine if the velocity of these rocks is anisotropic. A new system, including 100 remote recording units, was developed to facilitate the collection of data in this type of survey. The data collected in this study reveal a large variation in velocity and thickness of the uppermost crustal rocks that is probably typical of the High Cascade province. A regional structural pattern surrounding Mount Hood, where there is a marked thinning of low-velocity near-surface rocks, suggests that the present edifice of Mount Hood lies on top of a much larger structure, possibly the roof of a large batholith that was emplaced prior to the eruption of the volcanic rocks that form the modern mountain. A straightforward time term analysis of the data reveals this structure. When the time terms are used to correct the travel time observations for the variation in structure under each station, the remaining set of residuals indicates a variation of travel time with azimuth. While a systematic error in the time term reduction cannot be ruled out, this variation is probably the result of anisotropic velocity structure in the underlying refractor. The alignment of the direction of maximum velocity with the direction of maximum principal stress and the decrease in anisotropy with depth suggest that the anisotropic velocity structure is related to the opening and closing of fractures in the present stress field. The direction of maximum velocity is about N25°W, and the degree of anisotropy appears to vary with depth from about 3% near the surface to zero at depths below 8 km.

Journal of Geophysical Research Solid Earth↗
USGS1985

A study of the depth of weathering and its relationship to the mechanical properties of near-surface rocks in the Mojave Desert

Weathered granite extends 70 m deep at Hi Vista in the arid central Mojave Desert of southern California. The low strength of this granite is due to the alteration of biotite and chlorite montmorillonite. Deep weathering probably occurs in most granites, although we cannot rule out some anomalous mechanisms at Hi Vista. Geophysical instruments set in these slightly altered rocks are limited by the unstable behavior of the rocks. Thus, tectonic signals from instruments placed in shallow boreholes give vague results. Geophysical measurements of these weathered rocks resemble measurements of granitic rocks near major faults. The rheology of the rocks in which instruments are placed limits the useful sensitivity of the instruments. ?? 1985 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗
USGS2010

Using airborne geophysical surveys to improve groundwater resource management models

Increasingly, groundwater management requires more accurate hydrogeologic frameworks for groundwater models. These complex issues have created the demand for innovative approaches to data collection. In complicated terrains, groundwater modelers benefit from continuous high‐resolution geologic maps and their related hydrogeologic‐parameter estimates. The USGS and its partners have collaborated to use airborne geophysical surveys for near‐continuous coverage of areas of the North Platte River valley in western Nebraska. The survey objectives were to map the aquifers and bedrock topography of the area to help improve the understanding of groundwater‐surface‐water relationships, leading to improved water management decisions. Frequency‐domain heliborne electromagnetic surveys were completed, using a unique survey design to collect resistivity data that can be related to lithologic information to refine groundwater model inputs. To render the geophysical data useful to multidimensional groundwater models, numerical inversion is necessary to convert the measured data into a depth‐dependent subsurface resistivity model. This inverted model, in conjunction with sensitivity analysis, geological ground truth (boreholes and surface geology maps), and geological interpretation, is used to characterize hydrogeologic features. Interpreted two‐ and three‐dimensional data coverage provides the groundwater modeler with a high‐resolution hydrogeologic framework and a quantitative estimate of framework uncertainty. This method of creating hydrogeologic frameworks improved the understanding of flow path orientation by redefining the location of the paleochannels and associated bedrock highs. The improved models reflect actual hydrogeology at a level of accuracy not achievable using previous data sets.

Conference Paper↗
USGS1993

Application of borehole geophysics in defining the wellhead protection area for a fractured crystalline bedrock aquifer

Wellbore geophysical techniques were used to characterize fractures and flow in a bedrock aquifer at a site near Blackwater Brook in Dover, New Hampshire. The primary focus ofthis study was the development of a model to assist in evaluating the area surrounding a planned water supply well where contaminants introduced at the land surface might be induced to flow towards a pumping well. Well logs and geophysical surveys used in this study included lithologic logs based on examination of cuttings obtained during drilling; conventional caliper and natural gamma logs; video camera and acoustic televiewer surveys; highresolution vertical flow measurements under ambient conditions and during pumping; and borehole fluid conductivity logs obtained after the borehole fluid was replaced with deionized water. These surveys were used for several applications: 1) to define a conceptual model of aquifer structure to be used in groundwater exploration; 2) to estimate optimum locations for test and observation wells; and 3) to delineate a wellhead protection area (WHPA) for a planned water supply well. Integration of borehole data with surface geophysical and geological mapping data indicated that the study site lies along a northeast-trending intensely fractured contact zone between surface exposures of quartz monzonite and metasedimentary rocks. Four of five bedrock boreholes at the site were estimated to produce more than 150 gallons per minute (gpm) (568 L/min) of water during drilling. Aquifer testing and other investigations indicated that water flowed to the test well along fractures parallel to the northeast-trending contact zone and along other northeast and north-northwest-trending fractures. Statistical plots of fracture strikes showed frequency maxima in the same northeast and north-northwest directions, although additional maxima occurred in other directions. Flowmeter surveys and borehole fluid conductivity logging after fluid replacement were used to identify water-producing zones in the boreholes; fractures associated with inflow into boreholes showed a dominant northeast orientation. Borehole fluid conductivity logging after fluid replacement also gave profiles of such water-quality parameters as fluid electrical conductivity (FEC), pH, temperature, and oxidation-reduction potential, strengthening the interpretation of crossconnection of boreholes by certain fracture zones. The results of this study showed that the application of these borehole geophysical techniques at the Blackwater Brook site led to an improved understanding of such parameters as fracture location, attitude, flow direction and velocity, and water quality; all of which are important in the determination of a WHPA.

Log Analyst↗
USGS2005

Fractured-aquifer hydrogeology from geophysical logs; the passaic formation, New Jersey

The Passaic Formation consists of gradational sequences of mudstone, siltstone, and sandstone, and is a principal aquifer in central New Jersey. Ground‐water flow is primarily controlled by fractures interspersed throughout these sedimentary rocks and characterizing these fractures in terms of type, orientation, spatial distribution, frequency, and transmissivity is fundamental towards understanding local fluid‐transport processes. To obtain this information, a comprehensive suite of geophysical logs was collected in 10 wells roughly 46 m in depth and located within a .05 km 2 area in Hopewell Township, New Jersey. A seemingly complex, heterogeneous network of fractures identified with an acoustic televiewer was statistically reduced to two principal subsets corresponding to two distinct fracture types: (1) bedding‐plane partings and (2) high‐angle fractures. Bedding‐plane partings are the most numerous and have an average strike of N84°W and dip of 20° N. The high‐angle fractures are oriented subparallel to these features, with an average strike of N79° E and dip of 71° S, making the two fracture types roughly orthogonal. Their intersections form linear features that also retain this approximately east‐west strike. Inspection of fluid temperature and conductance logs in conjunction with flow meter measurements obtained during pumping allows the transmissive fractures to be distinguished from the general fracture population. These results show that, within the resolution capabilities of the logging tools, approximately 51 (or 18 percent) of the 280 total fractures are water producing. The bedding‐plane partings exhibit transmissivities that average roughly 5 m 2 /day and that generally diminish in magnitude and frequency with depth. The high‐angle fractures have average transmissivities that are about half those of the bedding‐plane partings and show no apparent dependence upon depth. The geophysical logging results allow us to infer a distinct hydrogeologic structure within this aquifer that is defined by fracture type and orientation. Fluid flow near the surface is controlled primarily by the highly transmissive, subhorizontal bedding‐plane partings. As depth increases, the high‐angle fractures apparently become more dominant hydrologically.

Ground Water↗
USGS2006

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗