Predicting offshore temperatures in Monterey Bay based on coastal observations using linear forecast models

Abstract The linear systems approach is used to forecast offshore near-surface and subsurface temperatures in Monterey Bay and further offshore based on coastal sea surface temperatures (SSTs) at Pacific Grove. SST from Pacific Grove provided the input to the system and the forecast parameters or outputs were temperature at 1 m and 100 m at the M1 buoy located 20 km from Pacific Grove near the center of Monterey Bay, and temperature at 1 m at the M2 buoy located 55 km from Pacific Grove. To forecast temperatures at the M1 and M2 buoys, Box-Jenkins, State-Space, ARX, and ARMAX models were employed. Model formulation, implementation, forecasting procedures, and methods of evaluation are presented. Seven and 30-day forecasts were routinely made for the daily observations although other forecast horizons were employed. For all models and variables, RMS differences between the forecasts and the observations increased rapidly between 1 and 15 days. Beyond about 30 days, RMS differences tended to remain almost constant with increasing forecast horizon. Overall, model forecasts were best for temperature at 100 m at the M1 buoy, due to the fact that temperature is well conserved at depth. Differences in performance between the models were small but the ARMAX model often produced forecasts that were slightly better than the rest, a result that we attribute to a more complete specification of the noise. Although the Box-Jenkins and State-Space models have the potential to produce better forecasts, because more terms must be specified to implement them, the opportunity to produce less-than-optimal results is also greater. Finally, because of seasonal changes in the circulation of Monterey Bay, it is possible that causality was violated, upon occasion, placing certain constraints on the results. Models based on the linear systems approach, where they can be implemented, could serve as a useful adjunct to hydrodynamic ocean circulation models by providing additional information for model initialization, evaluation, and data assimilation. Using the same approach, operational forecasts of the coastal circulation could be made by including forecast winds and the predicted tides as inputs, and CODAR-observed surface currents as the output. In a less glamorous but still useful role, they could be used to fill significant gaps in offshore records where data continuity and quality are important.

[1]  H. Lilliefors On the Kolmogorov-Smirnov Test for Normality with Mean and Variance Unknown , 1967 .

[2]  R. E. Kalman,et al.  On the general theory of control systems , 1959 .

[3]  Joseph D. Petruccelli,et al.  An Automatic Procedure for Identification, Estimation and Forecasting Univariate Self Exiting Threshold Autoregressive Models , 1988 .

[4]  John G. Proakis,et al.  Digital Signal Processing: Principles, Algorithms, and Applications , 1992 .

[5]  Lennart Ljung,et al.  System Identification: Theory for the User , 1987 .

[6]  H. Saunders,et al.  Linear Systems and Digital Signal Processing , 1985 .

[7]  Johannes Ledolter,et al.  Statistical methods for forecasting , 1983 .

[8]  L. C. Breaker,et al.  The circulation of Monterey Bay and related processes , 1989 .

[9]  Chris Chatfield,et al.  The Analysis of Time Series: An Introduction , 1981 .

[10]  P. Strub,et al.  Atmospheric conditions during the spring and fall transitions in the coastal ocean off western United States , 1988 .

[11]  Barbara M. Hickey,et al.  The California current system—hypotheses and facts☆ , 1979 .

[12]  L. Breaker What's happening in Monterey Bay on seasonal to interdecadal time scales , 2005 .

[13]  J. Bendat,et al.  Random Data: Analysis and Measurement Procedures , 1971 .

[14]  Laurence C. Breaker,et al.  Development of a Real-Time Regional Ocean Forecast System with Application to a Domain off the U.S. East Coast , 2004 .

[15]  L. Breaker Nonlinear aspects of sea surface temperature in Monterey Bay , 2006 .

[16]  M. A. Yoder,et al.  Dsp First: A Multimedia Approach , 1997 .

[17]  Andrew Harvey,et al.  Forecasting, Structural Time Series Models and the Kalman Filter , 1990 .

[18]  W. Cleveland Robust Locally Weighted Regression and Smoothing Scatterplots , 1979 .

[19]  R. Katz Use of cross correlations in the search for teleconnections , 1988 .

[20]  W. J. Conover,et al.  Practical Nonparametric Statistics , 1972 .

[21]  W. Brass,et al.  Perspectives in Population Prediction: Illustrated by the Statistics of England and Wales , 1974 .

[22]  Jeffrey D. Paduan,et al.  Remotely sensed surface currents in Monterey Bay from shore-based HF radar (Coastal Ocean Dynamics Application Radar) , 1996 .

[23]  H. Akaike A new look at the statistical model identification , 1974 .

[24]  Franklin B. Schwing,et al.  Bifurcated flow from an upwelling center: a cold water source for Monterey Bay , 1994 .

[25]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1972 .

[26]  P. Lewis,et al.  Analysis of a 12-year record of sea-surface temperatures off Pt. Sur, California , 1983 .

[27]  J. Sprott Chaos and time-series analysis , 2001 .

[28]  Mark D. Semon,et al.  POSTUSE REVIEW: An Introduction to Error Analysis: The Study of Uncertainties in Physical Measurements , 1982 .

[29]  Thomas Kailath,et al.  Linear Systems , 1980 .

[31]  William S. Cleveland The elements of graphing data , 1980 .