Geology Reports⌕ Search

Geology topics

Eric D. Swain

Publications and source records attributed to Eric D. Swain.

36 records · Page 2Linked to original sources
USGS2010

Temperature inverted haloclines provide winter warm-water refugia for manatees in southwest Florida

Florida manatees ( Trichechus manatus latirostris ) overwintering in the Ten Thousand Islands and western Everglades have no access to power plants or major artesian springs that provide warm-water refugia in other parts of Florida. Instead, hundreds of manatees aggregate at artificial canals, basins, and natural deep water sites that act as passive thermal refugia (PTR). Monitoring at two canal sites revealed temperature inverted haloclines, which provided warm salty bottom layers that generally remained above temperatures considered adverse for manatees. At the largest PTR, the warmer bottom layer disappeared unless significant salt stratification was maintained by upstream freshwater inflow over a persistent tidal wedge. A detailed three-dimensional hydrology model showed that salinity stratification inhibited vertical convection induced by atmospheric cooling. Management or creation of temperature inverted haloclines may be a feasible and desirable option for resource managers to provide passive thermal refugia for manatees and other temperature sensitive aquatic species.

Florida↗
USGS2010

Spatial and stage-structured population model of the American crocodile for comparison of comprehensive Everglades Restoration Plan (CERP) alternatives

As part of the U.S. Geological Survey Priority Ecosystems Science (PES) initiative to provide the ecological science required during Everglades restoration, we have integrated current regional hydrologic models with American crocodile (Crocodylus acutus) research and monitoring data to create a model that assesses the potential impact of Comprehensive Everglades Restoration Plan (CERP) efforts on the American crocodile. A list of indicators was created by the Restoration Coordination and Verification (RECOVER) component of CERP to help determine the success of interim restoration goals. The American crocodile was established as an indicator of the ecological condition of mangrove estuaries due to its reliance upon estuarine environments characterized by low salinity and adequate freshwater inflow. To gain a better understanding of the potential impact of CERP restoration efforts on the American crocodile, a spatially explicit crocodile population model has been created that has the ability to simulate the response of crocodiles to various management strategies for the South Florida ecosystem. The crocodile model uses output from the Tides and Inflows in the Mangroves of the Everglades (TIME) model, an application of the Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) simulator. TIME has the capability to link to the South Florida Water Management Model (SFWMM), which is the primary regional tool used to assess CERP restoration scenarios. A crocodile habitat suitability index and spatial parameter maps that reflect salinity, water depth, habitat, and nesting locations are used as driving functions to construct crocodile finite rate of increase maps under different management scenarios. Local stage-structured models are integrated with a spatial landscape grid to display crocodile movement behavior in response to changing environmental conditions. Restoration efforts are expected to affect salinity levels throughout the habitat of the American crocodile. This modeling effort examines how CERP restoration alternatives will affect growth and survival rates of hatchling and juvenile crocodiles, hatchling dispersal to suitable nursery habitat, and relative abundance and distribution in response to changing salinity and water depth for all stage classes of crocodiles. The response of the American crocodile to restoration efforts will provide a quantifiable measure of restoration success. By applying the crocodile model to proposed restoration alternatives and predicting population responses, we can choose alternatives that approximate historical conditions, enhance habitat for multiple species, and identify future research needs.

Open-File Report↗
USGS2009

Integrated Science: Florida Manatees and Everglades Hydrology

Predicting and monitoring restoration effects on Florida manatees, which are known to make extended movements, will be incomplete if modeling and monitoring are limited to the smaller areas defined by the various res-toration components. U.S. Geological Survey (USGS) efforts, thus far, have focused on (1) collecting manatee movement data throughout the Ten Thousand Islands (TTI) region, and (2) developing an individual-based model for manatees to illustrate manatee responses to changes in hydrology related to the Picayune Strand Restoration Project (PSRP). In 2006, new regional research was begun to extend an Everglades hydrology model into the TTI region; extend the manatee movement model into the southern estuaries of Everglades National Park (ENP); and integrate hydrology and manatee data, models, and monitoring across the TTI region and ENP. Currently (2008), three research tasks are underway to develop the necessary modeling components to assess restoration efforts across the Greater Everglades Ecosystem.

Fact Sheet↗
USGS2009

Development, Testing, and Application of a Coupled Hydrodynamic Surface-Water/Groundwater Model (FTLOADDS) with Heat and Salinity Transport in the Ten Thousand Islands/Picayune Strand Restoration Project Area, Florida

A numerical model application was developed for the coastal area inland of the Ten Thousand Islands (TTI) in southwestern Florida using the Flow and Transport in a Linked Overland/Aquifer Density-Dependent System (FTLOADDS) model. This model couples a two-dimensional dynamic surface-water model with a three-dimensional groundwater model, and has been applied to several locations in southern Florida. The model application solves equations for salt transport in groundwater and surface water, and also simulates surface-water temperature using a newly enhanced heat transport algorithm. One of the purposes of the TTI application is to simulate hydrologic factors that relate to habitat suitability for the West Indian Manatee. Both salinity and temperature have been shown to be important factors for manatee survival. The inland area of the TTI domain is the location of the Picayune Strand Restoration Project, which is designed to restore predevelopment hydrology through the filling and plugging of canals, construction of spreader channels, and the construction of levees and pump stations. The effects of these changes are simulated to determine their effects on manatee habitat. The TTI application utilizes a large amount of input data for both surface-water and groundwater flow simulations. These data include topography, frictional resistance, atmospheric data including rainfall and air temperature, aquifer properties, and boundary conditions for tidal levels, inflows, groundwater heads, and salinities. Calibration was achieved by adjusting the parameters having the largest uncertainty: surface-water inflows, the surface-water transport dispersion coefficient, and evapotranspiration. A sensitivity analysis did not indicate that further parameter changes would yield an overall improvement in simulation results. The agreement between field data from GPS-tracked manatees and TTI application results demonstrates that the model can predict the salinity and temperature fluctuations which affect manatee behavior. Comparison of the existing conditions simulation with the simulation incorporating restoration changes indicated that the restoration would increase the period of inundation for most of the coastal wetlands. Generally, surface-water salinity was lowered by restoration changes in most of the wetlands areas, especially during the early dry season. However, the opposite pattern was observed in the primary canal habitat for manatees, namely, the Port of the Islands. Salinities at this location tended to be moderately elevated during the dry season, and unchanged during the wet season. Water temperatures were in close agreement between the existing conditions and restoration simulations, although minimum temperatures at the Port of the Islands were slightly higher in the restoration simulation as a result of the additional surface-water ponding and warming that occurs in adjacent wetlands. The TTI application output was used to generate salinity and temperature time series for comparison to manatee field tracking data and an individually-based manatee-behavior model. Overlaying field data with salinity and temperature results from the TTI application reflects the effect of warm water availability and the periodic need for low-salinity drinking water on manatee movements. The manatee-behavior model uses the TTI application data at specific model nodes along the main manatee travel corridors to determine manatee migration patterns. The differences between the existing conditions and restoration scenarios can then be compared for manatee refugia. The TTI application can be used to test a variety of hydrologic conditions and their effect on important criteria.

Scientific Investigations Report↗
USGS2007

Application of FTLOADDS to Simulate Flow, Salinity, and Surface-Water Stage in the Southern Everglades, Florida

The Comprehensive Everglades Restoration Plan requires numerical modeling to achieve a sufficient understanding of coastal freshwater flows, nutrient sources, and the evaluation of management alternatives to restore the ecosystem of southern Florida. Numerical models include a regional water-management model to represent restoration changes to the hydrology of southern Florida and a hydrodynamic model to represent the southern and western offshore waters. The coastal interface between these two systems, however, has complex surface-water/ground-water and freshwater/saltwater interactions and requires a specialized modeling effort. The Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) code was developed to represent connected surface- and ground-water systems with variable-density flow. The first use of FTLOADDS is the Southern Inland and Coastal Systems (SICS) application to the southeastern part of the Everglades/Florida Bay coastal region. The need to (1) expand the domain of the numerical modeling into most of Everglades National Park and the western coastal area, and (2) better represent the effect of water-delivery control structures, led to the application of the FTLOADDS code to the Tides and Inflows in the Mangroves of the Everglades (TIME) domain. This application allows the model to address a broader range of hydrologic issues and incorporate new code modifications. The surface-water hydrology is of primary interest to water managers, and is the main focus of this study. The coupling to ground water, however, was necessary to accurately represent leakage exchange between the surface water and ground water, which transfers substantial volumes of water and salt. Initial calibration and analysis of the TIME application produced simulated results that compare well statistically with field-measured values. A comparison of TIME simulation results to previous SICS results shows improved capabilities, particularly in the representation of coastal flows. This improvement most likely is due to a more stable numerical representation of the coastal creek outlets. Sensitivity analyses were performed by varying frictional resistance, leakage, barriers to flow, and topography. Changing frictional resistance values in inland areas was shown to improve water-level representation locally, but to have a negligible effect on area-wide values. These changes have only local effects and are not physically based (as are the unchanged values), and thus have limited validity. Sensitivity tests indicate that the overall accuracy of the simulation is diminished if leakage between surface water and ground water is not simulated. The inclusion of a major road as a complete barrier to surface-water flow influenced the local distribution and timing of flow; however, the changes in total flow and individual creekflows were negligible. The model land-surface altitude was lowered by 0.1 meter to determine the sensitivity to topographic variation. This topographic sensitivity test produced mixed results in matching field data. Overall, the representation of stage did not improve definitively. A final calibration utilized the results of the sensitivity analysis to refine the TIME application. To accomplish this calibration, the friction coefficient was reduced at the northern boundary inflow and increased in the southwestern corner of the model, the evapotranspiration function was varied, additional data were used for the ground-water head boundary along the southeast, and the frictional resistance of the primary coastal creek outlet was increased. The calibration improved the match between measured and simulated total flows to Florida Bay and coastal salinities. Agreement also was improved at most of the water-level sites throughout the model domain.

Scientific Investigations Report↗
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↗
USGS2005

Implementation and use of direct-flow connections in a coupled ground-water and surface-water model

The U.S. Geological Survey's MODFLOW finite-difference ground-water flow model has been coupled with three surface-water packages - the MODBRANCH, River, and Stream packages - to simulate surface water and its interaction with ground water. Prior to the development of the coupling packages, the only interaction between these modeling packages was that leakage values could be passed between MODFLOW and the three surface-water packages. To facilitate wider and more flexible uses of the models, a computer program was developed and added to MODFLOW to allow direct flows or stages to be passed between any of the packages and MODFLOW. The flows or stages calculated in one package can be set as boundary discharges or stages to be used in another package. Several modeling packages can be used in the same simulation depending upon the level of sophistication needed in the various reaches being modeled. This computer program is especially useful when any of the River, Stream, or MODBRANCH packages are used to model a river flowing directly into or out of wetlands in direct connection with the aquifer and represented in the model as an aquifer block. A field case study is shown to illustrate an application.

Groundwater↗
USGS2005

A model for simulation of surface-water integrated flow and transport in two dimensions: user's guide for application to coastal wetlands

The computation of hydrodynamic flow in two dimensions is implemented in the Surface-Water Integrated Flow and Transport in Two-Dimensions (SWIFT2D) model using a numerical code that also includes computation of reactive constituent transport, density variation effects, drying and wetting of periodically inundated areas, and hydraulic structures. The model can be utilized in a variety of settings where velocity and concentration gradients can be assumed to have minimal vertical variations, and the representation of flow in two horizontal dimensions is sufficient. The finite-difference forms of the equations of mass continuity and momentum are solved in two dimensions by the use of a staggered timestep solution solved with an efficient alternating-direction implicit solution. The finite-difference forms of the constituent transport equations can be solved in conjunction with the flow equations. If salt transport is simulated, an equation of state relates the density terms in the momentum equation to salinity. The traditional formulation of SWIFT2D has been applied to numerous estuaries, bays, and harbors throughout the world. Modifications have been made to expand SWIFT2D for applicability to shallow coastal wetlands. These modifications include the representation of spatially and temporally varying rainfall and evapotranspiration, wind sheltering owing to effects of emergent vegetation, and changes in frictional resistance with depth. These modifications expand the versatility of the code?s applications to include open freshwater or saltwater conditions along coasts and within embayments and estuaries as well as associated fresh, brackish, and hypersaline wetlands and marshes linked to such water bodies. Inclusion of precipitation and evapotranspiration processes also permits long-term simulations.

Open-File Report↗
USGS2004

Assigning boundary conditions to the Southern Inland and Coastal Systems (SICS) model using results from the South Florida Water Management Model (SFWMM)

The Comprehensive Everglades Restoration Plan (CERP) requires the testing and evaluation of different water-management scenarios for southern Florida. As part of CERP, the South Florida Water Management District is using its regional hydrologic model, the South Florida Water Management Model (SFWMM), to evaluate different hydrologic scenarios. The SFWMM was designed specifically for the inland freshwater areas in southern Florida, and extends only slightly into Florida Bay. Thus, the U.S. Geological Survey developed the Southern Inland and Coastal Systems (SICS) model, which is an integrated surface-water and ground-water model designed to simulate flows, stages, and salinities in the southern Everglades and Florida Bay. Modifications to the SICS boundary conditions allow the local-scale SICS model to be linked to the regional-scale SFWMM. The linked model will be used to quantify the effects of restoration alternatives on flows, stages, and salinities in the SICS area. This report describes the procedure for linking the SICS model with the SFWMM. The linkage is shown to work by comparing the results of a linked 5-year simulation with the results from a simulation in which the model boundaries are assigned using field data. The surface-water module of the SICS model is driven by areal influences and lateral boundaries. The areal influences (wind, rainfall, and evapotranspiration) remain the same when the SICS model is modified to link to the SFWMM. Four types of lateral boundaries (discharge, water level, no flow, and salinity) are used in the SICS model. Two of three discharge boundaries (at Taylor Slough Bridge and C-111 Canal) in the current SICS model domain are converted to water-level boundaries to increase accuracy. The only change to the third discharge boundary (at Levee 31W) is that the flow data are derived from SFWMM model output instead of using measured field data flows. Three water-level boundaries are modified only by receiving their data from SFWMM model output data. Additionally, two marine water-level boundaries remain the same because the SFWMM does not include Florida Bay and, therefore, this model cannot provide input data for these boundaries. The SICS no-flow boundaries remain intact because no additional data, provided by the SFWMM, suggest that any significant flow occurs along these boundaries. The Florida Bay salinity boundary is not modified because the SFWMM does not contain any salinity data that can be used to modify the model. The ground-water module of the SICS model contains a general-head boundary and a no-flow boundary. The general-head boundary, which extends along the edges of the wetland part of the SICS model domain, is modified by acquiring stage values from SFWMM cells that correspond in location to the SICS model cells. Values from the SFWMM cells are bilinearly interpolated and assigned to the appropriate SICS general-head boundary cells in all layers of the ground-water model. The ground-water no-flow boundary in Florida Bay is unaltered because the SFWMM does not include this area. A 5-year simulation was developed to test the linkage of the SICS model with the SFWMM. Results from the linked model are similar to those obtained from the original SICS model in which boundaries are assigned using field data. The simulated discharges at the coastal creeks along Florida Bay are about 5 percent lower than the field data simulation; water levels in the wetlands are about 4 percent lower, and salinities at the various coastal creeks are slightly higher.

Florida↗
USGS2004

Simulation of Integrated Surface-Water/Ground-Water Flow and Salinity for a Coastal Wetland and Adjacent Estuary

The SWIFT2D surface-water flow and transport code, which solves the St. Venant equations in two dimensions, was coupled with the SEAWAT variable-density ground-water code to represent hydrologic processes in coastal wetlands and adjacent estuaries. The integrated code was applied to the southern Everglades of Florida to quantify flow and salinity patterns and to evaluate effects of hydrologic processes. Results indicate that most surface water within Taylor Slough flows through Joe Bay and into Florida Bay through Trout Creek. Overtopping of the Buttonwood Embankment, a narrow but continuous ridge that separates the coastal wetlands from Florida Bay, does occur in response to tropical storms, but the net overflow is only 1.5 percent of creek discharge. The net leakage rate for the coastal wetland is about zero with nearly equal upward (17.1 cm/yr) and downward (17.4 cm/yr) rates. During the dry season, the coastal wetland increases in salinity to 30-35 practical salinity units but is flushed each year with the onset of the wet season. Model results demonstrate that surface-water/ground-water interactions, density-dependent flow, and wind affect flow and salinity patterns.

Open-File Report↗
USGS2004

Two-dimensional hydrodynamic simulation of surface-water flow and transport to Florida Bay through the Southern Inland and Coastal Systems (SICS)

Successful restoration of the southern Florida ecosystem requires extensive knowledge of the physical characteristics and hydrologic processes controlling water flow and transport of constituents through extremely low-gradient freshwater marshes, shallow mangrove-fringed coastal creeks and tidal embayments, and near-shore marine waters. A sound, physically based numerical model can provide simulations of the differing hydrologic conditions that might result from various ecosystem restoration scenarios. Because hydrology and ecology are closely linked in southern Florida, hydrologic model results also can be used by ecologists to evaluate the degree of ecosystem restoration that could be achieved for various hydrologic conditions. A robust proven model, SWIFT2D, (Surface-Water Integrated Flow and Transport in Two Dimensions), was modified to simulate Southern Inland and Coastal Systems (SICS) hydrodynamics and transport conditions. Modifications include improvements to evapotranspiration and rainfall calculation and to the algorithms that describe flow through coastal creeks. Techniques used in this model should be applicable to other similar low-gradient marsh settings in southern Florida and elsewhere. Numerous investigations were conducted within the SICS area of southeastern Everglades National Park and northeastern Florida Bay to provide data and parameter values for model development and testing. The U.S. Geological Survey and the National Park Service supported investigations for quantification of evapotranspiration, vegetative resistance to flow, wind-induced flow, land elevations, vegetation classifications, salinity conditions, exchange of ground and surface waters, and flow and transport in coastal creeks and embayments. The good agreement that was achieved between measured and simulated water levels, flows, and salinities through minimal adjustment of empirical coefficients indicates that hydrologic processes within the SICS area are represented properly in the SWIFT2D model, and that the spatial and temporal resolution of these processes in the model is adequate. Sensitivity analyses were conducted to determine the effect of changes in boundary conditions and parameter values on simulation results, which aided in identifying areas of greatest uncertainty in the model. The parameter having the most uncertainty (most in need of further field study) was the flow coefficient for coastal creeks. Smaller uncertainties existed for wetlands frictional resistance and wind. Evapotranspiration and boundary inflows indicated the least uncertainty as determined by varying parameters used in their formulation and definition. Model results indicated that wind was important in reversing coastal creek flows. At Trout Creek (the major tributary connecting Taylor Slough wetlands with Florida Bay), flow in the landward direction was not simulated properly unless wind forcing was included in the simulation. Simulations also provided insight into the major influence that wind has on salinity mixing along the coast, the varying distribution of wetland flows at differing water levels, and the importance of topography in controlling flows to the coast. Slight topographic variations were shown to highly influence the routing of water. A multiple regression analysis was performed to relate inflows at the northern boundary of Taylor Slough bridge to a major pump station (S-332) north of the SICS model area. This analysis allows Taylor Slough bridge boundary conditions to be defined for the model from operating scenarios at S-332, which should facilitate use of the SICS model as an operational tool.

Florida↗
USGS2003

An analytical formulation of two‐dimensional groundwater dispersion induced by surficial recharge variability

A predominant cause of dispersion in groundwater is advective mixing due to variability in seepage rates. Hydraulic conductivity variations have been extensively researched as a cause of this seepage variability. In this paper the effect of variations in surface recharge to a shallow surficial aquifer is investigated as an important additional effect. An analytical formulation has been developed that relates aquifer parameters and the statistics of recharge variability to increases in the dispersivity. This is accomplished by solving Fourier transforms of the small perturbation forms of the groundwater flow equations. Two field studies are presented in this paper to determine the statistics of recharge variability for input to the analytical formulation. A time series of water levels at a continuous groundwater recorder is used to investigate the temporal statistics of hydraulic head caused by recharge, and a series of infiltrometer measurements are used to define the spatial variability in the recharge parameters. With these field statistics representing head fluctuations due to recharge, the analytical formulation can be used to compute the dispersivity without an explicit representation of the recharge boundary. Results from a series of numerical experiments are used to define the limits of this analytical formulation and to provide some comparison. A sophisticated model has been developed using a particle‐tracking algorithm (modified to account for temporal variations) to estimate groundwater dispersion. Dispersivity increases of 9 percent are indicated by the analytical formulation for the aquifer at the field site. A comparison with numerical model results indicates that the analytical results are reasonable for shallow surficial aquifers in which two‐dimensional flow can be assumed.

Water Resources Research↗
USGS1997

Review of selected features of the natural system model, and suggestions for applications in south Florida

A study was conducted to review selected features of the Natural System Model, version 4.3. The Natural System Model is a regional-scale model that uses recent climatic data and estimates of historic vegetation and topography to simulate pre-canal-drainage hydrologic response in south Florida. Equations used to represent the hydrologic system and the numerical solution of these equations in the model were documented and reviewed. Convergence testing was performed using 1965 input data, and selected other aspects of the model were evaluated. Some conclusions from the evaluation of the Natural System Model include the following observations. Simulations were generally insensitive to the temporal resolution used in the model. However, reduction of the computational cell size from 2-mile by 2-mile to 2/3-mile by 2/3-mile resulted in a decrease in spatial mean ponding depths for October of 0.35 foot for a 3-hour time step. Review of the computer code indicated that there is no limit on the amount of water that can be transferred from the river system to the overland flow system, on the amount of seepage from the river to the ground-water system, on evaporation from the river system, or on evapotranspiration from the overland-flow system. Oscillations of 0.2 foot or less in simulated river stage were identified and attributed to a volume limiting function which is applied in solution of the overland-flow equations. The computation of the resistance coefficient is not consistent with the computation of overland-flow velocity. Ground-water boundary conditions do not always ensure a no-flow condition at the boundary. These inconsistencies had varying degrees of effects on model simulations, and it is likely that simulations longer than 1 year are needed to fully identify effects. However, inconsistencies in model formulations should not be ignored, even if the effects of such errors on model results appear to be small or have not been clearly defined. The Natural System Model can be a very useful tool for estimating pre-drainage hydrologic response in south Florida. The model includes all of the important physical processes needed to simulate a water balance. With a few exceptions, these hydrologic processes are represented in a reasonable manner using empirical, semiempirical, and mechanistic relations. The data sets that have been assembled to represent physical features, and hydrologic and meteorological conditions are quite extensive in their scope. Some suggestions for model application were made. Simulation results from the Natural System Model need to be interpreted on a regional basis, rather than cell by cell. The available evidence suggests that simulated water levels should be interpreted with about a plus or minus 1 foot uncertainty. It is probably not appropriate to use the Natural System Model to estimate pre-drainage discharges (as opposed to hydroperiods and water levels) at a particular location or across a set of adjacent computational cells. All simulated results for computational cells within about 10 miles of the model boundaries have a higher degree of uncertainty than results for the interior of the model domain. It is most appropriate to interpret the Natural System Model simulation results in connection with other available information. Stronger linkages between hydrologic inputs to the Everglades and the ecological response of the system would enhance restoration efforts.

Florida↗
USGS1996

South Florida Ecosystem Program: Quantifying freshwater discharge for coastal hydraulic control structures in eastern Dade County, Florida

The South Florida Ecosystem Restoration Program is an intergovernmental effort, involving a number of agencies, to reestablish and maintain the ecosystem of south Florida. One element of the restoration effort is the development of a firm scientific basis for resource decision making. The U.S. Geological Survey (USGS), one of the agencies, provides scientific information as part of the South Florida Ecosystem Restoration Program. The USGS began their ow program, called the South Florida Ecosystem Program, in fiscal year 1995 for the purpose of gathering hydrologic, cartographic, and geologic data that relate to the mainland of south Florida, Florida Bay, and the Florida Keys and Reef ecosystems. As part of the South Florida Ecosystem Program, the USGS, in cooperation with the South Florida Water Management District (SFWMD), has conducted a study to determine discharge ratings for 16 coastal hydraulic control structures in eastern Dade County, Fla. Discharge data are needed to quantify water that can be made available for water supply and ecosystem restoration and to calibrate regional hydrologic models.

Florida↗
USGS1996

Review and evaluation of a model for simulating the natural hydrology of South Florida

The South Florida Ecosystem Program is an intergovernmental effort to re-establish and maintain the ecosystem of South Florida. One element of the restoration effort is the development of a firm scientific basis for making management decisions. The U.S. Geologcal Survey (USGS) is one of the agencies that provides this needed scientific information through the USGS South Florida Ecosystem Program.

Florida↗
USGS1996

A coupled surface-water and ground-water flow model (MODBRANCH) for simulation of stream-aquifer interaction

Ground-water and surface-water flow models traditionally have been developed separately, with interaction between subsurface flow and streamflow either not simulated at all or accounted for by simple formulations. In areas with dynamic and hydraulically well-connected ground-water and surface-water systems, stream-aquifer interaction should be simulated using deterministic responses of both systems coupled at the stream-aquifer interface. Accordingly, a new coupled ground-water and surface-water model was developed by combining the U.S. Geological Survey models MODFLOW and BRANCH; the interfacing code is referred to as MODBRANCH. MODFLOW is the widely used modular three-dimensional, finite-difference ground-water model, and BRANCH is a one-dimensional numerical model commonly used to simulate unsteady flow in open- channel networks. MODFLOW was originally written with the River package, which calculates leakage between the aquifer and stream, assuming that the stream's stage remains constant during one model stress period. A simple streamflow routing model has been added to MODFLOW, but is limited to steady flow in rectangular, prismatic channels. To overcome these limitations, the BRANCH model, which simulates unsteady, nonuniform flow by solving the St. Venant equations, was restructured and incorporated into MODFLOW. Terms that describe leakage between stream and aquifer as a function of streambed conductance and differences in aquifer and stream stage were added to the continuity equation in BRANCH. Thus, leakage between the aquifer and stream can be calculated separately in each model, or leakages calculated in BRANCH can be used in MODFLOW. Total mass in the coupled models is accounted for and conserved. The BRANCH model calculates new stream stages for each time interval in a transient simulation based on upstream boundary conditions, stream properties, and initial estimates of aquifer heads. Next, aquifer heads are calculated in MODFLOW based on stream stages calculated by BRANCH, aquifer properties, and stresses. This process is repeated until convergence criteria are met for head and stage. Because time steps used in ground-water modeling can be much longer than time intervals used in surface- water simulations, provision has been made for handling multiple BRANCH time intervals within one MODFLOW time step. An option was also added to BRANCH to allow the simulation of channel drying and rewetting. Testing of the coupled model was verified by using data from previous studies; by comparing results with output from a simpler, four-point implicit, open-channel flow model linked with MODFLOW; and by comparison to field studies of L-31N canal in southern Florida.

Techniques of Water-Resources Investigations↗
USGS1991

Coupled surface-water and ground-water model

In areas with dynamic and hydraulically well connected ground-water and surface-water systems, it is desirable that stream-aquifer interaction be simulated with models of equal sophistication and accuracy. Accordingly, a new, coupled ground-water and surface-water model was developed by combining the U.S. Geological Survey models MODFLOW and BRANCH. MODFLOW is the widely used modular three-dimensional, finite-difference, ground-water model and BRANCH is a one-dimensional numerical model commonly used to simulate flow in open-channel networks. Because time steps used in ground-water modeling commonly are much longer than those used in surface-water simulations, provision has been made for handling multiple BRANCH time steps within one MODFLOW time step. Verification testing of the coupled model was done using data from previous studies and by comparing results with output from a simpler four-point implicit open-channel flow model linked with MODFLOW.

Conference Paper↗