Geology topics

W. Barclay Shoemaker

Publications and source records attributed to W. Barclay Shoemaker.

13 recordsLinked to original sources

USGS2024

Evaporation from the interior of Lake Okeechobee—A large freshwater lake in Florida, 2013–16

In 2012, a platform at the approximate center of Lake Okeechobee in central Florida was instrumented to continuously measure evaporation with the Bowen-ratio energy-budget method as part of a long-term partnership between the South Florida Water Management District and the U.S. Geological Survey. The primary goal for the study was to quantify daily rates of open-water evaporation. A secondary goal was to assess differences in evaporation rates among alternate methods and determine if instrumentation and operational expenses associated with the Bowen-ratio method could be reduced. Mean annual evaporation from Lake Okeechobee for 2013–16 was about 1,825 millimeters per year. Annual evaporation from 2013 to 2016 was 1,760, 1,840, 1,810, and 1,890 millimeters per year, respectively. These evaporation rates are among the highest rates observed in Florida based on scientifically vetted methods such as evaporation pans, lysimeters, eddy-covariance, or Bowen-ratio methods. The high evaporation rates are largely a result of frequent clear-sky conditions over the interior of Lake Okeechobee, which allows solar radiation to reach the water surface and drive open-water evaporation. Cloud formation over the interior of Lake Okeechobee is suppressed because of a relatively large heat capacity for water that buffers convective fluxes of air that form clouds while rising and cooling. Estimated evaporation rates obtained using five alternative methods were compared to measured Bowen-ratio energy-budget daily, monthly, and annual evaporation: the Penman, Priestly-Taylor, Mass-Transfer, Simple, and Turc equations. All five methods performed relatively well (within 10 percent of the Bowen ratio annual totals). The Penman, Priestley-Taylor, and Mass-Transfer methods captured relatively large evaporation rates that occurred in the winter due to cold fronts, because these methods account for large wind speeds and vapor pressure deficits associated with the regional cold fronts. For operational implementation, the Simple, Mass-Transfer, or Turc methods are likely preferable because of their simplicity, limited data requirements, and improved accuracy for computing monthly and annual evaporation totals. The Turc equation computed monthly evaporation within 8 percent of the Bowen-ratio method, while requiring only air temperature and solar radiation data. The Simple equation achieved similar accuracy while requiring only solar radiation data.

Florida↗

USGS2021

Carbon fluxes and potential soil accumulation within Greater Everglades cypress and pine forested wetlands

In forested wetlands, accumulation of organic matter in soil is partly governed by carbon fluxes where photosynthesis, respiration, lateral advection of waterborne carbon, fire-derived carbon emissions, and methanogenesis are balanced by changes in stored carbon. Stored carbon can eventually accumulate as soil over time if net primary productivity exceeds biomass decomposition. For this study, potential soil accumulation was estimated based on four years of continuous daily carbon cycling data and a one-dimensional mass-balance model of landscape-atmospheric exchange for cypress and pine forested wetlands in the Greater Everglades of south Florida. The mass-balance model was driven by eddy-covariance estimates of vertical net ecosystem exchange of carbon dioxide and methane. Key findings include confirmation of a basic premise of the historic Everglades restoration project; specifically, more water either from rainfall or water management encourages soil carbon accumulation and thus conservation of soils that support biologic activity and ecosystem services. For example, an anomalous wet season for south Florida that flooded the forested wetlands through the traditional dry season was followed by the most productive year for photosynthetic carbon uptake and potential soil accumulation. On the other hand, methane emissions were enhanced by the anomalous wet season and extended flooding – which confirmed a complex tradeoff to consider if wetlands are managed for both soil conservation and reduction of greenhouse gas emissions. Potential soil accumulation rates were about 1.7, 2.8, and 18 millimeters per year at the Dwarf Cypress, Cypress Swamp, and Pine Upland ecosystems, assuming soil C density values of 0.07, 0.09, and 0.02 grams of carbon per cubic centimeter, respectively. For these values of soil C density, the accumulation rates are considered a “best-case” upper limit because the lateral export of carbon in the canals and creeks that drain the study area were assumed negligible.

Florida↗

USGS2015

Carbon exchange between the atmosphere and subtropical forested cypress and pine wetlands

Carbon dioxide exchange between the atmosphere and forested subtropical wetlands is largely unknown. Here we report a first step in characterizing this atmospheric–ecosystem carbon (C) exchange, for cypress strands and pine forests in the Greater Everglades of Florida as measured with eddy covariance methods at three locations (Cypress Swamp, Dwarf Cypress and Pine Upland) for 2 years. Links between water and C cycles are also examined at these three sites, as are methane emission measured only at the Dwarf Cypress site. Each forested wetland showed net C uptake from the atmosphere both monthly and annually, as indicated by the net ecosystem exchange (NEE) of carbon dioxide (CO 2 ). For this study, NEE is the difference between photosynthesis and respiration, with negative values representing uptake from the atmosphere that is retained in the ecosystem or transported laterally via overland flow (unmeasured for this study). Atmospheric C uptake (NEE) was greatest at the Cypress Swampp (−900 to −1000 g C m 2 yr −1 ), moderate at the Pine Upland (−650 to −700 g C m 2 yr −1 ) and least at the Dwarf Cypress (−400 to −450 g C m 2 yr −1 ). Changes in NEE were clearly a function of seasonality in solar insolation, air temperature and flooding, which suppressed heterotrophic soil respiration. We also note that changes in the satellite-derived enhanced vegetation index (EVI) served as a useful surrogate for changes in NEE at these forested wetland sites.

Florida↗

USGS2011

Modifications to the conduit flow process mode 2 for MODFLOW-2005

As a result of rock dissolution processes, karst aquifers exhibit highly conductive features such as caves and conduits. Within these structures, groundwater flow can become turbulent and therefore be described by nonlinear gradient functions. Some numerical groundwater flow models explicitly account for pipe hydraulics by coupling the continuum model with a pipe network that represents the conduit system. In contrast, the Conduit Flow Process Mode 2 (CFPM2) for MODFLOW‐2005 approximates turbulent flow by reducing the hydraulic conductivity within the existing linear head gradient of the MODFLOW continuum model. This approach reduces the practical as well as numerical efforts for simulating turbulence. The original formulation was for large pore aquifers where the onset of turbulence is at low Reynolds numbers (1 to 100) and not for conduits or pipes. In addition, the existing code requires multiple time steps for convergence due to iterative adjustment of the hydraulic conductivity. Modifications to the existing CFPM2 were made by implementing a generalized power function with a user‐defined exponent. This allows for matching turbulence in porous media or pipes and eliminates the time steps required for iterative adjustment of hydraulic conductivity. The modified CFPM2 successfully replicated simple benchmark test problems.

Ground Water↗

USGS2011

Evapotranspiration over spatially extensive plant communities in the Big Cypress National Preserve, southern Florida, 2007-2010

Evapotranspiration (ET) was quantified over plant communities within the Big Cypress National Preserve (BCNP) using the eddy covariance method for a period of 3 years from October 2007 to September 2010. Plant communities selected for study included Pine Upland, Wet Prairie, Marsh, Cypress Swamp, and Dwarf Cypress. These plant communities are spatially extensive in southern Florida, and thus, the ET measurements described herein can be applied to other humid subtropical locations such as the Everglades. The 3-year mean annual ET was about 1,000, 1,050, 1,100, 930, and 900 mm (millimeters) at the Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland, and Marsh sites, respectively. Spatial differences in annual ET were considerable due to the recovery of the Marsh site from extensive forest fire and drought conditions. Temporal variability in annual ET was relatively small at sites that were well watered (Dwarf Cypress, Wet Prairie, Cypress Swamp, Pine Upland) over the 3-year study. In other words, locations that were well watered appeared to have similar annual ET rates. Diurnal water-level variability was observed in response to ET and was less at flooded sites than at dry sites. For example, surface-water levels declined about 1.5 mm in response to ET at the flooded Cypress Swamp site during July 18-22, 2008 and declined about 10 mm in response to ET at the flooded Dwarf Cypress site from April 18-27, 2008. Specific yield was computed using ET estimates and diurnal water-level variability measured at the dry Pine Upland site as a check on the accuracy of the eddy covariance method. Water levels repeatedly dropped about 15 mm on average in response to ET at the Pine Upland site from April 27 to May 4, 2008. ET was about 3 mm on each of these days, resulting in a reasonable estimate for specific yield of 0.2 for the sandy soils at the Pine Upland site. Monthly ET estimates exhibited seasonal variation. ET was generally greatest between March to October when solar radiation was relatively large, and least from November to February when solar radiation was small. Monthly ET was greatest in the spring and summer at the Cypress Swamp site, reaching rates as large as 140 mm. The large ET rates at this site coincide with the most active period of cypress growth during late spring and early summer. Cypress trees begin to senesce in late summer reducing transpiration. Net radiation and available energy explained most of the variability in ET observed at all five sites. Mean annual and monthly net radiation varied among the sites in response to cloud cover and the albedo of the land surface and plant community. Net radiation was greatest at the Cypress Swamp site, averaging about 130 W/m 2 (watts per square meter) during the 3-year study. Net radiation was generally less at the Dwarf Cypress site, averaging about 115 W/m 2 over 3 years. The Dwarf Cypress site apparently has the largest albedo, which likely is due to the sparse canopy and a highly reflective, calcareous, periphyton-covered land surface. Furthermore, mean annual net radiation was least in the first year of the study, which likely was due to greater cloud cover during a relatively wet year. In contrast, net radiation was greatest in the second year of the study, which likely was due to less cloud cover during a relatively dry year.

Scientific Investigations Report↗

USGS2008

Permeameter data verify new turbulence process for MODFLOW

A sample of Key Largo Limestone from southern Florida exhibited turbulent flow behavior along three orthogonal axes as reported in recently published permeameter experiments. The limestone sample was a cube measuring 0.2 m on edge. The published nonlinear relation between hydraulic gradient and discharge was simulated using the turbulent flow approximation applied in the Conduit Flow Process (CFP) for MODFLOW-2005 mode 2, CFPM2. The good agreement between the experimental data and the simulated results verifies the utility of the approach used to simulate the effects of turbulent flow on head distributions and flux in the CFPM2 module of MODFLOW-2005.

Ground Water↗

USGS2008

Effects of turbulence on hydraulic heads and parameter sensitivities in preferential groundwater flow layers

A conduit flow process (CFP) for the Modular Finite Difference Ground‐Water Flow model, MODFLOW‐2005, has been created by the U.S. Geological Survey. An application of the CFP on a carbonate aquifer in southern Florida is described; this application examines (1) the potential for turbulent groundwater flow and (2) the effects of turbulent flow on hydraulic heads and parameter sensitivities. Turbulent flow components were spatially extensive in preferential groundwater flow layers, with horizontal hydraulic conductivities of about 5,000,000 m d −1 , mean void diameters equal to about 3.5 cm, groundwater temperature equal to about 25°C, and critical Reynolds numbers less than or equal to 400. Turbulence either increased or decreased simulated heads from their laminar elevations. Specifically, head differences from laminar elevations ranged from about −18 to +27 cm and were explained by the magnitude of net flow to the finite difference model cell. Turbulence also affected the sensitivities of model parameters. Specifically, the composite‐scaled sensitivities of horizontal hydraulic conductivities decreased by as much as 70% when turbulence was essentially removed. These hydraulic head and sensitivity differences due to turbulent groundwater flow highlight potential errors in models based on the equivalent porous media assumption, which assumes laminar flow in uniformly distributed void spaces.

Water Resources Research↗

USGS2007

Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005

This report documents the Conduit Flow Process (CFP) for the modular finite-difference ground-water flow model, MODFLOW-2005. The CFP has the ability to simulate turbulent ground-water flow conditions by: (1) coupling the traditional ground-water flow equation with formulations for a discrete network of cylindrical pipes (Mode 1), (2) inserting a high-conductivity flow layer that can switch between laminar and turbulent flow (Mode 2), or (3) simultaneously coupling a discrete pipe network while inserting a high-conductivity flow layer that can switch between laminar and turbulent flow (Mode 3). Conduit flow pipes (Mode 1) may represent dissolution or biological burrowing features in carbonate aquifers, voids in fractured rock, and (or) lava tubes in basaltic aquifers and can be fully or partially saturated under laminar or turbulent flow conditions. Preferential flow layers (Mode 2) may represent: (1) a porous media where turbulent flow is suspected to occur under the observed hydraulic gradients; (2) a single secondary porosity subsurface feature, such as a well-defined laterally extensive underground cave; or (3) a horizontal preferential flow layer consisting of many interconnected voids. In this second case, the input data are effective parameters, such as a very high hydraulic conductivity, representing multiple features. Data preparation is more complex for CFP Mode 1 (CFPM1) than for CFP Mode 2 (CFPM2). Specifically for CFPM1, conduit pipe locations, lengths, diameters, tortuosity, internal roughness, critical Reynolds numbers (NRe), and exchange conductances are required. CFPM1, however, solves the pipe network equations in a matrix that is independent of the porous media equation matrix, which may mitigate numerical instability associated with solution of dual flow components within the same matrix. CFPM2 requires less hydraulic information and knowledge about the specific location and hydraulic properties of conduits, and turbulent flow is approximated by modifying horizontal conductances assembled by the Block-Centered Flow (BCF), Layer-Property Flow (LPF), or Hydrogeologic-Unit Flow Packages (HUF) of MODFLOW-2005. For both conduit flow pipes (CFPM1) and preferential flow layers (CFPM2), critical Reynolds numbers are used to determine if flow is laminar or turbulent. Due to conservation of momentum, flow in a laminar state tends to remain laminar and flow in a turbulent state tends to remain turbulent. This delayed transition between laminar and turbulent flow is introduced in the CFP, which provides an additional benefit of facilitating convergence of the computer algorithm during iterations of transient simulations. Specifically, the user can specify a higher critical Reynolds number to determine when laminar flow within a pipe converts to turbulent flow, and a lower critical Reynolds number for determining when a pipe with turbulent flow switches to laminar flow. With CFPM1, the Hagen-Poiseuille equation is used for laminar flow conditions and the Darcy-Weisbach equation is applied to turbulent flow conditions. With CFPM2, turbulent flow is approximated by reducing the laminar hydraulic conductivity by a nonlinear function of the Reynolds number, once the critical head difference is exceeded. This adjustment approximates the reductions in mean velocity under turbulent ground-water flow conditions.

Techniques and Methods↗

USGS2006

Alternate corrections for estimating actual wetland evapotranspiration from potential evapotranspiration

Corrections can be used to estimate actual wetland evapotranspiration (AET) from potential evapotranspiration (PET) as a means to define the hydrology of wetland areas. Many alternate parameterizations for correction coefficients for three PET equations are presented, covering a wide range of possible data-availability scenarios. At nine sites in the wetland Everglades of south Florida, USA, the relatively complex PET Penman equation was corrected to daily total AET with smaller standard errors than the PET simple and Priestley-Taylor equations. The simpler equations, however, required less data (and thus less funding for instrumentation), with the possibility of being corrected to AET with slightly larger, comparable, or even smaller standard errors. Air temperature generally corrected PET simple most effectively to wetland AET, while wetland stage and humidity generally corrected PET Priestley-Taylor and Penman most effectively to wetland AET. Stage was identified for PET Priestley-Taylor and Penman as the data type with the most correction ability at sites that are dry part of each year or dry part of some years. Finally, although surface water generally was readily available at each monitoring site, AET was not occurring at potential rates, as conceptually expected under well-watered conditions. Apparently, factors other than water availability, such as atmospheric and stomata resistances to vapor transport, also were limiting the PET rate. ?? 2006, The Society of Wetland Scientists.

Wetlands↗

USGS2005

Estimating changes in heat energy stored within a column of wetland surface water and factors controlling their importance in the surface energy budget

Changes in heat energy stored within a column of wetland surface water can be a considerable component of the surface energy budget, an attribute that is demonstrated by comparing changes in stored heat energy to net radiation at seven sites in the wetland areas of southern Florida, including the Everglades. The magnitude of changes in stored heat energy approached the magnitude of net radiation more often during the winter dry season than during the summer wet season. Furthermore, the magnitude of changes in stored heat energy in wetland surface water generally decreased as surface energy budgets were upscaled temporally. A new method was developed to estimate changes in stored heat energy that overcomes an important data limitation, namely, the limited spatial and temporal availability of water temperature measurements. The new method is instead based on readily available air temperature measurements and relies on the convolution of air temperature changes with a regression‐defined transfer function to estimate changes in water temperature. The convolution‐computed water temperature changes are used with water depths and heat capacity to estimate changes in stored heat energy within the Everglades wetland areas. These results likely can be adapted to other humid subtropical wetlands characterized by open water, saw grass, and rush vegetation type communities.

Water Resources Research↗

USGS2004

Comparison of Estimated Areas Contributing Recharge to Selected Springs in North-Central Florida by Using Multiple Ground-Water Flow Models

Areas contributing recharge to springs are defined in this report as the land-surface area wherein water entering the ground-water system at the water table eventually discharges to a spring. These areas were delineated for Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs in north-central Florida using four regional ground-water flow models and particle tracking. As expected, different models predicted different areas contributing recharge. In general, the differences were due to different hydrologic stresses, subsurface permeability properties, and boundary conditions that were used to calibrate each model, all of which are considered to be equally feasible because each model matched its respective calibration data reasonably well. To evaluate the agreement of the models and to summarize results, areas contributing recharge to springs from each model were combined into composite areas. During 1993-98, the composite areas contributing recharge to Blue Spring, Silver Springs, Alexander Springs, and Silver Glen Springs were about 130, 730, 110, and 120 square miles, respectively. The composite areas for all springs remained about the same when using projected 2020 ground-water withdrawals.

Open-File Report↗

USGS2003

MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model--Documentation of the SEAWAT-2000 Version with the Variable-Density Flow Process (VDF) and the Integrated MT3DMS Transport Process (IMT)

SEAWAT-2000 is the latest release of the SEAWAT computer program for simulation of three-dimensional, variable-density, transient ground-water flow in porous media. SEAWAT-2000 was designed by combining a modified version of MODFLOW-2000 and MT3DMS into a single computer program. The code was developed using the MODFLOW-2000 concept of a process, which is defined as ?part of the code that solves a fundamental equation by a specified numerical method.? SEAWAT-2000 contains all of the processes distributed with MODFLOW-2000 and also includes the Variable-Density Flow Process (as an alternative to the constant-density Ground-Water Flow Process) and the Integrated MT3DMS Transport Process. Processes may be active or inactive, depending on simulation objectives; however, not all processes are compatible. For example, the Sensitivity and Parameter Estimation Processes are not compatible with the Variable-Density Flow and Integrated MT3DMS Transport Processes. The SEAWAT-2000 computer code was tested with the common variable-density benchmark problems and also with problems representing evaporation from a salt lake and rotation of immiscible fluids.

Open-File Report↗

USGS2003

Potential for saltwater intrusion into the lower Tamiami aquifer near Bonita Springs, southwestern Florida

A study was conducted to examine the potential for saltwater intrusion into the lower Tamiami aquifer beneath Bonita Springs in southwestern Florida. Field data were collected, and constant- and variable-density ground-water flow simulations were performed that: (1) spatially quantified modern and seasonal stresses, (2) identified potential mechanisms of saltwater intrusion, and (3) estimated the potential extent of saltwater intrusion for the area of concern. MODFLOW and the inverse modeling routine UCODE were used to spatially quantify modern and seasonal stresses by calibrating a constant-density ground-water flow model to field data collected in 1996. The model was calibrated by assuming hydraulic conductivity parameters were accurate and by estimating unmonitored ground-water pumpage and potential evapotranspiration with UCODE. Uncertainty in these estimated parameters was quantified with 95-percent confidence intervals. These confidence intervals indicate more uncertainty (or less reliability) in the estimates of unmonitored ground-water pumpage than estimates of pan-evaporation multipliers, because of the nature and distribution of observations used during calibration. Comparison of simulated water levels, streamflows, and net recharge with field data suggests the model is a good representation of field conditions. Potential mechanisms of saltwater intrusion into the lower Tamiami aquifer include: (1) lateral inland movement of the freshwater-saltwater interface from the southwestern coast of Florida; (2) upward leakage from deeper saline water-bearing zones through natural upwelling and upconing, both of which could occur as diffuse upward flow through semiconfining layers, conduit flow through karst features, or pipe flow through leaky artesian wells; (3) downward leakage of saltwater from surface-water channels; and (4) movement of unflushed pockets of relict seawater. Of the many potential mechanisms of saltwater intrusion, field data and variable-density ground-water flow simulations suggest that upconing is of utmost concern, and lateral encroachment is of second-most concern. This interpretation is uncertain, however, because the predominance of saltwater intrusion through leaky artesian wells with connection to deeper, more saline, and higher pressure aquifers was difficult to establish. Effective management of ground-water resources in southwestern Florida requires an understanding of the potential extent of saltwater intrusion in the lower Tamiami aquifer near Bonita Springs. Variable-density, ground-water flow simulations suggest that when saltwater is at dynamic equilibrium with 1996 seasonal stresses, the extent of saltwater intrusion is about 100 square kilometers areally and 70,000 hectare-meters volumetrically. The volumetric extent of saltwater intrusion was most sensitive to changes in recharge, ground-water pumpage, sea level, salinity of the Gulf of Mexico, and the potentiometric surface of the sandstone aquifer, respectively.

Florida↗