Geology topics

Brian R. Clark

Publications and source records attributed to Brian R. Clark.

At least 19 recordsLinked to original sources

USGS2021

Effects of climate and land-use change on thermal springs recharge—A system-based coupled surface-water and groundwater-flow model for Hot Springs National Park, Arkansas

A three-dimensional hydrogeologic framework of the Hot Springs anticlinorium beneath Hot Springs National Park, Arkansas, was constructed to represent the complex hydrogeology of the park and surrounding areas to depths exceeding 9,000 feet below ground surface. The framework, composed of 6 rock formations and 1 vertical fault emplaced beneath the thermal springs, was discretized into 19 layers, 429 rows, and 576 columns and incorporated into a 3-dimensional steady-state groundwater-flow model constructed in MODFLOW-2005. Historical daily mean thermal spring flows were simulated for one stress period of approximately 34 years (1980–2014), chosen to represent the period of record for historical climate data used in the quantification of the boundary conditions. The groundwater-flow model was manually calibrated to historical daily mean thermal spring flows of 88,000 cubic feet per day observed over a 12-year period of record (1990–1995 and 1998–2005) at the thermal springs collection system. Calibration was achieved by calculating starting heads and general head boundary conditions from the Bernoulli equation and then adjusting the horizontal and vertical hydraulic conductivities of the rock formations and vertical fault and the hydraulic conductance of head-dependent flux boundaries. The groundwater-flow model was coupled to a surface-water model developed in the Precipitation-Runoff Modeling System (PRMS) by using PRMS-simulated gravity drainage as a specified flux recharge boundary condition in the groundwater-flow model. Together, the coupled models were used to (1) locate the areas of groundwater recharge to the thermal springs in the discretized hydrogeologic framework by using forward and reverse particle-tracking capabilities of MODPATH, (2) simulate the effects of variable recharge rates on the spring flows at the thermal springs, and (3) assess possible effects of climate and land-use change on the long-term variability of spring flows at the thermal springs. Forward and backward particle-tracking maps indicated that the most prevalent areas of recharge in the discretized hydrogeologic framework used in this study were within about 0.6–0.9 mile of the thermal springs. Forward particle tracking indicated a recharge area southwest of the thermal springs that corresponded to a location where the predominant lithologies are the Arkansas Novaculite, Hot Springs Sandstone, and Bigfork Chert. Backward particle tracking indicated a second localized area of recharge to the northeast of the thermal springs that corresponded to a location where the dominant lithology is the Bigfork Chert. The groundwater-flow model indicated that the most probable recharge formations are the Arkansas Novaculite, Bigfork Chert, and Hot Springs Sandstone. The simulated effects of climate and land-use changes on the variability of the spring-flow rates at the thermal springs generally resulted in reductions of thermal spring flow attributed to urban development and more extreme climates characterized by elevated mean surface air temperatures. The groundwater-flow model predicted a linear relation between the thermal spring discharge and the cumulative recharge volume applied to the hydrogeologic framework, and the positive slope of the predicted relation between recharge and simulated thermal spring flow indicates that more extreme precipitation events that supply more recharge may in fact increase the thermal spring-flow rates.

Arkansas↗

USGS2020

Comparison of groundwater storage changes from GRACE satellites with monitoring and modeling of major U.S. aquifers

GRACE satellite data are widely used to estimate groundwater storage (GWS) changes in aquifers globally; however, comparisons with GW monitoring and modeling data are limited. Here we compared GWS changes from GRACE over 15 yr (2002–2017) in 14 major U.S. aquifers with groundwater-level (GWL) monitoring data in ~23,000 wells and with regional and global hydrologic and land surface models. Results show declining GWS trends from GRACE data in the six southwestern and south-central U.S. aquifers, totaling −90 km 3 over 15 yr, related to long-term (5–15 yr) droughts, and exceeding Lake Mead volume by ~2.5×. GWS trends in most remaining aquifers were stable or slightly rising. GRACE-derived GWS changes agree with GWL monitoring data in most aquifers (correlation coefficients, R = 0.52–0.95), showing that GRACE satellites capture groundwater (GW) dynamics. Regional GW models (eight models) generally show similar or greater GWS trends than those from GRACE. Large discrepancies in the Mississippi Embayment aquifer, with modeled GWS decline approximately four times that of GRACE, may reflect uncertainties in model storage parameters, stream capture, pumpage, and/or recharge rates. Global hydrologic models (2003–2014), which include GW pumping, generally overestimate GRACE GWS depletion (total: approximately −172 to −186 km 3 ) in heavily exploited aquifers in southwestern and south-central U.S. by ~2.4× (GRACE: −74 km 3 ), underscoring needed modeling improvements relative to anthropogenic impacts. Global land surface models tend to track GRACE GWS dynamics better than global hydrologic models. Intercomparing remote sensing, monitoring, and modeling data underscores the importance of considering all data sources to constrain GWS uncertainties.

Water Resources Research↗

USGS2020

Water priorities for the Nation—U.S. Geological Survey Integrated Water Availability Assessments

The United States faces growing challenges to its water supply, infrastructure, and aquatic ecosystems because of population growth, climate change, floods and droughts, and aging water delivery systems. To help address these challenges, the U.S. Geological Survey (USGS) Water Resources Mission Area has established new strategic priorities that capitalize on the operational and scientific strengths of the USGS to address these complex societal issues. The USGS Integrated Water Availability Assessment Program within the Water Resources Mission Area will provide nationally consistent assessments of water available for human and ecological needs and identify factors that limit water availability. Integrated Water Availability Assessments will assess current conditions, long-term trends, factors responsible for current conditions and observed trends, and provide forecasts of future changes at national and regional scales.

Conterminous United States↗

USGS2019

Groundwater availability in the Ozark Plateaus aquifer system

Executive Summary The study described in this report, initiated by the U.S. Geological Survey in 2014, was designed to evaluate fresh groundwater resources within the Ozark Plateaus, central United States, as an area within a broader national assessment of groundwater availability. The goals of the Ozark study were to evaluate historical effects of human activities on water levels and groundwater availability, quantify groundwater resources now and under probable future pumping and climate conditions, and evaluate existing monitoring networks for their value in making better predictions of future groundwater resources. Previous studies include simulation of local-scale groundwater flow under varying temporal scales, or simulation of the regional system under steady-state conditions. While these studies are useful, particularly for the problem for which they were designed, there is a need to look at the larger regional system under transient conditions to fully evaluate the water resource over time. This study focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water-level declines. The regional scale provides a broad view of the sources and demands on the system with time. The study area covers approximately 68,000 square miles in the central United States in parts of Missouri, Arkansas, Kansas, and Oklahoma and encompasses the Ozark Plateaus Physiographic Province (Ozark Plateaus), including the Salem Plateau, Springfield Plateau, and Boston Mountains. Groundwater is withdrawn from the Ozark Plateaus aquifer system (Ozark system) for public supply and for domestic, agriculture (including irrigation and aquaculture), livestock, and non-agricultural use (including industrial, thermoelectric power generation, mining, and commercial). The Ozark system provides an important drinking-water supply for people living in the Ozark Plateaus because public supply and domestic use combined constitute the largest groundwater use. Precipitation is the ultimate source of freshwater to the Ozark system; most rainfall occurs during April, May, and June, and precipitation increases generally from north to south across the study area. Groundwater use currently accounts for only 10 percent of the total water use in the areas overlying the Ozark system, but provides a critical drinking-water resource because public supply and domestic groundwater withdrawals are largely from groundwater resources. The 380 million gallons per day of groundwater withdrawn from the Ozark system in 2010 accounts for approximately 2 percent of recharge. Although groundwater use represents a small component of the hydrologic budget, because of low storage in aquifer units, cones of depression with steep water-level gradients can develop quickly around pumping centers. The amount of water entering and leaving the aquifer system from 1900 to about 1965 was relatively constant at a rate of about 13 billion gallons per day (Bgal/d). Much of this inflow of water is discharged through streams in the system to balance the hydrologic budget. Changes in storage over time (from outflows to inflows) reflect the large variability in recharge: if recharge decreases, water levels will decrease, resulting in less groundwater discharge to streams and more water released from aquifer storage. Conversely, when recharge increases, water levels increase, more groundwater discharges to streams, and aquifer storage is replenished. Although pumping generally increased from 1900 to 2016, it does not appear to correlate with the change in storage over the same time period. Regionally, simulated change in groundwater storage corresponds with changes in recharge, more so than with increases in pumping. Average recharge was 11.6 Bgal/d for the period 1900 to 2016. Recharge was generally above average from predevelopment to 1965, followed by a period of below-average recharge from 1965 to about 1980. Recharge remained consistently above average from 1980 to about 1988, after which there was a period of average or below-average recharge, reflected by a decline through the mid-2000s. The implications and potential effects of increased pumping and long-term climate change on the Ozark Plateaus hydrologic system and groundwater availability are a concern for communities and resource managers in the area. Pumping varies from year to year, but is generally expected to moderately increase with population, industrial, and agricultural needs. Most climate models predict warmer minimum and maximum air temperatures by midcentury in the Ozark Plateaus area, especially from midspring through early fall. Three scenarios were developed to simulate possible future conditions from 2016 to 2060 and assess the potential effects on the hydrologic system and availability of water resources. For each scenario, changes in water levels and hydrologic budget components were evaluated from predevelopment (1900) to present (2016) and 45 years into the future (2060). The baseline scenario represents an extension of the average (1996 to 2016) seasonal pumping and recharge values. The pumping scenario is an extension of the average (1996 to 2016) seasonal recharge values with increases in pumping following the historical trend for the period 2016–2060 of up to 120 percent of the 1996 to 2016 average seasonal pumping values. The general circulation model (GCM) scenario is an extension of the average (1996 to 2016) seasonal pumping values and variable recharge based on seasonal averages of soil water storage from a water-balance model using temperature and precipitation from multiple GCMs. The general patterns of water-level decline are similar for each scenario. The areas of water-level decline in southwest Missouri and northeast Oklahoma are only marginally different by 2060 from those of 2009. In one area south of Springfield, Mo., water-level declines are less in the baseline and GCM scenarios than in 2009. This may be the result of a transition from groundwater use to surface-water supplies for a larger percentage of the demand in the area. For all three scenarios, forecasted pumping, recharge, and aquifer properties play an important role in determining the uncertainty of water-level forecasts at 94 real-time observation wells. Simulated aquifer properties in the productive middle and lower Ozark aquifers and the St. Francois confining unit of the Ozark system contribute most to predictive uncertainty in water levels at approximately 35 percent of the real-time observation wells. Out of the 94 real-time observation wells, 82 are developed in the lower Ozark aquifer.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2018

Extraction and development of inset models in support of groundwater age calculations for glacial aquifers

The U.S. Geological Survey developed a regional model of Lake Michigan Basin (LMB). This report describes the construction of five MODFLOW inset models extracted from the LMB regional model and their application using the particle-tracking code MODPATH to simulate the groundwater age distribution of discharge to wells pumping from glacial deposits. The five study areas of the inset model correspond to 8-digit hydrologic unit code (HUC8) basins. Two of the basins are tributary to Lake Michigan from the east, two are tributary to the lake from the west, and one is just west of the western boundary of the Lake Michigan topographic basin. The inset models inherited many of the inputs to the parent LMB model, including the hydrostratigraphy and layering scheme, the hydraulic conductivity assigned to bedrock layers, recharge distribution, and water use in the form of pumping rates from glacial and bedrock wells. The construction of the inset models entailed modifying some inputs, most notably the grid spacing (reduced from cells 5,000 feet on a side in the parent LMB model to 500 feet on a side in the inset models). The refined grid spacing allowed for more precise location of pumped wells and more detailed simulation of groundwater/surface-water interactions. The glacial hydraulic conductivity values, the top bedrock surface elevation, and the surface-water network input to the inset models also were modified. The inset models are solved using the MODFLOW–NWT code, which allows for more robust handling of conditions in unconfined aquifers than previous versions of MODFLOW. Comparison of the MODFLOW inset models reveals that they incorporate a range of hydrogeologic conditions relative to the glacial part of the flow system, demonstrated by visualization and analysis of model inputs and outputs and reflected in the range of ages generated by MODPATH for existing and hypothetical glacial wells. Certain inputs and outputs are judged to be candidate predictors that, if treated statistically, may be capable of explaining much of the variance in the simulated age metrics. One example of a predictor that model results indicate strongly affects simulated age is the depth of the well open interval below the simulated water table. The strength of this example variable as an overall predictor of groundwater age and its relation to other predictors can be statistically tested through the metamodeling process. In this way the inset models are designed to serve as a training area for metamodels that estimate groundwater age in glacial wells, which in turn will contribute to ongoing studies, under the direction of the U.S. Geological Survey National Water Quality Assessment, of contaminant susceptibility of shallow groundwater across the glacial aquifer system.

Illinois, Indiana, Michigan, Wisconsin↗

USGS2018

The Ozark Plateaus Regional Aquifer Study—Documentation of a groundwater-flow model constructed to assess water availability in the Ozark Plateaus

Recent short-term drought conditions have emphasized the need to better understand the delicate balance between abundance, sustainability, and scarcity of groundwater in the Ozark Plateaus aquifer system. In 2014, the U.S. Geological Survey began construction of a groundwater-flow model as a tool for the assessment of groundwater availability in the Ozark Plateaus aquifer system. The model was developed to benefit concurrent and future investigations involving groundwater-pumping scenarios, optimization, particle transport, and groundwater-monitoring network analysis. The groundwater model simulates 116 years (1900–2015) of hydrologic conditions and the response of the groundwater system to changes in stress including changes in recharge and groundwater pumping for water supply. Semiseasonal stress periods were simulated from the later part of 1991 to 2015 and represent higher demand and lower recharge in the spring and summer months and lower demand and higher recharge in the fall and winter months. Groundwater pumping increases throughout the simulation period with a maximum rate of about 600 million gallons per day (Mgal/d). The process of matching historical hydrologic data for the Ozark Plateaus aquifer system model was accomplished by a combination of manual changes to parameter values and automated calibration methods. Observation data used in the development and evaluation of the model included 19,045 hydraulic-head observations from 6,683 wells within the model area. Observation data also included stream leakage estimates summed to calculate a net gain or net loss value for approximately 81 named streams. The majority (mean of over 95 percent) of the recharge component is discharged through streams simulated in the model. The total simulated discharge to streams fluctuates seasonally between 7,500 and 17,500 Mgal/d with a mean outflow of 11,500 Mgal/d. Much of the remaining balance between modeled recharge inflows and stream outflows is made up by water moving into or out of storage in the aquifer system resulting in changes in modeled groundwater levels. The goal of the model was to develop a model capable of suitable accuracy at regional scales. The intent was not to reproduce individual local-scale details, which are typically not possible given the uniform cell size of 1 square mile. Although the model may not represent each local-scale detail, the model can be applied for a better understanding of the regional flow system and to evaluate responses to changes in climate and groundwater pumping.

Ozark Plateaus aquifer system↗

USGS2017

Groundwater development stress: Global-scale indices compared to regional modeling

The increased availability of global datasets and technologies such as global hydrologic models and the Gravity Recovery and Climate Experiment (GRACE) satellites have resulted in a growing number of global-scale assessments of water availability using simple indices of water stress. Developed initially for surface water, such indices are increasingly used to evaluate global groundwater resources. We compare indices of groundwater development stress for three major agricultural areas of the United States to information available from regional water budgets developed from detailed groundwater modeling. These comparisons illustrate the potential value of regional-scale analyses to supplement global hydrological models and GRACE analyses of groundwater depletion. Regional-scale analyses allow assessments of water stress that better account for scale effects, the dynamics of groundwater flow systems, the complexities of irrigated agricultural systems, and the laws, regulations, engineering, and socioeconomic factors that govern groundwater use. Strategic use of regional-scale models with global-scale analyses would greatly enhance knowledge of the global groundwater depletion problem.

Groundwater↗

USGS2017

Challenges for creating a site-specific groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010

Hydrologic budgets to determine groundwater availability are important tools for water-resource managers. One challenging component for developing hydrologic budgets is quantifying water use through time because historical and site-specific water-use data can be sparse or poorly documented. This research developed a groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010 that related county-level aggregated water-use data to site-specific well locations and aquifer units. A simple population-based linear model, constrained to 0 million liters per day in 1900, provided the best means to extrapolate groundwater-withdrawal rates pre-1950s when there was a paucity of water-use data. To disaggregate county-level data to individual wells across a regional aquifer system, a programmatic hierarchical process was developed, based on the level of confidence that a well pumped groundwater for a specific use during a specific year. Statistical models tested on a subset of the best-available site-specific water-use data provided a mechanism to bracket historic groundwater use, such that groundwater-withdrawal rates ranged, on average, plus or minus 38% from modeled values. Groundwater withdrawn for public supply and domestic use accounted for between 48 and 74% of total groundwater use since 1901, highlighting that groundwater provides an important drinking-water resource. The compilation, analysis, and spatial and temporal extrapolation of water-use data remain a challenging task for water scientists, but is of paramount importance to better quantify groundwater use and availability.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2017

Simulation of groundwater flow in the glacial aquifer system of northeastern Wisconsin with variable model complexity

The U.S. Geological Survey, National Water-Quality Assessment seeks to map estimated intrinsic susceptibility of the glacial aquifer system of the conterminous United States. Improved understanding of the hydrogeologic characteristics that explain spatial patterns of intrinsic susceptibility, commonly inferred from estimates of groundwater age distributions, is sought so that methods used for the estimation process are properly equipped. An important step beyond identifying relevant hydrogeologic datasets, such as glacial geology maps, is to evaluate how incorporation of these resources into process-based models using differing levels of detail could affect resulting simulations of groundwater age distributions and, thus, estimates of intrinsic susceptibility. This report describes the construction and calibration of three groundwater-flow models of northeastern Wisconsin that were developed with differing levels of complexity to provide a framework for subsequent evaluations of the effects of process-based model complexity on estimations of groundwater age distributions for withdrawal wells and streams. Preliminary assessments, which focused on the effects of model complexity on simulated water levels and base flows in the glacial aquifer system, illustrate that simulation of vertical gradients using multiple model layers improves simulated heads more in low-permeability units than in high-permeability units. Moreover, simulation of heterogeneous hydraulic conductivity fields in coarse-grained and some fine-grained glacial materials produced a larger improvement in simulated water levels in the glacial aquifer system compared with simulation of uniform hydraulic conductivity within zones. The relation between base flows and model complexity was less clear; however, the relation generally seemed to follow a similar pattern as water levels. Although increased model complexity resulted in improved calibrations, future application of the models using simulated particle tracking is anticipated to evaluate if these model design considerations are similarly important for understanding the primary modeling objective - to simulate reasonable groundwater age distributions.

Wisconsin↗

USGS2016

Hydrogeology and hydrologic conditions of the Ozark Plateaus aquifer system

The hydrogeology and hydrologic characteristics of the Ozark Plateaus aquifer system were characterized as part of ongoing U.S. Geological Survey efforts to assess groundwater availability across the Nation. The need for such a study in the Ozark Plateaus physiographic province (Ozark Plateaus) is highlighted by increasing demand on groundwater resources by the 5.3 million people of the Ozark Plateaus, water-level declines in some areas, and potential impacts of climate change on groundwater availability. The subject study integrates knowledge gained through local investigation within a regional perspective to develop a regional conceptual model of groundwater flow in the Ozark Plateaus aquifer system (Ozark system), a key phase of groundwater availability assessment. The Ozark system extends across much of southern Missouri and northwestern and north-central Arkansas and smaller areas of southeastern Kansas and northeastern Oklahoma. The region is one of the major karst landscapes in the United States, and karst aquifers are predominant in the Ozark system. Groundwater flow is ultimately controlled by aquifer and confining unit lithologies and stratigraphic relations, geologic structure, karst development, and the character of surficial lithologies and regolith mantle. The regolith mantle is a defining element of Ozark Plateaus karst, affecting recharge, karst development, and vulnerability to surface-derived contaminants. Karst development is more advanced—as evidenced by larger springs, hydraulic characteristics, and higher well yields—in the Salem Plateau and in the northern part of the Springfield Plateau (generally north of the Arkansas-Missouri border) as compared with the southern part of the Springfield Plateau in Arkansas, largely due to thinner, less extensive regolith and purer carbonate lithology. Precipitation is the ultimate source of all water to the Ozark system, and the hydrologic budget for the Ozark system includes inputs from recharge, losing-stream sections, and groundwater inflows and losses of water to gaining-stream sections, groundwater withdrawals, and surface-water and groundwater outflows to neighboring systems. Groundwater recharge, estimated by a soil-water-balance model, represents about 24 percent, or 11 inches, of 43.9 inches annual precipitation. Recharge is spatially variable, being greater in the northern Springfield Plateau and Salem Plateau than in the southern Springfield Plateau (generally south of the Arkansas border) because of differences in regolith mantle extent and thickness and carbonate lithology and hydraulic properties. Increased precipitation and decreased agricultural land use during the period 1951 through 2011 increased recharge by approximately 5 percent. Although all Ozark streams have losing, neutral, and gaining sections, they are dominantly gaining and are a net sink for groundwater with nearly 90 percent of groundwater recharge returned to springs and streams. Groundwater pumping is a small but important loss of water in the Ozark system hydrologic budget; water-level declines and local cones of depression have been observed around pumping centers and strong concerns exist over potential effects on stream and spring flow. Data indicate that societal needs for freshwater resources in the Ozark Plateaus will continue to increase and will do so in the context of changing climate and hydrology. Groundwater will continue to be an important part of supporting these societal needs and also local ecosystems. The unique character and hydrogeologic variability across the Ozark system will control how the system responds to future stress. Groundwater of the Ozark system in the northern study area is more dynamic, has greater storage and larger flux, and has greater potential for further development than in the part of the study area south of the Arkansas-Missouri border. Further south in Arkansas, a line exists, roughly defined as 5 miles south of the Springfield Plateau-Boston Mountains boundary, beyond which further extensive municipal or commercial development appears unlikely under current economic and resource-need conditions. A small part of the Ozark system groundwater budget is currently drafted for use, leaving an apparently large component available for further development and use—particularly in the northern Springfield Plateau and Salem Plateau; however, the effects of increased pumping on groundwater’s role in maintaining ecosystems and ecosystem services are not quantitatively well understood, and the close relation between groundwater and surface water highlights the importance of further quantitative assessment.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2016

Altitudes and thicknesses of hydrogeologic units of the Ozark Plateaus aquifer system in Arkansas, Kansas, Missouri, and Oklahoma

A hydrogeologic framework was constructed to represent the altitudes and thicknesses of hydrogeologic units within the Ozark Plateaus aquifer system as part of a regional groundwater-flow model supported by the U.S. Geological Survey Water Availability and Use Science Program. The Ozark Plateaus aquifer system study area is nearly 70,000 square miles and includes parts of Arkansas, Kansas, Missouri, and Oklahoma. Nine hydrogeologic units were selected for delineation within the aquifer system and include the Western Interior Plains confining system, the Springfield Plateau aquifer, the Ozark confining unit, the Ozark aquifer, which was divided into the upper, middle, and lower Ozark aquifers to better capture the spatial variation in the hydrologic properties, the St. Francois confining unit, the St. Francois aquifer, and the basement confining unit. Geophysical and well-cutting logs, along with lithologic descriptions by well drillers, were compiled and interpreted to create hydrologic altitudes for each unit. The final compiled dataset included more than 23,000 individual altitude points (excluding synthetic points) representing the nine hydrogeologic units within the Ozark Plateaus aquifer system.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2015

Chemical variations in Yellowknife Bay formation sedimentary rocks analyzed by ChemCam on board the Curiosity rover on Mars

The Yellowknife Bay formation represents a ~5 m thick stratigraphic section of lithified fluvial and lacustrine sediments analyzed by the Curiosity rover in Gale crater, Mars. Previous works have mainly focused on the mudstones that were drilled by the rover at two locations. The present study focuses on the sedimentary rocks stratigraphically above the mudstones by studying their chemical variations in parallel with rock textures. Results show that differences in composition correlate with textures and both manifest subtle but significant variations through the stratigraphic column. Though the chemistry of the sediments does not vary much in the lower part of the stratigraphy, the variations in alkali elements indicate variations in the source material and/or physical sorting, as shown by the identification of alkali feldspars. The sandstones contain similar relative proportions of hydrogen to the mudstones below, suggesting the presence of hydrous minerals that may have contributed to their cementation. Slight variations in magnesium correlate with changes in textures suggesting that diagenesis through cementation and dissolution modified the initial rock composition and texture simultaneously. The upper part of the stratigraphy (~1 m thick) displays rocks with different compositions suggesting a strong change in the depositional system. The presence of float rocks with similar compositions found along the rover traverse suggests that some of these outcrops extend further away in the nearby hummocky plains.

Journal of Geophysical Research E: Planets↗

USGS2015

Simulations of potential future conditions in the cache critical groundwater area, Arkansas

A three-dimensional finite-difference model for part of the Mississippi River Valley alluvial aquifer in the Cache Critical Groundwater Area of eastern Arkansas was constructed to simulate potential future conditions of groundwater flow. The objectives of this study were to test different pilot point distributions to find reasonable estimates of aquifer properties for the alluvial aquifer, to simulate flux from rivers, and to demonstrate how changes in pumping rates for different scenarios affect areas of long-term water-level declines over time. The model was calibrated using the parameter estimation code. Additional calibration was achieved using pilot points with regularization and singular value decomposition. Pilot point parameter values were estimated at a number of discrete locations in the study area to obtain reasonable estimates of aquifer properties. Nine pumping scenarios for the years 2011 to 2020 were tested and compared to the simulated water-level heads from 2010. Hydraulic conductivity values from pilot point calibration ranged between 42 and 173 m/d. Specific yield values ranged between 0.19 and 0.337. Recharge rates ranged between 0.00009 and 0.0006 m/d. The model was calibrated using 2,322 hydraulic head measurements for the years 2000 to 2010 from 150 observation wells located in the study area. For all scenarios, the volume of water depleted ranged between 5.7 and 23.3 percent, except in Scenario 2 (minimum pumping rates), in which the volume increased by 2.5 percent.

Arkansas↗

USGS2014

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2—two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗

USGS2013

Simulation of the June 11, 2010, flood along the Little Missouri River near Langley, Arkansas, using a hydrologic model coupled to a hydraulic model

A substantial flood event occurred on June 11, 2010, causing the Little Missouri River to flow over much of the adjacent land area, resulting in catastrophic damages. Twenty fatalities occurred and numerous automobiles, cabins, and recreational vehicles were destroyed within the U.S. Department of Agriculture-Forest Service Albert Pike Recreation Area, at a dispersed campsite area in the surrounding Ouachita National Forest lands, and at a nearby privately owned camp. The Little Missouri River streamgage near Langley, Arkansas, reached a record streamflow of 70,800 cubic feet per second and a stage (water level) of 23.5 feet at 5:30 a.m., with a 10-foot rise occurring in slightly more than 1 hour. To better understand the flood event on June 11, 2010, the U.S. Geological Survey, in cooperation with the U.S. Department of Agriculture-Forest Service, developed a precipitation-runoff hydrologic model, U.S. Army Corps of Engineers Hydrologic Engineering Center Hydrologic Modeling System (HEC-HMS), coupled with a one-dimensional unsteady-state hydraulic model, U.S. Army Corps of Engineers Hydrologic Engineering Center River Analysis System (HEC-RAS), to simulate precipitation runoff and streamflow characteristics along the Little Missouri River and at various tributaries within the 68-square mile watershed upstream from the Langley streamgage. Within the proximity of two campgrounds, the Little Missouri River just downstream from the confluence of Brier Creek had a peak simulated streamflow of 49,300 cubic feet per second at 4:08 a.m.; the simulated streamflow stayed within 500 cubic feet per second of the peak for nearly 15 minutes. The simulated water surface increased an average of 0.5 feet every 5 minutes for a total of 2 hours, with a maximum rate of rise of 2 feet in 15 minutes. The Little Missouri River just downstream from the confluence of Brier Creek had a peak simulated water-surface elevation of 935.0 feet, a maximum water depth of 22.2 feet, and a maximum stream channel velocity of 12.6 feet per second at 4:15 a.m. The results from the precipitation-runoff hydrologic model, the one-dimensional unsteady-state hydraulic model, and a separate two-dimensional model developed as part of a coincident study, each complement the other in terms of streamflow timing, water-surface elevations, and velocities propagated by the June 11, 2010, flood event. The simulated grids for water depth and stream velocity from each model were directly compared by subtracting the one-dimensional hydraulic model grid from the two-dimensional model grid. The absolute mean difference for the simulated water depth was 0.9 foot. Additionally, the absolute mean difference for the simulated stream velocity was 1.9 feet per second.

Arkansas↗

USGS2013

Enhancements to the Mississippi Embayment Regional Aquifer Study (MERAS) groundwater-flow model and simulations of sustainable water-level scenarios

Arkansas continues to be one of the largest users of groundwater in the Nation. As such, long-term planning and management are essential to ensure continued availability of groundwater and surface water for years to come. The Mississippi Embayment Regional Aquifer Study (MERAS) model was developed previously as a tool to evaluate groundwater availability within the Mississippi embayment, which encompasses much of eastern Arkansas where the majority of groundwater is used. The Arkansas Water Plan is being updated for the first time since 1990 and serves as the State’s primary, comprehensive water-resources planning and guidance document. The MERAS model was selected as the best available tool for evaluation of specific water-use pumping scenarios that are currently being considered by the State of Arkansas. The model, developed as part of the U.S. Geological Survey Groundwater Resources Program’s assessment of the Nation’s groundwater availability, is proving to be invaluable to the State as it works toward development of a sustained yield pumping strategy. One aspect of this investigation was to evaluate multiple methods to improve the match of observed to simulated groundwater levels within the Mississippi River Valley alluvial and middle Claiborne (Sparta) aquifers in the MERAS model. Five primary methods were evaluated: (1) explicit simulation of evapotranspiration (ET), (2) upgrade of the Multi-Node Well (MNW2) Package, (3) geometry improvement within the Streamflow Routing (SFR) Package, (4) parameter estimation of select aquifer properties with pilot points, and (5) modification of water-use estimates. For the planning purposes of the Arkansas Water Plan, three scenarios were developed to evaluate potential future conditions: (1) simulation of previously optimized pumping values within the Mississippi River Valley alluvial and the middle Claiborne aquifers, (2) simulated prolonged effects of pumping at average recent (2000–5) rates, and (3) simulation of drawdown constraints on most pumping wells. The explicit simulation of ET indicated little, if any, improvement of model fit at the expense of much longer simulation time and was not included in further simulations. Numerous attempts to fully utilize the MNW2 Package were unsuccessful in achieving model stability, though modifications made to the water-use dataset remained intact. Final improvements in the residual statistics may be attributed to a single method, or a cumulative effect of all other methods (geometry improvement with the SFR Package, parameter estimation with pilot points, and modification of water-use estimates) attempted. The root mean squared error (RMSE) for all observations in the model is 22.65 feet (ft) over a range in observed hydraulic head of 741.66 ft. The RMSE for water-level observations in the Mississippi River Valley alluvial aquifer is 14.14 ft (an improvement of almost 3 ft) over a range in observed hydraulic head of 297.25 ft. The RMSE for the Sparta aquifer is 32.02 ft (an improvement of approximately 3 ft) over a range in observed hydraulic head of 634.94 ft. Three scenarios were developed to utilize a steady-state version of the MERAS model. Scenario 1 was developed to use pumping values resulting from the optimization of baseline rates (typically 1997 pumping rates) from previous optimization modeling of the alluvial aquifer and the Sparta aquifer. Scenario 2 was developed to evaluate the prolonged effects of pumping from the alluvial aquifer at recent pumping rates. Scenario 3A was designed to evaluate withdrawal limits from the alluvial aquifer by utilizing drawdown constraints equal to an altitude of approximately 50 percent of the predevelopment saturated thickness of the alluvial aquifer or 30 ft above the bottom of the alluvial aquifer, whichever was greater. The results of scenario 1 indicate large water-level declines throughout the area of the alluvial aquifer, regardless of the substitution of the optimized pumping values from earlier model simulations. The results of scenario 2 also indicate large areas of water-level decline, as compared to half of the saturated thickness, throughout the alluvial aquifer. The results of scenario 3A reveal some effects from the inclusion of multiple aquifers in a single simulation. The initial configuration of scenario 3A resulted in water levels well below the defined drawdown constraint, and some areas of depleted aquifer (water levels that are near or below the bottom of the aquifer) in east-central Arkansas. A fourth simulation (scenario 3B) was configured to apply the same drawdown constraints from the alluvial aquifer wells to the Sparta aquifer wells in the depleted area. These drawdown constraints reduce leakage from the alluvial aquifer to the underlying Sparta aquifer. This configuration did not produce depleted areas within the alluvial aquifer. Scenarios 3A and 3B indicate that even when pumping is limited in the alluvial aquifer, water levels in the alluvial aquifer may continue to decline in some areas because of pumping in the underlying Sparta aquifer.

Arkansas↗

USGS2011

A new tool to assess groundwater resources in the Mississippi embayment

What is the Mississippi Embayment? The Mississippi embayment study area encompasses approximately 78,000 square miles in eight States and includes large parts of Arkansas, Louisiana, Mississippi, and Tennessee, and smaller areas of Alabama, Illinois, Kentucky, and Missouri (fig. 1). The Mississippi embayment is essentially a basin that slopes toward the Gulf of Mexico and is filled with sediments of alternating sand, silt, and clay layers. There are two principal aquifers in the embayment-the Mississippi River Valley alluvial aquifer (alluvial aquifer) and the middle Claiborne aquifer (fig. 1). The shallow alluvial aquifer is the primary source of groundwater for irrigation in the largely agricultural region, while the deeper middle Claiborne aquifer is a primary source of drinking water for many of the 5.2 million people living in the embayment. The U.S. Geological Survey (USGS) is conducting large-scale multidisciplinary regional studies of groundwater availability for the Nation. Studies comprise individual assessments of regional groundwater-flow systems that encompass varied terrains and document a comprehensive regional and national perspective of groundwater resources. Collectively, these studies are the foundation for the national assessment of groundwater availability and are conducted in cooperation with other Federal, State, local governments, and the private sector. Numerical groundwater-flow models are used in these studies to document effects of human activities and climate variability on groundwater levels, changes in aquifer storage, and flow between groundwater and surface-water bodies. As part of the Mississippi Embayment Regional Aquifer Study (MERAS), a numerical model was constructed of 13 layers over 78,000 square miles representing multiple aquifers and confining units for the period of 1870 to 2007. The model is a tool that was used to assess and better understand groundwater resources.

Mississippi↗

USGS2011

Groundwater availability of the Mississippi embayment

Groundwater is an important resource for agricultural and municipal uses in the Mississippi embayment. Arkansas ranks first in the Nation for rice and third for cotton production, with both crops dependent on groundwater as a major source of irrigation requirements. Multiple municipalities rely on the groundwater resources to provide water for industrial and public use, which includes the city of Memphis, Tennessee. The demand for the groundwater resource has resulted in groundwater availability issues in the Mississippi embayment including: (1) declining groundwater levels of 50 feet or more in the Mississippi River Valley alluvial aquifer in parts of eastern Arkansas from agricultural pumping, (2) declining groundwater levels of over 360 feet over the last 90 years in the confined middle Claiborne aquifer in southern Arkansas and northern Louisiana from municipal pumping, and (3) litigation between the State of Mississippi and a Memphis water utility over water rights in the middle Claiborne aquifer. To provide information to stakeholders addressing the groundwater-availability issues, the U.S. Geological Survey Groundwater Resources Program supported a detailed assessment of groundwater availability through the Mississippi Embayment Regional Aquifer Study (MERAS). This assessment included (1) an evaluation of how these resources have changed over time through the use of groundwater budgets, (2) development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends, and (3) application of statistical tools to evaluate the importance of individual observations within a groundwater-monitoring network. An estimated 12 million acre-feet per year (11 billion gallons per day) of groundwater was pumped in 2005 from aquifers in the Mississippi embayment. Irrigation constitutes the largest groundwater use, accounting for approximately 10 million acre-feet per year (9 billion gallons per day) in 2000 from the Mississippi River Valley alluvial aquifer in Arkansas, Louisiana, Mississippi, and Missouri, and to a lesser extent in Illinois, Kentucky, and Tennessee. Predevelopment groundwater flow is represented in the MERAS model as a steady-state stress period, assumed to be prior to 1870. The simulated groundwater-flow budget indicates the largest predevelopment inflow to the system is net recharge to the alluvial aquifer. This inflow is balanced by outflow to gaining streams. Overall, water enters as net recharge to the alluvial aquifer or through outcrop areas of the various hydrogeologic units. Away from the outcrop areas, groundwater flow in the deeper formations is primarily upward into overlying units, ultimately discharging to streams through the alluvial aquifer. Total net recharge and discharge (sum of inflows or outflows) for the model ranged from about 0.66 million acre-feet per year during predevelopment to 20.16 million acre-feet per year by the end of the simulation (final simulated irrigation period in summer of 2006). This change in the model budget reflects increases in withdrawals compared to predevelopment conditions. Cumulative storage within aquifers simulated in the MERAS model indicates overall depletion of 140 million acre-feet (equivalent to 2.8 feet of water covering the entire study area). Postdevelopment inflow to the system is still through net recharge to the alluvial aquifer and the outcrop areas of the several hydrogeologic units, however, the flow between each unit is no longer upward to the alluvial aquifer. Groundwater flow during the summer of 2006 was primarily downward to offset demand from pumping. Early in the model simulation (1870-1920s), the primary components of the water budget were simulated as outflow from stream leakage and inflow from net recharge. As pumpage increased through time, water that would otherwise flow to streams reversed, and net stream leakage became an inflow to the system. The largest reversals began in the mid-1980s, but indications of the reversal began in the early 1960s with a trend in loss of streamflow leakage coupled with the first consistent inflow from storage. While groundwater pumped out of the alluvial aquifer was derived primarily from storage, pumpage out of the middle Claiborne aquifer was derived primarily from other aquifers (up to 15 percent from the alluvial aquifer), followed by flow from storage and net recharge. The potential consequences of climate change have been identified as a major concern facing the sustainability of the Nation's groundwater resources. To address this concern, two climate simulations were developed through the use of the MERAS model by extending the simulation period by 30 years to the year 2038 using extrapolated precipitation based on frequency analysis of historic climate cycles. There is little difference between the dry and wet scenarios in terms of percent water-level change. Both scenarios resulted in 14.6 to 13.9 percent of the area containing more than 100 feet of decline, 14.5 to 13.8 percent containing between 75 and 100 feet of decline, and 15.8 to 15.7 percent containing 51 to 75 feet of decline in the alluvial aquifer. The middle Claiborne aquifer water-level changes also were similar between the two scenarios. These scenarios indicate that even with a 25-percent increase in precipitation from that of the dry scenario, there is little difference in the resultant water levels. This is in large part because of the magnitude of differences between changes in net recharge and changes in pumping. When compared to the volume of water pumped out of the system, the effect of this change in net recharge is negligible. The groundwater-level monitoring network used to construct the 2007 middle Claiborne aquifer potentiometric surface was used as an example case to demonstrate statistical technique and to evaluate the importance of individual groundwater-level observations. To calculate the importance of each water-level observation to a prediction, predictions were specified as water-level altitudes near the end of the dry scenario simulation. These predictions were located near the center of cones of depression. Many of the observations that have a high importance are in close proximity to stressed areas of the aquifer.

Mississippi↗