Translator Disclaimer
1 July 2008 Water Quality Changes in the Guana Tolomato Matanzas National Estuarine Research Reserve, Florida, Associated with Four Tropical Storms
Author Affiliations +

The objective of this study was to document the effects of extreme wind and rainfall conditions associated with tropical storms on physiochemical variability in a tidal creek, Pellicer Creek, in northeast Florida. High-frequency salinity and meteorological data from the Guana Tolomato Matanzas National Estuarine Research Reserve were examined at a range of temporal scales, from 30-minute to annual intervals. Monthly measures of nutrient and water clarity parameters were compared to salinity variations. It was hypothesized that the four tropical storms impacting the region in 2004 (i.e., Charley, Frances, Ivan, and Jeanne) altered tidal regimes and watershed inputs to Pellicer Creek sufficiently to generate water column conditions that deviated significantly from nonstorm periods. The four tropical systems of 2004 suppressed tidally induced salinity variations. Strong northeasterly winds associated with the storm events initially prompted salinity spikes. However, high rainfall levels during the course of each event ultimately caused strong declines in salinity for extended periods of time. Nitrogen concentrations in Pellicer Creek were significantly elevated after storm events. Because primary production in many of the coastal environments along the east coast of Florida, as well as around the world, is nitrogen limited, increases in nitrogen input represent a potential for enhanced algal production and biomass. Given the major changes in watershed characteristics and global climate patterns expected in future years, the ability to predict the influences of these changes on the estuarine environment will be an essential part of designing, implementing, and justifying management efforts.


Salinity, nutrient concentrations, water clarity, hydrology, organismal abundance, and species diversity are among the key elements that define the overall structure and function of estuaries (Berquist et al., 2006; Hackney, Burbanck, and Hackney, 1976; Imberger et al., 1983; Kennish, 2004; McMillan, 1974; Montague and Ley, 2003; Orlando et al., 1994). Recent trends in cultural eutrophication and land use within estuarine watersheds have affected the first four parameters and in turn altered the latter two in many ecosystems (Cloern, 2001; Nixon, 1995). Defining the nature and magnitude of these effects, particularly on broad temporal scales, is difficult without incorporating the influence of variations in meteorological conditions, particularly as it relates to the hydrology of ecosystems (Cloern et al., 1985; Eyre and Pont, 2003; Hama and Handa, 1994; Moore et al., 2006; Sharler and Baird, 2003). Variability in meteorological conditions occurs over a range of temporal scales and intensities, from daily changes in rainfall and wind to storm events and multidecadal shifts in climate (Enfield, Mestas-Nuñez, and Trimble, 2001; Gerten and Adrian, 2000; Goodrich, 1988; Livingston et al., 1997; Scavia et al., 2002).

Storm events, because of their intensity, can alter the structure and function of ecosystems, despite their relatively short duration (Boesch, Diaz, and Virnstein, 1976; Blood et al., 1991; Cahoon, 2006; Greening, Doering, and Corbett, 2006; Mallin and Corbett, 2006; Tab and Jones, 1962; Tilmant et al., 1994; Van Dolah and Anderson, 1991). It is obvious from even cursory empirical observations that storms can have significant and multifaceted effects on estuarine ecosystems, although the magnitude of these alterations is dependent upon the intensity, location, and duration of particular storms. The extreme wind and rainfall conditions associated with storms alter tidal regimes, water circulation patterns, and watershed inputs to estuaries. Documenting these effects can be challenging because most monitoring programs do not deal with the rapid and dramatic changes that occur during storm events, but are limited to sampling under less extreme conditions.

Studies of multiple large storm events involving continuous water quality monitoring, or monitoring more than monthly, are rare. Paerl et al. (2001) described the impacts of three sequential hurricanes in 1999 (i.e., Dennis, Floyd, and Irene) on the Pamlico Sound estuary in North Carolina. Data used in the study included continuous (15-minute sampling intervals) and weekly monitoring in the Neuse and Pamlico River estuaries and monthly poststorm sampling in Pamlico Sound. In the Sound, floodwaters created a combination of vertical stratification and high organic matter loads that resulted in hypoxic conditions for approximately 3 weeks. In another example, Walker (2001) used time-series data from three stations in a southern Louisiana bay system to assess the impacts of two major storms on physical conditions in the water column. He found increases in water level and salinity with onshore winds, as well as decreases in water level and salinity with offshore winds. Large peaks in suspend solids were also associated with the storm events.

The environmental monitoring networks established by the National Estuarine Research Reserve (NERR) program provide an opportunity to examine the effects of storm events on a relevant time scale (Kennish, 2004). In this study, environmental monitoring data from the Guana Tolomato Matanzas (GTM) NERR in northeast Florida were used to investigate the effects of the 2004 Atlantic hurricane season on water quality in a tidal creek ecosystem, Pellicer Creek (Figure 1).

Four tropical storms affected northeast Florida in the summer of 2004. Hurricane Charley, the first major storm of 2004, struck the southwest coast of Florida in the middle of August as a Category 4 on the Saffir-Simpson hurricane scale (Pasch, Brown, and Blake, 2004), then traveled northeasterly across the state and exited near Daytona Beach as a tropical storm (Figure 2). The next storm, Hurricane Frances, made landfall on Florida's southeast coast as a Category 2 on September 5, 2004, then traveled into the Gulf of Mexico and north through the Big Bend region of Florida, affecting most of north Florida (Figure 2; Beven, 2004). The third storm of 2004, Hurricane Ivan, made landfall on the southern coast of Alabama on September 16, 2004, then crossed the southeastern United States before re-entering the Atlantic Ocean and returning south along the Florida coast (Figure 2; Stewart, 2004). Hurricane Jeanne, the final major storm of the 2004 season, made landfall on September 26, 2004, on Florida's east coast, as a Category 3 storm and traveled west-northwest, weakening to a tropical storm as it traveled up the west coast of the state (Figure 2; Lawrence and Cobb, 2004).

GTMNERR weather and water quality observations from 2004 were compared to 2003, which was not subject to any major land-falling storms in Florida. We hypothesized that the rain and wind associated with storms would significantly alter the typical diel patterns of marine and freshwater exchange within Pellicer Creek, thereby affecting a range of important response variables, including salinity, light transparency (i.e., turbidity and dissolved organic matter), nutrient concentrations, and phytoplankton abundances (i.e., chlorophyll a concentrations). One of the elements of consideration in the study was the longevity and magnitude of alterations associated with each storm, as well as potential carryover effects from one storm to the next.


Site Description

Pellicer Creek is a major tributary of the Matanzas River and is located within the GTMNERR. The GTMNERR lies within the southern temperate climate of northeast Florida (Chen and Gerber, 1990). The area has an average of 140 cm of rainfall annually, with a summer wet season from June to September (Figure 3). Pellicer Creek experiences semidiurnal tides with an average range of 0.6 m (NOAA, 2004). The majority of the creek is surrounded by public conservation lands. Salinity at the sampling site ranges from 0 ppt to 33 ppt.

Data Collection

The GTMNERR manages a meteorological station at the mouth of Pellicer Creek and a water quality monitoring station located approximately 4 km upstream from the weather station (Figure 1). At the weather station, data were recorded with a CR10X-2M data logger at 15-minute intervals. At the water quality station, a YSI 6600 continuous monitoring data sonde collected measurements on salinity and water depth at 30-minute intervals (NOAA, 2004). Weather and water quality data from 2003 and 2004 were downloaded from the Centralized Data Management Office website (; NOAA, 2004). As 2003 is the first complete year of data for the GTMNERR continuous monitoring program and was not subject to any major storms in Florida, it was used as a reference year. Even though Florida experienced an unusually active storm season in 2004, the annual rainfall total for the Pellicer Creek watershed was only slightly higher than the rainfall total in 2003 (Figure 3; NOAA, 2006). On an annual basis, both years were near the average rainfall total for 1994–2006. Due to data logger malfunctions, the weather station did not record data for a total of 32 days during portions of September, October, and November 2003. Malfunctions at the water quality station resulted in missing salinity and depth data for May 1 to May 12, 2004 and June 27 to July 7, 2004. These missing data may have caused slightly biased averages, but were assumed not to affect general interpretations because they represented only a small portion of the high frequency dataset.

United States Geological Survey (USGS) discharge measurements for 2003–2004 taken at a station just upstream (<500 m) from the Pellicer Creek water quality station were downloaded from The months of August, September, and December of 2003 and January, April, October, November, and December of 2004 were incomplete and therefore not used in data analysis.

In addition to continuous in situ monitoring data, water was collected on a monthly basis using grab sampling methods at the Pellicer Creek site in accordance with NERR guidelines (Kennish, 2004; NOAA, 2004). One integrated sample from the entire column of water and two samples from 1 m above the bottom were collected with a Polyvinyl Chloride (PVC) pole (Venrick, 1978). In addition, an ISCO automatic sampler was deployed monthly to collect water samples from 1 m above the bottom every 2.5 hours for one complete tidal cycle (25 hours). Samples from both collection methods were brought back to the lab on ice and processed in a similar manner for the determination of nutrient concentrations, including total nitrogen (TN), total dissolved nitrogen (TDN), total phosphorus (TP), total dissolved phosphorus (TDP), nitrate + nitrite (NO3 + NO2), ammonium (NH4), and soluble reactive phosphorus (SRP), as well as water clarity parameters such as chlorophyll a, colored dissolved organic matter (CDOM), and turbidity.

Water Chemistry

Whole water samples were used to determine concentrations of TN, TP, chlorophyll a, and turbidity. To determine TN and TP, samples were digested and measured colorometrically on a Bran-Luebbe autoanalyzer (TN) and a dual-beam scanning spectrophotometer (TP) (APHA, 1998). Analysis methods for TDN and TDP were the same as those for TN and TP, respectively, but samples were filtered first through PALL A/E glass-fiber filters (1-μm pore size) (APHA, 1998). Chlorophyll a was processed using the Sartory and Grobbelaar (1984) ethanol extraction method. Chlorophyll a concentrations (uncorrected) were determined spectrophotometrically according to Standard Methods (APHA, 1998). Turbidity was measured in Nephelometric Turbidity Units using a LaMotte Model 2020 Turbidimeter.

For soluble nutrient and CDOM analyses, whole water was filtered through glass-fiber filters (1- μm pore size). Concentrations of NH4 and NO3 + NO2 were determined colorometrically on a Bran-Luebbe autoanalyzer (APHA, 1998; Strickland and Parsons, 1972). Concentrations of SRP were measured on a dual-beam scanning spectrophotometer at 882 nm following standard methods (APHA, 1998). Values of CDOM were measured against a platinum-cobalt standard using a dual-beam scanning spectrophotometer (APHA, 1998).

Data Analysis

Descriptive statistics of the 2003–2004 time series were used to explore inter- and intraannual variability in salinity, rainfall, wind speed, wind direction, nutrient concentrations, and water clarity. Sanger et al. (2002) and Wenner et al. (2004) measured storm effects by looking at the time it takes for a water quality parameter to return to prestorm conditions. They defined “prestorm” as 2 weeks before the storm and “poststorm” as 4 weeks after. A similar approach was taken in this study. Descriptive statistics were calculated for salinity values before and after storms and for summer seasons of both years. Summer was defined as the period between June and September of 2003 and 2004. The SAS statistical software version 8 (SAS Institute, Inc., 1999) was used to calculate all test statistics (α = 0.05).

Nonparametric tests were used after determining that most parameters were not normally distributed according to the Kolmogorov-Smirnov test for goodness-of-fit. The Wilcox-on sign-rank test was used to explore interannual variation by testing for differences between median values during the summer of 2003 and summer of 2004. The Spearman's Rank Correlation was used to examine the relationship between the rate of freshwater flow and salinity. Discharge and salinity were strongly negatively correlated (ρ = −0.95, p <0.0001). Because the discharge dataset was sparse, correlations were calculated between salinity, rainfall, and all other parameters to determine the associations between water column conditions and the interchange of marine and freshwater masses. Simple regressions were performed to evaluate the dependency of environmental parameters on salinity. Multiple regressions were not useful because of signs of multicollinearity.

Graphs of 30-minute salinity and water depth data for each diel ISCO sampling event were used to determine the times at which high and low tides occurred. Water samples and in situ observations from those times were used in correlation and regression calculations. The rainfall data used in the correlation analysis, obtained from the GTMNERR weather station, consisted of daily totals averaged for the intervals between monthly ISCO sampling events. September, October, and November 2003 rainfall averages were not used in correlation calculations because of data logger malfunctions during those months.


Basic Physiochemical Conditions

The bimodal rainfall pattern described for the study area by Jordan (1984) and Chen and Gerber (1990) was reflected in the elevated precipitation in summer and late winter/early spring of both 2003 and 2004 (Figure 4). In 2003, the latter half of July and the first half of August experienced frequent days of moderate rainfall, with a corresponding drop in daily salinity averages at the end of July. By August 2, mean daily salinity was below 1 ppt. Salinity values did not rise above 1 ppt until August 25, after which salinity increased until a rainfall event occurred on October 7 and salinity values temporarily declined. Wind speed during this period was minimal. By contrast, the tropical storms of 2004 caused several instances of elevated mean daily wind speed, and corresponding spikes were observed in salinity and water depth (Figure 4). Several days of rain from July 15–20, 2004, resulted in a daily mean salinity below 1 ppt on July 19. Afterward, salinity values increased to approximately 15 ppt until the passage of the first tropical storm. The USGS monthly discharge data for 2003–2004 indicated that September 2004 experienced the highest mean discharge of all months for which data were available (Figure 5). Salinity, precipitation, water depth, and wind speed were all more variable during the summer of 2004 than during the summer of 2003 (Figure 4).

Median chlorophyll a concentrations were similar for both years of the study (Table 1). Median TN and TDN concentrations were higher in 2004 than 2003. Much of the TN was in soluble form (TDN) in both years. Median NO3 + NO2 and NH4 represented roughly 10% of TDN. The NH4 levels were about four times as high as NO3 + NO2 concentrations in both years. Median TP, TDP, and SRP concentrations were lower in 2004 than 2003, and more than 70% of TP was in the soluble form (TDP). Roughly 50% of mean TDP was present as SRP.

Salinity and Meteorological Conditions during Storm Events

The salinity and weather time series available for Pellicer Creek made it possible to examine immediate storm effects. During the passage of Hurricane Charley, the Pellicer Creek weather station recorded a maximum wind speed of 17.9 m/s and a total of 129.6 mm of rainfall (Table 2). The weather station was experiencing southeasterly winds before the passage of Charley, but shifted to the north-northeast and increased in strength during the storm. North-northeasterly winds blew up-creek, coinciding with a temporary increase in salinity. During the 2 weeks before Charley, mean daily salinity was 12.2 ppt. After the storm passed, heavy rainfall in the Pellicer Creek catchment caused salinities in the creek to drop to less than 1 ppt (Figure 6). Salinity stayed below 1 ppt for 11 days, after which it began to rise, showing more typical tidal oscillations. Salinity remained below the 2-week pre-storm average for a total of 20 days until the next storm event on September 5.

During Hurricane Frances, maximum wind speed at Pellicer Creek was 18.8 m/s and the storm dropped 204.5 mm of rain at the Pellicer Creek weather station (Table 2). Average salinity 2 weeks before the storm was 1.8 ppt due to freshwater inputs caused by Charley. On September 5, salinity peaked at 29.8 ppt during a wind-induced saltwater surge up the creek and then quickly dropped to zero after passage of the storm (Figure 6). Water remained fresh and exhibited little salinity variability until September 20, when salinities were elevated due to northeast winds associated with Tropical Storm Ivan (Figure 6). Salinities fell after the initial spike and remained under 1 ppt until September 25.

Salinity in Pellicer Creek increased to 8.9 ppt on September 26 during Hurricane Jeanne, then dropped on September 27 and remained below 1 ppt for 10 days (Figure 6). After Jeanne, salinity remained below the pre-Charley average of 12.2 ppt for more than 3 weeks (September 28–October 21). The passage of each storm resulted in spikes in water depth due to a combination of wind-induced movement of marine water upstream and subsequent peaks in freshwater outflow from the creek watershed due to high rainfall inputs to the Pellicer Creek catchment.

Temporal Trends in Water Column Characteristics

The influence of multiple storms in 2004 on salinity in Pellicer Creek was evident in a comparison of summer median values in 2003 (Table 3). Median salinity during the summer of 2003 was significantly higher than during the summer of 2004. Rainfall and wind speed medians were not different between years, but both showed larger ranges (higher maxima) during the summer of 2004 than during the summer of 2003.

Salinity was inversely related to rainfall amounts, reflecting the influence of discharge from the catchment (Figure 7; Table 4). Water clarity parameters (i.e., CDOM and turbidity) showed bimodal distributions, coinciding with the typical late winter and summer rainfall peaks described for this region by Jordan (1984) and Chen and Gerber (1990) (Figure 7). The fact that CDOM levels were closely tied to watershed discharge was reflected in the strong inverse correlation between salinity and CDOM and the strong positive relationship between rainfall and CDOM (Table 4). Simple regression analysis showed that salinity accounted for 56% of the variability in CDOM (Table 5). Turbidity showed weak positive relationships (nonsignificant) with both salinity and rainfall, as might be expected from the complex suite of factors that contribute to turbidity, including wind-induced sediment re-suspension (Figure 7, Table 4).

Chlorophyll a also showed a bimodal distribution, peaking in the early summer and late fall of both years, before the peaks in rainfall (Figure 7). August and September of 2004 exhibited lower chlorophyll a concentrations and suppressed diel variability compared with the same months in 2003 (Figure 7). Phytoplankton levels (in terms of chlorophyll a concentration) entering Pellicer Creek from the estuary were typically higher than those entering the creek from the catchment. This pattern was illustrated by the positive correlation between chlorophyll a and salinity (Table 4). Regression analysis indicated that salinity was responsible for 47% of the variability in chlorophyll a concentrations (Table 5). This pattern was further illustrated by the observation that chlorophyll a concentrations were lower during the poststorm high flow period of September 14–15, 2004, than during the period of September 23–24, 2003, which was not preceded by tropical storms (Figure 8).

Nitrogen concentrations generally followed rainfall levels (Figure 7). Low salinities corresponded with high TN values, indicating that freshwater flow from the creek watershed is enriched with nitrogen. This observation is further supported by the fact that TN, TDN, and NO3 + NO2 were all negatively correlated with salinity (Table 4). Also, TN and TDN were significantly positively correlated with rainfall. Regression equations indicated that salinity was responsible for 50% of the variability in TN and TDN (Table 5). The difference in nitrogen concentrations between flow regimes is illustrated by comparing the diel sampling period of September 23–24, 2003, and the poststorm period of September 14–15, 2004 (Figure 8). The TN concentrations in the creek in the latter period were twofold higher than in the former period.

A similar contrast between high and presumably low watershed outflow periods was not observed for phosphorus. The TDP was inversely correlated with salinity and positively correlated with rainfall, whereas TP showed no significant relationships (Table 4). This pattern was exemplified by the similarity in TP concentrations observed during the diel sampling event of September 23–24, 2003, and the high flow period of September 14–15, 2004 (Figure 8). The proportion of particulate phosphorus in TP concentrations peaked in June and July of both years (Figure 9). In contrast, TP was 100% dissolved phosphorus during the diel sampling events of August and September 2004 (Figure 9). The disparity in the relationships of TP and TDP to salinity suggests that particulate phosphorus concentrations were controlled by factors other than freshwater discharge (i.e., wind-induced sediment resuspension).


The results of this study provide a glimpse into the time scale and magnitude of water quality changes in a tidal creek ecosystem after multiple storm events. One of the most direct and obvious effects of storms is the alteration of salinity regimes as it relates to wind-induced tidal surge and rainfall-induced freshwater runoff from the watershed. It is difficult to model the relationship between salinity and rainfall with simple regression analyses. To correctly define the dependence of salinity on rainfall and describe tropical storm effects, it is necessary to employ a transfer function model with Autoregressive Moving Average (ARMA) errors, time lags, and multiple explanatory variables (e.g., wind speed, wind direction, and tidal effects; Niu, Edmiston, and Bailey, 1998). However, for the applications in this study, it was not necessary to model rainfall effects on salinity, as salinity was a reasonable measure of the proportion of marine and freshwater inputs to the creek at a particular location and time. Model relationships developed from the limited discharge data available for the Pellicer Creek confirm the existence of a strong correlation between discharge and salinity. Based on the latter observation, salinity values were used throughout the study to explore the relative influence of fresh and salt water on water column characteristics in the creek.

Salinity was influenced by storm events in several significant ways. Strong onshore and upstream winds during storms drove high salinity marine water well up into the tidal creek, beyond the normal range associated with tidal exchange, although the effect was relatively short-lived. Conversely, the high rainfall totals associated with the tropical storms resulted in significantly increased discharges of freshwater from the watershed, which drove down salinities for extended periods of time, overwhelming the normal diel tidal fluctuations. These findings are comparable to those of Sanger et al. (2002) and Wenner et al. (2004) who concluded that surges associated with storms that approached the east coast from the Atlantic generally cause short-term increases in salinity, while rainfall inputs result in longer-term salinity decreases. The observations in this study also support those of Steward et al. (2006) who studied the effects of the 2004 Atlantic hurricane season on seagrass communities of Florida's east coast. They measured water quality conditions in the Indian River Lagoon approximately weekly during and after the 2004 hurricane season and observed salinity and water clarity decreases after the storms.

The relative magnitude of the storm effects on Pellicer Creek are exemplified by comparing salinity regimes in the summer periods of 2003 and 2004. Median salinities for 2004 were nearly threefold lower than 2003, although both years had extended periods of low salinity. In 2003, the period was restricted to August, the month when rainfall totals peaked. In 2004, multiple extended periods of low salinity were observed from July through early October, in large part due to the multiple storm events in the late summer.

The timing of the tropical storm season has important ecological implications (Mallin and Corbett, 2006; Michener et al., 1997). Typically, the ground is not saturated in early fall and soils are dry enough to absorb initial rain events and recharge wetlands. However, excessive rainfall associated with storms can change the magnitude of river discharge, thereby altering the normal seasonal variability in salinity and nutrient concentrations.

In conjunction with alterations in salinity regimes, storm events affected other water column characteristics in Pellicer Creek related to the substantial rainfall inputs to the watershed. Storm-induced flooding of watersheds can alter nutrient inputs to estuaries (Paerl et al., 2006) and impact the potential for hypoxia (Tomasko, Anastasiou, and Kovach, 2006). As seen from the current data, nitrogen concentrations in Pellicer Creek outflows were significantly elevated after major rainfall events. This effect is illustrated by the more than twofold higher TN concentration at the Pellicer Creek site on comparable dates in 2003 and poststorm 2004. The multiple storm impact is also reflected in the higher median concentrations of both TN and TDN in 2004 compared to 2003. The strong negative correlation between salinity and TN further demonstrates the role of the creek watershed in nitrogen load to the estuary. Conversely, the lack of significant relationships between TP or SRP and rainfall or salinity suggests a lesser role for the creek as a source of phosphorus to the estuary.

Because primary production in many coastal environments along the east coast of Florida (Phlips et al., 2002), as well as around the world (Cloern, 2001), are nitrogen limited, increases in nitrogen input represent a potential for enhanced algal production and biomass (Mallin and Corbett, 2006; Paerl et al., 2006). Recent research in the Suwannee River estuary has shown that increases in nutrient load during flood periods can be correlated to increases in the standing crop of phytoplankton in the nearshore environment (Bledsoe et al., 2004; Bledsoe and Phlips, 2000). Even episodic nutrient pulses associated with storm events may have important consequences for estuarine ecology, such as the stimulation of harmful algal blooms. For example, it has been hypothesized that a major bloom of the toxic dinoflagellate Pyrodinium bahamense in the Indian River Lagoon of Florida may have been in part stimulated by the pulse in nutrient load resulting from the high rainfall El Niño event of 2001/2002 (Phlips et al., 2004a, 2004b, 2006). Increases in nutrient concentrations and organic matter can also have consequences for higher trophic levels. In Chesapeake Bay, annual differences in biomass and the abundance of macroinvertebrates were attributed to differences in rainfall patterns and consequent differences in nutrient inputs between years (Bilkovic et al., 2006).

The limited spatial extent of this study prevents any conclusions about how far nutrients are dispersed from the Pellicer Creek catchment into the adjacent estuary or the ultimate fate of the bioavailable forms of nutrients. The relative importance of dilution and transport out of the system in this highly flushed estuary may also contribute to the impact of nutrient loads on the response of primary production (Phlips et al., 2004b). However, at the Pellicer Creek water quality station, patterns in chlorophyll a concentrations, and the associations between chlorophyll a and salinity, indicate that phytoplankton biomass entering Pellicer Creek from the estuary on the flood tide is greater than the biomass coming out from the creek watershed. This observation is analogous to that observed in other similar ecosystems, such as the Suwannee River estuary, where phytoplankton biomass in the river is depressed by low light availability and short water residence time, and biomass peaks just outside the river mouth where the high nutrients from the river mix with clear oceanic water (Bledsoe et al., 2004; Bledsoe and Phlips, 2000).

The pronounced deviations in salinity regimes caused by storm events can result in changes in the structure and function of tidal creeks and estuaries (Frazer et al., 2006; Paerl et al., 2001, 2006; Switzer et al., 2006). For example, in the estuarine environment associated with Pellicer Creek, commercial blood ark (Anadara ovalis), ponderous ark (Noetia ponderosa), and hard clam (Mercenaria mercenaria) populations were devastated by the multiple freshwater pulses associated with the 2004 storms and did not begin to recover in the region until 2006 (Nuñez, personal communication). Situations of freshwater pulses affecting benthic community structure were also documented at Cape Fear, North Carolina, where six hurricanes in 4 years resulted in shifts of benthic infauna from marine to freshwater community structure (Mallin et al., 1999).

Estuarine communities are not only affected by changes in the relative magnitude of freshwater inputs, but also by the rate of change. In a Florida Bay study, Montague and Ley (1993) found that stations with large salinity fluctuations had relatively low biomass and unstable species composition. Ridler, Dent, and Arrington (2006) observed negative posthurricane impacts on seagrass communities in Southeast Florida during the 2004 storm season. By contrast, recent observations of seagrass communities in the Indian River Lagoon before and after the 2004 storm season show considerable resilience to the events (Steward et al., 2006). Similarly, research on the response of fish communities to the 2004 storms showed an initial short-term impact but rapid recovery (Paperno et al., 2006; Stevens, Blewett, and Casey, 2006). The results of the present study indicate that the tropical storms of 2004 greatly accelerated the rates and magnitudes of change in salinity within Pellicer Creek, but the effects on the benthic plant or fish communities remain undocumented.


The meteorological and water quality monitoring networks established by the NERR program provide opportunities to examine ecological processes that are stochastic in nature and subject to considerable short-term variability, such as storm events. The results of this study at the GTMNERR reveal that multiple storm events can have significant impacts on a range of key water quality parameters that play major roles in defining ecosystem structure and function. Recent advances in the development of hydrodynamic models of estuarine circulation, such as those recently described for the GTMNERR (Sheng et al., this issue), will ultimately provide an opportunity to predict how future storm patterns will affect the spatial and temporal distribution of water quality parameters. The NERR monitoring networks will provide essential data for calibration and validation of these models. Given the major changes in watershed characteristics and global climate patterns expected in future years, the ability to predict the impacts of these changes will be an essential part of designing, implementing, and justifying management efforts.


We would like to acknowledge all of the past and present Phlips Lab employees who helped with the sampling, processing, and chemistry needed to obtain the data for this project. In particular, we would like to thank Jean Lockwood, Ken Black, Lance Riley, and Katie O'Donnell. This project was made possible by the presence of NOAA's NERR monitoring program and was in part funded by an award from the Estuarine Reserves Division, Office of Ocean and Coastal Resource Management, National Ocean Service, NOAA.



APHA (American Public Health Association) 1998. Standard Methods for the Examination of Water and Wastewater, 20th Edition. Baltimore, Maryland: United Book Press, Inc. Google Scholar


D. C. Berquist, J. A. Hale, P. Baker, and S. M. Baker . 2006. Development of ecosystem indicators for the Suwannee River Estuary: oyster reef habitat quality along a salinity gradient. Estuaries and Coasts 29 3:353–360. Google Scholar


J. L. Beven 2004. Tropical cyclone report: Hurricane Frances. (accessed September 19, 2006). Google Scholar


D. M. Bilkovic, M. Roggero, C. H. Hershner, and K. H. Havens . 2006. Influence of land use on macrobenthic communities in nearshore estuarine habitats. Estuaries and Coasts 29:6B. 1185–1195. Google Scholar


E. Bledsoe and E. J. Phlips . 2000. Nutrient versus light limitation of phytoplankton in the Suwannee River and Estuary. Estuaries 23:458–473. Google Scholar


E. Bledsoe, E. J. Phlips, C. E. Jett, and K. A. Donnelly . 2004. The relationships among phytoplankton biomass, nutrient loading and hydrodynamics in an inner-shelf estuary, the Suwannee River estuary, Florida, USA. Ophelia 58:29–47. Google Scholar


E. R. Blood, P. Anderson, P. A. Smith, C. Nybro, and K. A. Ginsberg . 1991. Effects of Hurricane Hugo on coastal soil solution chemistry in South Carolina. Biotropica 23:4a. 348–355. Google Scholar


D. F. Boesch, R. J. Diaz, and R. W. Virnstein . 1976. Effects of Tropical Storm Agnes on soft-bottom macrobenthic communities of the James and York Estuaries and the Lower Chesapeake Bay. Chesapeake Science 17 4:246–259. Google Scholar


D. R. Cahoon 2006. A review of major storm impacts on coastal wetland elevations. Estuaries and Coasts 29:6A. 889–898. Google Scholar


E. Chen and J. F. Gerber . 1990. Climate. In: Meyers, R.L. and Ewel, J.J. (eds.), Ecosystems of Florida. Orlando, Florida: University of Central Florida Press, pp. 11–34. Google Scholar


J. E. Cloern 2001. Our evolving conceptual model of the coastal eutrophication problem. Marine Ecology Progress Series 210:223–253. Google Scholar


J. E. Cloern, B. E. Cole, R. L. Wong, and A. E. Alpine . 1985. Temporal dynamics of estuarine phytoplankton: a case study of San Francisco Bay. Hydrobiologia 129:153–176. Google Scholar


D. B. Enfield, A. M. Mestas-Nuñez, and P. J. Trimble . 2001. The Atlantic multidecadal oscillation and its relation to rainfall and river flows in the continental U.S. Geophysical Research Letters 28 10:2077–2080. Google Scholar


B. D. Eyre and D. Pont . 2003. Intra- and inter-annual variability in the different forms of diffuse nitrogen and phosphorus delivered to seven sub-tropical east Australian estuaries. Estuarine, Coastal and Shelf Science 57:137–148. Google Scholar


T. K. Frazer, S. K. Notestein, C. A. Jacoby, C. Jones Littles, S. R. Keller, and R. A. Swett . 2006. Effects of storm-induced salinity changes on submersed aquatic vegetation in Kings Bay, Florida. Estuaries and Coasts 29:6A. 943–953. Google Scholar


D. Gerten and R. Adrian . 2000. Climate-driven changes in spring phytoplankton dynamics and the sensitivity of shallow polymictic lakes to the North Atlantic Oscillation. Limnology and Oceanography 45 5:1058–1066. Google Scholar


D. M. Goodrich 1988. On meteorologically induced flushing in three U.S. east coast estuaries. Estuarine, Coastal and Shelf Science 26:111–121. Google Scholar


H. Greening, P. Doering, and C. Corbett . 2006. Hurricane impacts on coastal ecosystems. Estuaries and Coasts 29:6A. 877–879. Google Scholar


C. T. Hackney, W. D. Burbanck, and O. P. Hackney . 1976. Biological and physical dynamics of a Georgia tidal creek. Chesapeake Science 17 4:271–280. Google Scholar


J. Hama and N. Handa . 1994. Variability of the biomass, chemical composition and productivity of phytoplankton in Kinuura Bay, Japan during the rainy season. Estuarine, Coastal and Shelf Science 39:497–509. Google Scholar


J. Imberger, T. Berman, R. R. Christian, E. B. Sherr, D. E. Whitney, L. R. Pomeroy, R. G. Wiegert, and W. J. Wiebe . 1983. The influence of water motion on the distribution and transport of materials in a salt marsh estuary. Limnology and Oceanography 28 2:201–214. Google Scholar


C. L. Jordan 1984. Florida's weather and climate: implications for water. In: Fernald, E.A. and Patton, D.J. (eds.), Water Resources Atlas of Florida. Tallahassee, Florida: Florida State University, pp. 18–35. Google Scholar


M. J. Kennish 2004. NERRS research and monitoring initiatives. Journal of Coastal Research 45:1–8. Google Scholar


M. B. Lawrence and H. D. Cobb . 2004. Tropical cyclone report: Hurricane Jeanne. (accessed September 19, 2006). Google Scholar


R. J. Livingston, X. Niu, F. G. Lewis III, and G. C. Woodsum . 1997. Freshwater input to a gulf estuary: long-term control of trophic organization. Ecological Applications 7 1:277–299. Google Scholar


M. A. Mallin and C. A. Corbett . 2006. How hurricane attributes determine the extent of environmental effects: multiple hurricanes and different coastal systems. Estuaries and Coasts 29:6A. 1046–1061. Google Scholar


M. A. Mallin, M. H. Posey, G. C. Shank, M. R. McIver, S. H. Ensign, and T. D. Alphin . 1999. Hurricane effects on water quality and benthos in the cape fear watershed: natural and anthropogenic impacts. Ecological Applications 9 1:350–362. Google Scholar


C. A. McMillan 1974. Salt tolerance of mangroves and submerged aquatic plants. In: Reimold, R.J. and Queen, W.H. (eds.), The Ecology of Halophytes. New York: Academic Press, pp. 379–390. Google Scholar


W. K. Michener, E. R. Blood, K. L. Bildstein, M. M. Brinson, and L. R. Gardner . 1997. Climate change, hurricanes and tropical storms, and rising sea level in coastal wetlands. Ecological Applications 7 3:770–801. Google Scholar


C. L. Montague and J. A. Ley . 1993. A possible effect of salinity fluctuation on abundance of benthic vegetation and associated fauna in northeastern Florida Bay. Estuaries 16 4:703–717. Google Scholar


S. K. Moore, M. E. Baird, and I. M. Suthers . 2006. Relative effects of physical and biological processes on nutrient and phytoplankton dynamics in a shallow estuary after a storm event. Estuaries and Coasts 29 1:81–95. Google Scholar


National Oceanic and Atmospheric Administration Staff 1994–2006. Climatological Data Annual Summary, ISSN 0145-0484. Asheville, North Carolina: National Oceanic and Atmospheric Administration, National Climatic Data Center, (13), 98–110. Google Scholar


National Oceanic and Atmospheric Administration, Office of Ocean and Coastal Resource Management, National Estuarine Research Reserve System-Wide Monitoring Program Data Set 2004. Centralized Data Management Office, Baruch Marine Field Lab, University of South Carolina: Scholar


X-F. Niu, H. L. Edmiston, and G. O. Bailey . 1998. Time series models for salinity and other environmental factors in the Apalachicola estuarine system. Estuarine, Coastal and Shelf Science 46:549–563. Google Scholar


S. W. Nixon 1995. Coastal marine eutrophication: a definition, social consequences and future concerns. Ophelia 41:199–220. Google Scholar


S. P. Orlando Jr., P. H. Wendt, C. J. Klein, M. E. Pattillo, K. C. Dennis, and G. H. Ward . 1994. Salinity Characteristics of South Atlantic estuaries. Silver Spring, Maryland: National Oceanic and Atmospheric Administration, Office of Ocean Resources Conservation and Assessment. Google Scholar


H. W. Paerl, J. D. Bales, L. W. Ausley, C. P. Buzzelli, L. B. Crowder, L. A. Eby, J. M. Fear, M. Go, B. L. Peierls, T. L. Richardson, and J. S. Ramus . 2001. Ecosystem impacts of three sequential hurricanes (Dennis, Floyd, and Irene) on the United States' largest lagoonal estuary, Pamlico Sound, NC. Ecology 98 10:5655–5660. Google Scholar


H. W. Paerl, L. M. Valdes, A. R. Joyner, B. L. Peierls, M. F. Piehler, S. R. Riggs, R. R. Christian, L. A. Eby, L. B. Crowder, J. S. Ramus, E. J. Clesceri, C. P. Buzzelli, and R. A. Luettich Jr. . 2006. Ecological response to hurricane events in the Pamlico Sound System, North Carolina and implications for assessment and management in a regime of increased frequency. Estuaries and Coasts 29:6A. 1033–1045. Google Scholar


R. Paperno, D. M. Tremain, D. H. Adams, A. P. Sebastian, J. T. Sauer, and J. Dutka-Gianelli . 2006. The disruption and recovery of fish communities in the Indian River Lagoon, Florida following two hurricanes in 2004. Estuaries and Coasts 29:6A. 1004–1010. Google Scholar


R. J. Pasch, D. P. Brown, and E. S. Blake . 2004. Tropical cyclone report: Hurricane Charley. (accessed September 19, 2006). Google Scholar


E. J. Phlips, S. Badylak, E. L. Bledsoe, and M. Cichra . 2006. Factors influencing the distribution and abundance of Pyrodinium bahamense in coastal ecosystems of Florida. Marine Ecology Progress Series 322:99–115. Google Scholar


E. J. Phlips, S. Badylak, and T. Grosskopf . 2002. Factors affecting the abundance of phytoplankton in a restricted subtropical lagoon, the Indian River Lagoon, Florida, USA. Estuarine, Coastal and Shelf Science 55:385–402. Google Scholar


E. J. Phlips, S. Badylak, S. Youn, and K. Kelley . 2004a. The occurrence of potentially toxic dinoflagellates and diatoms in a subtropical lagoon, the Indian River Lagoon, Florida, USA. Harmful Algae 3:39–49. Google Scholar


E. J. Phlips, N. Love, S. Badylak, P. Hansen, C. V. John, and R. Gleeson . 2004b. A comparison of water quality and hydrodynamic characteristics of the Guana Tolomato Matanzas National Estuarine Research Reserve and the Indian River Lagoon in Florida. Journal of Coastal Research Special Issue No. 45:93–109. Google Scholar


M. S. Ridler, R. C. Dent, and D. A. Arrington . 2006. Effects of two hurricanes on Syringodium filiforme, manatee grass, within the Loxahatchee River Estuary, southeast Florida. Estuaries and Coasts 29:6A. 1019–1025. Google Scholar


D. M. Sanger, M. D. Arendt, Y. Chen, E. L. Wenner, A. F. Holland, D. Edwards, and J. Caffrey . 2002. A synthesis of water quality data: National Estuarine Research Reserve System-Wide Monitoring Program (1995–2000). Technical Report Series 2002: 3, No. 500. Silver Spring, Maryland: U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Ocean Service, Office of Ocean and Coastal Resource Management. 135. p. Google Scholar


D. P. Sartory and J. U. Grobbelaar . 1984. Extraction of chlorophyll a from freshwater phytoplankton for spectrophotometric analysis. Hydrobiologia 114:177–187. Google Scholar


SAS Institute, Inc 1999. SAS Version 8. Cary, North Carolina: SAS Institute, Inc. Google Scholar


D. Scavia, J. C. Field, D. F. Boesch, R. W. Buddemeier, V. Burkett, D. R. Cayan, M. Fogarty, M. A. Harwell, R. W. Howarth, C. Mason, D. J. Reed, T. C. Royer, A. H. Sallenger, and J. G. Titus . 2002. Climate change impacts on U.S. coastal and marine systems. Estuaries 25 2:149–164. Google Scholar


U. M. Sharler and D. Baird . 2003. The influence of catchment management on salinity, nutrient stochiometry and phytoplankton biomass of Eastern Cape estuaries, South Africa. Estuarine, Coastal and Shelf Science 56:735–748. Google Scholar


P. W. Stevens, D. A. Blewett, and J. P. Casey . 2006. Short-term effects of a low dissolved oxygen event on estuarine fish assemblages following the passage of Hurricane Charley. Estuaries and Coasts 29:6A. 997–1003. Google Scholar


J. S. Steward, R. W. Virnstein, M. A. Lasi, L. J. Morris, J. D. Miller, L. M. Hall, and W. A. Tweedale . 2006. The impacts of the 2004 hurricanes on hydrology, water quality, and seagrass in the central Indian River Lagoon, Florida. Estuaries and Coasts 29:6A. 954–965. Google Scholar


S. R. Stewart 2004. Tropical cyclone report: Hurricane Ivan. (accessed September 19, 2006). Google Scholar


J. D. H. Strickland and T. R. Parsons . 1972. A Practical Handbook of Seawater Analysis: Determination of Ammonia (Oxidation Method). Ottawa, Canada: Fisheries Research Board of Canada. Google Scholar


T. S. Switzer, B. L. Winner, N. M. Dunham, J. A. Whittington, and M. Thomas . 2006. Influence of sequential hurricanes on nekton communities in a southeast Florida estuary: short-term effects in the context of historical variations in freshwater inflow. Estuaries and Coasts 29:6A. 1011–1018. Google Scholar


D. C. Tabb and A. C. Jones . 1962. Effect of Hurricane Donna on aquatic fauna of North Florida Bay. Transactions of the American Fisheries Society 91:375–378. Google Scholar


J. T. Tilmant, R. W. Curry, R. Jones, A. Szmant, J. C. Zieman, M. Flora, M. B. Robblee, D. Smith, R. W. Snow, and H. Wanless . 1994. Hurricane Andrew's effects on marine resources. BioScience 44 4:230–237. Google Scholar


D. A. Tomasko, C. Anastasiou, and C. Kovach . 2006. Dissolved oxygen dynamics in Charlotte Harbor and its contributing watershed, in response to Hurricanes Charley, Frances, and Jeanne—impacts and recovery. Estuaries and Coasts 29:6A. 932–938. Google Scholar


R. F. Van Dolah and G. S. Anderson . 1991. Effects of Hurricane Hugo on salinity and dissolved oxygen conditions in the Charleston Harbor Estuary. Journal of Coastal Research 8:83–94. Google Scholar


E. L. Venrick 1978. Sampling Techniques. In: Sournia, A. (ed.), Phytoplankton Manual. Paris, France: United Nations Educational, Scientific and Cultural Organization, pp. 33–67. Google Scholar


N. D. Walker 2001. Tropical storm and hurricane wind effects on water level, salinity, and sediment transport in the river-influenced Atchafalaya-Vermilion Bay system, Louisiana, USA. Estuaries 24 4:498–508. Google Scholar


E. Wenner, D. Sanger, M. Arendt, A. F. Holland, and Y. Chen . 2004. Variability in dissolved oxygen and other water-quality variables within the National Estuarine Research Reserve System. Journal of Coastal Research Special Issue No. 45:17–38. Google Scholar


Figure 1.

Location map of Pellicer Creek and Guana Tolomato Matanzas National Estuarine Research Reserve monitoring stations.


Figure 2.

Storm tracks (Charley, France, Ivan, and Jeanne).


Figure 3.

Annual rainfall totals, Hastings, Florida (NOAA, 2006).


Figure 4.

Daily averages for entire time series (2003–2004). (A) Salinity; (B) total precipitation; (C) water depth; (D) wind speed.


Figure 5.

The USGS mean monthly discharge in Pellicer Creek (*data not collected).


Figure 6.

Mean daily salinity at the Pellicer Creek water quality station 2 weeks before Hurricane Charley and 4 weeks after Hurricane Jeanne.


Figure 7.

The 2003–2004 mean monthly rainfall (represented by gray bars) collected from the Pellicer Creek weather station, plotted relative to maximum (dashed top line) and minimum (solid bottom line) monthly nutrient concentrations, phytoplankton biomass, and water clarity parameters determined from diel sampling. (A) Salinity →?(F) TP.


Figure 8.

Diel patterns in salinity, chlorophyll, and nutrient concentrations observed during the period of September 23–24 of 2003 (A) and the poststorm high flow period of September 14–15, 2004 (B).


Figure 9.

Mean monthly particulate phosphorus (TP-TDP) concentrations.


Table 1.

Median and [range] of chlorophyll and nutrient values obtained from monthly grab and diel sampling in Pellicer Creek during 2003 and 2004 as well as p values obtained from Wilcoxon signed rank tests between years (H0: 2003 = 2004, N = 166)


Table 2.

Total rainfall and maximum wind speed experienced at Pellicer Creek during the passage of the 2004 tropical systems


Table 3.

Median and range (based on daily averages) describing the distribution of salinity, daily rainfall totals, and wind speed for the 2003–2004 wet seasons and p values obtained from Wilcoxon signed rank tests between summer seasons (H0: summer 2003 = summer 2004)


Table 4.

Spearman correlation coefficients between rainfall, salinity, and all other parameters


Table 5.

Simple regression coefficients of salinity vs. all other parameters at one high and one low tide per month (N = 48)

Nicole G. Dix, Edward J. Phlips, and Richard A. Gleeson "Water Quality Changes in the Guana Tolomato Matanzas National Estuarine Research Reserve, Florida, Associated with Four Tropical Storms," Journal of Coastal Research 2008(10055), 26-37, (1 July 2008).
Published: 1 July 2008

Back to Top