Search

Guha, Hillol and Sorab Panday, 2012. Impact of Sea Level Rise on Groundwater Salinity in a Coastal Community of South Florida. Journal of the American Water Resources Association (JAWRA) 48(3): 510-529. DOI: 10.1111/j.1752-1688.2011.00630.x Abstract: Freshwater resources of coastal communities in the United States and world over are threatened by the rate of sea level rise. According to recent estimates by various governmental agencies and climate researchers, the global sea level rise is likely to be between 0.6 and 2.1 m by the year 2100. South Florida is a coastal community and much of its coastline is subject to sea level rise and potential impacts to wetlands and the water resources of the area. To understand what the impact of sea level rise would cause to the groundwater level and salinity intrusion, an integrated groundwater and surface water model was developed for North Miami-Dade and Broward Counties of South Florida. The model was calibrated against daily groundwater heads, base flows in canals, and chloride concentrations for a period of one year and six months. Three separate sensitivity analyses were conducted by increasing the sea level by 0.6, 0.9, and 1.22 m. Results of the simulations shows increase of groundwater heads in some areas from 4 to 15%; whereas the average relative chloride concentrations increased significantly by 100-600% in some wells. The increase in groundwater elevations and chloride concentrations varies from location of the wells and its proximity to the coast. The model results indicate that even a 0.6 m increase in sea level (which is the conservative estimate) is likely to impair the vital freshwater resources in many parts of South Florida. [ABSTRACT FROM AUTHOR]

Various approaches have been proposed to manage the nonlinearities associated with the unconfined flow equation and to simulate perched groundwater conditions using the MODFLOW family of codes. The approaches comprise a variety of numerical techniques to prevent dry cells from becoming inactive and to achieve a stable solution focused on formulations of the unconfined, partially-saturated, groundwater flow equation. Keeping dry cells active avoids a discontinuous head solution which in turn improves the effectiveness of parameter estimation software that relies on continuous derivatives. Most approaches implement an upstream weighting of intercell conductance and Newton-Raphson linearization to obtain robust convergence. In this study, several published approaches were implemented in a stepwise manner into MODFLOW for comparative analysis. First, a comparative analysis of the methods is presented using synthetic examples that create convergence issues or difficulty in handling perched conditions with the more common dry-cell simulation capabilities of MODFLOW. Next, a field-scale three-dimensional simulation is presented to examine the stability and performance of the discussed approaches in larger, practical, simulation settings. [ABSTRACT FROM AUTHOR]

A diverse group of problems requires quantification of the entire hydrologic cycle by the integrated simulation of water flow in the surface and subsurface regimes. In a transient integrated simulation of the water cycle, the time step size is a key factor in controlling the solution accuracy and the simulation efficiency for a given spatial discretization. In general, if the time step size is sufficiently small, the resulting solution becomes more accurate but with higher computational cost. Thus, to maintain an acceptable level of solution accuracy in the entire simulation domain, the time step size is restricted by the relatively rapid responses in the surface flow regime. As the relatively rapid responses are typically limited to a small portion of the surface domain compared with the groundwater system, a large portion of the domain tends to be temporally overdiscretized. The implicit subtime stepping approach described here can apply smaller subtime steps only to the subdomain where the accuracy requirements are needed. In this work, generalized formulations for implicit subtime stepping in the numerical solution of the nonlinear coupled surface-subsurface equations were derived and implemented into the integrated model HydroGeoSphere. Application to several problems showed that implicit subtime stepping can significantly improve the simulation efficiency with minimal loss in accuracy. The methodology was successfully applied to enhance the computational efficiency of an integrated flow simulation in the San Joaquin Valley, California, where the characteristic response time near surface drainage streams is orders of magnitude shorter than in the groundwater regime. [ABSTRACT FROM AUTHOR]