Florida Administrative Code and Florida Statues require water management districts in Florida to identify “specific geographical areas that have water resource problems which have become critical or are anticipated to become critical within the next 20 years.” SJRWMD is assessing water supply needs and sources using regional numerical ground water flow models and local analytical ground water flow models to project which areas are expected to have inadequate water resources to meet the anticipated water supply demand in the year 2010.
To ensure adequate placement of the locations of the regional and local ground water flow models and to identify areas were future work may need to be performed, SJRWMD completed a study, using a geographic information system (GIS) methodology, parallel to the ground water modeling that identifies areas having the potential for public water supply problems. The study discussed in this paper details the methodology and results using a GIS approach to identify areas having a high potential for (1) vegetative impacts or changes, (2) salt water intrusion, and/or (3) an increase in public water supply demand.
The study area covers the entire SJRWMD. SJRWMD includes all or part of 19 counties in northeast Florida, extending from the Florida-Georgia state line in the north to Vero Beach in the south. The eastern border of SJRWMD is the Atlantic Ocean, and the area of jurisdiction extends as far west as central Florida. Wetlands, lakes, and springs are important features of the landscape; however, they are not persistent and are dispersed throughout SJRWMD. Land surface elevations change from 0 feet above mean sea level (msl) to approximately 310 feet msl.
The hydrogeology and other related parameters that would affect public water supply problems vary significantly throughout SJRWMD. The hydrogeology typically consists of a surficial aquifer, a semi-confining unit, and the Floridan aquifer, which is the primary source of water supply. In some areas the surficial aquifer is productive enough to be a source for public water supply, and in other areas the surficial aquifer is non-existent. The hydrogeology also varies from regions where no confining unit is present to regions with a confining unit many hundreds of feet thick. The Floridan aquifer, in some locations, is near land surface and in other locations it is at depths of 500 to 1,000 feet. Chloride concentrations vary from no potable water through the entire Floridan aquifer to very good water quality. The decline in the elevation of the potentiometric surface of the Floridan aquifer varies in SJRWMD based upon the hydrogeology and amount of new withdrawals in the last 10 years and the aggregation of withdrawals over the last 50 years.
The methodology essentially was an overlay procedure in the GRID module of ARC/INFO. Boniol et al. (1993) determined that a 380-foot cell size provided a good calibration between the actual and simulated 1990 potentiometric surfaces; therefore, all the grids were created with a grid cell size of 380 feet.
Eight surfaces were created, each relating in some way to (1) vegetative impacts, (2) salt water intrusion, and/or (3) increases in public water supply demand. Based upon the distribution of the data and the hydrology, the data within each grid were weighted from 1 to 5, where 1 represented a factor having the least affect on the impacts being evaluated and 5 represented a factor having the greatest affect on the impacts being evaluated. The grids also were weighted against each other on the same scaling system, from 1 to 5, based upon the interrelationships among the grids which would influence the impacts being evaluated in this study. Each grid was multiplied by the designated weight, then all of the grids were added together to produce a final grid of scores.
The predevelopment surface, created primarily from 1936 data, was a U.S. Geological Survey (USGS) contour map (Johnston et al., 1980) where the contours represented the potentiometric surface of the Floridan aquifer in feet above mean sea level. The map was digitized into the GIS. The coverage was processed through the triangulated irregular network (TIN) module, then through the GRID module to create the final grid representing the predevelopment potentiometric surface of the Floridan aquifer.
The 1990 potentiometric surface of the Floridan aquifer (Boniol et. al., 1993) was similarly created and used previously as part of the SJRWMD recharge project. This surface began originally as a point coverage of monitoring well water levels from May 1990. The point coverage was processed through the TIN module and then the GRID module to create a grid of the 1990 potentiometric surface of the Floridan aquifer.
The final surface was calculated by subtracting the predevelopment potentiometric surface grid from the 1990 potentiometric surface grid. This procedure produced a surface showing the amount of long-term decline in the potentiometric surface.
The 1980 potentiometric surface was created using monitoring well water levels from May of 1980. The coordinates and water levels of the well points were used to create a point coverage. A TIN surface was created from the point coverage and a grid surface was created from the TIN surface using the ARC/INFO GRID module. The 1990 potentiometric surface, previously described, was again used in the analysis. The 1980 potentiometric surface grid was subtracted from the 1990 potentiometric surface grid to create a grid showing the short-term term decline in the potentiometric surface change of the Floridan aquifer.
The thickness of the confining unit surface also was used as part of the SJRWMD recharge project (Boniol et al., 1993). The initial coverages were compiled by evaluating geophysical logs throughout SJRWMD to determine the top of the confining unit and the top of the Floridan aquifer. Both of these coverages were processed through the TIN module and through the GRID module. Then the top of the confining unit was subtracted from the top of the Floridan aquifer to obtain the thickness of the confining unit above the Floridan aquifer.
The vertical head difference grid was created originally for use in the SJRWMD recharge project (Boniol et al., 1993). Two grids were used to create the vertical head difference grid. The first grid was the 1990 potentiometric surface grid already discussed above. The second grid was a grid of surficial aquifer water levels. To create the surficial aquifer water level grid a regression equation (Y= 0.901x - 1.61) was used showing the relationship between land surface elevation and the depth to the water table. A grid of hypsography was then processed through the regression equation to create a depth to the water table surface. The depth to the water table was subtracted from hypsography to yield elevations of the surficial aquifer. The 1990 potentiometric surface grid was subtracted from the surficial aquifer water level grid to produce the vertical head difference between the surficial and Floridan aquifers grid.
To produce this surface, first a coverage was produced showing the thickness of potable water in the Floridan aquifer. Causey and Leve (1977), of the USGS, prepared a water quality map that was used as a base map for this surface. The map was digitized into the GIS and the contours were then edited based upon current information from ground water flow models, monitoring wells, and time-domain electromagnetic (TDEM) studies. A TIN surface was created from the contours, then a grid surface was created from the TIN surface. The grid surface was evaluated to ensure it matched the edited contour coverage.
A map showing the total thickness of the Floridan aquifer was digitized into the GIS from Miller (1986). This contour coverage was processed through the TIN module, then through the GRID module to produce a grid representing the total thickness of the Floridan aquifer.
To calculate the percentage of potable water in the Floridan aquifer, the following equation was used:
(thickness of potable water MINUS total thickness of Floridan aquifer) MULTIPLIED BY 100 DIVIDED BY total thickness of Floridan aquifer
The area of no potable water, identified as the area where the chlorides were greater than 250 mg/L throughout the entire aquifer, was selected out of the thickness of potable water coverage. The ARC/INFO NEAR function, which calculates the distance from any point to the nearest point, node, or arc, was used to calculate the distance of any potential wellfield from the area of no potable water. The NEAR function in ARC/INFO was run on the point coverage to determine the distance at each point, representing a potential wellfield, to the areas delineated as having no potable water.
The wetland surface was created from the hydrography layer in the USGS Digital Line Graph (DLG) coverage. The layers designated as wetlands were reselected from the coverage and placed in a separate coverage. The layers reselected were hydrography, intermittent water bodies, deep water wetlands, marsh or wetlands, and mangroves. The polygon coverage of wetlands was processed through the GRID module to create a grid representing wetlands.
The projected 2010 population density surface was created using a coverage of 1990 census tract population data. The 1990 population data were projected to 2010 by multiplying the census tracts within each county by a percentage representing the projected population change between 1990 and 2010 on a county level. After calculating the projected 2010 population figures, the population value for each census tract was divided by the area of the census tract to obtain the projected 2010 population density in persons per square mile. Dividing the 2010 population figures by the area of each census tract normalizes the data by accounting for varying population figures over varying census tract sizes.
A histogram distribution of the data was created for each grid. In some cases the distribution of the values led to the assigned weights, as was the case with the potentiometric surface grids. For other surfaces, weights were determined based upon knowledge of hydrology and geology and the relative influence of the factor or surface on (1) vegetative impacts, (2) salt water intrusion, and/or (3) an increase in water supply demand.
The vertical head difference between the surficial and Floridan aquifers grid was difficult to categorize into weights because at one extreme vegetative impacts occur and at the other extreme saltwater impacts occur. It was assumed that saltwater intrusion impacts could be more accurately accounted for, therefore, areas with low values having the potential for saltwater intrusion were weighted the highest. Specifically, in the wetlands surface, it was assumed that there would be little or no affect in locations where there were no wetlands; consequently, wetlands were weighted the highest.
Table 1. Grid weights
SURFACE...................................WEIGHT long-term change in the......................1 potentiometric surface of the Floridan aquifer short-term change in the................ ....3 potentiometric surface of the Floridan aquifer confining unit thickness above the...........3 Floridan aquifer vertical head difference between the........ 3 surficial and Floridan aquifers percentage of potable water in the...........1 Floridan aquifer horizontal distance to areas of no...........1 potable water in the Floridan aquifer wetlands.....................................3 projected 2010 population density............5
The weighted grids were overlaid and weighted against each other before being added together to create a cumulative final score grid. The horizontal distance to no potable water grid and the percentage of fresh water grid were weighted using a value of 1 because they are conceptually and physically similar to each other. This weighting procedure causes the general distribution of the data to automatically appear twice in the overall scores. The long-term change in the potentiometric surface grid was weighted using a value of 1, and the short-term change in the potentiometric surface grid was weighted using a value of 3 because the current and near-future changes were more important to this study. The confining unit thickness above the Floridan aquifer grid was weighted with a value of 3 because it is an important factor affecting vegetative impacts. The vertical head difference between the surficial and Floridan aquifers grid was also weighted with a value of 3 because it affects both vegetative impacts and saltwater upconing. The wetlands grid was weighted with a value of 3 because of the high sensitivity of wetlands to changes in the water table and a general emphasis in Florida to preserve and/or enhance wetland functionality. The projected population density grid was weighted the highest, with a value of 5, because this is the only surface that accounts for projected future changes.
The possible range of scores on the final grid was from 20 to 100, where 20 indicated the least potential for impacts and 100 indicated the greatest potential for impacts. The distribution fell in to a classic bell-shaped histogram. The low, medium, and high potential categories were selected based upon the distribution. The public water supply service areas and regional ground water models were overlaid on top of the grid.
The low potential category was selected to have scores less than 40, the medium potential category, where the bell of the histogram occurred, was selected to be between the scores of 40 and 60, and the high potential category was selected to have scores greater than 60.
The locations of the water supply service areas further defined those areas that need special attention, although it should be noted that some wellfields do occur outside of the water supply service areas. The cumulative effect of all the wellfields in an area can cause a general decline in the elevation of the potentiometric surface outside the water supply service areas.
The locations of the regional ground water flow models were overlaid on the final scored surface and the distribution of the regional ground water flow models shows that virtually all of the areas that have high scores fall within the flow model boundaries. It may be necessary to develop more site-specific subregional ground water flow or transport models to fully understand the effects associated with future water supply development. Several cities that fall outside any of the regional ground water flow models, including Ocala and Vero Beach, have been evaluated with local-scale analytical models.
The methodology used in this study provides a tool which could be used to identify areas in need of water resource assessment and as an approach to evaluate the need for and location of new numerical and analytical ground water flow models.
Boniol, D., M. Williams, and D. Munch. 1993. Mapping recharge to the Floridan aquifer using a geographic information system. Technical Publication SJ93-5. Palatka, Fla.: St. Johns River Water Management District.
Causey, L.V., and G.W. Leve. 1977. Thickness of the potable-water zone in the Floridan aquifer. Map Series No. 74. Tallahassee, Fla.: U.S. Geological Survey,
[ESRI]. Environmental Systems Research Institute. 1992.. ARC/INFO Release 6.1.1, Redlands, Calif.
Johnston, R.H., R.E. Krause, F.W. Meyer, P.D. Ryder, C.H. Tibbals, and J.D. Hunn. 1980. Estimated potentiometric surface for the Tertiary limestone aquifer system, south-eastern United States prior to development. Open File Report 80-406. Atlanta, Ga.: U.S. Geological Survey.
Miller, J.A. 1986. Hydrogeologic framework of the Floridan aquifer system in Florida and parts of Georgia, Alabama and South Carolina. Professional Paper 1403-B. Tallahassee, Fla.: U.S. Geological Survey.
home page