(right) Input equations and parameters for running the wave simulation program (Bostater et al. [11, 12]) for simulating the wave patch shown in Figure 4 and the measured and simulation distributions (left image). Comparison of the calculated Weibull cumulative curve (solid line) and the cumulative curve of the wave energy spectrum extracted from the measured wave displacements with estimated Weibull parameters of α = 0.55 (Hz) and β = 1.31 (dimensionless) as shown in the table.
This paper discusses a shallow-water wave height measurement method that uses high definition video cameras to image a water surface wave patch. Wave height time series are extracted from water surface video sequences. Wave features such as the wavelength distribution and energy contained in a wave patch (W/m2) were obtained by analyzing the extracted wave height time series and expressing the wind-driven wave energy as a wave energy spectrum. A Weibull probability distribution was used as the mathematical form of the energy spectrum. Wave spectra are used as input to a wave patch simulation model that generates simulated wind-driven wave images. The measurement protocol is inexpensive, easy to implement, and useful to calibrate and validate wind-driven wave models. The protocol is used to understand resuspension of bottom muds due to wind waves in shallow waters. Scaled staff gauges made of polyvinyl chloride (PVC) materials are placed in shallow water and imaged at 30 Hz followed by frame based image analysis to extract wave height time series. Wave spectra calculated using the fast Fourier transform (FFT) results in a Weibull probability distribution function (WPDF) energy spectrum. The estimated wave spectrum is used to estimate wave energy in W/m2 followed by generation of wave patch simulations of the water surface. Simulated wave patches are compared with the sensor-based wave patch video measurements. Sensitivity analysis of coefficients α and β in the model are used to adjust the synthetic wave images to measured wave patch images. The approach allows one to obtain an estimate of the energy (W/m2) transferred from the local wind field to a water surface gravity wave patch.
- water surface waves
- gravity waves
- wind driven waves
- wave patches
- wave patch energy
- remote sensing
- staff gauges
- video imaging
- video analysis
- energy spectra
- Weibull distributions
- synthetic imaging
- wind energy
Differential heating of the earth causes large-scale atmospheric pressure gradients. These gradients and the acceleration due to rotation of the earth result in the major wind fields around the globe. The winds cause friction on the water surface and result in the generation of wind-driven water surface gravity waves. These wind-driven gravity waves are a major source of transfer of energy from the sun and atmosphere to the earth’s water system on a global scale as well as in small water regions. Measurements of these wind-driven shallow-water surface gravity wave characteristics such as wave height, period, and direction of propagation are key to understanding the magnitude of energy in a small gravity wave field or “wave patch.” These surface wind-driven waves cause orbital motion of water parcels . In water the wind stress or friction causes the downward movement or transport of momentum from the water surface to the water column. The downward transport of momentum causes internal friction in the water in the form of turbulent friction or viscosity, and the associated circular eddy-type motions influence the bottom lutocline within the bottom boundary layer [2, 3]. The resulting shearing forces at the bottom in shallow water cause resuspension of and fluidization of bottom sediments and associated nutrients and/or trace metals. The fluidization can result in fluidized mud, muck, high-bottom water turbidity, and nephelometric wave motions (internal wavelike motions that are observed in acoustic imaging of the bottom water column). Thus, water quality variables are influenced by the surface water wave field driven by surface winds. Bottom orbital or elliptical velocities due to these wind-driven gravity waves in shallow water can be calculated using surface wind-driven gravity wave information, such as significant wave height and wave spectra [3, 4].
Wind-driven gravity waves in water can be measured using numerous techniques and descriptions, including satellite altimeters, video cameras, wave buoys, and many types of wave gauges . Altimeters flown on satellites such as Seasat in 1978, Geosat from 1985 to 1988, ERS-1 and ERS-2 from 1991, TOPEX/Poseidon from 1992, Jason from 2001, and Envisat can provide data for monthly mean maps of wave heights and the variability of wave energy density in time and space . Wave height estimations based on binocular or trinocular imaging also allow acquisition of both spatial and temporal wave information for surface wave patches [6, 7, 14]. Many types of wave gauges including resistance-type, capacitance-type, and wave pressure gauges can also provide data for quantifiable wind-driven gravity wave height estimations . However, there is a need to develop measurement system methods [17, 18] for shallow-water areas since most systems developed to date can only be used in deep waters (>3–5 m).
In this research, a technique or protocol that combines gravity wave staff gauges, HD video cameras , and image analyses is used to measure wind-driven gravity waves in shallow-water environments that receive energy from the atmosphere. Wave height time series are extracted from video sequence taken using high-definition (HD) video cameras at a 30 Hz frame rate. Significant wave height, wave spectrum, and wave energy (W m2) are estimated from the video-based wave height time series.
In this chapter, techniques to describe the protocol are demonstrated for shallow water near the Atlantic Ocean and Space Coast Region of the Banana River estuary and watershed, near Melbourne, Florida. Simulated wave patch synthetic images are then generated using estimated wave spectrum derived from video imagery and compared to the video images.
2. Techniques and methods
2.1. Video imaging and staff gauge system
Staff gauges and HD video cameras (JVC Everio HD cameras) are used to measure the surface wave heights in shallow water (1–2 m) as shown in Figure 1(a). Four staff gauges are mounted in the water by inserting them into the bottom sediment. Staff gauges are constructed from white PVC pipes. Waterproof paper with 0.5 spaced horizontal lines is taped to each gauge pipe.
Images recorded at 1920 × 1080 pixels (30 per second of recording time) obtained from the video sequence yield spatial resolution (using lens zooming) is 0.04 cm/pixel. Each line is either 1 mm (≈2 pixels) or 0.5 mm (≈1 pixel) as shown in Figure 1(b). The HD cameras are mounted on a tripod on the shore or from a fixed platform or mount in the water. The video staff gauges allow one to define a wave patch of ~4.5 × 4.5 m (Wx = Wy = 4.5 m).
An inverted T-shaped (Figure 1(c)) aluminum and steel camera optical mount is used to mount one, two, or three HD cameras. Cameras can be time synchronized, using a known or fixed “digital zoom” mode. During measurement experiments, one camera is focused upon and zoomed (18×) on a staff gauge, and the other two cameras image the entire wave patch and record at 30 Hz or frames per second (fps) (Figure 2).
2.2. Image processing
In order to determine the water surface level based on video sequences, three steps are needed as summarized below:
Step 1: Determine the region of interest (ROI)
The water surface level at the staff gauges was clearly observed from the 18× zoomed video frame images. A representative image is shown in Figure 3(a). A ROI was selected to be the area between the top horizontal scale line on the gauge and the bottom margin of the image and was slightly less in width than the actual gauge as defined by the rectangular area in Figure 3(b). A thinner ROI than the actual gauge prevents pixels outside of the staff gauge influencing the results.
Step 2: Average the digital counts within the ROI
Within the ROI, the mean pixel value across each horizontal line is calculated to compress the ROI into a “representative line” shown in Figure 3(b). The average of the three color bands (R, G, B) is calculated in this step as well. The oscillating line in Figure 3(c) is the average pixel value of the example frame, and the air-water interface is obvious.
Step 3: Locate the air-water level on the gauge
The shape of the pixel value curve changes dramatically at the interface of water and gauge as shown in Figure 3(c), which is the water surface level. A threshold value is selected in order to automatically locate the interface within each video frame.
The threshold value was determined based on the average of part of the “representative line” (first 500 pixels from the top in the example shown in Figure 3 (b)) and it slightly varies from frame to frame. The key is to semiautomatically select a threshold value and apply this value to each frame. The location of the start of values lower than the threshold value is considered as the water surface. Figure 4 shows the result of 30 Hz time series of the water wave surface oscillations calculated from the example video sequence.
3. Wave spectrum development
A wave energy density spectrum of a wind-driven surface gravity wave patch can be expressed in units of and is related to the wave height or the variance of the wave displacement, given as , where is the wave displacement and the wave height (H) is twice the wave displacement or . Given a measured wave time series, the wave spectrum can be calculated using the fast Fourier transform (FFT) :
where is the wave spectrum density in the unit of cm2/Hz, is temporal frequency in the unit of Hz (1/s), is the Fourier transform of the original signal (dimensionless), is the conjugate of (dimensionless), is the number of the data point in the measured sequence (dimensionless), is the instant time corresponding to the data point in the unit of second (s), and is the length of measurement time in the unit of second (s). A two-parameter Weibull distribution function was used as a test model for mathematical expression of the spectrum . The Weibull probability density function (PDF) and cumulative density function (CDF) are given by :
where is a variable that follows the Weibull distribution. In this research Y has the unit of 1/s (temporal frequency) or 1/m (wavenumber (spatial frequency); see Figure 5(a)), and Table 1 on x-axis, = scale parameter has the units of the variable Y, which is 1/s or 1/m, = nondimensional shape parameter.
4. Experimental results
4.1. Image-based wave height time series and statistical characteristics
By using the three steps described in Section 2.2, the water surface-level change throughout the video sequence was obtained as shown in Table 1 and Figure 4. Mean wave height of the measured series was 7.9 cm, and significant wave height Hs was 13 cm. The significant wave height Hs is defined as average of the highest one-third of the wave heights . The speed of this wave series extracting process is about two frames per second.
4.2. Wave spectrum, wave energy, and simulation for the Banana River site
The wave spectrum coefficients estimated using the video-based wave series were and . The calculated wave energy spectrum of the measured wave height time series on temporal frequency domain (Hz) is calculated from Eq. (4):
where is the wave energy spectral density (m2/Hz), is the temporal frequency (Hz), = 0.55 is the Weibull scale parameter (Hz), and = 1.31 is the Weibull shape parameter (dimensionless). The estimated wave spectrum can then be used in the spectral wave patch model [11, 12]. Since the wave simulation program by Bostater et al. [11, 12] was designed to generate a shallow-water gravity wave patch that was independent of time, a dispersion relation was applied to transfer the energy spectrum from temporal frequency domain (Hz) to spatial frequency domain (1/m):
where is the temporal frequency (Hz), is the spatial frequency wavenumber (1/m), is the gravity acceleration (m/s2), and is the water depth (m).
The wave energy in the units of can now be calculated using the estimated wave spectrum, based on the area under the spectrum curve :
where is the wave energy in units of (), is the wave energy spectrum in units of , is the temporal frequency in units of Hz, is the density of the water used in this research, and is the acceleration due to gravity. The discrete approximation of wave energy in this research is calculated by summing the area under the spectral curve:
where is the total number of discrete frequencies integral, is the magnitude of the energy spectrum at each frequency ωn in units of , and is the discrete frequency interval in units of Hz. A representative shallow-water wind-driven gravity wave can be simulated using the estimated wave spectrum from measurements as shown in Figure 4(a). A wave patch simulation model specific for shallow water such as Banana River and Indian River Lagoon was used . The simulated mean wave height is 7.6 cm, and significant wave height is 0.1 m, which is consistent with the video-based measurements. Sensitivity analysis of the simulation model can also be conducted in order to determine the effects of the model coefficients and on the simulated wave patches. The results suggest that the model coefficients and can affect wave height scale and wave pattern of the simulated wave patches. Therefore, it is reasonable to assume that and are related to physical or environmental variables that affect the wave conditions, such as wind speed, wind duration, wind direction, fetch, water depth, bottom slope, etc. When the wave patch model coefficients are adjusted, realistic random water waves are produced as shown in Figure 4(b). The wave energy spectrum (Eq. (4)) can also be introduced or multiplied by a nondimensional scaling coefficient A = 15 in this example for sensitivity analysis related to wave patch simulations shown Figures 4 and 5 (Table 1).
A sensitivity analysis of the wave model parameter α (scale parameter) and the β (nondimensional shape parameter) allows one to clearly see how the water surface wave field would look for different parameter values. The sensitivity simulation results are shown in Figure 5 with selected wave patch synthetic images for selected model parameters. Figure 6 shows example wave height displacements, wavenumbers, and significant wave height, and the resulting wave patch energy is simulated wave patches. Figure 7 shows the associated (a–f) wave energy and significant wave heights obtained from synthetic image model runs. The results demonstrate the utility of making observations, followed by sensitivity analysis and resulting wave patch energy, and associated significant wave heights derived from the synthetic images which are representative measured wind-driven surface gravity wave field patches.
This research developed an imaging methodology to measure wind-driven gravity waves in shallow water using high-definition (HD) video cameras and specially constructed staff gauges. Wave spectrum, wave energy, and significant wave height are estimated from video-based wave height measurement, which can be used to estimate bottom velocity and sediment resuspension. Simulated wave field images are generated using estimated wave spectrum; wave characteristics (significant wave heights) based on the simulation images agree with the in situ video-based measurement protocol and methods described in Section 3 above.
Coefficients and in the spectral wave simulation model can affect the wave height and pattern of the simulated wave patches, and they are very likely to be related with some physical and environmental variables. Ongoing research is investigating these variables for these small amplitude water surface waves. The approach developed and described above utilizes the conceptualization shown in Figure 8. The energy stored in a surface water gravity wave field drives the sediment resuspension process in the bottom boundary layer and lutocline consisting of fluid mud and muck. The flowchart shown in Figure 9 summarizes the protocol and methods described above.
Each time series of water surface displacements is linearly detrended before spectral analysis and calculation of the wind wave energy in the wave field. The mono imaging used in the protocol is the simplest method because water surface elevations can be easily obtained using in situ staff gauges. Algorithms can then be used to calculate wave displacement and wave heights using image processing techniques that analyze each video frame as an instantaneous realization of the water surface elevation or displacement at 0.033 s time intervals. The key for successful mono imaging is the bottom mounted staff gauge used for automatic water surface elevation detection. The staff gauge used in this research was scaled to the appropriate measured wave height field and the gauge and inserted into the bottom sediments.
Simulation of wind-driven water wave patch results shown in Figure 5 indicates that as increases, the wave displacement increases and smaller wavelength waves appear on top of the larger waves. The effect of is just the opposite of . Increasing causes decrease in wave displacements with coincident disappearance of smaller waves on top of the larger waves.
The wave energy spectrum graphs (obtained from the synthetic images) shown in Figures 6 and 7 also provide a visual demonstration of the effect of and on wind-driven waves using the Weibull model. Increasing causes the energy distribution to spread into larger wavenumber in the spatial frequency domain (1/m). This reveals that a simulated wave patch with larger is consistent with waves with a greater variety of wavelengths. However, increasing causes the narrowing of the energy distribution to smaller wavenumber (1/m) and suppresses the frequency of larger wavelength waves in a simulated wave patch.
From an energy point of view, wave energy in units of W/m2 increases with increasing and decreasing . The wave energy can be calculated within the simulated wave patch area (4.5 × 4.5 m) using the energy distribution from the sum of the variance in the temporal frequency domain . The difference of energy stored in the simulated surface wave patches results in differences in the surface wave heights. In this research, significant wave height was used to represent the simulated surface state.
The significant wave height is defined as the average of the highest one-third of the wave heights . It can be used as a representation of wave height of an observed or simulated wave patch. The tendency of simulated significant wave heights changes with and as shown in Figures 5–7 and suggests that simulated water surfaces are smooth with a small and large using the Weibull wave spectrum. A larger and smaller tend to produce a rougher surface with more energy. However, considering the shallow-water circumstances of this research, and should be limited to a range appropriate to the experimental conditions from field measurements. The authors believe that the techniques, method, and the resulting protocol for the measurement of wind-driven surface water gravity waves described and summarized in Figures 9 and 10 will allow other scientists, engineers, and meteorologists access to an inexpensive and rapidly deployable instrumental approach for shallow waters.
Practical use and applications of water surface video imaging include studying of coastal nearshore processes . The use of fixed sensor video platforms such as bridges, piers, and building which view the littoral zone for water quality constituent determinations has also been demonstrated  including the mathematical methods to correct video hyperspectral push broom imagery taken at high oblique angles. Video imaging of bottom features such as seagrasses and sand bottom features has also been demonstrated using practical imaging techniques . The imaging technique described above is being used to assist in assessing moving fluid mud and muck in Indian River Lagoon and Banana River estuarine areas and tributaries . In many shallow-water estuarine regions, wind-driven gravity waves are responsible for resuspension and transport of muds and decaying organic matter within the lutocline [17, 18]. Figure 11 shows the results of the vertical profile of horizontal mass transport (fluxes) using the water wave video imaging technique in conjunction with simultaneous deployment of a vertical array of six passive Sondes  during October 2017. Wind speed during the video wave gauge measurements were 16-17 knots.
Sonde array  deployment was 16–17 knots, with east winds in a shallow ≈0.7 m water column and 20–30 cm amplitude wind-driven gravity waves at the location in Banana River near Pineda Causeway [2, 16] near the Atlantic Ocean in Florida (28°12′33.02″N, 80°38′14.72″W). The use of the staff gauge video imaging and Sondes shown in Figure 11 showed that the shallow-water gravity waves increase the bottom water fine particulate transport or flux of 15,071 g m−2 day−1 dry weight—nearly 450 times greater than the surface values of 13.1 g m−2 day−1. The vertical profile of sediment and particulate mass flux or movement under the waves follows the same profile found in previous studies in tributaries during high flow conditions . The above recent result describes the application of the novel approach to water wave imaging in shallow littoral zones. In essence, the novel approach  requires no expensive equipment , can be easily constructed, and  is self-calibrated to video imagery. The method can be used in conjunction with sophisticated wave patch imaging models and when used as a system can be used to improve scientific and engineering understanding of sediment transport due to wind-driven waves. The imaging techniques are currently being used in the field of transportation construction engineering concerned with pile driving rebound imaging and related soil engineering problems.
The work presented in this paper has been supported in part by the Northrop Grumman Corporation, the NASA, the Kennedy Space Center, the KB Science, the National Science Foundation, the US-Canadian Fulbright Program, and the US Department of Education, the FIPSE and European Union’s grant Atlantis STARS (Sensing Technology and Robotic Systems) to Florida Institute of Technology, the Budapest University of Technology and Economics (BME), and the Belgium Royal Military Academy, Brussels, in order to support the involvement of students in obtaining international dual US-EU undergraduate engineering degrees.