Search USGS⌕ Search

Geology topics

Eve L. Kuniansky

Publications and source records attributed to Eve L. Kuniansky.

At least 19 recordsLinked to original sources
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2021

National-scale reservoir thermal energy storage pre-assessment for the United States

The U.S. Geological Survey is performing a pre-assessment of the cooling potential for reservoir thermal energy storage (RTES) in five generalized geologic regions (Basin and Range, Coastal Plains, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Reservoir models are developed for the metropolitan areas of eight cities (Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon) so that computed metrics can be compared to evaluate RTES potential across diverse climates, geologic settings, and physiography. Permeable, semi-confined/confined units that underlie more-utilized aquifers and contain low-quality groundwater are selected for each city. Energy storage metrics are computed for the anticipated total thickness of stratigraphy for which RTES might be feasible, including estimated required well spacing, thermal storage capacity, and thermal recovery efficiency over time. Falta et al. (2016) showed that for a modern 25,000 square-foot (2,323 square-meter), two-story office building, cooling needs exceed heating demand for almost every region of the country. We therefore use Falta et al.’s cooling demand for each city as the representative RTES stress condition for metric computation, allowing comparisons across regions. Results indicate that favorable RTES conditions exist in each region, particularly in the Illinois Basin, Coastal Plains, and Basin and Range. Thermal recovery efficiencies are very high in all regions and increase over time. The thermal storage capacity metric is most informative in the pre-assessment and underscores the importance of mapping reservoir thicknesses and porosities to permit detailed mapping of thermal storage capacity per unit area as a key RTES resource classification standard. This assessment provides a basic understanding of the RTES potential in several metropolitan areas and geologic regions throughout the United States and will aid further evaluation of national RTES efficacy.

Conference Paper↗
USGS2019

Interpretation of dye tracing data collected November 13–December 2, 2017, at the Savoy Experimental Watershed as part of the Advanced Groundwater Field Techniques in Karst Terrains course, Savoy, Arkansas

The first course on the use of advanced groundwater field techniques for karst aquifers was conducted November 13–17, 2017, at the University of Arkansas Savoy Experimental Watershed (SEW), which is located on pastures for beef livestock research conducted by the Department of Animal Sciences at the University of Arkansas at Savoy, Arkansas. The SEW is an interdisciplinary, collaborative, long-term research site for the study of animal-waste management in a mantled karst setting. The course focused on advanced field activities appropriate for karst aquifer studies: dye tracing, groundwater/surface-water interactions, geophysical methods, and geochemistry. This report summarizes the data collected and interpreted from the dye tracing part of the November 2017 course, other USGS field courses, and past dye tracing investigations conducted by University of Arkansas students.

Arkansas↗
USGS2018

Hydrogeologic setting, conceptual groundwater flow system, and hydrologic conditions 1995–2010 in Florida and parts of Georgia, Alabama, and South Carolina

The hydrogeologic setting and groundwater flow system in Florida and parts of Georgia, Alabama, and South Carolina is dominated by the highly transmissive Floridan aquifer system. This principal aquifer is a vital source of freshwater for public and domestic supply, as well as for industrial and agricultural uses throughout the southeastern United States. Population growth, increased tourism, and increased agricultural production have led to increased demand on groundwater from the Floridan aquifer system, particularly since 1950. The response of the Floridan aquifer system to these stresses often poses regional challenges for water-resource management that commonly transcend political or jurisdictional boundaries. To help water-resource managers address these regional challenges, the U.S. Geological Survey (USGS) Water Availability and Use Science Program began assessing groundwater availability of the Floridan aquifer system in 2009. The current conceptual groundwater flow system was developed for the Floridan aquifer system and adjacent systems partly on the basis of previously published USGS Regional Aquifer-System Analysis (RASA) studies, specifically many of the potentiometric maps and the modeling efforts in these studies. The Floridan aquifer system extent was divided into eight hydrogeologically distinct subregional groundwater basins delineated on the basis of the estimated predevelopment (circa 1880s) potentiometric surface: (1) Panhandle, (2) Dougherty Plain-Apalachicola, (3) Thomasville-Tallahassee, (4) Southeast Georgia-Northeast Florida-South South Carolina, (5) Suwannee, (6) West-central Florida, (7) East-central Florida, and (8) South Florida. The use of these subregions allows for a more detailed analysis of the individual basins and the groundwater flow system as a whole. The hydrologic conditions and associated groundwater budget were updated relative to previous RASA studies to include additional data collected since the 1980s and to reflect the entire groundwater flow system, including the surficial, intermediate, and Floridan aquifer systems for a contemporary period (1995–2010). Inflow to the groundwater flow system of 33,700 million gallons per day (Mgal/d) was assumed to be exclusively from net recharge (precipitation minus evapotranspiration and surface runoff). Outflow from the groundwater flow system included spring discharge (7,700 Mgal/d) and groundwater withdrawals (5,200 Mgal/d). Estimates for all components of the groundwater system were not possible because of large uncertainties associated with internal leakage, coastal discharge, and discharge to streams and lakes. A numerical modeling analysis is required to improve this hydrologic budget calculation and to forecast future changes in groundwater levels and aquifer storage caused by groundwater withdrawals, land-use change, and the effects of climate variability and change.

Alabama, Florida, Georgia, South Carolina↗
USGS2016

Simulating groundwater flow in karst aquifers with distributed parameter models—Comparison of porous-equivalent media and hybrid flow approaches

Understanding karst aquifers, for purposes of their management and protection, poses unique challenges. Karst aquifers are characterized by groundwater flow through conduits (tertiary porosity), and (or) layers with interconnected pores (secondary porosity) and through intergranular porosity (primary or matrix porosity). Since the late 1960s, advances have been made in the development of numerical computer codes and the use of mathematical model applications towards the understanding of dual (primary [matrix] and secondary [fractures and conduits]) porosity groundwater flow processes, as well as characterization and management of karst aquifers. The Floridan aquifer system (FAS) in Florida and parts of Alabama, Georgia, and South Carolina is composed of a thick sequence of predominantly carbonate rocks. Karst features are present over much of its area, especially in Florida where more than 30 first-magnitude springs occur, numerous sinkholes and submerged conduits have been mapped, and numerous circular lakes within sinkhole depressions are present. Different types of mathematical models have been applied for simulation of the FAS. Most of these models are distributed parameter models based on the assumption that, like a sponge, water flows through connected pores within the aquifer system and can be simulated with the same mathematical methods applied to flow through sand and gravel aquifers; these models are usually referred to as porous-equivalent media models. The partial differential equation solved for groundwater flow is the potential flow equation of fluid mechanics, which is used when flow is dominated by potential energy and has been applied for many fluid problems in which kinetic energy terms are dropped from the differential equation solved. In many groundwater model codes (basic MODFLOW), it is assumed that the water has a constant temperature and density and that flow is laminar, such that kinetic energy has minimal impact on flow. Some models have been developed that incorporate the submerged conduits as a one-dimensional pipe network within the aquifer rather than as discrete, extremely transmissive features in a porous-equivalent medium; these submerged conduit models are usually referred to as hybrid models and may include the capability to simulate both laminar and turbulent flow in the one-dimensional pipe network. Comparisons of the application of a porous-equivalent media model with and without turbulence (MODFLOW-Conduit Flow Process mode 2 and basic MODFLOW, respectively) and a hybrid (MODFLOW-Conduit Flow Process mode 1) model to the Woodville Karst Plain near Tallahassee, Florida, indicated that for annual, monthly, or seasonal average hydrologic conditions, all methods met calibration criteria (matched observed groundwater levels and average flows). Thus, the increased effort required, such as the collection of data on conduit location, to develop a hybrid model and its increased computational burden, is not necessary for simulation of average hydrologic conditions (non-laminar flow effects on simulated head and spring discharge were minimal). However, simulation of a large storm event in the Woodville Karst Plain with daily stress periods indicated that turbulence is important for matching daily springflow hydrographs. Thus, if matching streamflow hydrographs over a storm event is required, the simulation of non-laminar flow and the location of conduits are required. The main challenge in application of the methods and approaches for developing hybrid models relates to the difficulty of mapping conduit networks or having high-quality datasets to calibrate these models. Additionally, hybrid models have long simulation times, which can preclude the use of parameter estimation for calibration. Simulation of contaminant transport that does not account for preferential flow through conduits or extremely permeable zones in any approach is ill-advised. Simulation results in other karst aquifers or other parts of the FAS may differ from the comparison demonstrated herein.

Florida↗
USGS2016

Custom map projections for regional groundwater models

For regional groundwater flow models (areas greater than 100,000 km 2 ), improper choice of map projection parameters can result in model error for boundary conditions dependent on area (recharge or evapotranspiration simulated by application of a rate using cell area from model discretization) and length (rivers simulated with head-dependent flux boundary). Smaller model areas can use local map coordinates, such as State Plane (United States) or Universal Transverse Mercator (correct zone) without introducing large errors. Map projections vary in order to preserve one or more of the following properties: area, shape, distance (length), or direction. Numerous map projections are developed for different purposes as all four properties cannot be preserved simultaneously. Preservation of area and length are most critical for groundwater models. The Albers equal-area conic projection with custom standard parallels, selected by dividing the length north to south by 6 and selecting standard parallels 1/6th above or below the southern and northern extent, preserves both area and length for continental areas in mid latitudes oriented east-west. Custom map projection parameters can also minimize area and length error in non-ideal projections. Additionally, one must also use consistent vertical and horizontal datums for all geographic data. The generalized polygon for the Floridan aquifer system study area (306,247.59 km 2 ) is used to provide quantitative examples of the effect of map projections on length and area with different projections and parameter choices. Use of improper map projection is one model construction problem easily avoided.

Groundwater↗
USGS2015

The current status of mapping karst areas and availability of public sinkhole-risk resources in karst terrains of the United States

Subsidence from sinkhole collapse is a common occurrence in areas underlain by water-soluble rocks such as carbonate and evaporite rocks, typical of karst terrain. Almost all 50 States within the United States (excluding Delaware and Rhode Island) have karst areas, with sinkhole damage highest in Florida, Texas, Alabama, Missouri, Kentucky, Tennessee, and Pennsylvania. A conservative estimate of losses to all types of ground subsidence was $125 million per year in 1997. This estimate may now be low, as review of cost reports from the last 15 years indicates that the cost of karst collapses in the United States averages more than $300 million per year. Knowing when a catastrophic event will occur is not possible; however, understanding where such occurrences are likely is possible. The US Geological Survey has developed and maintains national-scale maps of karst areas and areas prone to sinkhole formation. Several States provide additional resources for their citizens; Alabama, Colorado, Florida, Indiana, Iowa, Kentucky, Minnesota, Missouri, Ohio, and Pennsylvania maintain databases of sinkholes or karst features, with Florida, Kentucky, Missouri, and Ohio providing sinkhole reporting mechanisms for the public.

Hydrogeology Journal↗
USGS2015

Revised hydrogeologic framework of the Floridan aquifer system in Florida and parts of Georgia, Alabama, and South Carolina

The hydrogeologic framework for the Floridan aquifer system has been revised throughout its extent in Florida and parts of Georgia, Alabama, and South Carolina. The updated framework generally conforms to the original framework established by the U.S. Geological Survey in the 1980s, except for adjustments made to the internal boundaries of the Upper and Lower Floridan aquifers and the individual higher and contrasting lower permeability zones within these aquifers. The system behaves as one aquifer over much of its extent; although subdivided vertically into two aquifer units, the Upper and Lower Floridan aquifers. In the previous framework, discontinuous numbered middle confining units (MCUI–VII) were used to subdivide the system. In areas where less-permeable rocks do not occur within the middle part of the system, the system was previously considered one aquifer and named the Upper Floridan aquifer. In intervening years, more detailed data have been collected in local areas, resulting in some of the same lithostratigraphic units in the Floridan aquifer system being assigned to the Upper or Lower Floridan aquifer in different parts of the State of Florida. Additionally, some of the numbered middle confining units are found to have hydraulic properties within the same order of magnitude as the aquifers. A new term “composite unit” is introduced for lithostratigraphic units that cannot be defined as either a confining or aquifer unit over their entire extent. This naming convention is a departure from the previous framework, in that stratigraphy is used to consistently subdivide the aquifer system into upper and lower aquifers across the State of Florida. This lithostratigraphic mapping approach does not change the concept of flow within the system. The revised boundaries of the Floridan aquifer system were mapped by considering results from local studies and regional correlations of lithostratigraphic and hydrogeologic units or zones. Additional zones within the aquifers have been incorporated into the framework to allow finer delineation of permeability variations within the aquifer system. These additional zones can be used to progressively divide the system for assessing groundwater and surface-water interaction, saltwater intrusion, and offshore movement of groundwater at greater detail if necessary. The lateral extent of the updip boundary of the Floridan aquifer system is modified from previous work based on newer data and inclusion of parts of the updip clastic facies. The carbonate and clastic facies form a gradational sequence, generally characterized by limestone of successively younger units that extend progressively farther updip. Because of the gradational nature of the carbonate-clastic sequence, some of the updip clastic aquifers have been included in the Floridan aquifer system, the Southeastern Coastal Plain aquifer system, or both. Thus, the revised updip limit includes some of these clastic facies. Additionally, the updip limit of the most productive part of the Floridan aquifer system was revised and indicates the approximate updip limit of the carbonate facies. The extent and altitude of the freshwater-saltwater interface in the aquifer system has been mapped to define the freshwater part of the flow system.

Alabama, Florida, Georgia, South Carolina↗
USGS2014

Simulation of groundwater flow and saltwater movement in the Onslow County area, North Carolina: predevelopment-2010

Onslow County, North Carolina, is located within the designated Central Coastal Plain Capacity Use Area (CCPCUA). The CCPCUA was designated by law as a result of groundwater level declines of as much as 200 feet during the past four decades within aquifers in rocks of Cretaceous age in the central Coastal Plain of North Carolina and a depletion of water in storage from increased groundwater withdrawals in the area. The declines and depletion of water in storage within the Cretaceous aquifers increase the potential for saltwater migration—both lateral encroachment and upward leakage of brackish water. Within the CCPCUA, a reduction in groundwater withdrawals over a period of 16 years from 2003 to 2018 is mandated. Under the CCPCUA rules, withdrawals in excess of 100,000 gallons per day from any of the Cretaceous aquifer well systems are subject to water-use reductions of as much as 75 percent. To assess the effects of the CCPCUA rules and to assist with groundwater-management decisions, a numerical model was developed to simulate the groundwater flow and chloride concentrations in the surficial Castle Hayne, Beaufort, Peedee, and Black Creek aquifers in the Onslow County area. The model was used to (1) simulate groundwater flow from 1900 to 2010; (2) assess chloride movement throughout the aquifer system; and (3) create hypothetical scenarios of future groundwater development. After calibration of a groundwater flow model and conversion to a variable-density model, five scenarios were created to simulate future groundwater conditions in the Onslow County area: (1) full implementation of the CCPCUA rules with three phases of withdrawal reductions simulated through 2028; (2) implementation of only phase 1 withdrawal reductions of the CCPCUA rules and simulated through 2028; (3) implementation of only phases 1 and 2 withdrawal reductions of the CCPCUA rules and simulated through 2028; (4) full implementation of the CCPCUA rules with the addition of withdrawals from the Castle Hayne aquifer in Onslow County at the fully permitted amount in the final stress period and simulated through 2028; and (5) full implementation of the CCPCUA rules as in scenario 1 except simulated through 2100. Results from the scenarios give an indication of the water-level recovery in the Black Creek aquifer throughout each phase of the CCPCUA rules in Onslow County. Furthermore, as development of the Castle Hayne aquifers was increased in the scenarios, cones of depression were created around pumping centers. Additionally, the scenarios indicated little to no change in chloride concentrations for the time periods simulated.

North Carolina↗
USGS2014

U.S. Geological Survey Karst Interest Group Proceedings, Carlsbad, New Mexico, April 29-May 2, 2014

Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex, and involve biological, chemical, and physical changes. These factors, combined with the diverse climatic regimes under which karst development in these rocks has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs has resulted in the most exceptional karst areas of the United States being designated as national or state parks; commercial caverns and known privately owned caves number in the tens of thousands. Both public and private properties provide access for scientists to study the flow of groundwater in situ. Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and as unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas can impact the ecosystem and water quality of these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Ground-Water Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as other federal and state agencies, and with researchers from academia and institutes. The KIG also encourages younger scientists by participation of students in the poster and oral sessions. To accomplish its mission, the KIG has organized a series of workshops that are held near nationally important karst areas. To date (2014) six KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst features. Proceedings of the workshops are published by the USGS and are available online at http://water.usgs.gov/ogw/karst/kig . The first KIG workshop was held in St. Petersburg, Florida, February 13–16, 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held August 20–22, 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held September 12–15, 2005, in Rapid City, South Dakota, nearby to karst features in evaporites and limestones of the Madison Group in the Black Hills of South Dakota, including Wind Cave National Park and Jewel Cave National Monument. The workshop also included a featured presentation by Thomas Casadevall, Central Region Director, USGS, on the status of earth science at the USGS and evening trips to Jewel Cave led by Mike Wiles, National Park Service (NPS) and Wind Cave led by Rod Horrocks, NPS. The fourth KIG workshop was held May 27–29, 2008, and hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave with Rickard Toomey and Rick Olsen, NPS. The fifth workshop was held April 26–29, 2011, and was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River of the southern Ozarks and a keynote presentation on paleokarst in the United States by Art and Peggy Palmer. This sixth and current 2014 KIG workshop is hosted by the National Cave and Karst Research Institute (NCKRI) in Carlsbad, New Mexico, with Director of NCKRI, George Veni, serving as co-chair of the workshop with Eve Kuniansky, USGS. The session planning committee for this sixth workshop includes Van Brahana, USGS retired and University of Arkansas Professor Emeritus; Tom Byl, USGS and Tennessee State University; Zelda Bailey, former Director of NCKRI and retired Director, National Institute of Standards and Technology, Boulder Laboratory, Colorado; Patrick Tucci, USGS retired; and Mike Bradley, Allan Clark, Geoff Delin, Daniel Doctor, James Kaufmann, Eve Kuniansky, Randy Orndorff, Larry Spangler, and Dave Weary of the USGS. The karst hydrology field trip on Thursday will be led by Lewis Land (NCKRI karst hydrologist) and the optional Friday field trip on the geology of Carlsbad Caverns National Park will be led by George Veni. The keynote speaker is Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at NCKRI, who will address the future of karst research. Additionally, there is a featured presentation “Irish karst and its management,” by Caoimhe Hickey, The Geological Survey of Ireland, preceding a panel discussion on “Collaboration During Times of Limited Resources.” The extended abstracts of USGS authors were peer reviewed and approved for publication by the U.S. Geological Survey. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to our editorial standards or stratigraphic nomenclature and is not research conducted or data collected by the USGS. However, all articles had at a minimum of two peer reviews, and all articles were edited for consistency of appearance in the published Proceedings. The use of trade, firm or product names in any article is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS, Office of Groundwater, provides technical support for the Karst Interest Group website and public availability of the Proceedings from these workshops, and the USGS Groundwater Resources Program funds the publication costs. Finally, the cover illustration is the work of Ann Tihansky, USGS, used since the first KIG workshop in 2000.

Scientific Investigations Report↗
USGS2014

Taking the mystery out of mathematical model applications to karst aquifers—A primer

Advances in mathematical model applications toward the understanding of the complex flow, characterization, and water-supply management issues for karst aquifers have occurred in recent years. Different types of mathematical models can be applied successfully if appropriate information is available and the problems are adequately identified. The mathematical approaches discussed in this paper are divided into three major categories: 1) distributed parameter models, 2) lumped parameter models, and 3) fitting models. The modeling approaches are described conceptually with examples (but without equations) to help non-mathematicians understand the applications.

Conference Paper↗
USGS2012

Tabulated Transmissivity and Storage Properties of the Floridan Aquifer System in Florida and Parts of Georgia, South Carolina, and Alabama

A goal of the U.S. Geological Survey Groundwater Resources Program is to assess the availability of fresh water within each of the principal aquifers in the United States with the greatest groundwater withdrawals. The Floridan aquifer system (FAS), which covers an area of approximately 100,000 square miles in Florida and parts of Georgia, Alabama, Mississippi, and South Carolina, is one such principal aquifer, having the fifth largest groundwater withdrawals in the Nation, totaling 3.64 billion gallons per day in 2000. Compilation of FAS hydraulic properties is critical to the development and calibration of groundwater flow models that can be used to develop water budgets spatially and temporally, as well as to evaluate resource changes over time. Wells with aquifer test data were identified as Upper Floridan aquifer (UFA), Lower Floridan aquifer (LFA), Floridan aquifer system (FAS, Upper Floridan with some middle and/or Lower Floridan), or middle Floridan confining unit (MCU), based on the identification from the original database or report description, or comparison of the open interval of the well with previously published maps. This report consolidates aquifer hydraulic property data obtained from multiple databases and reports of the U.S. Geological Survey, various State agencies, and the Water Management Districts of Florida, that are compiled into tables to provide a single information source for transmissivity and storage properties of the FAS as of October 2011. Transmissivity calculated from aquifer pumping tests and specific-capacity data are included. Values for transmissivity and storage coefficients are intended for use in regional or sub regional groundwater flow models; thus, any tests (aquifer pumping tests and specific capacity data) that were conducted with packers or for open intervals less than 30 feet in length are excluded from the summary statistics and tables of this report, but are included in the database. The transmissivity distribution from the aquifer pumping tests is highly variable. The transmissivity based on aquifer pumping tests (from 1,045 values for the UFA and FAS) ranges from 8 to about 9,300,000 square feet per day (ft 2 /d) and values of storage coefficient (646 reported) range from 3x10 -9 to 0.41. The 64 transmissivity values for the LFA range from about 130 to 4,500,000 ft 2 /d, and the 17 storage coefficient values range from 7x10 -8 to 0.03. The 14 transmissivity values for the MCU range from 1 to about 600,000 ft 2 /d and the 10 storage coefficient values range from 8x10 -8 to 0.03. Transmissivity estimates for the UFA and FAS for 442 specific capacity tests range from approximately 200 to 1,000,000 ft 2 /d.

Alabama, Florida, Georgia, South Carolina↗
USGS2012

Transmissivity of the Upper Floridan aquifer in Florida and parts of Georgia, South Carolina, and Alabama

The Floridan aquifer system (FAS) covers an area of approximately 100,000 square miles in Florida and parts of Georgia, South Carolina, Alabama, and Mississippi. Groundwater wells for water supply were first drilled in the late 1800s and by the year 2000, the FAS was the primary source of drinking water for about 10 million people. One of the methods for assessing groundwater availability is the development of regional or subregional groundwater flow models of the aquifer system that can be used to develop water budgets spatially and temporally, as well as evaluate the groundwater resource change over time. Understanding the distribution of transmissivity within the FAS is critical to the development of groundwater flow models. The map presented herein differs from previously published maps of the FAS in that it is based on interpolation of 1,487 values of transmissivity. The transmissivity values in the dataset range from 8 to 9,000,000 feet squared per day (ft 2 /d) with the majority of the values ranging from 10,000 to 100,000 ft 2 /d. The wide range in transmissivity (6 orders of magnitude) is typical of carbonate rock aquifers, which are characterized by a wide range in karstification. Commonly, the range in transmissivity is greatest in areas where groundwater flow creates conduits in facies that dissolve more readily or areas of high porosity units that have interconnected vugs, with diameters greater than 0.1 foot. These are also areas where transmissivity is largest. Additionally, first magnitude springsheds and springs are shown because in these springshed areas, the estimates of transmissivity from interpolation may underestimate the actual range in transmissivity. Also shown is an area within the Gulf Trough in Georgia where high yielding wells are unlikely to be developed in the Upper Floridan aquifer. The interpolated transmissivity ranges shown on this map reflect the geologic structure and karstified areas. Transmissivity is large in the areas where the system is unconfined, such as west-central Florida and southwest Georgia just northwest of the Gulf Trough. Transmissivity is small along the Gulf Trough and Southwest Georgia Embayment (referred to as Apalachicola Embayment in some reports). Transmissivity is also small in the thin, updip part of the system near its northern boundary. Another area of large transmissivity coincides with the Southeast Georgia Embayment.

Alabama, Florida, Georgia, South Carolina↗
USGS2011

Effects of model layer simplification using composite hydraulic properties

Groundwater provides much of the fresh drinking water to more than 1.5 billion people in the world (Clarke et al., 1996) and in the United States more that 50 percent of citizens rely on groundwater for drinking water (Solley et al., 1998). As aquifer systems are developed for water supply, the hydrologic system is changed. Water pumped from the aquifer system initially can come from some combination of inducing more recharge, water permanently removed from storage, and decreased groundwater discharge. Once a new equilibrium is achieved, all of the pumpage must come from induced recharge and decreased discharge (Alley et al., 1999). Further development of groundwater resources may result in reductions of surface water runoff and base flows. Competing demands for groundwater resources require good management. Adequate data to characterize the aquifers and confining units of the system, like hydrologic boundaries, groundwater levels, streamflow, and groundwater pumping and climatic data for recharge estimation are to be collected in order to quantify the effects of groundwater withdrawals on wetlands, streams, and lakes. Once collected, three-dimensional (3D) groundwater flow models can be developed and calibrated and used as a tool for groundwater management. The main hydraulic parameters that comprise a regional or subregional model of an aquifer system are the hydraulic conductivity and storage properties of the aquifers and confining units (hydrogeologic units) that confine the system. Many 3D groundwater flow models used to help assess groundwater/surface-water interactions require calculating ?effective? or composite hydraulic properties of multilayered lithologic units within a hydrogeologic unit. The calculation of composite hydraulic properties stems from the need to characterize groundwater flow using coarse model layering in order to reduce simulation times while still representing the flow through the system accurately. The accuracy of flow models with simplified layering and hydraulic properties will depend on the effectiveness of the methods used to determine composite hydraulic properties from a number of lithologic units.

Florida↗
USGS2010

Effects of Changes in Irrigation Practices and Aquifer Development on Groundwater Discharge to the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico

Since 1990, about 75 acres of black mangroves have died in the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico. Although many factors can contribute to the mortality of mangroves, changes in irrigation practices, rainfall, and water use resulted in as much as 25 feet of drawdown in the potentiometric surface of the aquifer in the vicinity of the reserve between 1986 and 2002. To clarify the issue, the U.S. Geological Survey, in cooperation with the Puerto Rico Department of Natural and Environmental Resources, conducted a study to ascertain how aquifer development and changes in irrigation practices have affected groundwater levels and groundwater flow to the Mar Negro area of the reserve. Changes in groundwater flow to the mangrove swamp and bay from 1986 to 2004 were estimated in this study by developing and calibrating a numerical groundwater flow model. The transient simulations indicate that prior to 1994, high irrigation return flows more than offset the effect of reduced groundwater withdrawals. In this case, the simulated discharge to the coast in the modeled area was 19 million gallons per day. From 1994 through 2004, furrow irrigation was completely replaced by micro-drip irrigation, thus eliminating return flows and the simulated average coastal discharge was 7 million gallons per day, a reduction of 63 percent. The simulated average groundwater discharge to the coastal mangrove swamps in the reserve from 1986 to 1993 was 2 million gallons per day, compared to an average simulated discharge of 0.2 million gallons per day from 1994 to 2004. The average annual rainfall for each of these periods was 38 inches. The groundwater discharge to the coastal mangrove swamps in the Jobos Bay National Estuarine Research Reserve was estimated at about 0.5 million gallons per day for 2003-2004 because of higher than average annual rainfall during these 2 years. The groundwater flow model was used to test five alternatives for increasing groundwater discharge to the coastal mangrove swamps to approximately 1.4 million gallons per day: (1) artificially recharging the aquifer with injection wells or (2) by increasing irrigation return flow by going back to furrow irrigation; (3) termination of groundwater withdrawals near the mangroves; (4) reduction of groundwater withdrawals at irrigation wells by 50 percent; and (5) a combination of alternatives 2 and 4 increasing irrigation return flows and decreasing irrigation withdrawals. Each alternative assumed average climatic conditions and groundwater withdrawals at 2004 rates. Alternative 1 required 1.5 million gallons per day of injected water. Alternative 2 required flooding 958 acres with a rate of 1.84 million gallons per day if no crops are grown. Alternative 3 required the termination of 2.44 million gallons per day of withdrawals to achieve 1.34 million gallons per day of discharge to the mangroves. Alternative 4 did not achieve the objective with only 0.80 million gallons per day simulated discharge to the mangroves, while requiring a 1.26 million gallon per day reduction in groundwater withdrawals. Alternative 5 required flooding fields with additional 1.13 million gallons of day and the same reduction in groundwater withdrawals, but did achieve the objective of about 1.4 million gallons per day discharge to the mangroves. Alternative 1, incorporating injection wells near the reserve required the least amount of water to raise groundwater levels and maintain discharge of 1.4 million gallons per day through the mangroves.

Scientific Investigations Report↗
USGS2009

Effects of model layer simplification using composite hydraulic properties

The effects of simplifying hydraulic property layering within an unconfined aquifer and the underlying confining unit were assessed. The hydraulic properties of lithologic units within the unconfined aquifer and confining unit were computed by analyzing the aquifer-test data using radial, axisymmetric two-dimensional (2D) flow. Time-varying recharge to the unconfined aquifer and pumping from the confined Upper Floridan aquifer (USA) were simulated using 3D flow. Conceptual flow models were developed by gradually reducing the number of lithologic units in the unconfined aquifer and confining unit by calculating composite hydraulic properties for the simplified lithologic units. Composite hydraulic properties were calculated using either thickness-weighted averages or inverse modeling using regression-based parameter estimation. No significant residuals were simulated when all lithologic units comprising the unconfined aquifer were simulated as one layer. The largest residuals occurred when the unconfined aquifer and confining unit were aggregated into a single layer (quasi-3D), with residuals over 100% for the leakage rates to the confined aquifer and the heads in the confining unit. Residuals increased with contrasts in vertical hydraulic conductivity between the unconfined aquifer and confining unit. Residuals increased when the constant-head boundary at the bottom of the Upper Floridan aquifer was replaced with a no-flow boundary.

Florida↗