Water level variations in the Neuse and Pamlico Estuaries, North Carolina due to local and remote forcing

Abstract

Water level time series records from the Neuse and Pamlico River Estuaries were statistically compared to local and distant wind field data, water level records within the Pamlico Sound and also coastal ocean sites to determine the relative contribution of each time series to water levels in the Neuse and Pamlico Estuaries. The objectives of this study were to examine these time series data using various statistical methods (i.e. autoregressive, empirical orthogonal function analysis (EOF), exploratory data analysis (EDA)) to determine short- and long-time-scale variability, and to develop predictive statistical models that can be used to estimate past water level fluctuations in both the Neuse Estuary (NE) and Pamlico Estuary (PE). Short- and long-time-scale similarities were observed in all time series of estuarine, Pamlico Sound and subtidal coastal ocean water level and wind component data, due to events (nor'easters, fronts and tropical systems) and seasonality. Empirical orthogonal function analyses revealed a strong coastal ocean and wind field contribution to water level in the NE and PE. Approximately 95% of the variation was captured in the first two EOF components for water level data from the NE, sound and coastal ocean, and 70% for the PE, sound and coastal ocean. Spectral density plots revealed strong diurnal signals in both wind and water level data, and a strong cross correlation and coherency between the NE water level and the North/South wind component. There was good agreement between data and predictions using autoregressive statistical models for the NE ( R 2 = 0.92) and PE ( R 2 = 0.76). These methods also revealed significant autoregressive lags for the NE (days 1 and 3) and for the PE (days 1, 2 and 3). Significant departures from predictions are attributed to local meteorological and hydrological events. The autoregressive techniques showed significant predictive improvement over ordinary least squares methods. The results are considered within the context of providing long time-scale hindcast data for the two estuaries, and the importance of these data for multidisciplinary researchers and managers.