9. Conclusion
The present work presents a description of the propagation of waves in the interior of a submerged canopy with anisotropic flow resistance; diffusion was neglected. The theory is based on the assumptions that boundary layer effects on top of the canopy can be neglected, that the waves are (almost) linear and there is no wave breaking. Furthermore, it is assumed that the outer flow is irrotational, while the inner flow is subject to friction. The second order wave properties such as wave energy density, wave energy flux, the corresponding vegetated group velocity, the radiation stress tensor and the Eulerian and Lagrangian Stokes drifts were derived. It was shown that the vegetated group velocity does not equal the classical group velocity from non-dissipative linear wave theory and there is an non-equal splitting between the potential, Ep, and kinetic, Ek, wave energy densities. Both of these properties have a direct implication for the large scale modelling of waves in vegetation, since they can cause shoaling, refraction and conservative decay in the wave height due to changes in the ratio Ek/Ep. The derivations of the Stokes drift have shown that the Eulerian and Lagrangian formulations are identical even in the presence of vegetation, but it is important to include the Lagrangian flux at the top of the canopy due to the discontinuity in the horizontal velocity. The expressions are not easily comparable due to their complexity, but through an evaluation of the finite Stokes velocity at the free surface, it was seen that the horizontal Eulerian and Lagrangian volume fluxes are identical. The Stokes drift increases in the presence of vegetation. The finite vertical Stokes drift must be balanced by a gradient in the Eulerian, mean velocity field (return flow) in the limit of a steady solution. It was hypothesised that the Stokes velocity is established faster than the return flow, thus the vertical Stokes velocity is (partly) responsible for the mass flux needed to establish the mean waveinduced setup. This mean wave-induced setup is required to balance the gradient in the radiation stress tensor. Further validation of the theory is needed. The validation material should preferably include data on (i) spatial variation in the surface elevation and (ii) orbital velocities over the height of the canopy for rigid stems. This will allow for a direct evaluation of the drag coefficients in a reduced velocity field.