Search USGS⌕ Search

Geology topics

Steven P. Phillips

Publications and source records attributed to Steven P. Phillips.

At least 19 recordsLinked to original sources
USGS2017

Electrical resistivity investigation of fluvial geomorphology to evaluate potential seepage conduits to agricultural lands along the San Joaquin River, Merced County, California, 2012–13

Increased flows in the San Joaquin River, part of the San Joaquin River Restoration Program, are designed to help restore fish populations. However, increased seepage losses could result from these higher restoration flows, which could exacerbate existing drainage problems in neighboring agricultural lands and potentially damage crops. Channel deposits of abandoned river meanders that are hydraulically connected to the river could act as seepage conduits, allowing rapid and widespread water-table rise during restoration flows. There is a need to identify the geometry and properties of these channel deposits to assess their role in potential increased seepage effects and to evaluate management alternatives for reducing seepage. Electrical and electromagnetic surface geophysical methods have provided a reliable proxy for lithology in studies of fluvial and hyporheic systems where a sufficient electrical contrast exists between deposits of differing grain size. In this study, direct-current (DC) resistivity was used to measure subsurface resistivity to identify channel deposits and to map their subsurface geometry. The efficacy of this method was assessed by using DC resistivity surveys collected along a reach of the San Joaquin River in Merced County, California, during the summers of 2012 and 2013, in conjunction with borings and associated measurements from a hydraulic profiling tool. Modeled DC resistivity data corresponded with data from cores, hand-auger samples, a hydraulic profiling tool, and aerial photographs, confirming that DC resistivity is effective for differentiating between silt and sand deposits in this setting. Modeled DC resistivity data provided detailed two-dimensional cross-sectional resistivity profiles to a depth of about 20 meters. The distribution of high-resistivity units in these profiles was used as a proxy for identifying areas of high hydraulic conductivity. These data were used subsequently to guide the location and depth of wells installed onsite for monitoring flow in the channel deposits. Estimates of the cross-sectional area of channel deposits from DC resistivity pseudosections can provide critical input for groundwater-flow models designed to simulate river seepage and evaluate seepage-management alternatives.

California↗
USGS2015

Sustainable groundwater management in California

The U.S. Geological Survey (USGS) uses data collection, modeling tools, and scientific analysis to help water managers plan for, and assess, hydrologic issues that can cause “undesirable results” associated with groundwater use. This information helps managers understand trends and investigate and predict effects of different groundwater-management strategies.

California↗
USGS2015

Hydrologic model of the Modesto Region, California, 1960-2004

Strategies for managing water supplies and groundwater quality in the Modesto region of the eastern San Joaquin Valley, California, are being formulated and evaluated by the Stanislaus and Tuolumne Rivers Groundwater Basin Association. Management issues and goals in the basin include an area in the lower part of the basin that requires drainage of the shallow water table to sustain agriculture, intra- and inter-basin migration of poor-quality groundwater, and efficient management of surface and groundwater supplies. To aid in the evaluation of water-management strategies, the U.S. Geological Survey and the Stanislaus and Tuolumne Rivers Groundwater Basin Association have developed a hydrologic model that simulates monthly groundwater and surface-water flow as governed by aquifer-system properties, annual and seasonal variations in climate, surface-water flow and availability, water use, and land use. The model was constructed by using the U.S. Geological Survey groundwater-modeling software MODFLOW-OWHM with the Farm Process. Available measurements of groundwater pumped for municipal, irrigation, and drainage purposes are specified in the model, as are deliveries of surface water. Private irrigation pumping and recharge associated with agricultural land use were estimated by using the Farm Process in MODFLOW-OWHM, which simulates landscape processes associated with irrigated agriculture and other land uses. The distribution of hydraulic conductivity in the aquifer system was constrained by using data from more than 3,500 drillers' logs. The model was calibrated to 4,061 measured groundwater levels in 109 wells and 2,739 mean monthly surface-water flows measured at 6 streamgages during 1960-2004 by using a semi-automated method of parameter estimation. The model fit to groundwater levels was good, with an absolute mean residual of 0.8 feet; 74 percent of simulated heads were within 10 feet of those observed. The model fit to streamflow was biased low, but reasonable overall; the absolute mean residual of streamflow was 780 cubic feet per second, and 68 percent of simulated streamflows were within 500 cubic feet per second of observed. Hydrographs both of groundwater levels and streamflow indicated overall an acceptable fit to observed trends. Simulated private agricultural pumpage ranged from about 780,000 to 1,380,000 acre-feet per year and averaged about 1,000,000 acre-feet per year from 1960 to 2004. Simulated deep percolation, or groundwater recharge from precipitation and irrigation, varied with climate and land use from about 1,100,000 to 1,700,000 acre-feet per year, averaging 1,360,000 acre-feet per year. Key limitations of the model with respect to estimating these large components of the water budget are the uncertainty associated with actual irrigation deliveries and irrigation efficiencies and the lack of metered data for private agricultural groundwater pumping. Different assumptions with respect to irrigation deliveries and efficiencies, and other model input, would result in different estimates of private agricultural groundwater use. The simulated exchange between groundwater and surface water was a small percentage of streamflow, typically ranging within a loss or gain of about 2 cubic feet per second per mile. The simulated exchange compared reasonably with limited independent estimates available, but substantial uncertainty is associated with these estimates.

California↗
USGS2015

Hindcast of water availability in regional aquifer systems using MODFLOW Farm Process

Coupled groundwater and surface-water components of the hydrologic cycle can be simulated by the Farm Process for MODFLOW (MF-FMP) in both irrigated and non-irrigated areas and aquifer-storage and recovery systems. MF-FMP is being applied to three productive agricultural regions of different scale in the State of California, USA, to assess the availability of water and the impacts of alternative management decisions. Hindcast simulations are conducted for similar periods from the 1960s to near recent times. Historical groundwater pumpage is mostly unknown in one region (Central Valley) and is estimated by MF-FMP. In another region (Pajaro Valley), recorded pumpage is used to calibrate model-estimated pumpage. Multiple types of observations are used to estimate uncertain parameters, such as hydraulic, land-use, and farm properties. MF-FMP simulates how climate variability and water-import availability affect water demand and supply. MF-FMP can be used to predict water availability based on anticipated changes in anthropogenic or natural water demands. Keywords groundwater; surface-water; irrigation; water availability; response to climate variability/change

Conference Paper↗
USGS2014

Documentation of a groundwater flow model (SJRRPGW) for the San Joaquin River Restoration Program study area, California

To better understand the potential effects of restoration flows on existing drainage problems, anticipated as a result of the San Joaquin River Restoration Program (SJRRP), the U.S. Geological Survey (USGS), in cooperation with the U.S. Bureau of Reclamation (Reclamation), developed a groundwater flow model (SJRRPGW) of the SJRRP study area that is within 5 miles of the San Joaquin River and adjacent bypass system from Friant Dam to the Merced River. The primary goal of the SJRRP is to reestablish the natural ecology of the river to a degree that restores salmon and other fish populations. Increased flows in the river, particularly during the spring salmon run, are a key component of the restoration effort. A potential consequence of these increased river flows is the exacerbation of existing irrigation drainage problems along a section of the river between Mendota and the confluence with the Merced River. Historically, this reach typically was underlain by a water table within 10 feet of the land surface, thus requiring careful irrigation management and (or) artificial drainage to maintain crop health. The SJRRPGW is designed to meet the short-term needs of the SJRRP; future versions of the model may incorporate potential enhancements, several of which are identified in this report. The SJRRPGW was constructed using the USGS groundwater flow model MODFLOW and was built on the framework of the USGS Central Valley Hydrologic Model (CVHM) within which the SJRRPGW model domain is embedded. The Farm Process (FMP2) was used to simulate the supply and demand components of irrigated agriculture. The Streamflow-Routing Package (SFR2) was used to simulate the streams and bypasses and their interaction with the aquifer system. The 1,300-square mile study area was subdivided into 0.25-mile by 0.25-mile cells. The sediment texture of the aquifer system, which was used to distribute hydraulic properties by model cell, was refined from that used in the CVHM to better represent the natural heterogeneity of aquifer-system materials within the model domain. In addition, the stream properties were updated from the CVHM to better simulate stream-aquifer interactions, and water-budget subregions were refined to better simulate agricultural water supply and demand. External boundary conditions were derived from the CVHM. The SJRRPGW was calibrated for April 1961 to September 2003 by using groundwater-level observations from 133 wells and streamflow observations from 19 streamgages. The model was calibrated using public-domain parameter estimation software (PEST) in a semi-automated manner. The simulated groundwater-level elevations and trends (including seasonal fluctuations) and surface-water flow magnitudes and trends reasonably matched observed data. The calibrated model is planned to be used to assess the potential effects of restoration flows on agricultural lands and the relative capabilities of proposed SJRRP actions to reduce these effects.

California↗
USGS2014

Groundwater-flow and land-subsidence model of Antelope Valley, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley groundwater basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, groundwater provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most groundwater pumping in the valley occurs in the Antelope Valley groundwater basin, which includes the rapidly growing cities of Lancaster and Palmdale. Groundwater-level declines of more than 270 feet in some parts of the groundwater basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may increase reliance on groundwater. In 2011, the Los Angeles County Superior Court of California ruled that the Antelope Valley groundwater basin is in overdraft—groundwater extractions are in excess of the Court-defined safe yield of the groundwater basin. The Court determined that the safe yield of the adjudicated area of the basin was 110,000 acre-feet per year (acre-ft/yr). Natural recharge is an important component of total groundwater recharge in Antelope Valley; however, the exact quantity and distribution of natural recharge, primarily in the form of mountain-front recharge, is uncertain, with total estimates ranging from 30,000 to 160,000 acre-ft/yr. Technical experts, retained by parties to the adjudication, used 60,000 acre-ft/yr to estimate the sustainable yield of the basin, and this value was used in this study. In order to better understand the uncertainty associated with natural recharge and to provide a tool to aid in groundwater management, a numerical model of groundwater flow and land subsidence in the Antelope Valley groundwater basin was developed using old and new geohydrologic information. The groundwater-flow system consists of three aquifers: the upper, middle, and lower aquifers. The three aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to groundwater development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Groundwater flowed from the recharge areas to discharge areas around the playas where it discharged from the aquifer system as either evapotranspiration or from springs. Partial barriers to horizontal groundwater flow, such as faults, have been identified in the groundwater basin. Water-level declines owing to groundwater development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the groundwater system. Infiltration of return flow from agricultural irrigation has become an important source of recharge to the aquifer system. The groundwater-flow model of the basin was discretized horizontally into a grid of 130 rows and 118 columns of square cells 1 kilometer (0.621 mile) on a side, and vertically into four layers representing the upper (two layers), middle (one layer), and lower (one layer) aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915–95, by using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley groundwater basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of groundwater-flow and land-subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Groundwater pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual agricultural pumpage peaked at 395,000 acre-feet (acre-ft) in 1951 and then declined because of declining agricultural production. Recharge from irrigation return flows was assumed to be 30 percent of agricultural pumpage; delays associated with return flow moving through the unsaturated zone were also simulated. The annual quantity of mountain-front recharge initially was based on estimates from previous studies. The model was calibrated using the PEST software suite; prior information from the area was incorporated through the use of Tikhonov regularization. During model calibration, the estimated mountain-front recharge was reduced from the previous estimate of 30,300 acre-ft/yr to 29,150 acre-ft/yr. Results of the simulations using the calibrated model indicate that simulated groundwater pumpage exceeded recharge in most years, resulting in an estimated cumulative depletion in groundwater storage of 8,700,000 acre-ft during the transient-simulation period (1915–2005). About 15,000,000 acre-ft of cumulative groundwater pumpage was simulated during the transient-simulation period (1915–2005), reaching a maximum rate of about 400,000 acre-ft/yr in 1951. Groundwater pumpage resulted in simulated hydraulic heads declining by more than 150 feet (ft) compared to 1915 conditions in agricultural areas. The decline in hydraulic head in the groundwater basin is the result of this depletion of groundwater storage. In turn, the simulated decline in hydraulic head in the groundwater basin has resulted in the decrease in natural discharge from the basin and has caused compaction of aquitards, resulting in land subsidence. The areal distribution of total simulated land subsidence for 2005, after about 90 years of groundwater development, indicates that land subsidence occurred throughout almost the entire Lancaster subbasin, with a maximum of about 9.4 ft in the central and eastern parts of the subbasin. An important objective of this study was to systematically address the uncertainty in estimates of natural recharge and related aquifer parameters by using the groundwater-flow and land-subsidence model with observational data and expert knowledge. After the model was calibrated to the observations and a reasonable parameter set obtained, the parameter null space—parameter values that do not appreciably affect the model calibration but may have importance for prediction—was identified. The effect of parameter uncertainty on the estimation of mountain-front recharge was addressed using the Null-Space Monte Carlo method. The Pareto trade-off method of visualizing uncertainty was also used to portray the reasonableness of larger natural-recharge rates. Results indicate that the total mountain-front recharge likely ranges between 28,000 and 44,000 acre-ft/yr, which is appreciably less than published estimates of 60,000 acre-ft/yr. Additionally, expected errors associated with agricultural pumpage estimates used in this study were found to have relatively little effect on the estimates of mountain-front recharge, reflecting the difficulty in increasing recharge through manipulation of other components of the water budget. The calibrated model was used to simulate the response of the aquifer to potential future pumping scenarios: (1) no change in the distribution of pumpage, or status quo; (2) redistribution of pumpage; and (3) artificial recharge. All three of these scenarios specify a total pumpage throughout the Antelope Valley of 110,000 acre-ft/yr according to the safe yield value ruled by the Los Angeles County Superior Court of California. This reduction in groundwater pumpage is assumed uniform throughout the basin, based on a 10-percent reduction of the total pumpage in 2005 to achieve the 110,000 acre-ft/yr level. The calibrated Antelope Valley groundwater-flow and land-subsidence model was used to simulate the hydrologic effects of the three groundwater-management scenarios during a 50-year period by using the reduced, temporally constant, pumpage distribution. Results from the first scenario indicated that the total drawdown observed since predevelopment would continue, with values exceeding 325 ft near Palmdale; consequently, land subsidence would also continue, with additional subsidence (since 2005) exceeding 3 ft in the central part of the Lancaster subbasin. The second scenario evaluated redistributing pumpage from areas in the Lancaster subbasin where simulated hydraulic-head declines were the greatest to areas where declines were smallest. Neither a formal optimization algorithm nor water-rights allocations were considered when redistributing the pumpage. Results indicated that hydraulic heads near Palmdale, where the pumpage was reduced, would recover by about 200 ft compared to 2005 conditions, with only 30 ft of additional drawdown in the northwestern part of the Lancaster subbasin, where the pumpage was increased. The magnitude of the simulated additional land subsidence decreased slightly compared to the first, status quo, scenario but land subsidence continued to be simulated throughout most of the northern part of the Lancaster subbasin. The third scenario consisted of two artificial-recharge simulations along the Upper Amargosa Creek channel and at a site located north of Antelope Buttes. Results indicate that applying artificial recharge at these sites would yield continued drawdowns and associated land subsidence. However, the magnitudes of drawdown and subsidence would be smaller than those simulated in the status quo scenario, indicating that artificial-recharge operations in the Antelope Valley could be expected to reduce the magnitude and extent of continued water-level declines and associated land subsidence.

California↗
USGS2013

Groundwater contributions of flow, nitrate, and dissolved organic carbon to the lower San Joaquin River, California, 2006-08

The influence of groundwater on surface-water quality in the San Joaquin River, California, was examined for a 59-mile reach from the confluence with Salt Slough to Vernalis. The primary objective of this study was to quantify the rate of groundwater discharged to the lower San Joaquin River and the contribution of nitrate and dissolved organic carbon concentrations to the river. Multiple lines of evidence from four independent approaches were used to characterize groundwater contributions of nitrogen and dissolved organic carbon. Monitoring wells (in-stream and bank wells), streambed synoptic surveys (stream water and shallow groundwater), longitudinal profile surveys by boat (continuous water-quality parameters in the stream), and modeling (MODFLOW and VS2DH) provided a combination of temporal, spatial, quantitative, and qualitative evidence of groundwater contributions to the river and the associated quality. Monitoring wells in nested clusters in the streambed (in-stream wells) and on both banks (bank wells) along the river were monitored monthly from September 2006 to January 2009. Nitrate concentrations in the bank wells ranged from less than detection—that is, less than 0.01 milligrams per liter (mg/L) as nitrogen (N)—to approximately 13 mg/L as N. Nitrate was not detected at 17 of 26 monitoring wells during the study period. Dissolved organic carbon concentrations among monitoring wells were highly variable, but they generally ranged from 1 to 4 mg/L. In a previous study, 14 bank wells were sampled once in 1988 following their original installation. With few exceptions, specific conductivity and nitrate concentrations measured in this study were virtually identical to those measured 20 years ago. Streambed synoptic measurements were made by using a temporarily installed drive-point piezometer at 113 distinct transects across the stream during 4 sampling events. Nitrate concentrations exceeded the detection limit of 0.01 mg/L as N in 5 percent of groundwater samples collected from the in-stream wells as part of the synoptic surveys. Only 7 of the 113 cross-sectional transects had nitrate concentrations greater than 1 mg/L as N. In contrast, surface waters in the San Joaquin River tended to have nitrate concentrations in the 1–3 mg/L as N range. A zone of lower oxygen (less than 2 mg/L) in the streambed could limit nitrate contributions from regional groundwater flow because nitrate can be converted to nitrogen gas within this zone. Appreciable concentrations of ammonium (average concentration was 1.92 mg/L as N, and 95th percentile was 10.34 mg/L as N) in the shallow groundwater, believed to originate from anoxic mineralization of streambed sediments, could contribute nitrogen to the overlying stream as nitrate following in-stream nitrification, however. Dissolved organic carbon concentrations were highly variable in the shallow groundwater below the river (1 to 6 ft below streambed) and generally ranged between 1 and 5 mg/L, but had maximum concentrations in the 15–25 mg/L range. The longitudinal profile surveys were not particularly useful in identifying groundwater discharge areas. However, the longitudinal approach described in this report was useful as a baseline survey of measured water-quality parameters and for identifying tributary inflows that affect surface-water concentrations of nitrate. Results of the calibrated MODFLOW model indicated that the simulated groundwater discharge rate was approximately 1.0 cubic foot per second per mile (cfs/mi), and the predominant horizontal groundwater flow direction between the deep bank wells was westward beneath the river. The modeled (VS2DH) flux values (river gain versus river loss) were calculated for the irrigation and non-irrigation season, and these fluxes were an order of magnitude less than those from MODFLOW. During the irrigation season, the average river gain was 0.11 cfs/mi, and the average river loss was −0.05 cfs/mi. During the non-irrigation season, the average river gain was 0.10 cfs/mi, and the average river loss was -0.08 cfs/mi. Information on groundwater interactions and water quality collected for this study was used to estimate loads of nitrate and dissolved organic carbon from the groundwater to the San Joaquin River. Estimated loads of dissolved inorganic nitrogen and dissolved organic carbon were calculated by using concentrations measured during four streambed synoptic surveys and the estimated groundwater discharge rate to the San Joaquin River from MODFLOW of 1 cfs/mi. The estimated groundwater loads to the San Joaquin River for dissolved inorganic nitrogen and dissolved organic carbon were 300 and 350 kilograms per day, respectively. These loads represent 9 and 7 percent, respectively, of the estimated instantaneous surface-water loads for dissolved inorganic nitrogen and dissolved organic carbon at the most downstream site, Vernalis, measured during the four streambed synoptic surveys.

California↗
USGS2010

Mixing effects on apparent reaction rates and isotope fractionation during denitrification in a heterogeneous aquifer

Gradients in contaminant concentrations and isotopic compositions commonly are used to derive reaction parameters for natural attenuation in aquifers. Differences between field‐scale (apparent) estimated reaction rates and isotopic fractionations and local‐scale (intrinsic) effects are poorly understood for complex natural systems. For a heterogeneous alluvial fan aquifer, numerical models and field observations were used to study the effects of physical heterogeneity on reaction parameter estimates. Field measurements included major ions, age tracers, stable isotopes, and dissolved gases. Parameters were estimated for the O 2 reduction rate, denitrification rate, O 2 threshold for denitrification, and stable N isotope fractionation during denitrification. For multiple geostatistical realizations of the aquifer, inverse modeling was used to establish reactive transport simulations that were consistent with field observations and served as a basis for numerical experiments to compare sample‐based estimates of “apparent” parameters with “true“ (intrinsic) values. For this aquifer, non‐Gaussian dispersion reduced the magnitudes of apparent reaction rates and isotope fractionations to a greater extent than Gaussian mixing alone. Apparent and true rate constants and fractionation parameters can differ by an order of magnitude or more, especially for samples subject to slow transport, long travel times, or rapid reactions. The effect of mixing on apparent N isotope fractionation potentially explains differences between previous laboratory and field estimates. Similarly, predicted effects on apparent O 2 threshold values for denitrification are consistent with previous reports of higher values in aquifers than in the laboratory. These results show that hydrogeological complexity substantially influences the interpretation and prediction of reactive transport.

Water Resources Research↗
USGS2008

Simulations of ground-water flow and particle pathline analysis in the zone of contribution of a public-supply well in Modesto, eastern San Joaquin Valley, California

Shallow ground water in the eastern San Joaquin Valley is affected by high nitrate and uranium concentrations and frequent detections of pesticides and volatile organic compounds (VOC), as a result of ground-water development and intensive agricultural and urban land use. A single public-supply well was selected for intensive study to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Modesto area. A network of 23 monitoring wells was installed, and water and sediment samples were collected within the approximate zone of contribution of the public-supply well, to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state local ground-water-flow and transport model was developed to evaluate the age of ground water reaching the well and to evaluate the vulnerability of the well to nonpoint source input of nitrate and uranium. Particle tracking was used to compute pathlines and advective travel times in the ground-water flow model. The simulated ages of particles reaching the public-supply well ranged from 9 to 30,000 years, with a median of 54 years. The age of the ground water contributed to the public-supply well increased with depth below the water table. Measured nitrate concentrations, derived primarily from agricultural fertilizer, were highest (17 milligrams per liter) in shallow ground water and decreased with depth to background concentrations of less than 2 milligrams per liter in the deepest wells. Because the movement of water is predominantly downward as a result of ground-water development, and because geochemical conditions are generally oxic, high nitrate concentrations in shallow ground water are expected to continue moving downward without significant attenuation. Simulated long-term nitrate concentrations indicate that concentrations have peaked and will decrease in the public-supply well during the next 100 years because of the low nitrate concentrations in recharge beneath the urban area and the increasing proportion of urban-derived ground water reaching the well. The apparent lag time between peak input concentrations and peak concentrations in the well is about 20 to 30 years. Measured uranium concentrations were also highest (45 micrograms per liter) in shallow ground water, and decreased with depth to background concentrations of about 0.5 microgram per liter. Naturally-occurring uranium adsorbed to aquifer sediments is mobilized by oxygen-rich, high-alkalinity water. Alkalinity increased in shallow ground water in response to agricultural development. As ground-water pumping increased in the 1940s and 1950s, this alkaline water moved downward through the ground-water flow system, mobilizing the uranium adsorbed to aquifer sediments. Ground water with high alkalinity and high uranium concentrations is expected to continue to move deeper in the system, resulting in increased uranium concentrations with depth in ground water. Because alkalinity (and correspondingly uranium) concentrations were high in shallow ground water beneath both the urban and the agricultural land, long-term uranium concentrations in the public-supply well are expected to increase as the proportion of uranium-affected water contributed to the well increases. Assuming that the alkalinity near the water table remains the same, the simulation of long-term alkalinity in the public-supply well indicates that uranium concentrations in the public-supply well will likely approach the maximum contaminant level; however, the time to reach this level is more than 100 years because of the significant proportion of old, unaffected water at depth that is contributed to the public-supply well.

California↗
USGS2007

Simulation of Multiscale Ground-Water Flow in Part of the Northeastern San Joaquin Valley, California

The transport and fate of agricultural chemicals in a variety of environmental settings is being evaluated as part of the U.S. Geological Survey (USGS) National Water-Quality Assessment Program. One of the locations being evaluated is a 2,700-km2 (square kilometer) regional study area in the northeastern San Joaquin Valley surrounding the city of Modesto, an area dominated by irrigated agriculture in a semi-arid climate. Ground water is a key source of water for irrigation and public supply, and exploitation of this resource has altered the natural flow system. The aquifer system is predominantly alluvial, and an unconfined to semiconfined aquifer overlies a confined aquifer in the southwestern part of the study area; these aquifers are separated by the lacustrine Corcoran Clay. A regional-scale 16-layer steady-state model of ground-water flow in the aquifer system in the regional study area was developed to provide boundary conditions for an embedded 110-layer steady-state local-scale model of part of the aquifer system overlying the Corcoran Clay along the Merced River. The purpose of the local-scale model was to develop a better understanding of the aquifer system and to provide a basis for simulation of reactive transport of agricultural chemicals. The heterogeneity of aquifer materials was explicitly incorporated into the regional and local models using information from geologic and drillers? logs of boreholes. Aquifer materials were differentiated in the regional model by the percentage of coarse-grained sediments in a cell, and in the local model by four hydrofacies (sand, silty sand, silt, and clay). The calibrated horizontal hydraulic conductivity values of the coarse-grained materials in the zone above the Corcoran Clay in the regional model and of the sand hydrofacies used in the local model were about equal (30?80 m/d [meter per day]), and the vertical hydraulic conductivity values in the same zone of the regional model (median of 0.012 m/d), which is dominated by the finer-grained materials, were about an order of magnitude less than that for the clay hydrofacies in the local model. Data used for calibrating both models included long-term hourly water-level measurements in 20 short-screened wells installed by the USGS in the Modesto and Merced River areas. Additional calibration data for the regional model included water-level measurements in 11 wells upslope and 17 wells downslope from these areas. The root mean square error was 2.3 m (meter) for all wells in the regional model and 0.8 m for only the USGS wells; the associated average errors were 0.9 m and 0.3 m, respectively. The root mean square error for the 12 USGS wells along a transect in the local model area was 0.08 m; the average error was 0.0 m. Particle tracking was used with the local model to estimate the concentration of an environmental tracer, sulfur hexafluoride, in 10 USGS transect wells near the Merced River that were sampled for this constituent. Measured and estimated concentrations in the mid-depth and deepest wells, which would be most sensitive to errors in hydraulic conductivity estimates, were consistent. The combined results of particle tracking and sulfur hexafluoride analysis suggest that most water sampled from the transect wells was recharged less that 25 years ago.

Scientific Investigations Report↗
USGS2005

Calibration of a texture-based model of a ground-water flow system, western San Joaquin Valley, California

The occurrence of selenium in agricultural drain water from the western San Joaquin Valley, California, has focused concern on the semiconfined ground-water flow system, which is underlain by the Corcoran Clay Member of the Tulare Formation. A two-step procedure is used to calibrate a preliminary model of the system for the purpose of determining the steady-state hydraulic properties. Horizontal and vertical hydraulic conductivities are modeled as functions of the percentage of coarse sediment, hydraulic conductivities of coarse-textured (K coarse ) and fine-textured (K fine ) end members, and averaging methods used to calculate equivalent hydraulic conductivities. The vertical conductivity of the Corcoran (K corc ) is an additional parameter to be evaluated. In the first step of the calibration procedure, the model is run by systematically varying the following variables: (1) K coarse /K fine , (2) K coarse /K corc , and (3) choice of averaging methods in the horizontal and vertical directions. Root mean square error and bias values calculated from the model results are functions of these variables. These measures of error provide a means for evaluating model sensitivity and for selecting values of K coarse , K fine , and K corc for use in the second step of the calibration procedure. In the second step, recharge rates are evaluated as functions of K coarse , K corc , and a combination of averaging methods. The associated K fine values are selected so that the root mean square error is minimized on the basis of the results from the first step. The results of the two-step procedure indicate that the spatial distribution of hydraulic conductivity that best produces the measured hydraulic head distribution is created through the use of arithmetic averaging in the horizontal direction and either geometric or harmonic averaging in the vertical direction. The equivalent hydraulic conductivities resulting from either combination of averaging methods compare favorably to field- and laboratory-based values.

California↗
USGS2004

Water-Quality Assessment of the San Joaquin-Tulare Basins: Entering a new decade

In 1991, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey began to address the need for consistent and scientifically sound information for managing the Nation's water resources. The long-term goals of this program are to assess the status of the quality of freshwater streams and aquifers, to describe trends or changes in water quality over time, and to provide a sound understanding of the natural and human factors that affect the quality of these resources (Hirsch and others, 1988). Investigations are being conducted within major river basins and aquifer systems, or 'study units,' throughout the Nation to provide a framework for national and regional water-quality assessments. In 2001, the NAWQA Program began its second decade of intensive water-quality assessments. Forty-two of the original 59 study units (reduced by elimination or combination) are being revisited (Gilliom and others, 2001). The San Joaquin-Tulare Basins study unit (fig. 1), located in central California, was a part of the first decadal cycle of the Program investigations and remains in the second cycle.

California↗
USGS2004

Estimation of a water budget for 1972-2000 for the Grasslands Area, central part of the Western San Joaquin Valley, California

Equitable implementation of regulations restricting discharges from agricultural drains into the San Joaquin River requires a greater understanding of the influence of extreme precipitation events on the ground-water flow system. As part of a larger investigation, this study estimated ground-water recharge and ground-water pumpage, two important components of the water budget in the Grasslands drainage area in the central part of the western San Joaquin Valley, California, for the water years 1972 through 2000. These estimates will be used as inputs to a numerical simulation model of the regional ground-water flow system in the continuing investigation. Crop-acreage and surface-water delivery data were compiled for 14 water districts and 6 other areas comprising approximately 97 percent of the 600-square-mile study area. Little ground-water pumpage data exists for the study area. A climate-based approach was employed to estimate annual water-table recharge flux and ground-water pumpage for 11 water-budget areas. Ground-water pumpage was estimated from the residual irrigation demand after crop consumption of surface water. Estimated recharge flux to the water table for the entire study area averaged 0.8 ft/yr, and estimated ground-water pumpage per unit area for the entire study area averaged 0.5 ft/yr. Increased discharges from agricultural drains in the late 1990s may have been due partly to 4 years of high recharge from precipitation over the 6-year period from 1993 to 1998. Knowledge of the ratio of annual crop water demand to annual potential evapotranspiration, expressed as an aggregate crop coefficient, Kd, will facilitate estimation of annual water-budget components in future studies. Annual aggregate crop coefficients, calculated each year for the entire study area, were nearly constant at 0.59 from 1983 to 2000, and reasonably constant at 0.53 prior to 1983. The overall trend suggests continuous reductions in recharge from irrigation over time. This reduction is most likely due to gradual improvements in irrigation management. The recharge and pumpage estimates are both sensitive to the consumption-distribution ratio, and the pumpage estimate is sensitive to the cropped acreage.

Scientific Investigations Report↗
USGS2003

Simulation of ground-water flow and land subsidence in the Antelope Valley ground-water basin, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley ground-water basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, ground water provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most ground-water pumping in the valley occurs in the Antelope Valley ground-water basin, which includes the rapidly growing cities of Lancaster and Palmdale. Ground-water-level declines of more than 200 feet in some parts of the ground-water basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may continue to increase reliance on ground water. To better understand the ground-water flow system and to develop a tool to aid in effectively managing the water resources, a numerical model of ground-water flow and land subsidence in the Antelope Valley ground-water basin was developed using old and new geohydrologic information. The ground-water flow system consists of three aquifers: the upper, middle, and lower aquifers. The aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to ground-water development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Ground water flowed from the recharge areas to discharge areas around the playas where it discharged either from the aquifer system as evapotranspiration or from springs. Partial barriers to horizontal ground-water flow, such as faults, have been identified in the ground-water basin. Water-level declines owing to ground-water development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the ground-water system. Infiltration of return flows from agricultural irrigation has become an important source of recharge to the aquifer system. The ground-water flow model of the basin was discretized horizontally into a grid of 43 rows and 60 columns of square cells 1 mile on a side, and vertically into three layers representing the upper, middle, and lower aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915-95 using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley ground-water basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of ground-water flow and land subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Ground-water pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual pumpage peaked at 395,000 acre-feet (acre-ft) in 1952 and then declined because of declining agricultural production. Recharge from irrigation-return flows was estimated to be 30 percent of agricultural pumpage; the irrigation-return flows were simulated as recharge to the regional water table 10 years following application at land surface. The annual quantity of natural recharge initially was based on estimates from previous studies. During model calibration, natural recharge was reduced from the initial

Water-Resources Investigations Report↗
USGS2003

Determination of specific yield and water-table changes using temporal microgravity surveys collected during the second injection, storage, and recovery test at Lancaster, Antelope Valley, California, November 1996 through April 1997

To evaluate the feasibility of artificially recharging the ground-water system in the Lancaster area of the Antelope Valley, California, the U.S. Geological Survey, in cooperation with the Los Angeles County Department of Public Works and the Antelope Valley-East Kern Water Agency, conducted a series of injection, storage, and recovery tests between September 1995 and September 1998. A key component of this study was to measure the response of the water table to injection, which was difficult because the water table averaged 300 feet below land surface. Rather than install many expensive piezometers, microgravity surveys were conducted to determine specific yield and to measure the development of a ground-water mound during the injection of about 1,050 acre-feet of fresh water into an alluvial-aquifer system. The surveys were done prior to, during, and near the end of a 5-month injection period (November 12, 1996, to April 17, 1997). Results of the surveys indicate increases in gravity of as much as 66 microgals between a bedrock reference station and 20 gravity stations within a 1-square-mile area surrounding the injection site. The changes were assumed to have been caused by changes in the ground-water elevation. Gravity and ground-water levels were measured simultaneously at an existing well (7N/12W-34B1). The coupled measurements were used to calculate a specific yield of 0.13 for the alluvial aquifer near the well. To determine the gravitational effect of the injection mound on the gravity measurements made near well 7N/12W-34B1, a two-dimensional gravity model was used. Results of the model simulation show that the effect on gravity associated with the mass of the injection mound was minor and thus had a negligible effect on the calculation of specific yield. The specific yield of 0.13, therefore, was used to infer water-level changes at other gravity stations within the study area. The gravity-derived water-level changes were compared with simulated water-table changes.

Water-Resources Investigations Report↗
USGS2003

Analysis of tests of subsurface injection, storage, and recovery of freshwater in Lancaster, Antelope Valley, California

Ground-water levels in Lancaster, California, declined more than 200 feet during the 20th century, resulting in reduced ground-water supplies and more than 6 feet of land subsidence. Facing continuing population growth, water managers are seeking solutions to these problems. Injection of imported, treated fresh water into the aquifer system when it is most available and least expensive, for later use during high-demand periods, is being evaluated as part of a management solution. The U.S. Geological Survey, in cooperation with the Los Angeles County Department of Public Works and the Antelope Valley-East Kern Water Agency, monitored a pilot injection program, analyzed the hydraulic and subsidence-related effects of injection, and developed a simulation/optimization model to help evaluate the effectiveness of using existing and proposed wells in an injection program for halting the decline of ground-water levels and avoiding future land subsidence while meeting increasing ground-water demand. A variety of methods were used to measure aquifer-system response to injection. Water levels were measured continuously in nested (multi-depth) piezometers and monitoring wells and periodically in other wells that were within several miles of the injection site. Microgravity surveys were done to estimate changes in the elevation of the water table in the absence of wells and to estimate specific yield. Aquifer-system deformation was measured directly and continuously using a dual borehole extensometer and indirectly using continuous Global Positioning System (GPS), first-order spirit leveling, and an array of tiltmeters. The injected water and extracted water were sampled periodically and analyzed for constituents, including chloride and trihalomethanes. Measured injection rates of about 750 gallons per minute (gal/min) per well at the injection site during a 5-month period showed that injection at or above the average extraction rates at that site (about 800 gal/min) was hydraulically feasible. Analyses of these data took many forms. Coupled measurements of gravity and water-level change were used to estimate the specific yield near the injection wells, which, in turn, was used to estimate areal water-table changes from distributed measurements of gravity change. Values of the skeletal components of aquifer-system storage, which are key subsidence-related characteristics of the system, were derived from continuous measurements of water levels and aquifer-system deformation. A numerical model of ground-water flow was developed for the area surrounding Lancaster and used to estimate horizontal and vertical hydraulic conductivities. A chemical mass balance was done to estimate the recovery of injected water. The ground-water-flow model was used to project changes in ground-water levels for 10 years into the future, assuming no injection, no change in pumping distribution, and forecasted increases in ground-water demand. Simulated ground-water levels decreased throughout the Lancaster area, suggesting that land subsidence would continue as would the depletion of ground-water supplies and an associated loss of well production capacity. A simulation/optimization model was developed to help identify optimal injection and extraction rates for 16 existing and 13 proposed wells to avoid future land subsidence and to minimize loss of well production capacity while meeting increasing ground-water demands. Results of model simulations suggest that these objectives can be met with phased installation of the proposed wells during the 10-year period. Water quality was not considered in the optimization, but chemical-mass-balance results indicate that a sustained injection program likely would have residual effects on the chemistry of ground water.

Califronia↗