Wave run-up is the height to which waves run up the slope of a revetment, bank or dike, regardless of whether the waves are breaking or not. Conversely, wave run-down is the height to which waves recede. These heights are always measured vertically (and not along the slope). The wave run-up height, denoted by , , or , is a very important parameter in coastal engineering as, together with the design highest still water level, it determines the required crest height of a dike or revetment.
The first scientific measurements of wave run-up were carried out by the Lorentz Committee in preparation for the works to close off the Zuiderzee. [1] The Committee measured the wave height and wave run-up at various locations in 1920, but established that state of the art methods for measuring waves in the field during storms were inadequate. As a result, scale tests were also undertaken, but these also proved to be of very limited efficacy due to the fact that only regular waves (idealised, periodic waves with constant amplitude and a fixed time period between successive wave crests, following a sinusoidal pattern) could be modelled at the time.
The methods and technology available to the committee at the time did not permit model testing of the more realistic and complex irregular waves (consisting of varying heights, periods and directions), which provide a more accurate representation of the actual conditions faced by coastal structures and shorelines.
It was found, however, that the depth in front of the dike is very important for wave run-up and that, at least for the range of observations in the committee's measurements, the slope ratio does not play a major role. Nearly all dikes in the Netherlands at that time had a slope of 1:3.
Current knowledge indicates that during storms and on gentle coastal slopes, the significant wave height is approximately half the water depth. This relationship appears to be accurate, and the observation is more pronounced for slopes around 1:3.
This research was continued during the Zuiderzee Works, and eventually led to the (old) Delft formula for wave run-up: [2]
in which:
This formula proved to be generally applicable for smooth slopes and relatively steep (storm) waves. Subsequently, it was discovered that longer (swell) waves resulted in higher run-up. To account for this, the wave period was incorporated into the formula using the Iribarren number, , leading to the development of Hunt's Formula: [3]
This formula was also valid for regular waves. The Old Delft Formula and Hunt's Formula are identical for waves with a steepness of 1/64, or about 2%. For higher values of , Hunt's formula has a limit value: .
In 1988, van der Meer provided formulae for wave run-up on rubble mound breakwaters, based on tests with rock-armoured straight slopes. He also introduced a notional permeability factor for the structure. [4] This factor also accounts for the effect of the pore volume. Defining at the run-up level of exceedance probability , the formula, valid for and head-on waves, is:
The term reaches a constant maximum value equal to in the case of permeable structures, i.e., . This corresponds to the region of surging waves, where there is no real wave breaking and where wave steepness and slope angle do not influence the run-up. [5] The coefficients and are presented in the table below:
Run-up exceedance probability level, p% | a | b | c | d |
---|---|---|---|---|
0.1 | 1.12 | 1.34 | 0.55 | 2.58 |
2 | 0.96 | 1.17 | 0.46 | 1.97 |
5 | 0.86 | 1.05 | 0.44 | 1.68 |
10 | 0.77 | 0.94 | 0.42 | 1.45 |
50 | 0.47 | 0.60 | 0.34 | 0.82 |
The values of the coefficients highlight the considerable variability of the run-up level from one wave to another in irregular seas. [6] For run-up levels on smooth slopes, work in the Netherlands by the Technische Adviescommissie voor de Waterkeringen (English: Technical Advisory Committee on water Defences) in 1974 discussed the reduction in run-up due to different types of surface roughness. [7]
In practical scenarios, waves are irregular, consisting of a combination of waves with varying heights, periods, and directions. These waves are typically analysed using statistical methods and spectral analysis, providing a more accurate representation of the actual conditions faced by coastal structures and shorelines. Consequently, it is not possible to define a single wave run-up value. Instead, a wave run-up with a specific probability of exceedance is used, typically set at 2%. This wave run-up represents the height exceeded by 2% of the waves in a wave field.
Research indicates that wave run-up follows a Rayleigh distribution, similar to the waves themselves. A probability of exceedance value has been chosen that is small enough to prevent overtopping waves from causing damage to the inner slope. The 2% value has been adopted internationally and was arbitrarily selected by the Dutch Waterloopkundig Laboratorium shortly before 1940. Considering the function, 1% or 5% could have also been possible. The choice of 2% was based on the duration of experimental designs, as a complete trial could be conducted in half a day.
In 1972, Jurjen Battjes, commissioned by the Dutch Technical Advisory Committee for Flood Defences, summarised the available research and provided a solid theoretical foundation. [8] [9] This work led to an improved version of Hunt's Formula, which explicitly included parameters for the angle of incidence of the waves, the effect of a berm, and the slope's roughness. However, the available experimental data on roughness and the berm were insufficient to establish a definitive formula.
Subsequent research was conducted in the following years, with an emphasis on wave overtopping as a more indicative factor for dike height than wave run-up. This research ultimately resulted in a Technical Report in 2002 by the Dutch organisation, TAW. [10] The wave run-up formula mentioned in this report remains in use, and the EurOtop manual has adopted it. The scope of validity has been further expanded in the EurOtop manual, featuring modified formulas. [11]
The EurOtop manual provides a general formula (Formula 1.4 in the manual) [11] for wave run-up:
with a maximum value around 3. The Iribarren number is then used based on the period determined using the first negative moment of the wave spectrum. Additionally, is the reduction coefficient for the factors described below.
The following equation is valid:
in which:
Reference type | [10] |
---|---|
Concrete | 1.0 |
Asphalt | 1.0 |
Concrete blocks | 1.0 |
Grass | 1.0 |
Basaltic columns | 0.90 |
Basalton | 0.90 |
Hydroblock, Haringmanblok | 0.90 |
Open stone asphalt | 0.90 |
Hillblock, Ronataille | 0.70 - 0.80 |
Single layer of armourstone | 0.70 |
Double layer of armourstone | 0.55 |
A range is provided for Hillblock and Ronataille materials, as their reduction coefficient is dependent on wave height. A similar phenomenon occurs with grass. When subjected to high waves, natural grass becomes very smooth, resulting in a reduction coefficient . However, for smaller waves — approximately 25 centimetres (9.8 in) or less — natural grass tends to be much rougher. In such cases, one may opt for a reduction coefficient of below 1.0.
When dealing with short-crested waves, the highest run-up is caused by head-on waves (those with an angle of incidence, ) and equates to the run-up for long-crested waves. For smooth slopes, the run-up decreases slightly with . The run-down typically ranges between a third and a half of the run-up. For breakwaters and revetments constructed with rock armour, the maximum run-down level may indicate the minimum downward extension of the primary armour, and a potential upper level for introducing a berm with a smaller armour size. [5]
For wave run-down there is a similar formula: [12]
Following storm events, a layer of floating debris, known as the flood mark or flotsam, often remains on the slope. This tide mark indicates the maximum wave run-up during the preceding storm. As the flood mark is situated near the height of maximum wave run-up and water levels are generally well-documented by nearby tide stations, it is straightforward to calculate the Ru2% of the storm by subtracting the observed storm surge level from the flood mark level.
In the past, authorities in the Netherlands systematically recorded these observations for most dikes, resulting in a collection of flood mark heights for each dike section. The statistics of flood mark heights can be utilised to determine dike height, which should comprise the design water level plus a safety height (freeboard). The freeboard at the design water level must be equal to the maximum permissible wave run-up. [13]
For a dike with an acceptable load exceedance probability per year, such as 1/500 (as with the temporary dike reinforcement in the Oosterschelde), it is necessary to determine the 1/500 wave run-up. This can be calculated if the 1/500 wave height at the toe of the dike is known. However, this value is rarely measured and must be determined using a computational model, such as SWAN. [14] In many instances, this process can be challenging and prone to errors.
By analysing flood mark heights, which involves simply plotting the data on logarithmic paper, it is possible to directly obtain values such as the 1/500 wave run-up, and consequently the required safety height. An example of this can be observed in the accompanying photo of the run-up and flood mark lines at a dike along the Bathpolder in Zeeland. The photo shows two flood mark lines, which represent the wave run-up of two subsequent storms (on October 12, 2009, with water levels at 2.2 metres (7 ft 3 in) and 1.9 metres (6 ft 3 in) above mean sea level) in the Bathpolder. In the foreground, there is a slope with Haringman blocks, while the background features a slope of Elastocoast. The wave height during these storms was approximately 0.5 metres (1 ft 7.7 in). The wave run-up was 0.8 and 0.85 metres (2 ft 7 in and 2 ft 9 in) above the storm surge level on the Elastocoast, and 1.05 and 1.1 metres (3 ft 5 in and 3 ft 7 in) above the storm surge level on the Haringman blocks. The slope gradient here is 1:4.2. As a Haringman block measures exactly 50 centimetres (19.7 in), the run-up can be assessed in this photo. Subsequent analysis reveals that the reduction coefficient γf for Haringman blocks here is 1.0, and for Elastocoast, it is 0.8.
To assess the safety of a dike and the durability of its grass cover, particularly on the sea or river side, a wave run-up simulator can be employed. The wave conditions for which a dike is designed are infrequent, and the strength of grass coverings varies. These dike conditions can be replicated in-situ using a wave run-up simulator, allowing the manager of the relevant flood defence system to determine if the grass cover is strong enough to withstand expected waves under extreme conditions. [11]
During these tests, the wave run-up simulator is placed on the outer slope and continuously filled with water at a constant flow rate. The flaps at the bottom of the simulator can be opened to varying extents, enabling the simulation of different wave run-up volumes. [15]
The wave run-up simulator is one method for assessing the strength of the grass cover. Another approach involves using a sod puller, which can determine the tensile strength of a sod and allows conversion of this tensile strength by an engineer into a strength under load caused by wave run-up. [16] In addition to simulating wave run-up, the simulation of wave impacts and wave overtopping can be achieved using specifically designed generators and simulators.
Wave run-up should not be confused with wave set-up (an increase in water level due to known waves) or with wind setup (storm surge, an increase in water level due to the driving force of wind).
The Afsluitdijk is a major dam and causeway in the Netherlands. It was constructed between 1927 and 1932 and runs from Den Oever in North Holland province to the village of Zurich in Friesland province, over a length of 32 kilometres (20 mi) and a width of 90 metres (300 ft), at an initial height above Amsterdam Ordnance Datum of between 6.7 metres (22 ft) along the section at Friesland, and 7.4 metres (24 ft) where it crosses the deep channel of the Vlieter. The height at the greater sea depths west of Friesland was required to be a minimum of 7 metres everywhere when originally constructed.
In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.
In differential geometry, a spray is a vector field H on the tangent bundle TM that encodes a quasilinear second order system of ordinary differential equations on the base manifold M. Usually a spray is required to be homogeneous in the sense that its integral curves t→ΦHt(ξ)∈TM obey the rule ΦHt(λξ)=ΦHλt(ξ) in positive re-parameterizations. If this requirement is dropped, H is called a semi-spray.
In mathematics, the parabolic cylinder functions are special functions defined as solutions to the differential equation
In spectroscopy, the Autler–Townes effect, is a dynamical Stark effect corresponding to the case when an oscillating electric field is tuned in resonance to the transition frequency of a given spectral line, and resulting in a change of the shape of the absorption/emission spectra of that spectral line. The AC Stark effect was discovered in 1955 by American physicists Stanley Autler and Charles Townes.
Expected shortfall (ES) is a risk measure—a concept used in the field of financial risk measurement to evaluate the market risk or credit risk of a portfolio. The "expected shortfall at q% level" is the expected return on the portfolio in the worst of cases. ES is an alternative to value at risk that is more sensitive to the shape of the tail of the loss distribution.
The lateral earth pressure is the pressure that soil exerts in the horizontal direction. It is important because it affects the consolidation behavior and strength of the soil and because it is considered in the design of geotechnical engineering structures such as retaining walls, basements, tunnels, deep foundations and braced excavations.
In financial mathematics, tail value at risk (TVaR), also known as tail conditional expectation (TCE) or conditional tail expectation (CTE), is a risk measure associated with the more general value at risk. It quantifies the expected value of the loss given that an event outside a given probability level has occurred.
In probability theory and statistics, the skew normal distribution is a continuous probability distribution that generalises the normal distribution to allow for non-zero skewness.
In fluid dynamics, a cnoidal wave is a nonlinear and exact periodic wave solution of the Korteweg–de Vries equation. These solutions are in terms of the Jacobi elliptic function cn, which is why they are coined cnoidal waves. They are used to describe surface gravity waves of fairly long wavelength, as compared to the water depth.
The uncertainty theory invented by Baoding Liu is a branch of mathematics based on normality, monotonicity, self-duality, countable subadditivity, and product measure axioms.
In fluid dynamics, wave setup is the increase in mean water level due to the presence of breaking waves. Similarly, wave setdown is a wave-induced decrease of the mean water level before the waves break. For short, the whole phenomenon is often denoted as wave setup, including both increase and decrease of mean elevation. This setup is primarily present in and near the coastal surf zone. Besides a spatial variation in the (mean) wave setup, also a variation in time may be present – known as surf beat – causing infragravity wave radiation.
In mathematical physics, the Wu–Sprung potential, named after Hua Wu and Donald Sprung, is a potential function in one dimension inside a Hamiltonian with the potential defined by solving a non-linear integral equation defined by the Bohr–Sommerfeld quantization conditions involving the spectral staircase, the energies and the potential .
In fluid dynamics, the Iribarren number or Iribarren parameter – also known as the surf similarity parameter and breaker parameter – is a dimensionless parameter used to model several effects of (breaking) surface gravity waves on beaches and coastal structures. The parameter is named after the Spanish engineer Ramón Iribarren Cavanilles (1900–1967), who introduced it to describe the occurrence of wave breaking on sloping beaches. The parameter used to describe breaking wave types on beaches; or wave run-up on – and reflection by – beaches, breakwaters and dikes.
Taylor–von Neumann–Sedov blast wave refers to a blast wave induced by a strong explosion. The blast wave was described by a self-similar solution independently by G. I. Taylor, John von Neumann and Leonid Sedov during World War II.
The Van der Meer formula is a formula for calculating the required stone weight for armourstone under the influence of (wind) waves. This is necessary for the design of breakwaters and shoreline protection. Around 1985 it was found that the Hudson formula in use at that time had considerable limitations. That is why the Dutch government agency Rijkswaterstaat commissioned Deltares to start research for a more complete formula. This research, conducted by Jentsje van der Meer, resulted in the Van der Meer formula in 1988, as described in his dissertation. This formula reads
Guderley–Landau–Stanyukovich problem describes the time evolution of converging shock waves. The problem was discussed by G. Guderley in 1942 and independently by Lev Landau and K. P. Stanyukovich in 1944, where the later authors' analysis was published in 1955.
Wave overtopping is the time-averaged amount of water that is discharged per structure length by waves over a structure such as a breakwater, revetment or dike which has a crest height above still water level.
Taylor–Maccoll flow refers to the steady flow behind a conical shock wave that is attached to a solid cone. The flow is named after G. I. Taylor and J. W. Maccoll, whom described the flow in 1933, guided by an earlier work of Theodore von Kármán.
General reference