In this paper an attempt is initiated to analyze long-term time series of wave data and to model them as a nonstationary stochastic process with yearly periodic mean value and standard deviation (periodically correlated or cyclostationary stochastic process). First, an analysis of annual mean values is performed in order to identify overyear trends. It turns out that it is very likely that an increasing trend is present in the examined hindcast data. The detrended time series Y(&tgr;) is then decomposed, using an appropriate seasonal standardization procedure, to a periodic mean value μ(&tgr;) and a residual time series W(&tgr;) multiplied by a periodic standard deviation &sgr;(&tgr;) of Y(&tgr;)=μ(&tgr;)+&sgr;(&tgr;)W(&tgr;). The periodic components μ(&tgr;) and &sgr;(&tgr;) are estimated and represented by means of low-order Fourier series, and the residual time series W(&tgr;) is examined for stationarity. For this purpose, spectral densities of W(&tgr;), obtained from different-season segments, are calculated and compared with each other. It is shown that W(&tgr;) can indeed be considered stationary, and thus Y(&tgr;) can be considered periodically correlated. This analysis has been applied to hindcast wave data from five locations in the North Atlantic Ocean. It turns out that the spectrum of W(&tgr;) is very weakly dependent on the site, a fact that might be useful for the geographic parameterization of wave climate. Finally applications of this modeling to simulation and extreme-value prediction are discussed. ¿ American Geophysical Union 1995