SWAN predictions of waves observed in shallow water onshore of complex bathymetry

SWAN model predictions, initialized with directional wave buoy observations in 550-m water depth offshore of a steep, submarine canyon, are compared with wave observations in 5.0-, 2.5-, and 1.0-m water depths. Although the model assumptions include small bottom slopes, the alongshore variations of the nearshore wave field caused by refraction over the steep canyon are predicted well over the 50 days of observations. For example, in 2.5-m water depth, the observed and predicted wave heights vary by up to a factor of 4 over about 1000 m alongshore, and wave directions vary by up to about 10°, sometimes changing from south to north of shore normal. Root-mean-square errors of the predicted wave heights, mean directions, periods, and radiation stresses (less than 0.13 m, 5°, 1 s, and 0.05 m³/s² respectively) are similar near and far from the canyon. Squared correlations between the observed and predicted wave heights usually are greater than 0.8 in all water depths. However, the correlations for mean directions and radiation stresses decrease with decreasing water depth as waves refract and become normally incident. Although mean wave properties observed in shallow water are predicted accurately, nonlinear energy transfers from near-resonant triads are not modeled well, and the observed and predicted wave energy spectra can differ significantly at frequencies greater than the spectral peak, especially for narrow-band swell.     

Modelling infragravity motions on a rip-channel beach

A non-linear shallow water wave model operating on the time-scale of wave groups is compared with measurements of infragravity motions on a rip-channel beach to verify the model concepts and assess the model performance. The measurements were obtained during the RIP-current EXperiment (RIPEX) in concert with the Steep Beach Experiment (SBE) performed at Sand City, Monterey Bay, CA, during the spring of 2001. The nearshore bathymetry was made up of shore-connected shoals incised by relatively narrow rip-channels spaced approximately 125 m apart. The comparison considers a 20-day period during which significant changes in both the offshore wave climate and nearshore bathymetry occurred. The temporal variation in infragravity conditions during the experiment is strong, with computational results typically explaining 70% to 80% of the observed infragravity motions within the nearshore. In contrast to the temporal variation, the alongshore spatial variation in infragravity intensity during the experiment is generally weak, even though the underlying bathymetry shows strong depth variations. Model computations suggest preferential coupling between the computed edge wave motions and the quasi-periodic bathymetry is present, a prerequisite for strong spatial variability. However, the infragravity field is dominated by cross-shore infragravity motions, which are only weakly coupled to the quasi-periodic bathymetry, resulting in a weak alongshore variability of the total infragravity motions.     

A ten-year data set for fetch- and depth-limited wave growth

This paper presents the key results from a ten-year data set for Lake IJssel and Lake Sloten in The Netherlands, containing information on wind, storm surges and waves, supplemented with SWAN 40.51 wave model results. The wind speeds U₁₀, effective fetches x and water depths d for the data set ranged from 0–24 m s^− 1, 0.8–25 km and 1.2–6 m respectively. For locations with non-sloping bottoms, the range in non-dimensional fetch x? ( = gxU₁₀^− 2) was about 25–80,000, while the range in dimensionless depth d? ( = g d U₁₀^− 2) was about 0.03–1.7. Land–water wind speed differences were much smaller than the roughness differences would suggest. Part of this seems due to thermal stability effects, which even play a role during near-gale force winds. For storm surges, a spectral response analysis showed that Lake IJssel has several resonant peaks at time scales of order 1 h. As for the waves, wave steepnesses and dimensionless wave heights H? ( = gH_m0U₁₀^− 2) agreed reasonably well with parametric growth curves, although there is no single curve to which the present data fit best for all cases. For strongly depth-limited waves, the extreme values of d? (0.03) and H_m0 / d (0.44) at the 1.7 m deep Lake Sloten were very close to the extremes found in Lake George, Australia. For the 5 m deep Lake IJssel, values of H_m0 / d were higher than the depth-limited asymptotes of parametric wave growth curves. The wave model test cases of this study demonstrated that SWAN underestimates H_m0 for depth-limited waves and that spectral details (enhanced peak, secondary humps) were not well reproduced from H_m0 / d = 0.2–0.3 on. SWAN also underestimated the quick wave response (within 0.3–1 h) to sudden wind increases. For the remaining cases, the new [Van der Westhuysen, A.J., Zijlema, M., and Battjes, J.A., 2007. Nonlinear saturation-based whitecapping dissipation in SWAN for deep and shallow water, Coast. Eng., 54, 151–170] SWAN physics yielded better results than the standard physics of Komen, G.J., Hasselmann, S., Hasselmann, K., 1984. On the existence of a fully developed wind-sea spectrum. J. Phys. Oceanogr. 14, 1271–1285, except for persistent overestimations that were found for short fetches. The present data set contains many interesting cases for detailed model validation and for further studies into the evolution of wind waves in shallow lakes.     

An empirical model for infragravity swash on barred beaches

Idealized computational simulations with the nearshore model XBeach were carried out for a series of barred beach configurations in order to quantify the impact of nearshore bars on infragravity swash. Results show that nearshore bar systems reduce infragravity swash energy at the shoreline. The amount of swash reduction was found to correlate with both bar depth and rip width, when a rip channel is present. In order to develop a generalized empirical model for significant infragravity swash for barred beaches, the simulations were used to extend the empirical swash model of Stockdon et al. (2006) to include bar characteristics. The developed empirical model relates significant infragravity swash to incident wave conditions and nearshore bar depth. With respect to Stockdon et al. (2006), this new model improves predictive skill by reducing root-mean-square error by 50% for the computational simulations and by 15% when applied to a range of field data.     

Potential impact of sea level rise on coastal surges in southeast Louisiana

Potential impacts of 0.5 and 1.0 m of relative sea level rise (RSLR) on hurricane surge and waves in southeast Louisiana are investigated using the numerical storm surge model ADCIRC and the nearshore spectral wave model STWAVE. The models were applied for six hypothetic hurricanes that produce approximately 100 yr water levels in southeastern Louisiana. In areas of maximum surge, the impact of RSLR on surge was generally linear (equal to the RSLR). In wetland or wetland-fronted areas of moderate peak surges (2-3 m), the surge levels were increased by as much as 1-3 m (in addition to the RSLR). The surge increase is as much as double and triple the RSLR over broad areas and as much as five times the RSLR in isolated areas. Waves increase significantly in shallow areas due to the combined increases in water depth due to RSLR and surge increases. Maximum increases in wave height for the modeled storms were 1-1.5 m. Surge propagation over broad, shallow, wetland areas is highly sensitive to RSLR. Wave heights also generally increased for all RSLR cases. These increases were significant (0.5-1.5 m for 1 m RSLR), but less dramatic than the surge increases.     

Directional spread parameter at intermediate water depth

The characteristics of directional spread parameters at intermediate water depth are investigated based on a cosine power ‘2s' directional spreading model. This is based on wave measurements carried out using a Datawell directional waverider buoy in 23 m water depth. An empirical equation for the frequency dependent directional spreading parameter is presented. Directional spreading function estimated based on the Maximum Entropy Method is compared with those obtained using a cosine power ‘2s' parameter model. A set of empirical equations relating the directional spreading parameter corresponding to the peak of wave spectrum to other wave parameters like significant wave height and period are obtained. It shows that the wave directional spreading at peak wave frequency can be related to the non-linearity parameter, which allows estimation of directional spreading without reference to wind information.     

Forecasting ocean waves: Comparing a physics-based model with statistical models

The literature on ocean wave forecasting falls into two categories, physics-based models and statistical methods. Since these two approaches have evolved independently, it is of interest to determine which approach can predict more accurately, and over what time horizons. This paper runs a comparative analysis of a well-known physics-based model for simulating waves near shore, SWAN, and two statistical techniques, time-varying parameter regression and a frequency domain algorithm. Forecasts are run for the significant wave height, over horizons ranging from the current period (i.e., the analysis time) to 15 h. Seven data sets, four from the Pacific Ocean and three from the Gulf of Mexico, are used to evaluate the forecasts. The statistical models do extremely well at short horizons, producing more accurate forecasts in the 1–5 hour range. The SWAN model is superior at longer horizons. The crossover point, at which the forecast error from the two methods converges, is in the area of 6 h. Based on these results, the choice of statistical versus physics-based models will depend on the uses to which the forecasts will be put. Utilities operating wave farms, which need to forecast at very short horizons, may prefer statistical techniques. Navies or shipping companies interested in oceanic conditions over longer horizons will prefer physics-based models.     

Modelling storm impacts on beaches, dunes and barrier islands

A new nearshore numerical model approach to assess the natural coastal response during time-varying storm and hurricane conditions, including dune erosion, overwash and breaching, is validated with a series of analytical, laboratory and field test cases. Innovations include a non-stationary wave driver with directional spreading to account for wave-group generated surf and swash motions and an avalanching mechanism providing a smooth and robust solution for slumping of sand during dune erosion. The model performs well in different situations including dune erosion, overwash and breaching with specific emphasis on swash dynamics, avalanching and 2DH effects; these situations are all modelled using a standard set of parameter settings. The results show the importance of infragravity waves in extending the reach of the resolved processes to the dune front. The simple approach to account for slumping of the dune face by avalanching makes the model easily applicable in two dimensions and applying the same settings good results are obtained both for dune erosion and breaching.     

变浅作用下的长周期重力波的实验研究

利用大型水槽设计了在由深水到近岸不同坡度处海浪在变浅作用下诱导产生的长周期重力波的实验。正态随机海浪在深水生成并沿斜坡向浅水传播,记录了不同水深处波面高度随时间的变化过程并进行统计分析和谱分析。实验数据分析结果表明,长周期重力波的能量随着水深的变浅而增高,其谱锋频率位于0.2~0.3fp附近,这里fp是深水正态海浪过程的谱峰频率。长周期重力波的能量与入射波的能量比与波面高度分布的偏度密切相关。进一步分析了两种波动的能量谱峰值比和波面高度分布偏度的相关关系,获得了经验关系,为预测近岸浅水长周期重力波提供了科学依据。     

Megacusps on rip channel bathymetry: Observations and modeling

The formation of beach megacusps along the shoreline of southern Monterey Bay, CA, is investigated using time-averaged video and simulated with XBeach, a recently developed coastal sediment transport model. Investigations focus on the hydrodynamic role played by the bay's ever-present rip channels. A review of four years of video and wave data from Sand City, CA, indicates that megacusps most often form shoreward of rip channels under larger waves (significant wave height (H_s) = 1.5–2.0 m). However, they also occasionally appear shoreward of shoals when waves are smaller (H_s ~ 1 m) and the mean water level is higher on the beach. After calibration to the Sand City site, XBeach is shown to hindcast measured shoreline change moderately well (skill = 0.41) but to overpredict the erosion of the swash region and beach face. Simulations with small to moderate waves (H_s = 0.5–1.2 m) suggest, similar to field data, that megacusps will form shoreward of either rip channels or shoals, depending on mean daily water level and pre-existing beach shape. A frequency-based analysis of sediment transport forcing is performed, decomposing transport processes to the mean, infragravity, and very-low-frequency (VLF) contributions for two highlighted cases. Results indicate that the mean flow plays the dominant role in both types of megacusp formation, but that VLF oscillations in sediment concentration and advective flow are also significant.     

Wave prediction in a port using a fully nonlinear Boussinesq wave model

A wave forecasting system using FUNWAVE-TVD which is based on the fully nonlinear Boussinesq equations by Chen(2006) was developed to provide an accurate wave prediction in the Port of Busan, South Korea. This system is linked to the Korea Operational Oceanographic System(KOOS) developed by Park et al.(2015). The computational domain covers a region of 9.6 km×7.0 km with a grid size of 2 m in both directions, which is sufficient to resolve short waves and dominant sea states. The total number of grid points exceeds 16 millions,making the model computational expensive. To provide real-time forecasting, an interpolation method, which is based on pre-calculated results of FUNWAVE-TVD and SWAN forecasting results at the FUNWAVE-TVD offshore boundary, was used. A total of 45 cases were pre-calculated, which took 71 days on 924 computational cores of a Linux cluster system. Wind wave generation and propagation from the deep water were computed using the SWAN in KOOS. SWAN results provided a boundary condition for the FUNWAVE-TVD forecasting system. To verify the model, wave observations were conducted at three locations inside the port in a time period of more than 7 months. A model/model comparison between FUNWAVE-TVD and SWAN was also carried out. It is found that, FUNWAVE-TVD improves the forecasting results significantly compared to SWAN which underestimates wave heights in sheltered areas due to incorrect physical mechanism of wave diffraction, as well as large wave heights caused by wave reflections inside the port.     

近岸较大区域波浪数值模型的比较

将适用于近岸较大区域波浪传播变形的三种模型,即基于抛物型缓坡方程的不规则波模型、引入浅水波浪谱 TMA 谱的 SWAN(simulating waves nearshore)模型以及采用默认 JONSWAP 谱的 SWAN模型应用于特拉华大学(University of Delaware)圆形浅滩实验进行比较.结果显示,抛物型缓坡方程和SWAN 的模拟结果与实验所测数据符合都比较好; SWAN 在非线性作用较强的浅滩中心及靠后部效果更佳,而抛物型缓坡方程由于没有考虑非线性作用,模拟得到的最大波高较实测值偏高,且波高变化较为剧烈.     

水动力时空变化对近岸风浪演化的影响——以渤海湾西南岸为例

采用SWAN模型和ADCIRC模型建立了风浪、潮汐和水流联合作用的耦合数值模式,并通过渤海湾西南岸实测资料对该模式进行了验证。利用该模式分析了近岸区水位和流场时空变化对风浪模拟结果的影响,计算结果表明水位变化对近岸区风浪模拟结果有显著影响,特别是中等大风过程高潮位时波高受水位影响的变化幅值可达0.5m以上,且水深越浅影响越大。但在岸滩平缓的近岸海域由于流速、流向的时空变化不太剧烈,流场作用和波浪辐射应力作用对波浪场的影响都基本可以忽略。在模拟近岸风浪过程时,应选用耦合模式。     

Evaluation of the wave transformation in an open bay with two spectral models

An evaluation of two state of the art phase averaged wave models for the transformation scale, SWAN and STWAVE, is carried out in the present work. The target area is the Obidos Bay located in the central part of the Portuguese continental nearshore. The wave input for the two models is provided by an offshore buoy. In order to compare the nearshore outputs of the wave models against in-situ measurements, a directional buoy and an ADCP, operating in intermediate water depth, are used. The wave parameters considered for comparisons are significant wave height, peak period and wave direction. Sensitivity analyses studies and evaluations in the spectral and geographical spaces concerning the results of the two models are also carried out in both intermediate and shallow water. The present study provides some information on the performances of the two wave models in different forcing conditions as well as on their sensitivity in relationship with various input parameters and some physical processes. STWAVE appears to be faster and more robust than SWAN, which on the other hand has more options and flexibility. In statistical terms the results are comparable.     

The effect of temporal wave averaging on the performance of an empirical shoreline evolution model

The effect of using time-averaged wave statistics in a simple empirical model for shoreline change is investigated. The model was first calibrated with a six-year time series of hourly wave conditions and weekly shoreline position at the Gold Coast, Australia. The model was then recalibrated with the hourly waves averaged over intervals up to 1 year. With wave averaging up to 2 days, model performance was approximately constant (squared correlation r² ~ 0.61–0.62), with only small changes in the values of empirical model parameters (e.g. the beach response coefficient c varied by less than 4%). With between 2 and 40 day averaging, individual storms are not resolved; model skill decreased only modestly (r² ~ 0.55), but c varied erratically by up to 40% of the original value. That is, optimal model coefficients depend on wave averaging, an undesirable result. With increased averaging (> 40 days) seasonal variability in the wave field is not resolved well and model skill declined markedly. Thus, temporal averaging of wave conditions increases numerical efficiency, but over-averaging degrades model performance and distorts best-fit values of model free parameters.     

Modeling of storm-induced coastal flooding for emergency management

This paper describes a model package that simulates coastal flooding resulting from storm surge and waves generated by tropical cyclones. The package consists of four component models implemented at three levels of nested geographic regions, namely, ocean, coastal, and nearshore. The operation is automated through a preprocessor that prepares the computational grids and input atmospheric conditions and manages the data transfer between components. The third generation spectral wave model WAM and a nonlinear long-wave model calculate respectively the wave conditions and storm surge over the ocean region. The simulation results define the water levels and boundary conditions for the model SWAN to transform the storm waves in coastal regions. The storm surge and local tides define the water level in each nearshore region, where a Boussinesq model uses the wave spectra output from SWAN to simulate the surf-zone processes and runup along the coastline. The package is applied to hindcast the coastal flooding caused by Hurricanes Iwa and Iniki, which hit the Hawaiian Island of Kauai in 1982 and 1992, respectively. The model results indicate good agreement with the storm-water levels and overwash debris lines recorded during and after the events, demonstrating the capability of the model package as a forecast tool for emergency management.     

Field observations of waves generated by passing ships: A note

This short contribution reports the results of a field study on the nearshore characteristics of waves generated by both conventional and high speed passenger ferries. The field observations took place in the late summer of 2005, at a beach close to the port of Mytilene (Island of Lesbos, Greece), and involved the visual observation of ship waves, using digital video recordings and image processing techniques. The results showed that passage of the fast ferry was associated with a longer, more complex and energetic nearshore event; this event not only did include higher nearshore waves (up to 0.74 m) and was organised in different wave packets, but it was also an order of magnitude longer (∼ 680 s) than the conventional ferry event. Regarding the effects on beach sediment dynamics, the fast ferry waves were estimated to be very efficient in mobilising the nearshore sediments in contrast to those of the conventional ferry. The fast ferry service appears to generate daily prolonged nearshore events, which contain waves with higher energy than those expected from the normal summer wind wave regime of the area; these events also include some high and very steep waves, which can be particularly erosive. Therefore, fast ferry wakes may have considerable impacts on the seasonal beach sediment dynamics/morphodynamics and the nearshore benthic ecology, as well as they may pose significant risks to bathers, affecting the recreational use of the beaches exposed to fast ferry traffic. Finally, the study has shown that satisfactory field observations of the nearshore characteristics of ship-generated (and wind) waves can be obtained using inshore deployments of calibrated poles, digital video cameras and appropriate image processing algorithms.     

Influence of bathymetric fluctuations on coastal storm surge

Natural events constantly alter nearshore bathymetric properties. Hurricanes particularly affect bathymetry as they pass over a body of water. To compute an accurate forecast or recreate a hurricane's effects through hindcasting techniques, an operational bathymetry data set must be known in advance. However, obtaining and maintaining current and accurate bathymetric data can be costly and difficult to manage. In this paper we examine the extent to which variations in nearshore bathymetry affect the storm surge at the coast. A common question for wave and surge modeling is, “how good is the bathymetric data?” If we can allow for a range of fluctuations in the bathymetry without significantly adjusting the results of the surge predictions, we can potentially save months of field work and millions of dollars. A one-dimensional (1D) analytical solution for waves and water level is developed for initial testing. In the 1D case we find that as long as the amplitudes of the bathymetric fluctuations are less than 60% of the original depth, the surge at the coast is within ± 10% of the surge generated on the initial bottom slope. If the fluctuation produces a hole, a deepening of the local bathymetry, within 80% of the local water depth, the coastal storm surge calculated is still within 10% of the unperturbed value computed for bottom slopes shallower than 1:20. In addition, we find there is an optimum distance offshore for each sloped profile that corresponds to a depth between 25 and 40 m, beyond which the effects of bathymetric fluctuations begin to decrease. A coupled 2D modeling system is implemented to test our hypothesis along a realistic coastline. After selecting three study sites, we vary the bathymetry at the selected locations by ± 20%. Consistent with the 1D tests, the storm surge at the shoreline varies by less than 5%.     

Surface wave simulation during winter with sea ice in the Bohai Sea

Sea ice can attenuate wave energy significantly when waves propagate through ice covers.In this study,a third-generation wave model called simulating wave nearshore(SWAN)was advanced to include damping of wave energy due to friction in the boundary layer below the ice.With the addition of an eddy viscosity wave-ice model,the resulting new SWAN model was applied to simulate wave height in the Bohai Sea during the freezing winter.Its performance was validated with available buoy data near the ice edge,and the new model showed an improvement in accuracy because it considered the ice effect on waves.We then performed a wave hindcast for the Bohai Sea during a freezing period in the winter of 2016 that had the severest ice conditions in recent years and found that the mean significant wave height changed by approximately 16.52%.In the Liaodong Bay,where sea ice concentration is highest,the change reached 32.57%,compared with the most recent SWAN model version.The average influence of sea ice on wave height simulation was also evaluated over a five-year(2013-2017)hindcast during January and February.We found that the wave height decrease was more significant in storm conditions even the eddy viscosity wave-ice model itself showed no advantage on damping stronger waves.     

Wave dissipation by vegetation with layer schematization in SWAN

The energy of waves propagating through vegetation is dissipated due to the work done by the waves on the vegetation. Dalrymple et al. (1984) estimated wave dissipation by integrating the force on a cylinder over its vertical extent. This was extended by Mendez and Losada (2004) to include varying depths and the effects of wave damping due to vegetation and wave breaking for narrow-banded random waves. This paper describes the wave dissipation over a vegetation field by the implementation of the Mendez and Losada formulation in a full spectrum model SWAN, with an extension to include a vertical layer schematization for the vegetation. The present model is validated with the original equation and results from Mendez and Losada (2004). The sensitivity of the model to the shape of the frequency spectrum, directional spreading and layer schematization are investigated. The model is then applied to field measurements by using a vegetation factor. This model has the ability to calculate two-dimensional wave dissipation over a vegetation field including some important aspects such as breaking and diffraction as used in SWAN model.