Inverse modeling of surface-water discharge to achieve restoration salinity performance measures in Florida Bay, Florida

Summary The use of numerical modeling to evaluate regional water-management practices involves the simulation of various alternative water-delivery scenarios, which typically are designed intuitively rather than analytically. These scenario simulations are used to analyze how specific water-management practices affect factors such as water levels, flows, and salinities. In lieu of testing a variety of scenario simulations in a trial-and-error manner, an optimization technique may be used to more precisely and directly define good water-management alternatives. A numerical model application in the coastal regions of Florida Bay and Everglades National Park (ENP), representing the surface- and ground-water hydrology for the region, is a good example of a tool used to evaluate restoration scenarios. The Southern Inland and Coastal System (SICS) model simulates this area with a two-dimensional hydrodynamic surface-water model and a three-dimensional ground-water model, linked to represent the interaction of the two systems with salinity transport. This coastal wetland environment is of great interest in restoration efforts, and the SICS model is used to analyze the effects of alternative water-management scenarios. The SICS model is run within an inverse modeling program called UCODE. In this application, UCODE adjusts the regulated inflows to ENP while SICS is run iteratively. UCODE creates parameters that define inflow within an allowable range for the SICS model based on SICS model output statistics, with the objective of matching user-defined target salinities that meet ecosystem restoration criteria. Preliminary results obtained using two different parameterization methods illustrate the ability of the model to achieve the goals of adjusting the range and reducing the variance of salinity values in the target area. The salinity variance in the primary zone of interest was reduced from an original value of 0.509 psu2 to values 0.418 psu2 and 0.342 psu2 using different methods. Simulations with one, two, and three target areas indicate that optimization is limited near model boundaries and the target location nearest the tidal boundary may not be improved. These experiments indicate that this method can be useful for designing water-delivery schemes to achieve certain water-quality objectives. Additionally, this approach avoids much of the intuitive type of experimentation with different flow schemes that has often been used to develop restoration scenarios.