Steven P. Phillips

DETECTION OF AQUIFER SYSTEM COMPACTION AND LAND SUBSIDENCE USING INTERFEROMETRIC SYNTHETIC APERTURE RADAR, ANTELOPE VALLEY, MOJAVE DESERT, CALIFORNIA

Interferometric synthetic aperture radar (InSAR) has great potential to detect and quantify land subsidence caused by aquifer system compaction. InSAR maps with high spatial detail and resolution of range displacement (±10 mm in change of land surface elevation) were developed for a groundwater basin (∼103 km2) in Antelope Valley, California, using radar data collected from the ERS-1 satellite. These data allow comprehensive comparison between recent (1993–1995) subsidence patterns and those detected historically (1926–1992) by more traditional methods. The changed subsidence patterns are generally compatible with recent shifts in land and water use. The InSAR-detected patterns are generally consistent with predictions based on a coupled model of groundwater flow and aquifer system compaction. The minor inconsistencies may reflect our imperfect knowledge of the distribution and properties of compressible sediments. When used in conjunction with coincident measurements of groundwater levels and other geologic information, InSAR data may be useful for constraining parameter estimates in simulations of aquifer system compaction.

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

[1] 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 O2 reduction rate, denitrification rate, O2 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 O2 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.

Limited Occurrence of Denitrification in Four Shallow Aquifers in Agricultural Areas of the United States

The ability of natural attenuation to mitigate agricultural nitrate contamination in recharging aquifers was investigated in four important agricultural settings in the United States. The study used laboratory analyses, field measurements, and flow and transport modeling for monitoring well transects (0.5 to 2.5 km in length) in the San Joaquin watershed, California, the Elkhorn watershed, Nebraska, the Yakima watershed, Washington, and the Chester watershed, Maryland. Ground water analyses included major ion chemistry, dissolved gases, nitrogen and oxygen stable isotopes, and estimates of recharge date. Sediment analyses included potential electron donors and stable nitrogen and carbon isotopes. Within each site and among aquifer-based medians, dissolved oxygen decreases with ground water age, and excess N(2) from denitrification increases with age. Stable isotopes and excess N(2) imply minimal denitrifying activity at the Maryland and Washington sites, partial denitrification at the California site, and total denitrification across portions of the Nebraska site. At all sites, recharging electron donor concentrations are not sufficient to account for the losses of dissolved oxygen and nitrate, implying that relict, solid phase electron donors drive redox reactions. Zero-order rates of denitrification range from 0 to 0.14 micromol N L(-1)d(-1), comparable to observations of other studies using the same methods. Many values reported in the literature are, however, orders of magnitude higher, which is attributed to a combination of method limitations and bias for selection of sites with rapid denitrification. In the shallow aquifers below these agricultural fields, denitrification is limited in extent and will require residence times of decades or longer to mitigate modern nitrate contamination.

Alternative to Agricultural Drains in California's San Joaquin Valley: Results of a Regional‐Scale Hydrogeologic Approach

The occurrence of selenium in agricultural drainage water derived from the central part of the western San Joaquin Valley has focused concern on alternatives to agricultural drains for managing shallow, poor-quality groundwater. A transient, three-dimensional simulation model was developed to evaluate the response of the water table to alternatives that affect recharge to or discharge from the groundwater flow system. The modeled area is 551 mi2 (1 mi2 = 2.59 km2) and includes both the semiconfined and confined zones above and below the Corcoran Clay Member of the Tulare Formation of Pleistocene age. The simulation model was calibrated using hydrologic data from 1972 to 1988, and was extended to the year 2040 to forecast for various management alternatives, including maintenance of present practices, land retirement, reduced recharge, increased groundwater pumping, and combinations of these alternatives. Maintenance of present practices results in a worsening of the situation: the total area subject to bare-soil evaporation increases from 224 mi2 in 1990 to 344 mi2 in 2040, and drain flow increases from 25,000 ac ft/yr (1 ac ft = 1234 m3) to 28,000 ac ft/yr. Although land retirement results in elimination of bare-soil evaporation and drain flow in the areas retired, it has little to no effect in adjacent areas. In contrast, regional-scale changes in recharge and pumping are effective for regional management. The area subject to bare-soil evaporation can be reduced to 78 mi2, and drain flow to 8000 ac ft/yr if (1) recharge is reduced by 15% (26,000 ac ft/yr) in areas that currently use surface and groundwater (362 mi2); (2) recharge is reduced by 40% (28,000 ac ft/yr) in areas that currently use only surface water (137 mi2); and (3) pumping rates are uniformly incremented by 0.5 ft/yr (160,000 ac ft/yr) in both areas. If these water budget changes were to be implemented in the study area, and in adjacent areas with similiar Hydrogeologic characteristics, then approximately 400,000 ac ft/yr of surface water would be made available. Thus a shift in the hydrologic budget in the central part of the western San Joaquin Valley improves the prospects for sustaining agriculture in the area, and could provide substantial water resources for other uses.

Distribution of wells in the central part of the western San Joaquin Valley, California

Information from 5,860 wells in the central part of the western San Joaquin Valley, California, was collected from several sources and compiled into a common data base. Only 2,547 wells had sufficient information for classification into four categories based on the hydrogeology: wells perforated in the semiconfined zone at depths less than or equal to 50 feet, wells perforated in the semiconfined zone at depths greater than 50 feet, wells perforated in the semiconfined and confined zones, and wells perforated only in the confined zone. Additionally, wells perforated in the semiconfined zone at depths greater than 50 feet were classified by the type of deposits in which they were perforated (Coast Range alluvium or Sierran sand). A computerized data base system was developed to manage well information and to facilitate characterizing the nature and distribution of the wells. Wells perforated in the semiconfined zone at depths less than or equal to 50 feet are evenly distributed over part of the study area underlain by shallow ground water. These wells generally are used as observation wells. Most wells perforated in the semiconfined zone at depths greater than 50 feet are perforated in the Sierran sand. This concentration of wells perforated in the Sierran sand indicates a tendency for using the Sierran sand, where it exists, as a source of water. There are 533 wells perforated in both the semiconfined and confined zones and 410 wells perforated only in the confined zone. Most of these wells are upslope of the valley trough in areas where the Sierran sand is not present. Wells perforated only in the confined zone are concentrated near the creeks.

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

.........................................................................................................................................................