Search USGSSearch

USGS · 70019911

Quantifying macropore recharge: Examples from a semi-arid area

Abstract

The purpose of this paper is to illustrate the significantly increased resolution of determining macropore recharge by combining physical, chemical, and isotopic methods of analysis. Techniques for quantifying macropore recharge were developed for both small-scale (1 to 10 km2) and regional-scale areas in and semi-arid areas. The Southern High Plains region of Texas and New Mexico was used as a representative field site to test these methods. Macropore recharge in small-scale areas is considered to be the difference between total recharge through floors of topographically dosed basins and interstitial recharge through the same area. On the regional scale, macropore recharge was considered to be the difference between regional average annual recharge and interstitial recharge measured in the unsaturated zone. Stable isotopic composition of ground water and precipitation was used us an independent estimate of macropore recharge on the regional scale. Results of this analysis suggest that in the Southern High Plains recharge flux through macropores is between 60 and 80 percent of the total 11 mm/y. Between 15 and 35 percent of the recharge occurs by interstitial recharge through the basin floors. Approximately 5 percent of the total recharge occurs as either interstitial or matrix recharge between the basin floors, representing approximately 95 percent of the area. The approach is applicable to other arid and semi-arid areas that focus rainfall into depressions or valleys.The purpose of this paper is to illustrate the significantly increased resolution of determining macropore recharge by combining physical, chemical, and isotopic methods of analysis. Techniques for quantifying macropore recharge were developed for both small-scale (1 to 10 km2) and regional-scale areas in arid and semi-arid areas. The Southern High Plains region of Texas and New Mexico was used as a representative field site to test these methods. Macropore recharge in small-scale areas is considered to be the difference between total recharge through floors of topographically closed basins and interstitial recharge through the same area. On the regional scale, macropore recharge was considered to be the difference between regional average annual recharge and interstitial recharge measured in the unsaturated zone. Stable isotopic composition of ground water and precipitation was used as an independent estimate of macropore recharge on the regional scale. Results of this analysis suggest that in the Southern High Plains recharge flux through macropores is between 60 and 80 percent of the total 11 mm/y. Between 15 and 35 percent of the recharge occurs by interstitial recharge through the basin floors. Approximately 5 percent of the total recharge occurs as either interstitial or matrix recharge between the basin floors, representing approximately 95 percent of the area. The approach is applicable to other arid and semi-arid areas that focus rainfall into depressions or valleys.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

W.W. Wood, Ken A. Rainwater, D.B. Thompson. 2005-08-04. Quantifying macropore recharge: Examples from a semi-arid area. https://doi.org/10.1111/j.1745-6584.1997.tb00182.x

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A “semi-structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water