Geology ReportsSearch

SEARCH · Geology Reports

Results for “Florida Geological Survey Report of Investigations”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13Linked to original sources

Drainage infrastructure and groundwater system response to changes in sea level and precipitation, Broward County, Florida

Executive Summary Southeast Florida is highly susceptible to flooding because of its low topography and porous, highly permeable Biscayne aquifer. Rising seas will likely result in increased groundwater levels in parts of Broward County, Florida, that will reduce available soil storage and therefore increase the likelihood of inundation and flooding from precipitation events. Moreover, rising seas may also reduce the capacity of the coastal water-control structures to discharge inland waters to tidal areas, thereby increasing surface-water stage and nearby groundwater levels. Increased rainfall intensity will likely further increase peak surface-water stages and groundwater levels, more quickly fill the reduced soil storage capacity, and increase the likelihood for inundation. Managers and planners in Broward County, Florida, face the challenge of understanding and preparing for the consequent risk to residents, businesses, and critical infrastructure posed by increased sea level and precipitation. The U.S. Geological Survey, in cooperation with the Broward County Environmental Planning and Resilience Division, has developed a groundwater/surface-water model to evaluate the response of the drainage infrastructure and groundwater system in Broward County to projected increases in sea level and potential changes in precipitation. The model was constructed using Modular Finite-Difference Groundwater Flow Model Newton formulation, with the surface-water system represented using the Surface-Water Routing process and the Urban Runoff process. The aquifer layering and flow parameters rely heavily on existing hydrologic flow models developed by the U.S. Geological Survey for the same model area. The surface-water drainage system within this newly developed model actively simulates the extensive canal network using level-pool routing and active structures representing gates, weirs, culverts, and pumps. Steady-state and transient simulation results represented historical conditions (2013–17). The simulated historical groundwater levels and upstream stage and flow at the primary structures generally captured the behavior of the actual hydrologic system. Simulation results incorporating increased sea level and precipitation were used to evaluate the effects of these projected changes on the surface-water drainage system and wet season groundwater levels. Four future sea-level scenarios were simulated by modifying the historical inputs for both steady-state and the transient versions of the model to represent mean sea levels of 0.5, 2.0, 2.5, and 3.0 feet (ft) above the North American Vertical Datum of 1988. These mean sea levels correspond to sea-level rises of 1.05, 2.55, 3.05, and 3.55 ft, respectively, above the 2013–17 mean measured tidal stage. Additional simulations represented a 15-percent increase in rainfall rates using the transient model and a 15-percent increase in rainfall recharge using the steady-state model. The simulated results indicated that (1) the effects of increased sea level were more evident in the easternmost, coastal areas of the county where increases in groundwater levels are nearly equivalent to sea-level rise; (2) groundwater levels west of the coastal water-control structures only changed slightly in response to increased sea level for most scenarios; (3) when the control elevations of the gravity-controlled coastal water-control structures were surpassed by sea-level rise, the resulting increases in upstream stage in the connected primary canal resulted in increased groundwater levels that can propagate into the western parts of the county; (4) the historical west-to-east downward gradient in groundwater levels decreased with increased sea level, and groundwater levels were lower in central parts of the county than areas west and east for the higher sea-level scenarios; (5) simulated upstream stage for most of the primary coastal water-control structures increased with increased sea level, with the largest increases occurring at gravity-controlled structures having the lowest control elevations; (6) total flow through the primary structures increased as sea level increased because of additional groundwater leakage into the surface-water network; (7) the 3.0-ft mean sea-level rise scenario resulted in an increase of 37.10 square miles in area having a wet season average depth to groundwater of less than 2 ft, and an increase of 22.84 square miles in newly inundated areas compared to historical simulation results; and (8) a 15-percent increase in rainfall rate for the entire simulation period produced little increase in upstream stages at the primary structures and an increase in total flow through the primary structures proportional to the increase in rainfall.

Florida

Water use in the Apalachicola-Chattahoochee-Flint River Basin, Alabama, Florida, and Georgia, 2010, and water-use trends, 1985-2010

The Apalachicola-Chattahoochee-Flint (ACF) River Basin encompasses about 20,230 square miles in parts of Alabama, Florida, and Georgia. Increasing population growth and agricultural production from the 1970s to 2010 has prompted increases in water-resources development and substantially increased water demand in the basin. Since the 1980s, Alabama, Florida, Georgia, and the U.S. Army Corps of Engineers are parties to litigation concerning water management in the ACF River Basin. Estimating the 2010 water use in the ACF River Basin is one aspect of a multipart water resources study on the ACF River Basin that began in 2011. This ACF River Basin study is one focus area of the U.S. Geological Survey’s National Water Census program. The 2010 water-use estimates for the ACF River Basin are presented in this report. These estimates include an inventory of the quantity and sources of water withdrawn by category of use and location (State and river basin), and the surface-water returns in the ACF River Basin during 2010. Water-use trends from 1985 to 2010 in the basin also are presented. Offstream water-withdrawal data in the ACF River Basin are presented for each of the following categories: public supply, self-supplied domestic, self-supplied commercial, industrial, mining, agricultural (including crop irrigation, livestock, and aquaculture uses), and thermoelectric-power generation. Water-use data are compiled for the 14 subbasins in the ACF River Basin. For the counties in Alabama, Florida, and Georgia that are partially within the ACF River Basin, data are presented for only that part of the county that lies within the basin. A variety of Federal, State, local, private, and online sources in Alabama, Florida, and Georgia were used to gather surface-water and groundwater withdrawal, surface-water discharges (return flows), and water-use data for the ACF River Basin in 2010. The population in the ACF River Basin was 3.835 million in 2010, a 45-percent increase from the 1990 population of nearly 2.636 million. About 92 percent of the 2010 ACF population resided in Georgia with nearly 75 percent living in the Atlanta metropolitan area. In 2010, 1,645 million gallons per day (Mgal/d) of water were withdrawn from groundwater (576 Mgal/d) and surface-water (1,069 Mgal/d) sources in the ACF River Basin. About 89 percent of the groundwater and 83 percent of the surface-water withdrawals were from Georgia. About 5.6 percent of the total groundwater and nearly 4 percent of the total surface-water withdrawals in the ACF River Basin were from Florida, whereas about 5.3 percent of groundwater and nearly 16 percent of surface water were withdrawn in Alabama. Total water use (withdrawals plus public-supplied deliveries) in the ACF River Basin was 1,593 Mgal/d in 2010. About 56 Mgal/d of water withdrawn in the ACF River Basin was delivered (interbasin transfer) to basins beyond the ACF River Basin. About 564 Mgal/d of water was returned to surface-water bodies in the ACF River Basin. Most of that amount, 63 percent, was treated wastewater discharged by public wastewater-treatment facilities. Water used for once-through cooling by thermoelectric-power facilities accounted for nearly 24 percent of the surface-water returns in the basin. About 70 percent of all water withdrawals in the ACF River Basin were by self-supplied agricultural water users and public water suppliers. Agricultural withdrawals were greatest in the Flint River Basin (501 Mgal/d) with ground-water representing 84 percent of the withdrawals from that basin. Within the Flint River Basin, agricultural withdrawals were greatest in the Lower Flint River and Spring Creek subbasins. About 3.52 million people were served by public water suppliers in the ACF River Basin during 2010, and 88 percent of that population used surface water. Georgia had the largest public-supplied population, representing nearly 93 percent (3.17 million) of the public-supplied population in the ACF River Basin. Public water suppliers served 193,700 people (5.7 percent) in Alabama and 31,880 people in Florida (1.3 percent). Public-supply losses were estimated at 101 Mgal/d. Withdrawals for public supply (483 Mgal/d) and self-supplied industry (141 Mgal/d) were greatest in the Chattahoochee River Basin. Surface water accounted for 96 percent of all withdrawals in the Chattahoochee River Basin. Withdrawals for public supply were greatest in the Upper Chattahoochee River subbasin (366 Mgal/d), whereas self-supplied industrial withdrawals were greatest in the Lower Chattahoochee River subbasin (110 Mgal/d). Water-use trends in the ACF River Basin have varied during the 25 years between 1985 and 2010. Surface-water withdrawals declined between 1985 and 2000, sharply increased in 2000, and declined again between 2000 and 2010. In contrast, groundwater withdrawals increased between 1985 and 2000, declined in 2005, and increased between 2005 and 2010.

Alabama, Florida, Georgia

Ground-cover vegetation in wetland forests of the lower Suwannee River floodplain, Florida, and potential impacts of flow reductions

Ground-cover vegetation was surveyed in wetland forests in the lower Suwannee River floodplain, Florida, in a study conducted by the U.S. Geological Survey in cooperation with the Suwannee River Water Management District from 1996 to 1999. Increased water use in the basin, supplied primarily from ground water, could reduce ground-water discharge to the river and flows in the lower Suwannee River. Many of the 282 ground-cover species found in wetland forests of the floodplain have distributions that are related to flow-dependent hydrologic characteristics of forest types, and their distributions would change if flows were reduced. Overall species diversity in the floodplain might decrease, and the composition of ground-cover vegetation in all forest types might change with flow reductions. The study area included forests within the 10-year floodplain of the lower Suwannee River from its confluence with the Santa Fe River to the lower limit of forests near the Gulf of Mexico. The floodplain is divided into three reaches (riverine, upper tidal, and lower tidal) due to variations in hydrology, vegetation, and soils with proximity to the coast. The riverine (non-tidal) reach had the greatest number of total species (203) and species unique to that reach (81). Mitchella repens, Toxicodendron radicans, and Axonopus furcatus were the most frequently dominant species in riverine bottomland hardwoods. Free-floating aquatic species, such as Spirodela punctata and Lemna valdiviana, were the dominant species in the wettest riverine swamps. The upper tidal reach had the lowest number of total species (116), only two species unique to that reach, and the lowest density of ground cover (26 percent). Panicum commutatum and Crinum americanum were frequent dominant species in upper tidal forests. The lower tidal reach had the highest ground-cover density (43 percent) and the second highest number of total species (183) and number of species unique to that reach (55). Saururus cernuus and species of Carex were frequently dominant in lower tidal swamps. Lower tidal hammocks, the most elevated lower tidal forests, were dominated by Osmunda cinnamomea and Chasmanthium laxum. Flow reductions in the lower Suwannee River could change the flow-dependent hydrologic characteristics of wetland forests. Decreases in inundation and saturation in riverine forests could result in a decrease in the number and extent of semi-permanently inundated ponds. As a result, several species of free-floating, aquatic plants that grow only in riverine floodplain ponds might decrease in abundance or disappear if flows were reduced. Decreases in inundation and saturation could also result in a shift to more upland species in all riverine forests and upper tidal bottomland hardwoods. Upland species and some exotic species might increase in abundance in the floodplain, invading forests where hydrologic conditions have been altered by flow reductions. Depth and duration of inundation due to river flooding could decrease in all riverine and upper tidal forests, probably resulting in a shift of species to those that are typically found in forests with shallower, shorter-duration floods. Salinity in the lower tidal reach and adjacent areas of the upper tidal reach might increase with flow reductions, and the distribution of species might change due to varying tolerances of salinity among species. Species with low salt-tolerance unique to the lower tidal reach might disappear from the floodplain, and species with high salinity tolerance could increase in abundance, replacing less salt-tolerant species.

Water-Resources Investigations Report

Hydrology, vegetation, and soils of four north Florida River flood plains with an evaluation of state and federal wetland determinations

A study of hydrologic conditions, vegetation, and soils was made in wetland forests of four north Florida streams from 1987 to 1990. The study was conducted by the U.S. Geological Survey in cooperation with the Florida Department of Environmental Regulation to support State and Federal efforts to improve wetland delineation methodology in flood plains. Plant communities and soils were described and related to topographic position and long-term hydrologic conditions at 10 study plots located on 4 streams. Detailed appendixes give average duration, frequency, and depth of flooding; canopy, subcanopy, and ground-cover vegetation; and taxonomic classification, series, and profile descriptions of soils for each plot. Topographic relief, range in stage, and depth of flooding were greatest on the alluvial flood plain of the Ochlockonee River, the largest of the four streams. Soils were silty in the lower elevations of the flood plain, and tree communities were distinctly different in each topographic zone. The Aucilla River flood plain was dominated by levees and terraces with very few depressions or low backwater areas. Oaks dominated the canopy of both lower and upper terraces of the Aucilla flood plain. Telogia Creek is a blackwater stream that is a major tributary of the Ochlockonee River. Its low, wet flood plain was dominated by Wyssa ogeche (Ogeechee tupelo) trees, had soils with mucky horizons, and was inundated by frequent floods of very short duration. The St. Marks River, a spring-fed stream with high base flow, had the least topographic relief and lowest range in stage of the four streams. St. Marks soils had a higher clay content than the other streams, and limestone bedrock was relatively close to the surface. Wetland determinations of the study plots based on State and Federal regulatory criteria were evaluated. Most State and Federal wetland determinations are based primarily on vegetation and soil characteristics because hydrologic records are usually not available. In this study, plots were located near long-term gaging stations, thus wetland determinations based on plant and soil characteristics could be evaluated at sites where long-term hydrologic conditions were known. Inconsistencies among hydrology, vegetation, and soil determinations were greatest on levee communities of the Ochlockonee and Aucilla River flood plains. Duration of average annual longest flood was almost 2 weeks for both plots. The wetland species list currently used (1991) by the State lacks many ground-cover species common to forested flood plains of north Florida rivers. There were 102 ground-cover species considered upland plants by the State that were present on the nine annually flooded plots of this study. Among them were 34 species that grew in areas continuously flooded for an average of 5 weeks or more each year. Common flood-plain species considered upland plants by the State were: Hypoxis leptocarpa (yellow star-grass), and two woody vines, Brunnichia ovata (ladies' eardrops) and Campsis radican s (trumpet-creeper), which were common in areas flooded continuously for 6 to 9 weeks a year; Sebastiania fruticosa (Sebastian-bush), Chasmanthium laxum (spikegrass), and Panicum dichotomum (panic grass), which typically grew in areas flooded an average of 2 to 3 weeks or more per year; Vitis rotundifolia (muscadine) and Toxicodendron radicans (poison-ivy), usually occurring in areas flooded an average of 1 to 2 weeks a year; and Quercus virginiana (live oak) present most often in areas flooded approximately 1 week a year. Federal wetland regulations (1989) limited wetland jurisdiction to only those areas that are inundated or saturated during the growing season. However, year-round hydrologic records were chosen in this report to describe the influence of hydrology on vegetation, because saturation, inundation, or flowing water can have a variety of both beneficial and adverse effects on flood-plain vegetation at any time of the year. These effects can occur because: (1) Soil temperatures in north Florida are probably high enough in winter for anaerobic conditions to develop in saturated soils. (2) Many plants in the flood plains of north Florida are active in the nongrowing season and might be adapted to anaerobic conditions in winter. (3) Other effects of standing water, such as decreased light penetration, prevention of seed germination, and protection of seeds from herbivores; and effects of flowing water, such as scouring and deposition of sediments, seed dispersal, and mechanical injury, can occur any time of year.

Florida

Hydrology of Polk County, Florida

Local water managers usually rely on information produced at the State and regional scale to make water-resource management decisions. Current assessments of hydrologic and water-quality conditions in Polk County, Florida, commonly end at the boundaries of two water management districts (South Florida Water Management District and the Southwest Florida Water Management District), which makes it difficult for managers to determine conditions throughout the county. The last comprehensive water-resources assessment of Polk County was published almost 40 years ago. To address the need for current countywide information, the U.S. Geological Survey began a 3?-year study in 2002 to update information about hydrologic and water-quality conditions in Polk County and identify changes that have occurred. Ground-water use in Polk County has decreased substantially since 1965. In 1965, total ground-water withdrawals in the county were about 350 million gallons per day. In 2002, withdrawals totaled about 285 million gallons per day, of which nearly 95 percent was from the Floridan aquifer system. Water-conservation practices mainly related to the phosphate-mining industry as well as the decrease in the number of mines in operation in Polk County have reduced total water use by about 65 million gallons per day since 1965. Polk County is underlain by three principal hydrogeologic units. The uppermost water-bearing unit is the surficial aquifer system, which is unconfined and composed primarily of clastic deposits. The surficial aquifer system is underlain by the intermediate confining unit, which grades into the intermediate aquifer system and consists of up to two water-bearing zones composed of interbedded clastic and carbonate rocks. The lowermost hydrogeologic unit is the Floridan aquifer system. The Floridan aquifer system, a thick sequence of permeable limestone and dolostone, consists of the Upper Floridan aquifer, a middle semiconfining unit, a middle confining unit, and the Lower Floridan aquifer. The Upper Floridan aquifer provides most of the water required to meet demand in Polk County. Data from about 300 geophysical and geologic logs were used to construct hydrogeologic maps showing the tops and thicknesses of the aquifers and confining units within Polk County. Thickness of the surficial aquifer system ranges from several feet thick or less in the extreme northwestern part of the county and along parts of the Peace River south of Bartow to more than 200 feet along the southern part of the Lake Wales Ridge in eastern Polk County. Thickness of the intermediate aquifer system/intermediate confining unit is highly variable throughout the county because of past erosional processes and sinkhole formation. Thickness of the unit ranges from less than 25 feet in the extreme northwestern part of the county to more than 300 feet in southwestern Polk County. The altitude of the top of the Upper Floridan aquifer in the county ranges from about 50 feet above National Geodetic Vertical Datum of 1929 (NGVD 29) in the northwestern part to more than 250 feet below NGVD 29 in the southern part. Water levels in the Upper Floridan aquifer fluctuate seasonally, increasing during the wet season (June through September) and decreasing during the rest of the year. Water levels in the Upper Floridan aquifer also can change from year to year, depending on such factors as pumpage and climatic variations. In the southwestern part of the county, fluctuations in water use related to phosphate mining have had a major impact on ground-water levels. Hydrographs of selected wells in southwestern Polk County show a general decline in water levels that ended in the mid-1970s. This water-level decline coincides with an increase in water use associated with phosphate mining. A substantial increase in water levels that began in the mid-1970s coincides with a period of decreasing water use in the county. Despite reductions in water use since 1970, howev

Scientific Investigations Report

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida

Drier forest composition associated with hydrologic change in the Apalachicola River floodplain, Florida

Forests of the Apalachicola River floodplain had shorter flood durations, were drier in composition, and had 17 percent fewer trees in 2004 than in 1976. The change to drier forest composition is expected to continue for at least 80 more years. Floodplain drying was caused by large declines in river levels resulting from erosion of the river channel after 1954 and from decreased flows in spring and summer months since the 1970s. Water-level declines have been greatest at low and medium flows, which are the most common flows (occurring about 80 percent of the time). Water levels have remained relatively unchanged during large floods which continue to occur about three times per decade. A study conducted by the U.S. Geological Survey compared temporal changes in hydrologic conditions, forest composition, forest characteristics, and individual species of trees, as well as estimated the potential for change in composition of floodplain forests in the nontidal reach of the Apalachicola River. The study was conducted with the cooperation of the Florida Department of Environmental Protection and the Northwest Florida Water Management District. Forest composition and field observations from studies conducted in 1976-1984 (termed “1976 data”) were used as baseline data for comparison with data from plots sampled in 2004-2006 (“2004 data”). Flood durations were shorter in all periods subsequent to 1923-1976. The periods of record used to calculate flood durations for forest data were subsets of the complete record available (1923-2004). At sampled plots in all forest types and reaches combined, flood durations changed an average of more than 70 percent toward the baseline flood duration of the next drier forest type. For all forest types, changes in flood durations toward the next drier type were greatest in the upper reach (95.9 percent) and least in the lower reach (42.0 percent). All forests are expected to be 38.2 percent drier in species composition by 2085, the year when the median age of surviving 2004 subcanopy trees will reach the median age (99 years) of the 2004 large canopy trees. The change will be greatest for forests in the upper reach (45.0 percent). Forest composition changes from pre-1954 to 2085 were calculated using Floodplain Indices from 1976 and 2004 tree-size classes and replicate plots. Species composition in high bottomland hardwood forests is expected to continue to change, and some low bottomland hardwood forests are expected to become high bottomland hardwood forests. Organisms associated with floodplain forests will be affected by the changes in tree species, which will alter the timing of leaf-out, fruiting, and leaf-drop, the types of fruit and debris produced, and soil chemistry. Swamps will contain more bottomland hardwood species, but will also have an overall loss of tree density. The density of trees in swamps significantly decreased by 37 percent from 1976 to 2004. Of the estimated 4.3 million (17 percent) fewer trees that existed in the nontidal floodplain in 2004 than in 1976, 3.3 million trees belonged to four swamp species: popash, Ogeechee tupelo, water tupelo, and bald cypress. Water tupelo, the most important tree in the nontidal floodplain in terms of basal area and density, has declined in number of trees by nearly 20 percent since 1976. Ogeechee tupelo, the species valuable to the tupelo honey industry, has declined in number of trees by at least 44 percent. Greater hydrologic variability in recent years may be the reason swamps have had a large decrease in tree density. Drier conditions are detrimental for the growth of swamp species, and periodic large floods kill invading bottomland hardwood trees. The loss of canopy density in swamps may result in the swamp floor being exposed to more light with an increase in the amount of ground cover present, which in turn, would reduce tree replacement. The microclimate of the swamp floor would become warmer due to the decrease in shade and inundation. Soils would become dehydrated more quickly in dry periods and debris would decompose more quickly. A loss of tree density in swamps would lead to a decrease in tree and leaf litter biomass, which would have additional effects on swamp organisms. The loss of litter would result in a loss of substrate for benthic organisms in the floodplain and, ultimately, in the downstream waters of the river and estuary.

Florida

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The East-Central Florida Transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration (ET), runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into ET, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average ET over the 1995 to 2006 period was 34.47 inches per year, compared to the calculated average ET rate of 36.39 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.58 inches per year, compared with the calculated average of 3.39 inches per year from the model-independent water-budget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 920 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida

Assessment of saltwater intrusion in southern coastal Broward County, Florida

Of the counties in southeastern Florida, Broward County has experienced some of the most severe effects of saltwater intrusion into the surficial Biscayne aquifer because, before 1950, most public water-supply well fields in the county were constructed near the principal early population centers located less than 5 miles from the Atlantic Ocean. The construction of major regional drainage canals in the early 20th century caused a lowering of the water table and a gradual inland movement of the saltwater front toward the well fields. The U.S. Geological Survey began field investigations of saltwater intrusion in the Biscayne aquifer of southeastern Broward County in 1939. As part of the present study, the positions of the saltwater front in 1945, 1969, and 1993 were estimated using chloride concentrations of water samples collected between 1939 and 1994 from various monitoring and exploratory wells. The data indicate that, between 1945 and 1993, the saltwater front has moved as much as 0.5 mile inland in parts of the study area. The position and movement of the saltwater front were simulated numerically to help determine which of the various hydrologic factors and water-management features characterizing the coastal subsurface environment and its alteration by man are of significance in increasing or decreasing the degree of saltwater intrusion. Two representational methods were applied by the selection and use of appropriate model codes. The SHARP code simulates the position of the saltwater front as a sharp interface, which implies that no transition zone (a zone in which a gradational change between freshwater and saltwater occurs) separates freshwater and saltwater. The Subsurface Waste Injection Program (SWIP) code simulates a two-fluid, variable-density system using a convective-diffusion approach that includes a representation of the transition zone that occurs between the freshwater and saltwater bodies. The models were applied to: (1) approximately replicate predevelopment and current positions of the interface in the study area; and (2) study the relative importance of various factors affecting the interface position. The model analyses assumed a conceptual model of uniform easterly flow in the aquifer toward points of offshore discharge to tidewater. Measurements of water-table altitude and the depth to the interface in the study area exhibit an interrelation that differes substantially from the classical Ghyben-Herzberg relation. However, both model codes simulated water-table altitudes and interface positions that were generally consistent with the Ghyben-Herzberg relation but differed substantially from observed data. The simulate interface positions were inland of the known positions, and simulate water-table altitudes were higher than measured ones. The SHARP and SWIP simulations were in general agreement with each other when a low value of longitudinal dispersivity was specified in the SWIP simulation and also for higher values of longitudinal dispersivity when modified dispersion algorithms were used in SWIP that greatly reduced the simulated degree of vertical dispersion. Sensitivity analyses performed using the SHARP code indicated simulation results to be relatively insensitive to a substantial change in the specified slope of the base of the aquifer and moderately sensitive to a 150-percent change in net atmospheric recharge to the aquifer (rainfall minus evapotranspiration). Representing well-field pumping by the City of hallandale had only a minor, localized influence on the simulated regional interface position. Using various cross-sectional grid designs in applications of the SWIP code, near convergence of all lines of equal concentrations in the transition zone was achieved within a simulation time of 10 years. The simulated equilibrium interface location was sensitive to substantial spatial variations in the specified hydraulic conductivity values, but was relatively insensitive to seasonal varying

Florida

Potential for increased inundation in flood-prone regions of southeast Florida in response to climate and sea-level changes in Broward County, Florida, 2060–69

The U.S. Geological Survey, in cooperation with Broward County Environmental Planning and Resilience Division, has developed county-scale and local-scale groundwater/surface-water models to study the potential for increased inundation and flooding in eastern Broward County that are due to changes in future climate and sea-level rise. These models were constructed by using MODFLOW 2005, with the surface-water system represented by using the Surface-Water Routing process and a new Urban Runoff process. The local-scale model allowed the use of finer grid resolution in a selected area of the county, whereas the county-scale model provided boundary conditions for the local-scale model and insight into the hydrologic behavior of the larger system. The aquifer layering, properties, and boundaries relied heavily on a previous three-dimensional variable-density solute-transport model of the same area developed by the U.S. Geological Survey. The surface-water system within these new models actively simulates a part of the extensive canal network by using level-pool routing and active structure operations within the Surface-Water Routing process. These models were used to simulate a historical base-case period (1990–99) by using measured data and regional climate model rainfall and potential evapotranspiration output. The simulated flow and water-level results generally captured the behavior of the hydrologic system. A future period (2060–69) was simulated by using regional climate model rainfall and potential evapotranspiration output representing a wetter and drier future and low, intermediate, and high sea-level rise projections. The results were used to evaluate the potential effects on the surface-water drainage system, coastal-structure operation, and wet-season groundwater levels. Future period simulations using the county-scale model indicate that (1) the effects of the changing climate and sea level are much more evident in eastern and coastal areas of Broward County compared to western areas, with increases in groundwater level nearly equivalent to sea-level rise; (2) coastal groundwater-level increases are distributed farther inland in the wetter future scenarios than in the drier future scenarios; (3) water levels at the westernmost groundwater station locations exhibited little change caused by sea-level rise and showed more dependence on changes in precipitation; (4) there was a reduced west-to-east groundwater gradient with increasing sea-level rise; and (5) increased downstream tidal stage at the S–13 structure resulted in increased reliance on the pump to control upstream inland canal stages. Future simulations using the local-scale model indicate similar behavior as seen in the county-scale model: (1) the coastal areas exhibited the largest impacts in groundwater levels in the future scenarios; (2) the westernmost, interior areas exhibited little change during the future scenarios; and (3) there was an increased reliance on the pump at the S–13 coastal structure but to a lesser extent than indicated in the county-scale model because of the reduced temporal scale of the local-scale model. Possible adaptation and mitigation strategies were simulated to evaluate the county-scale and local-scale models’ ability to simulate hydrologic changes. Alterations to S–13 pump operations within the county-scale model were tested, and results indicate a reduced effect of sea-level rise inland of the control structure, but the affected area is spatially limited. The concept of using pumps to reduce the local groundwater levels in two neighborhood-sized areas was tested by using the local-scale model. The MODFLOW 2005 Drain package was used to remove groundwater by using drainage elevations set to zero, 1 foot, and 2 feet above average wet-season groundwater levels. Area 1 was well connected to coastal boundaries, and a high rate of groundwater removal was required, whereas the rate of groundwater removal required was greatly reduced in Area 2, which is less connected to tidal boundaries. Water for these scenarios was assumed to be pumped to tide with no downstream effects.

Florida

Chemical Characteristics, Water Sources and Pathways, and Age Distribution of Ground Water in the Contributing Recharge Area of a Public-Supply Well near Tampa, Florida, 2002-05

In 2001, the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey began a series of studies on the transport of anthropogenic and natural contaminants (TANC) to public-supply wells. The main goal of the TANC program was to better understand the source, transport, and receptor factors that control contaminant movement to public-supply wells in representative aquifers of the United States. Studies were first conducted at regional scales at four of the eight TANC study areas during 2002-03 and at small (local) scales during 2003-05 in California, Nebraska, Connecticut, and Florida. In the Temple Terrace study area near Tampa, Florida, multiple chemical indicators and geochemical and ground-water flow modeling techniques were used to assess the vulnerability of a public-supply well in the karstic Upper Floridan aquifer to contamination from anthropogenic and naturally occurring contaminants. During 2003-05, water samples were collected from the public-supply well and 13 surrounding monitoring wells that all tap the Upper Floridan aquifer, and from 15 monitoring wells in the overlying surficial aquifer system and the intermediate confining unit that are located within the modeled ground-water contributing recharge area of the public-supply well. Six volatile organic compounds and four pesticides were detected in trace concentrations (well below drinking-water standards) in water from the public-supply well, which had an open interval from 36 to 53 meters below land surface. These contaminants were detected more frequently in water samples from monitoring wells in the overlying clastic surficial aquifer system than in water from monitoring wells in the Upper Floridan aquifer in the study area. Likewise, nitrate-N concentrations in the public-supply well (0.72-1.4 milligrams per liter) were more similar to median concentrations in the oxic surficial aquifer system (2.1 milligrams per liter) than to median nitrate-N concentrations in the anoxic Upper Floridan aquifer (0.06 milligram per liter) under sulfate-reducing conditions. High concentrations of radon-222 and uranium in the public-supply well compared to those in monitoring wells in the Upper Floridan aquifer appear to originate from water moving downward through sands and discontinuous clay lenses that overlie the aquifer. Water samples also were collected from three overlapping depth intervals (38-53, 43-53, and 49-53 meters below land surface) in the public-supply well. The 49- to 53-meter interval was identified as a high-flow zone during geophysical logging of the wellbore. Water samples were collected from these depth intervals at a low pumping rate by placing a low-capacity submersible pump (less than 0.02 cubic meter per minute) at the top of each interval. To represent higher pumping conditions, a large-capacity portable submersible pump (1.6 cubic meters per minute) was placed near the top of the open interval; water-chemistry samples were collected using the low-capacity submersible pump. The 49- to 53-meter depth interval had distinctly different chemistry than the other two sampled intervals. Higher concentrations of nitrate-N, atrazine, radon, trichloromethane (chloroform), and arsenic (and high arsenic (V)/arsenic (III) ratios); lower concentrations of dissolved solids, strontium, iron, manganese, and lower nitrogen and sulfur isotope ratios were found in this highly transmissive zone in the limestone than in water from the two other depth intervals. Movement of water likely occurs from the overlying sands and clays of the oxic surficial aquifer system and intermediate confining unit (that contains high radon-222 and nitrate-N concentrations) into the anoxic Upper Floridan aquifer (that contains low radon-222 and nitrate-N concentrations). Differences in arsenic concentrations in water from the various depth intervals in the public-supply well (3.2-19.0 micrograms per liter) were related to pumping conditions. The high arsenic

Scientific Investigations Report

Relations between well-field pumping and induced canal leakage in east-central Miami-Dade County, Florida, 2010-2011

An extensive canal and water management system exists in south Florida to prevent flooding, replenish groundwater, and impede saltwater intrusion. The unconfined Biscayne aquifer, which underlies southeast Florida and provides water for millions of residents, interacts with the canal system. The Biscayne aquifer is composed of a highly transmissive karst limestone; therefore, canal stage and flow may be affected by production well pumping, especially in locations where production wells and canals are in proximity. The U.S. Geological Survey developed a local-scale, transient, numerical groundwater flow model of a well field in Miami-Dade County to (1) quantify relations between well-field pumping and C-2 Canal (herein referred to as the Snapper Creek Canal) leakage, (2) determine primary controls on canal leakage variability, and (3) summarize results that could simplify characterization of canal/well-field interactions in other locations. In addition to the groundwater model development, stable isotope data from water-quality samples were used to characterize the relations between production well pumping and canal leakage. The results from the groundwater model and the isotope data were used to refine the conceptual flow model of the study area. The groundwater flow model MODFLOW-NWT was used for simulating groundwater flow and quantifying interactions between pumping from the well field and Snapper Creek Canal leakage. Input data for the groundwater model included precipitation, evapotranspiration, pumping, canal stage, and regional groundwater elevation. The inverse modeling tool UCODE and groundwater data from June 2010 to July 2011 were used to calibrate the model. Parameter sensitivity analyses were performed with UCODE. Model sensitivities to geologic heterogeneity, non-laminar flow, and changes in the regional flow boundary were evaluated. The groundwater model generally fits the calibration criteria well within estimated error ranges for groundwater elevations and canal leakage values. The mean average error for heads simulated with the model was 0.19 meter, and head residuals were generally randomly distributed. The model simulated groundwater flow under ambient conditions without production well pumping to establish background leakage. Groundwater flow also was simulated with production well pumping to estimate induced leakage from the Snapper Creek Canal that occurs in response to pumping. Canal leakage was quantified as a percentage of total canal leakage. The percentage of leakage during pumping increased non-linearly with pumping rate, indicating a decreasing sensitivity of canal leakage to pumping at relatively large pumping magnitudes. The results for Snapper Creek Canal may serve as an upper limit for well-field interaction with surface-water features in Miami-Dade County, given the proximity (about 50 meters) of the pumping wells in this study to the Snapper Creek Canal. The isotopic compositions of hydrogen (H) and oxygen (O) in groundwater samples were used to distinguish sources for groundwater within the study area and to assess the extent of natural mixing and pumping-induced mixing with water in the Snapper Creek Canal. Water-level data and water-quality samples were collected from monitoring well clusters, production wells, and the Snapper Creek Canal during discrete sampling events under ambient and pumping conditions. Trends in the isotope data generally follow the regional west-to-east hydraulic gradient across the study area. Data collected within the monitoring-well clusters in closest proximity to the canal indicate that groundwater/surface-water interactions are greatest within the shallow flow zone of the aquifer, especially during pumping conditions. The isotopic composition of samples collected within the study area indicates that the shallow, highly transmissive preferential flow zone receives substantial recharge from the canal. The isotope data from the production wells which are open to the deeper flow zone within the aquifer, indicate only traces of mixing with a 2 H- and 18 O-enriched source, suggesting little canal admixture with waters of the deeper flow zone. Results from the groundwater model and the stable isotope data analysis indicate the importance of considering geologic heterogeneity when investigating the relations between pumping and canal leakage, not only at this site, but also at other sites with similar heterogeneous geology. The model results were consistently sensitive to the hydrogeologic framework and changes in hydraulic conductivities. The model and the isotope data indicate that the majority of the groundwater/surface-water interactions occurred within the shallow flow zone. A relatively lower-permeability geologic layer occurring between the shallowest and deep preferential flow zones lessens the interactions between the production wells and the canal.

Florida

Assessment of the hydraulic connection between ground water and the Peace River, west-central Florida

The hydraulic connection between the Peace River and the underlying aquifers along the length of the Peace River from Bartow to Arcadia was assessed to evaluate flow exchanges between these hydrologic systems. Methods included an evaluation of hydrologic and geologic records and seismic-reflection profiles, seepage investigations, and thermal infrared imagery interpretation. Along the upper Peace River, a progressive long-term decline in streamflow has occurred since 1931 due to a lowering of the potentiometric surface of the Upper Floridan aquifer by as much as 60 feet because of intensive ground-water withdrawals for phosphate mining and agriculture. Another effect from lowering the potentiometric surface has been the cessation of flow at several springs located near and within the Peace River channel, including Kissengen Spring, that once averaged a flow of about 19 million gallons a day. The lowering of ground-water head resulted in flow reversals at locations where streamflow enters sinkholes along the streambed and floodplain. Hydrogeologic conditions along the Peace River vary from Bartow to Arcadia. Three distinctive hydrogeologic areas along the Peace River were delineated: (1) the upper Peace River near Bartow, where ground-water recharge occurs; (2) the middle Peace River near Bowling Green, where reversals of hydraulic gradients occur; and (3) the lower Peace River near Arcadia, where ground-water discharge occurs. Seismic-reflection data were used to identify geologic features that could serve as potential conduits for surface-water and ground-water exchange. Depending on the hydrologic regime, this exchange could be recharge of surface water into the aquifer system or discharge of ground water into the stream channel. Geologic features that would provide pathways for water movement were identified in the seismic record; they varied from buried irregular surfaces to large-scale subsidence flexures and vertical fractures or enlarged solution conduits. Generally, the upper Peace River is characterized by a shallow, buried irregular top of rock, numerous observed sinkholes, and subsidence depressions. The downward head gradient provides potential for the Peace River to lose water to the ground-water system. Along the middle Peace River area, head gradients alternate between downward and upward, creating both recharging and discharging ground-water conditions. Seismic records show that buried, laterally continuous reflectors in the lower Peace River pinch out in the middle Peace River streambed. Small springs have been observed along the streambed where these units pinch out. This area corresponds to the region where highest ground-water seepage volumes were measured during this study. Further south, along the lower Peace River, upward head gradients provide conditions for ground-water discharge into the Peace River. Generally, confinement between the surficial aquifer and the confined ground-water systems in this area is better than to the north. However, localized avenues for surface-water and ground-water interactions may exist along discontinuities observed in seismic reflectors associated with large-scale flexures or subsidence features. Ground-water seepage gains or losses along the Peace River were quantified by making three seepage runs during periods of: (1) low base flow, (2) high base flow, and (3) high flow. Low and high base-flow seepage runs were performed along a 74-mile length of the Peace River, between Bartow and Nocatee. Maximum losses of 17.3 cubic feet per second (11.2 million gallons per day) were measured along a 3.2-mile reach of the upper Peace River. The high-flow seepage run was conducted to quantify losses in the Peace River channel and floodplain between Bartow and Fort Meade. Seepage losses calculated during high-flow along a 7.2-mile reach of the Peace River, from the Clear Springs Mine bridge to the Mobil Mine bridge, were approximately 10 percent of the river flow, or 118 cubic feet per second. Calculated seepages along the Peace River in Hardee and De Soto Counties were inconclusive, because most seepages were within the range of discharge measurement error. Two continuous aerial thermal infrared imagery surveys were conducted to locate sites of ground-water discharge along the Peace River. Although temperature and hydrologic conditions were ideal to observe spring flow using thermal infrared imaging techniques, no sources of ground-water discharge were identified using this method. Diffuse ground-water seepage may, however, provide significant ground-water discharge.

Florida

Evaluation of effects of changes in canal management and precipitation patterns on salinity in Biscayne Bay, Florida, using an integrated surface-water/groundwater model

Biscayne National Park, located in Biscayne Bay in southeast Florida, is one of the largest marine parks in the country and sustains a large natural marine fishery where numerous threatened and endangered species reproduce. In recent years, the bay has experienced hypersaline conditions (salinity greater than 35 practical salinity units) of increasing magnitude and duration. Hypersalinity events were particularly pronounced during April to August 2004 in nearshore areas along the southern and middle parts of the bay. Prolonged hypersaline conditions can cause degradation of water quality and permanent damage to, or loss of, brackish nursery habitats for multiple species of fish and crustaceans as well as damage to certain types of seagrasses that are not tolerant of extreme changes in salinity. To evaluate the factors that contribute to hypersalinity events and to test the effects of possible changes in precipitation patterns and canal flows into Biscayne Bay on salinity in the bay, the U.S. Geological Survey constructed a coupled surface-water/groundwater numerical flow model. The model is designed to account for freshwater flows into Biscayne Bay through the canal system, leakage of salty bay water into the underlying Biscayne aquifer, discharge of fresh and salty groundwater from the Biscayne aquifer into the bay, direct effects of precipitation on bay salinity, indirect effects of precipitation on recharge to the Biscayne aquifer, direct effects of evapotranspiration (ET) on bay salinity, indirect effects of ET on recharge to the Biscayne aquifer, and maintenance of mass balance of both water and solute. The model was constructed using the Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) simulator, version 3.3, which couples the two-dimensional, surface-water flow and solute-transport simulator SWIFT2D with the density-dependent, groundwater flow an solute-transport simulator SEAWAT. The model was calibrated by a trial-and-error method to fit observed groundwater heads, estimated base flow, and measured bay salinity and temperatures from 1996 to 2004, as well as the location of the freshwater-saltwater interface in the aquifer, by adjusting ET rate limiters, canal vertical hydraulic conductance, leakage rate coefficients (transition-layer thickness and hydraulic conductivity), Manning's n value, and delineation of rainfall zones. Although flow budget calculations indicate that precipitation, ET, and groundwater flux into the bay represent a small portion of the overall budget, these factors may be important in controlling salinity in some parts of the bay, for example the southern parts of the bay where the canal system is not extensively developed or controlled. The balance of precipitation and ET during the wet season generally results in a reduction of bay salinity, whereas the balance of precipitation and ET during the dry season generally results in an increase in bay salinity. During years when wet season precipitation is lower than average, for example less than 70 percent total precipitation for an average year, ET could outweigh precipitation over the bay for essentially the entire year. Hypersaline conditions are prone to occur near the end of the dry season because precipitation rates are generally lower, canal discharge rates (which are strongly correlated to precipitation rates) are also generally lower, and ET rates are higher than during the rest of the year. The hypersalinity event of 2004 followed several years of relatively low precipitation and correspondingly reduced canal structure releases and was unusually extensive, continuing into July. Thus, hypersalinity is ultimately the result of a cumulative deficit of precipitation. The model was used to test the effects of possible changes in canal flux and precipitation. Simulation results showed that by increasing, reducing, or modifying canal discharge rates, the effects on salinity in the bay were more pronounced in the northern part of the bay where there are more canals and canal-control structures. By doubling and halving precipitation, the effects on bay salinity were more pronounced in the southern part of the bay than in the northern part of the bay where there are fewer canals and canal-control structures. The model is designed to quantify factors that contribute to hypersaline conditions in Biscayne Bay and may be less appropriate for addressing other issues or examining conditions substantially different from those described in this report. Model results must be interpreted in light of model limitations, which include representation of the system and conceptual model, uncertainty in physical properties used to describe the system or processes, the scale and discretization of the system, and representation of the boundary conditions.

Florida

Ground-water hydrology and simulation of ground-water flow at Operable Unit 3 and surrounding region, U.S. Naval Air Station, Jacksonville, Florida

The Naval Air Station, Jacksonville (herein referred to as the Station), occupies 3,800 acres adjacent to the St. Johns River in Duval County, Florida. Operable Unit 3 (OU3) occupies 134 acres on the eastern side of the Station and has been used for industrial and commercial purposes since World War II. Ground water contaminated by chlorinated organic compounds has been detected in the surficial aquifer at OU3. The U.S. Navy and U.S. Geological Survey (USGS) conducted a cooperative hydrologic study to evaluate the potential for ground water discharge to the neighboring St. Johns River. A ground-water flow model, previously developed for the area, was recalibrated for use in this study. At the Station, the surficial aquifer is exposed at land surface and forms the uppermost permeable unit. The aquifer ranges in thickness from 30 to 100 feet and consists of unconsolidated silty sands interbedded with local beds of clay. The low-permeability clays of the Hawthorn Group form the base of the aquifer. The USGS previously conducted a ground-water investigation at the Station that included the development and calibration of a 1-layer regional ground-water flow model. For this investigation, the regional model was recalibrated using additional data collected after the original calibration. The recalibrated model was then used to establish the boundaries for a smaller subregional model roughly centered on OU3. Within the subregional model, the surficial aquifer is composed of distinct upper and intermediate layers. The upper layer extends from land surface to a depth of approximately 15 feet below sea level; the intermediate layer extends from the upper layer down to the top of the Hawthorn Group. In the northern and central parts of OU3, the upper and intermediate layers are separated by a low-permeability clay layer. Horizontal hydraulic conductivities in the upper layer, determined from aquifer tests, range from 0.19 to 3.8 feet per day. The horizontal hydraulic conductivity in the intermediate layer, determined from one aquifer test, is 20 feet per day. An extensive stormwater drainage system is present at OU3 and the surrounding area. Some of the stormwater drains have been documented to be draining ground water from the upper layer of the surficial aquifer, whereas other drains are only suspected to be draining ground water. The subregional model contained 78 rows and 148 columns of square model cells that were 100 feet on each side. Vertically, the surficial aquifer was divided into two layers; layer 1 represented the upper layer and layer 2 represented the intermediate layer. Steady-state ground-water flow conditions were assumed. The model was calibrated to head data collected on October 29 and 30, 1996. After calibration, the model matched all 67 measured heads to within the calibration criterion of 1 foot; and 48 of 67 simulated heads (72 percent) were within 0.5 foot. Model simulated recharge rates ranged from 0.4 inch per year in areas that were largely paved to 13.0 inches per year in irrigated areas. Simulated hydraulic conductivities in the upper layer at OU3 ranged from 0.5 foot per day in the north to 1.0 foot per day in the south. Simulated vertical leakance between the upper and intermediate layers ranged from 1.0x10-6 per day in an area with low-permeability clays to 4.3x10-2 per day in an area that had been dredged. Simulated transmissivities in the intermediate layer ranged from 25 feet squared per day in an area of low-permeability channel-fill deposits to a high of 1,200 feet squared per day in areas covering most of OU3. Simulated riverbed conductances ranged from 4 to 60 feet squared per day and simulated bottom conductances of leaking stormwater drains ranged from 5 to 20 feet squared per day. The direction and velocity of ground-water flow was determined using particle-tracking techniques. Ground-water flow in the upper layer was generally eastward toward the St. Johns River. However, leaking stormwater drains locally modified the flow system to create small areas with flow that was diverted to the drains. The flow velocities in the upper layer at OU3 were slow, averaging about 2 feet per year. The slow velocities were primarily the result of the low horizontal hydraulic conductivity and, secondarily, the result of the low recharge rate. The simulated rate at which ground water leaked into the stormwater drains was low, averaging about 0.0011 cubic feet per second per 100 feet of stormwater drainage conduit. Ground-water flow in the intermediate layer moved eastward toward and discharged into the St. Johns River. Flow velocities were significantly higher in this layer than in the upper layer. The velocity was about 35 and 12 feet per year in the northern and southern parts of OU3, respectively.

Florida

Anthropogenic Organic Compounds in Ground Water and Finished Water of Community Water Systems in the Northern Tampa Bay Area, Florida, 2002-04

As part of the U.S. Geological Survey's (USGS's) National Water-Quality Assessment (NAWQA) Program, a Source Water-Quality Assessment (SWQA) was conducted in the unconfined and semiconfined portions of the Upper Floridan aquifer system during 2002-04. SWQAs are two-phased sampling activities, wherein phase 1 was designed to evaluate the occurrence of 258 anthropogenic organic compounds (AOCs) in ground water used as source water for 30 of the largest-producing community water system (CWS) wells in the northern Tampa Bay area, Florida. The 258 AOCs included volatile organic compounds (VOCs), pesticides, and other anthropogenic organic compounds (OAOCs). Phase 2 was designed to monitor concentrations in the source water and also the finished water of CWSs for compounds most frequently detected during phase 1. During phase 1 of the SWQA study, 31 of the 258 AOCs were detected in source-water samples collected from CWS wells at low concentrations (less than 1.0 microgram per liter (ug/L)). Twelve AOCs were detected in at least 10 percent of samples. Concentrations from 16 of the 31 detected AOCs were about 2 to 5 orders of magnitude below human-health benchmarks indicating that concentrations were unlikely to be of potential human-health concern. The potential human-health relevance for the remaining 15 detected unregulated AOCs could not be evaluated because no human-health benchmarks were available for these compounds. Hydrogeology, population, and land use were examined to evaluate the effects of these variables on the source water monitored. Approximately three times as many detections of VOCs (27) and pesticides (34) occurred in unconfined areas than in the semiconfined areas (8 VOCs, 14 pesticides). In contrast, 1 OAOC was detected in unconfined areas, and 13 OAOCs were detected in semiconfined areas with 9 of the OAOC detections occurring in samples from two wells located near septic systems. Analyses of population and land use indicated that the number of compounds detected increased as the population surrounding each well increased. Detection frequencies and concentrations for VOCs (particularly chloroform) and pesticides were highest in residential land-use areas. The results of source-water samples from the 30 CWS wells monitored during phase 1 of this SWQA study were compared to four locally conducted studies. These general comparisons indicate that the occurrence of VOCs in other studies is similar to their occurrence in source water of CWSs monitored as part of this SWQA. However, pesticide compounds, especially atrazine and its breakdown products, occurred more frequently in the SWQA study than in the other four studies. Phase 2 of the SWQA assessed AOCs in samples from 11 of the 30 CWS wells and the associated finished water. Overall, 42 AOCs were detected in either source water or finished water and more compounds were detected in finished water than in source water. Specifically, 22 individual AOCs were detected in source water and 27 AOCs were detected in finished water. The total number of detections was greater in the finished water (80) than in the source water (49); however, this was largely due to the creation of disinfection by-products (DBPs) during water treatment. Excluding DBPs, about the same number of total detections was observed in source water (40) and finished water (44). During phase 2, AOC detected concentrations ranged from E0.003 (estimated) to 1,140 ug/L in the source water and from E0.003 to 36.3 ug/L in the finished water. Concentrations of 24 of the 42 compounds were compared to human-health benchmarks and were about 1 to 5 orders of magnitude below their human-health benchmarks indicating that concentrations are unlikely to be of potential human-health concern, excluding DBPs. Concentrations of carbon tetrachloride, however, were within 10 percent of its human-health benchmark, which is considered a level that may warrant inclusion of the compound in a low-concentration, t

Scientific Investigations Report

Relations between precipitation, groundwater withdrawals, and changes in hydrologic conditions at selected monitoring sites in Volusia County, Florida, 1995--2010

A study to examine the influences of climatic and anthropogenic stressors on groundwater levels, lake stages, and surface-water discharge at selected sites in northern Volusia County, Florida, was conducted in 2009 by the U.S. Geological Survey. Water-level data collected at 20 monitoring sites (17 groundwater and 3 lake sites) in the vicinity of a wetland area were analyzed with multiple linear regression to examine the relative influences of precipitation and groundwater withdrawals on changes in groundwater levels and lake stage. Analyses were conducted across varying periods of record between 1995 and 2010 and included the effects of groundwater withdrawals aggregated from municipal water-supply wells located within 12 miles of the project sites. Surface-water discharge data at the U.S. Geological Survey Tiger Bay canal site were analyzed for changes in flow between 1978 and 2001. As expected, water-level changes in monitoring wells located closer to areas of concentrated groundwater withdrawals were more highly correlated with withdrawals than were water-level changes measured in wells further removed from municipal well fields. Similarly, water-level changes in wells tapping the Upper Floridan aquifer, the source of municipal supply, were more highly correlated with groundwater withdrawals than were water-level changes in wells tapping the shallower surficial aquifer system. Water-level changes predicted by the regression models over precipitation-averaged periods of record were underestimated for observations having large positive monthly changes (generally greater than 1.0 foot). Such observations are associated with high precipitation and were identified as points in the regression analyses that produced large standardized residuals and/or observations of high influence. Thus, regression models produced by multiple linear regression analyses may have better predictive capability in wetland environments when applied to periods of average or below average precipitation conditions than during wetter than average conditions. For precipitation-averaged hydrologic conditions, water-level changes in the surficial aquifer system were statistically correlated solely with precipitation or were more highly correlated with precipitation than with groundwater withdrawals. Changes in Upper Floridan aquifer water levels and in water-surface stage (stage) at Indian and Scoggin Lakes tended to be highly correlated with both precipitation and withdrawals. The greater influence of withdrawals on stage changes, relative to changes in nearby surficial aquifer system water levels, indicates that these karstic lakes may be better connected hydraulically with the underlying Upper Floridan aquifer than is the surficial aquifer system at the other monitoring sites. At most sites, and for both aquifers, the 2-month moving average of precipitation or groundwater withdrawals included as an explanatory variable in the regression models indicates that water-level changes are not only influenced by stressor conditions across the current month, but also by those of the previous month. The relations between changes in water levels, precipitation, and groundwater withdrawals varied seasonally and in response to a period of drought. Water-level changes tended to be most highly correlated with withdrawals during the spring, when relatively large increases contributed to water-level declines, and during the fall when reduced withdrawal rates contributed to water-level recovery. Water-level changes tended to be most highly (or solely) correlated with precipitation in the winter, when withdrawals are minimal, and in the summer when precipitation is greatest. Water-level changes measured during the drought of October 2005 to June 2008 tended to be more highly correlated with groundwater withdrawals at Upper Floridan aquifer sites than at surficial aquifer system sites, results that were similar to those for precipitation-averaged conditions. Also, changes in stage at Indian and Scoggin Lakes were highly correlated with precipitation and groundwater withdrawals during the drought. Groundwater-withdrawal rates during the drought were, on average, greater than those for precipitation-averaged conditions. Accounting only for withdrawals aggregated from pumping wells located within varying radial distances of less than 12 miles of each site produced essentially the same relation between water-level changes and groundwater withdrawals as that determined for withdrawals aggregated within 12 miles of the site. Similarly, increases in withdrawals aggregated over distances of 1 to 12 miles of the sites had little effect on adjusted R-squared values. Analyses of streamflow measurements collected between 1978 and 2001 at the U.S. Geological Survey Tiger Bay canal site indicate that significant changes occurred during base-flow conditions during that period. Hypothesis and trend testing, together with analyses of flow duration, the number of zero-flow days, and double-mass curves indicate that, after 1988, when a municipal well field began production, base flow was statistically lower than the period before 1988. This decrease in base flow could not be explained by variations in precipitation between these two periods.

Florida

Estimating nitrate concentrations in groundwater at selected wells and springs in the surficial aquifer system and Upper Floridan aquifer, Dougherty Plain and Marianna Lowlands, Georgia, Florida, and Alabama, 2002-50

Groundwater from the surficial aquifer system and Upper Floridan aquifer in the Dougherty Plain and Marianna Lowlands in southwestern Georgia, northwestern Florida, and southeastern Alabama is affected by elevated nitrate concentrations as a result of the vulnerability of the aquifer, irrigation water-supply development, and intensive agricultural land use. The region relies primarily on groundwater from the Upper Floridan aquifer for drinking-water and irrigation supply. Elevated nitrate concentrations in drinking water are a concern because infants under 6 months of age who drink water containing nitrate concentrations above the U.S. Environmental Protection Agency maximum contaminant level of 10 milligrams per liter as nitrogen can become seriously ill with blue baby syndrome. In response to concerns about water quality in domestic wells and in springs in the lower Apalachicola–Chattahoochee–Flint River Basin, the Florida Department of Environmental Protection funded a study in cooperation with the U.S. Geological Survey to examine water quality in groundwater and springs that provide base flow to the Chipola River. A three-dimensional, steady-state, regional-scale groundwater-flow model and two local-scale models were used in conjunction with particle tracking to identify travel times and areas contributing recharge to six groundwater sites—three long-term monitor wells (CP-18A, CP-21A, and RF-41) and three springs (Jackson Blue Spring, Baltzell Springs Group, and Sandbag Spring) in the lower Apalachicola–Chattahoochee–Flint River Basin. Estimated nitrate input to groundwater at land surface, based on previous studies of nitrogen fertilizer sales and atmospheric nitrate deposition data, were used in the advective transport models for the period 2002 to 2050. Nitrate concentrations in groundwater samples collected from the six sites during 1993 to 2007 and groundwater age tracer data were used to calibrate the transport aspect of the simulations. Measured nitrate concentrations (as nitrogen) in wells and springs sampled during the study ranged from 0.37 to 12.73 milligrams per liter. Average apparent ages of groundwater calculated from measurements of chlorofluorocarbon, sulfur hexafluoride, and tritium from wells CP-18A, CP-21A,and RF-41 were about 23, 29, and 32 years, respectively. Average apparent ages of groundwater from Baltzell Springs Group, Sandbag Spring, and Jackson Blue Spring were about 16, 18, and 19 years, respectively. Simulated travel times of particles from the six selected sites ranged from less than 1 day to 511 years; both the minimum and maximum particle travel times were estimated for water from Jackson Blue Spring. Median simulated travel times of particles were about 30, 38, and 62 years for Jackson Blue Spring, Sandbag Spring, and Baltzell Springs Group, respectively. Study results indicated that travel times for approximately 50 percent of the particles from all spring sites were less than 50 years. The median simulated travel times of particles arriving at receptor wells CP-18A, CP-21A, and RF-41 were about 50, 35, and 36 years, respectively. All particle travel times were within the same order of magnitude as the tracer-derived average apparent ages for water, although slightly older than the measured ages. Travel time estimates were substantially greater than the measured age for groundwater reaching well CP-18A, as confirmed by the average apparent age of water determined from tracers. Local-scale particle-tracking models were used to predict nitrate concentrations in the three monitor wells and three springs from 2002 to 2050 for three nitrogen management scenarios: (1) fixed input of nitrate at the 2001 level, (2) reduction of nitrate inputs of 4 percent per year (from the previous year) from 2002 to 2050, and (3) elimination of nitrate input after 2001. Simulated nitrate concentrations in well CP-21A peaked at 7.82 milligrams per liter in 2030, and concentrations in background well RF-41 peaked at 1.10 milligrams per liter in 2020. The simulated particle travel times were longer than indicated by age dating analysis for groundwater in well CP-18A; to account for the poor calibration fit at this well, nitrate concentrations were shifted 21 years. With the shift, simulated nitrate concentrations in groundwater at CP-18A peaked at 13.76 milligrams per liter in 2026. For groundwater in Baltzell Springs Group, Jackson Blue Spring, and Sandbag Spring, simulated nitrate concentrations peaked at 3.77 milligrams per liter in 2006, 3.51 milligrams per liter in 2011, and 0.81 milligram per liter in 2018, respectively, under the three management scenarios. In management scenario 3 (elimination of nitrate input after 2001), simulated nitrate concentrations in Baltzell Springs Group decreased to less than background concentrations (0.10 milligram per liter) by 2033, and in Sandbag Spring concentrations decreased to less than background by 2041. Simulations using nitrate management scenarios 1 (fixed input of nitrate at 2001 levels) and 2 (reduction of 4.0 percent per year) indicate that nitrate concentrations in groundwater may remain above background concentrations through 2050 at all sites.

Alabama, Florida, Georgia