Search USGS⌕ Search

Geology topics

L. E. Flint

Publications and source records attributed to L. E. Flint.

17 recordsLinked to original sources
USGS2012

A basin-scale approach for assessing water resources in a semiarid environment: San Diego region, California and Mexico

Many basins throughout the world have sparse hydrologic and geologic data, but have increasing demands for water and a commensurate need for integrated understanding of surface and groundwater resources. This paper demonstrates a methodology for using a distributed parameter water-balance model, gaged surface-water flow, and a reconnaissance-level groundwater flow model to develop a first-order water balance. Flow amounts are rounded to the nearest 5 million cubic meters per year. The San Diego River basin is 1 of 5 major drainage basins that drain to the San Diego coastal plain, the source of public water supply for the San Diego area. The distributed parameter water-balance model (Basin Characterization Model) was run at a monthly timestep for 1940–2009 to determine a median annual total water inflow of 120 million cubic meters per year for the San Diego region. The model was also run specifically for the San Diego River basin for 1982–2009 to provide constraints to model calibration and to evaluate the proportion of inflow that becomes groundwater discharge, resulting in a median annual total water inflow of 50 million cubic meters per year. On the basis of flow records for the San Diego River at Fashion Valley (US Geological Survey gaging station 11023000), when corrected for upper basin reservoir storage and imported water, the total is 30 million cubic meters per year. The difference between these two flow quantities defines the annual groundwater outflow from the San Diego River basin at 20 million cubic meters per year. These three flow components constitute a first-order water budget estimate for the San Diego River basin. The ratio of surface-water outflow and groundwater outflow to total water inflow are 0.6 and 0.4, respectively. Using total water inflow determined using the Basin Characterization Model for the entire San Diego region and the 0.4 partitioning factor, groundwater outflow from the San Diego region, through the coastal plain aquifer to the Pacific Ocean, is calculated to be approximately 50 million cubic meters per year. The area-scale assessment of water resources highlights several hydrologic features of the San Diego region. Groundwater recharge is episodic; the Basin Characterization Model output shows that 90 percent of simulated recharge occurred during 3 percent of the 1982–2009 period. The groundwater aquifer may also be quite permeable. A reconnaissance-level groundwater flow model for the San Diego River basin was used to check the water budget estimates, and the basic interaction of the surface-water and groundwater system, and the flow values, were found to be reasonable. Horizontal hydraulic conductivity values of the volcanic and metavolcanic bedrock in San Diego region range from 1 to 10 m per day. Overall, results establish an initial hydrologic assessment formulated on the basis of sparse hydrologic data. The described flow variability, extrapolation, and unique characteristics represent a realistic view of current (2012) hydrologic understanding for the San Diego region.

California↗
USGS2012

A method for physically based model analysis of conjunctive use in response to potential climate changes

Potential climate change effects on aspects of conjunctive management of water resources can be evaluated by linking climate models with fully integrated groundwater-surface water models. The objective of this study is to develop a modeling system that links global climate models with regional hydrologic models, using the California Central Valley as a case study. The new method is a supply and demand modeling framework that can be used to simulate and analyze potential climate change and conjunctive use. Supply-constrained and demand-driven linkages in the water system in the Central Valley are represented with the linked climate models, precipitation-runoff models, agricultural and native vegetation water use, and hydrologic flow models to demonstrate the feasibility of this method. Simulated precipitation and temperature were used from the GFDL-A2 climate change scenario through the 21st century to drive a regional water balance mountain hydrologic watershed model (MHWM) for the surrounding watersheds in combination with a regional integrated hydrologic model of the Central Valley (CVHM). Application of this method demonstrates the potential transition from predominantly surface water to groundwater supply for agriculture with secondary effects that may limit this transition of conjunctive use. The particular scenario considered includes intermittent climatic droughts in the first half of the 21st century followed by severe persistent droughts in the second half of the 21st century. These climatic droughts do not yield a valley-wide operational drought but do cause reduced surface water deliveries and increased groundwater abstractions that may cause additional land subsidence, reduced water for riparian habitat, or changes in flows at the Sacramento-San Joaquin River Delta. The method developed here can be used to explore conjunctive use adaptation options and hydrologic risk assessments in regional hydrologic systems throughout the world.

Water Resources Research↗
USGS2008

A basin-scale approach to estimating stream temperatures of tributaries to the lower Klamath River, California

Stream temperature is an important component of salmonid habitat and is often above levels suitable for fish survival in the Lower Klamath River in northern California. The objective of this study was to provide boundary conditions for models that are assessing stream temperature on the main stem for the purpose of developing strategies to manage stream conditions using Total Maximum Daily Loads. For model input, hourly stream temperatures for 36 tributaries were estimated for 1 Jan. 2001 through 31 Oct. 2004. A basin-scale approach incorporating spatially distributed energy balance data was used to estimate the stream temperatures with measured air temperature and relative humidity data and simulated solar radiation, including topographic shading and corrections for cloudiness. Regression models were developed on the basis of available stream temperature data to predict temperatures for unmeasured periods of time and for unmeasured streams. The most significant factor in matching measured minimum and maximum stream temperatures was the seasonality of the estimate. Adding minimum and maximum air temperature to the regression model improved the estimate, and air temperature data over the region are available and easily distributed spatially. The addition of simulated solar radiation and vapor saturation deficit to the regression model significantly improved predictions of maximum stream temperature but was not required to predict minimum stream temperature. The average SE in estimated maximum daily stream temperature for the individual basins was 0.9 ?? 0.6??C at the 95% confidence interval. Copyright ?? 2008 by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America. All rights reserved.

Journal of Environmental Quality↗
USGS2008

Modeling soil moisture processes and recharge under a melting snowpack

Recharge into granitic bedrock under a melting snowpack is being investigated as part of a study designed to understand hydrologic processes involving snow at Yosemite National Park in the Sierra Nevada Mountains of California. Snowpack measurements, accompanied by water content and matric potential measurements of the soil under the snowpack, allowed for estimates of infiltration into the soil during snowmelt and percolation into the bedrock. During portions of the snowmelt period, infiltration rates into the soil exceeded the permeability of the bedrock and caused ponding to be sustained at the soil-bedrock interface. During a 5-d period with little measured snowmelt, drainage of the ponded water into the underlying fractured granitic bedrock was estimated to be 1.6 cm d?1, which is used as an estimate of bedrock permeability. The numerical simulator TOUGH2 was used to reproduce the field data and evaluate the potential for vertical flow into the fractured bedrock or lateral flow at the bedrock-soil interface. During most of the snowmelt season, the snowmelt rates were near or below the bedrock permeability. The field data and model results support the notion that snowmelt on the shallow soil overlying low permeability bedrock becomes direct infiltration unless the snowmelt rate greatly exceeds the bedrock permeability. Late in the season, melt rates are double that of the bedrock permeability (although only for a few days) and may tend to move laterally at the soil-bedrock interface downgradient and contribute directly to streamflow. ?? Soil Science Society of America.

Conference Paper↗
USGS2006

Reconstructing depositional processes and history from reservoir stratigraphy: Englebright Lake, Yuba River, northern California

Reservoirs provide the opportunity to link watershed history with its stratigraphic record. We analyze sediment cores from a northern California reservoir in the context of hydrologic history, watershed management, and depositional processes. Observations of recent depositional patterns, sediment‐transport calculations, and 137 Cs geochronology support a conceptual model in which the reservoir delta progrades during floods of short duration (days) and is modified during prolonged (weeks to months) drawdowns that rework topset beds and transport sand from topsets to foresets. Sediment coarser than 0.25–0.5 mm deposits in foresets and topsets, and finer material falls out of suspension as bottomset beds. Simple hydraulic calculations indicate that fine sand (0.063–0.5 mm) is transported into the distal bottomset area only during floods. The overall stratigraphy suggests that two phases of delta building occurred in the reservoir. The first, from dam construction in 1940 to 1970, was heavily influenced by annual, prolonged >20 m drawdowns of the water level. The second, built on top of the first, reflects sedimentation from 1970 to 2002 when the influence of drawdowns was less. Sedimentation rates in the central part of the reservoir have declined ∼25% since 1970, likely reflecting a combination of fewer large floods, changes in watershed management, and winnowing of stored hydraulic mining sediment.

Journal of Geophysical Research F: Earth Surface↗
USGS2006

Global synthesis of groundwater recharge in semiarid and arid regions

Global synthesis of the findings from ∼140 recharge study areas in semiarid and arid regions provides important information on recharge rates, controls, and processes, which are critical for sustainable water development. Water resource evaluation, dryland salinity assessment (Australia), and radioactive waste disposal (US) are among the primary goals of many of these recharge studies. The chloride mass balance (CMB) technique is widely used to estimate recharge. Average recharge rates estimated over large areas (40–374 000 km 2 ) range from 0·2 to 35 mm year −1 , representing 0·1–5% of long-term average annual precipitation. Extreme local variability in recharge, with rates up to ∼720 m year −1 , results from focussed recharge beneath ephemeral streams and lakes and preferential flow mostly in fractured systems. System response to climate variability and land use/land cover (LU/LC) changes is archived in unsaturated zone tracer profiles and in groundwater level fluctuations. Inter-annual climate variability related to El Niño Southern Oscillation (ENSO) results in up to three times higher recharge in regions within the SW US during periods of frequent El Niños (1977–1998) relative to periods dominated by La Niñas (1941–1957). Enhanced recharge related to ENSO is also documented in Argentina. Climate variability at decadal to century scales recorded in chloride profiles in Africa results in recharge rates of 30 mm year −1 during the Sahel drought (1970–1986) to 150 mm year −1 during non-drought periods. Variations in climate at millennial scales in the SW US changed systems from recharge during the Pleistocene glacial period (≥10 000 years ago) to discharge during the Holocene semiarid period. LU/LC changes such as deforestation in Australia increased recharge up to about 2 orders of magnitude. Changes from natural grassland and shrublands to dryland (rain-fed) agriculture altered systems from discharge (evapotranspiration, ET) to recharge in the SW US. The impact of LU change was much greater than climate variability in Niger (Africa), where replacement of savanna by crops increased recharge by about an order of magnitude even during severe droughts. Sensitivity of recharge to LU/LC changes suggests that recharge may be controlled through management of LU. In irrigated areas, recharge varies from 10 to 485 mm year −1 , representing 1–25% of irrigation plus precipitation. However, irrigation pumpage in groundwater-fed irrigated areas greatly exceeds recharge rates, resulting in groundwater mining. Increased recharge related to cultivation has mobilized salts that accumulated in the unsaturated zone over millennia, resulting in widespread groundwater and surface water contamination, particularly in Australia. The synthesis of recharge rates provided in this study contains valuable information for developing sustainable groundwater resource programmes within the context of climate variability and LU/LC change.

Hydrological Processes↗
USGS2003

Use of porosity to estimate hydraulic properties of volcanic tuffs

Correlations of hydraulic properties with easily measured physical properties are useful for purposes of site characterization in heterogeneous sites. Approximately 600 samples of volcanic rocks from Yucca Mountain, Nevada, representing lithologies with a large range of hydraulic properties, were analyzed to develop correlations of effective porosity with saturated hydraulic conductivity and moisture-retention curve-fit parameters that relate to lithologies of varying depositional history and alteration processes. Effective porosity, ??e, defined as the porosity calculated using drying at a relative humidity of -70 MPa, is used in a generalized Kozeny-Carman equation to predict saturated hydraulic conductivity, Ks = b??en, where b and n are constants. The entire dataset has an R2 of 0.36. When samples are grouped according to general lithology, correlations result in an R2 of 0.71 for the crystallized/vitric samples, 0.24 for samples with mineral alteration, and 0.34 for samples with microfractures, thus increasing the predictive capability over that of the total dataset. Published by Elsevier Science Ltd.

Advances in Water Resources↗
USGS2002

Estimating recharge at Yucca Mountain, Nevada, USA: Comparison of methods

Obtaining values of net infiltration, groundwater travel time, and recharge is necessary at the Yucca Mountain site, Nevada, USA, in order to evaluate the expected performance of a potential repository as a containment system for high-level radioactive waste. However, the geologic complexities of this site, its low precipitation and net infiltration, with numerous mechanisms operating simultaneously to move water through the system, provide many challenges for the estimation of the spatial distribution of recharge. A variety of methods appropriate for arid environments has been applied, including water-balance techniques, calculations using Darcy's law in the unsaturated zone, a soil-physics method applied to neutron-hole water-content data, inverse modeling of thermal profiles in boreholes extending through the thick unsaturated zone, chloride mass balance, atmospheric radionuclides, and empirical approaches. These methods indicate that near-surface infiltration rates at Yucca Mountain are highly variable in time and space, with local (point) values ranging from zero to several hundred millimeters per year. Spatially distributed net-infiltration values average 5 mm/year, with the highest values approaching 20 mm/year near Yucca Crest. Site-scale recharge estimates range from less than 1 to about 12 mm/year. These results have been incorporated into a site-scale model that has been calibrated using these data sets that reflect infiltration processes acting on highly variable temporal and spatial scales. The modeling study predicts highly non-uniform recharge at the water table, distributed significantly differently from the non-uniform infiltration pattern at the surface.

Hydrogeology Journal↗
USGS2002

Preliminary estimates of spatially distributed net infiltration and recharge for the Death Valley region, Nevada-California

A three-dimensional ground-water flow model has been developed to evaluate the Death Valley regional flow system, which includes ground water beneath the Nevada Test Site. Estimates of spatially distributed net infiltration and recharge are needed to define upper boundary conditions. This study presents a preliminary application of a conceptual and numerical model of net infiltration. The model was developed in studies at Yucca Mountain, Nevada, which is located in the approximate center of the Death Valley ground-water flow system. The conceptual model describes the effects of precipitation, runoff, evapotranspiration, and redistribution of water in the shallow unsaturated zone on predicted rates of net infiltration; precipitation and soil depth are the two most significant variables. The conceptual model was tested using a preliminary numerical model based on energy- and water-balance calculations. Daily precipitation for 1980 through 1995, averaging 202 millimeters per year over the 39,556 square kilometers area of the ground-water flow model, was input to the numerical model to simulate net infiltration ranging from zero for a soil thickness greater than 6 meters to over 350 millimeters per year for thin soils at high elevations in the Spring Mountains overlying permeable bedrock. Estimated average net infiltration over the entire ground-water flow model domain is 7.8 millimeters per year. To evaluate the application of the net-infiltration model developed on a local scale at Yucca Mountain, to net-infiltration estimates representing the magnitude and distribution of recharge on a regional scale, the net-infiltration results were compared with recharge estimates obtained using empirical methods. Comparison of model results with previous estimates of basinwide recharge suggests that the net-infiltration estimates obtained using this model may overestimate recharge because of uncertainty in modeled precipitation, bedrock permeability, and soil properties for locations such as the Spring Mountains. Although this model is preliminary and uncalibrated, it provides a first approximation of the spatial distribution of net infiltration for the Death Valley region under current climatic conditions.

California, Nevada↗
USGS2001

Hydrology of Yucca Mountain, Nevada

Yucca Mountain, located in southern Nevada in the Mojave Desert, is being considered as a geologic repository for high-level radioactive waste. Although the site is arid, previous studies indicate net infiltration rates of 5-10 mm yr-1 under current climate conditions. Unsaturated flow of water through the mountain generally is vertical and rapid through the fractures of the welded tuffs and slow through the matrix of the nonwelded tuffs. The vitric-zeolitic boundary of the nonwelded tuffs below the potential repository, where it exists, causes perching and substantial lateral flow that eventually flows through faults near the eastern edge of the potential repository and recharges the underlying groundwater system. Fast pathways are located where water flows relatively quickly through the unsaturated zone to the water table. For the bulk of the water a large part of the travel time from land surface to the potential repository horizon (~300 m below land surface) is through the interlayered, low fracture density, nonwelded tuff where flow is predominately through the matrix. The unsaturated zone at Yucca Mountain is being modeled using a three-dimensional, dual-continuum numerical model to predict the results of measurements and observations in new boreholes and excavations. The interaction between experimentalists and modelers is providing confidence in the conceptual model and the numerical model and is providing researchers with the ability to plan further testing and to evaluate the usefulness or necessity of further data collection.

Reviews of Geophysics↗
USGS1998

Characterization of hydrogeologic units using matrix properties, Yucca Mountain, Nevada

Determination of the suitability of Yucca Mountain, in southern Nevada, as a geologic repository for high-level radioactive waste requires the use of numerical flow and transport models. Input for these models includes parameters that describe hydrologic properties and the initial and boundary conditions for all rock materials within the unsaturated zone, as well as some of the upper rocks in the saturated zone. There are 30 hydrogeologic units in the unsaturated zone, and each unit is defined by limited ranges where a discrete volume of rock contains similar hydrogeologic properties. These hydrogeologic units can be easily located in space by using three-dimensional lithostratigraphic models based on relation- ships of the properties with the lithostratigraphy. Physical properties of bulk density, porosity, and particle density; flow properties of saturated hydraulic conductivity and moisture-retention characteristics; and the state variables (variables describing the current state of field conditions) of saturation and water potential were determined for each unit. Units were defined using (1) a data base developed from 4,892 rock samples collected from the coring of 23 shallow and 8 deep boreholes, (2) described lithostratigraphic boundaries and corresponding relations to porosity, (3) recognition of transition zones with pronounced changes in properties over short vertical distances, (4) characterization of the influence of mineral alteration on hydrologic properties such as permeability and moisture-retention characteristics, and (5) a statistical analysis to evaluate where boundaries should be adjusted to minimize the variance within layers. This study describes the correlation of hydrologic properties to porosity, a property that is well related to the lithostratigraphy and depositional and cooling history of the volcanic deposits and can, therefore, be modeled to be distributed laterally. Parameters of the hydrogeologic units developed in this study and the relation of flow properties to porosity that are described can be used to produce detailed and accurate representations of the core-scale hydrologic processes ongoing at Yucca Mountain.

Nevada↗
USGS1996

Improved method for measuring water imbibition rates on low-permeability porous media

Existing methods for measuring water imbibition rates are inadequate when imbibition rates are small (e.g., clay soils and many igneous rocks). We developed an improved laboratory method for performing imbibition measurements on soil or rock cores with a wide range of hydraulic properties. Core specimens are suspended from an electronic strain gauge (load cell) in a closed chamber while maintaining the lower end of the core in contact with a free water surface in a constant water level reservoir. The upper end of the core is open to the atmosphere. During imbibition, mass increase of the core is recorded continuously by a datalogger that converts the load cell voltage signal into mass units using a calibration curve. Computer automation allows imbibition rate measurement on as many as eight cores simultaneously and independently. Performance of each component of the imbibition apparatus was evaluated using a set of rock cores (2.5 cm in diameter and 2–5 cm in length) from a signle lithostratigraphic unit composed of non-to-moderately welded ash-flow tuff (a glass-rich pyroclastic rock partially fused by heat and pressure) with porosities ranging from 0.094 to 0.533 m 3 m -3 . Reproducibility of sample handling and testing procedures was demonstrated using replicate measurements. Precision and accuracy of load cell measurements were assessed using mass balance calculations and indicated agreement within a few tenths of a percent of total mass. Computed values of sorptivity, S , ranged from 8.83 × 10 -6 to 4.55 × 10 -4 m s -0.5 . The developed method should prove useful for measuring imbibition rates on a wide range of porous materials.

Soil Science Society of America Journal↗
USGS1995

Physical and hydrologic properties of outcrop samples from a nonwelded to welded tuff transition, Yucca Mountain, Nevada

Quantitative material-property data are needed to describe lateral and vertical spatial variability of physical and hydrologic properties and to model ground-water flow and radionuclide transport at the potential Yucca Mountain nuclear-waste repository site in Nevada. As part of ongoing site characterization studies of Yucca Mountain directed toward this understanding of spatial variability, laboratory measurements of porosity, bull* and particle density, saturated hydraulic conductivity, and sorptivity have been obtained for a set of outcrop samples that form a systematic,two dimensional grid that covers a large exposure of the basal Tiva Canyon Tuff of the Paintbrush Group of Miocene age at Yucca Mountain. The samples form a detailed vertical grid roughly parallel to the transport direction of the parent ash flows, and they exhibit material-property varia- tions in an interval of major lithologic change overlying a potential nuclear-waste repository at Yucca Mountain. The observed changes in hydrologic properties were systematic and consistent with the changes expected for the nonwelded to welded transition at the base of a major ash-flow sequence. Porosity, saturated hydraulic conductivity, and sorptivity decreased upward from the base of the Tiva Canyon Tuff, indicating the progressive compaction of ash- rich volcanic debris and the onset of welding with increased overburden pressure from the accumulating ash-flow sheet. The rate of decrease in the values of these material properties varied with vertical position within the transition interval. In contrast, bulk-density values increased upward, a change that also is consistent with progressive compaction and the onset of welding. Particle-density values remained almost constant throughout the transition interval, probably indicating compositional (chemical) homogeneity.

Water-Resources Investigations Report↗