Impulsive waves caused by subaerial landslides

This paper presents the experimental results of impulsive waves caused by subaerial landslides. A wide range of effective parameters are considered and studied by performing 120 laboratory tests. Considered slide masses are both rigid and deformable. The effects of bed slope angle, water depth, slide impact velocity, geometry, shape and deformation on impulse wave characteristics have been inspected. The impulse wave features such as amplitude, period and also energy conversation are studied. The effects of slide Froude number and deformation on energy conversation from slide into wave are also investigated. Based on laboratory measured data an empirical equation for impulse wave amplitude and period have been presented and successfully verified using available data of previous laboratory works.     

Hydrodynamic and ecological assessment of nearshore restoration: A modeling study

Along the Pacific Northwest coast, much of the estuarine habitat has been lost over the last century to agricultural land use, residential and commercial development, and transportation corridors. As a result, many of the ecological processes and functions have been disrupted. To protect and improve these coastal habitats that are vital to aquatic species, many projects are currently underway to restore estuarine and coastal ecosystems through dike breaches, setbacks, and removals. Understanding site-specific information on physical processes is critical for improving the success of such restoration actions. In this study, a three-dimensional hydrodynamic model was developed to simulate estuarine processes in the Stillaguamish River estuary, where restoration of a 160-acre parcel through dike setback has been proposed. The model was calibrated to observed tide, current, and salinity data for existing conditions and applied to simulate the hydrodynamic responses to two restoration alternatives. Model results were then combined with biophysical data to predict habitat responses within the restoration footprint. Results showed that the proposed dike removal would result in desired tidal flushing and conditions that would support four habitat types on the restoration footprint. At the estuary scale, restoration would substantially increase the proportion of area flushed with freshwater (<5 ppt) at flood tide. Potential implications of predicted changes in salinity and flow dynamics are discussed relative to the distribution of tidal marsh habitat.     

Trophic food web and ecosystem attributes of a water reservoir of the Ria Formosa (south Portugal)

Using a top-down modelling approach, a balanced trophic food web model of a water reservoir of the Ria Formosa lagoonal system was constructed. Few adjustments of input data were necessary to run the model since most of the parameters were estimated at the study site and the characteristics of the water reservoir allowed a high degree of control and precision in terms of sampling and data collection. Trophic levels of the 14 compartments included in the model varied between 1.0 for primary producers and detritus to 3.4 for carnivorous fish and the 14 groups were aggregated in a food chain with six trophic levels. The water reservoir has a detritus based food chain, with the majority of the biomass concentrated in the first two levels, the producers level and the herbivore/detritivore level (97.6%). The transfer efficiencies were low, and decreased with increasing level number, varying between 6.0 and 0.2%. The degree of “ecosystem maturity” was difficult to establish, but several parameters indicate that it could be high. The water reservoir studied, which has similar environmental and ecological characteristics as the Ria Formosa lagoon, is near its carrying capacity.     

Modelling the motion of an internal solitary wave over a bottom ridge in a stratified fluid

A series of laboratory experiments has been carried out to investigate the passage of an internal solitary wave of depression over a bottom ridge, in a two-layer fluid system for which the upper and lower layer is linearly-stratified and homogeneous respectively. Density, velocity and vorticity fields induced by the wave propagation over the ridge have been measured simultaneously at three locations, namely upstream, downstream and over the ridge crest, for a wide range of model parameters. Results are presented to show that wave breaking may occur for a sufficiently large wave amplitude and a strong ridge blockage factor, with accompanying mixing and overturning. Density field data are presented (i) to illustrate the overturning and mixing processes that accompany the wave breaking and (ii) to quantify the degree of mixing in terms of the wave and ridge parameters. For weak encounters, good agreement is obtained between the laboratory experimental results (velocity and vorticity fields induced by the wave propagation) and the predictions of a recently-developed fully nonlinear theory. Discrepancies between theory and experiment are discussed for cases in which breaking and mixing occur.     

Experiments on mixing and dissipation in internal solitary waves over two triangular obstacles

A series of laboratory experiments was undertaken in a stratified two-layer fluid to investigate the energetics of the interaction between an internal solitary wave (ISW) and triangular obstacles, as well as to determine the partitioning of ISW energy and its subsequent dynamics. The ISW energy was dissipated as a result of internal breaking and turbulent mixing induced by wave instability. Tests involving different combinations of triangular obstacles in various heights and intervals and ISW of different amplitudes were performed. The wave features resulting from the interaction of an ISW and double obstacles were found to differ from those of single obstacle. The incident energy of an ISW was either reflecting back from the obstacles, dissipated through turbulent mixing, or transmitted over the double obstacles. Reduction in wave energy increased as the intervals between obstacles reduced. For two obstacles in different heights, energy dissipation was greater in the case with a higher obstacle ahead of a lower one. However, the overall performance was dependent on the relative height of the obstacles, relative water depth of the upper and bottom layer, in addition to the intervals between the obstacles.     

Where is the best site for wave energy exploitation? Case study along the coast of northern Latium (ITALY)

The present work proposes a siting process for the detection of a suitable site for wave energy exploitation. The choice of a suitable site is based on the good agreement between energy availability, environmental sustainability, and equipped facilities to exploit wave energy. The case study in the northern Latium coast is explicative because in this area there are several activities that affect marine ecosystems, and the introduction of renewable energies promote the reduction of anthropic pressures. The nearshore wave power is studied through the numerical wave model (CMS-Wave), and available wave buoys along the coast were used to compare numerical results. In correspondence with Civitavecchia harbour, the largest nearshore wave energy was found; the large depth in front of the breakwater allows to conserve a great part of offshore wave power, with an average dissipation rate of 10 % less than offshore, with mean annual available wave energy of 25.4 MWh/m and seasonal fluctuation of 5.4 MWh/m. This area appears to be an optimal site for nearshore and shoreline wave energy device tests and installations, for energy availability (intermediate level respect Mediterreanean Sea), low potential environmental impact, easier accessibility, and policy oriented towards a larger sustainability of harbour activities.     

Salt-marsh erosion and restoration in relation to flood protection on the Wadden Sea barrier island Terschelling

This paper explores the impact of erosion and restoration measures on habitat development and on wave damping by a small salt marsh nestled alongside a dike on the Wadden island of Terschelling. The aim is to advance knowledge about the benefits and possible side-effects of salt-marsh restoration. Analysis of a time series of aerial photographs from 1944 to 2010 indicates that the salt marsh decreased steadily in size after maintenance of accretion works was terminated. In the western part of the marsh, which is accessible to sheep, vegetation is low (5–15 cm) and dominated by Salicornia europaea and by Spartina anglica. In the most intensively grazed parts, vegetation is very scarce. The eastern, inaccessible part of the salt marsh is covered by dense patches of the shrubby perennial Atriplex portulacoides and Spartina anglica (15–25 cm in height). SWAN wave models show that wave height at this location is significantly affected by the areal extent of the salt marsh as well as by the vegetation. High or dense vegetation are in the models nearly as effective in damping waves (with an initial height of 0.15 and 0.5 m) as widening the salt-marsh area by 350 m. A low density of low plants, as observed in the grazed part of the marsh, has almost no wave-damping effect. Even under conditions of sea level rise, a broader salt marsh vegetated with high plants significantly affects modelled wave height. Therefore, salt-marsh restoration is an adaptation measure worth exploring, though an array of effect types must be considered.     

Ecology of sandy beaches in Oman

The benthic macrofauna and physical features of 10 sandy beaches along the coast of Oman were surveyed quantitatively. This is a mesotidal regime mostly subject to low to moderate wave energy but more exposed in the south. Five northern beaches are tide-dominated, with low wave energy, and their profiles consist of a berm, a steep, swash-dominated upper shore and a broad tide-dominated terrace from mid-shore downwards. They are composed of moderately sorted fine to medium sands. Southern beaches experience greater wave energy, particularly during the summer southwest monsoon, and exhibit smoother, concave profiles with fine, fairly well sorted carbonate sand. 58 species and species groups were recorded, with crustanceans, polychaetes and molluscs dominant. In general species richness was high, at least 19–25 species per beach, but dry biomass moderate to low at 26–90g/m shoreline, with one high value of 450g/m. Total abundance was moderate at 3–73×10³ organisms/m of beach. Some zonation was evident with ocypodid crabs andTylos in the supralittoral, cirolanid isopods on the upper shore and a variety of species on the lower shore. The coast of Oman appears to constitute a single zoogeographic region, but with some regional differentiation between north and south due to varying physical conditions. Thus, OmanÕs beaches are characterized by tide-dominated morphodynamics and exceptionally high species richness.     

Impact of rubble mound groyne structural interventions in restoration of Koggala lagoon,Sri Lanka; numerical modelling approach

Physical processes of the lagoon are influenced by structural interventions. Understanding the complex reality of physical processes sometimes difficult with field observations thus a model provides a simplified abstract view. Two dimensional hydrodynamic model is used to describe, restoration efforts to Koggala lagoon, a combined freshwater and estuarine complex of rich ecosystem on the southern coast of Sri Lanka. The lagoon mouth was naturally closed by a sand bar which controlled the seawater intrusion. Due to large-scale sand removal at lagoon mouth, formation of the sandbar shifted towards the lagoon. After the removal of natural sand barrier, rubble mound groyne structures were built to avoid sand deposition in the lagoon and to protect the highway bridge from the sea wave attack. Construction of the groyne resulted in the lagoon mouth being permanently open which in turn led to many environmental problems with saline intrusion. The aim of this study is to evaluate the current situation of the lagoon and propose alternative structural interventions for minimization of seawater intrusion and subsequently improve lagoon ecosystem. Hydrological parameters were investigated and mathematical models for hydrodynamic behavior of the lagoon were applied in order to describe the lagoon physical processes and flow characteristics. Existing rubble mound structures were redesigned in order to minimize the seawater intrusion. Numerical simulations were carried out for two different mouth widths (40 m and 20 m) with appropriate structural interventions. Existing salting factor for the lagoon is 0.68 and numerical simulation results showed salting factor for 40 m and 20 m openings are 0.61 and 0.54 respectively. This shows the mouth width can be reduced up to 20 m in order to obtain a slating factor close to 0.5, which indicates the predominant influence of fresh water which in turn leads lagoon to a fresh water ecosystem.     

Effects of climate change on surface water management of Seyhan basin,Turkey

The goal of this study was to set-up the basis for climate change adaptation of water resources management policies in Seyhan River basin. The first priority was to identify the balances between water resources and water users with respect to existing and planned projects. In this respect various aspects of Seyhan basin were evaluated, including evaluation of existing water resources, determination of water demand of existing and planned projects, and water resources supply-demand characteristics. The global climate change model was downscaled to the basin scale, the results were associated with hydrometeorological monitoring network and finally the impact of climate change on surface water resources and demands were determined for specific projection years. Water resources management scenarios were developed to evaluate adaptation alternatives to climate change scenarios at the basin level. It was determined that even though there was no water stress in Seyhan basin in 2010, many parts of the basin were expected to suffer significant shortages over the coming years.     

Tsunami generation by submarine landslides: comparison of physical and numerical models

An idealised two-dimensional laboratory model of tsunamis generated by submarine landslides is described. The experimental configuration corresponds to the benchmark configuration suggested by other researchers in the international tsunami community. It comprises a semi-elliptical rigid landslide with a height to length ratio of 0.052 sliding down a 15° slope. The initial landslide submergence and specific gravity are varied, the second of which primarily determines the initial landslide acceleration. In these experiments the landslide motion is generally well approximated as consisting of two periods of constant acceleration. The first phase of positive acceleration finishes as the landslide reaches the base of the slope, while the second period of a slower deceleration continues until the landslide comes to rest along the horizontal base of the tank. A novel experimental technique, which utilises laser-induced fluorescence (LIF), is employed to measure the free surface displacement over the entire space and time domains. This enables the wave potential energy field to be computed directly and provides a vivid picture of the wave generation and development process. Particle tracking velocimetry provides detailed information on the landslide motion and also some data on the sub-surface velocity field. Experimental runs require multiple repeats (typically 35–50) of the same setup in order to capture the entire wave field with the desired resolution. Thus high level experimental repeatability is required, and this is demonstrated. A range of parameters relevant to hazard management are presented and discussed. Maximum crest and trough amplitudes of the offshore propagating waves are shown to be approximately proportional to the initial landslide acceleration and somewhat less strongly dependent on the initial landslide submergence. The maximum wave run-up experienced at the shoreline is shown to depend almost linearly on the magnitude of a high deceleration that occurs for a short period when the landslide nears the toe of the slope. The initial submergence and initial acceleration do not directly determine the maximum wave run-up, although for these experiments they impact indirectly on the magnitude of the deceleration. The efficiency of the energy transfer from the landslide potential energy to the wave field potential energy reaches values of up to 6% and is found to be strongly dependent on the initial submergence. However because of the link between the landslide mass and its acceleration, this efficiency is almost completely independent of the initial acceleration. The results from a numerical model based on linear, inviscid and irrotational wave theory, and solved with the boundary element method, are compared with the data from the experimental program. The numerical model accurately produces the generated sequence of wave crests and troughs, but slightly overpredicts their phase speed by between 2 and 4%. For all other parameters the numerical model predictions are within 25% of the experimental values, although this includes both under- and overprediction for the range of independent parameters covered.     

Waves of intermediate length through an array of vertical cylinders

We report a semi-analytical theory of wave propagation through a vegetated water. Our aim is to construct a mathematical model for waves propagating through a lattice-like array of vertical cylinders, where the macro-scale variation of waves is derived from the dynamics in the micro-scale cells. Assuming infinitesimal waves, periodic lattice configuration, and strong contrast between the lattice spacing and the typical wavelength, the perturbation theory of homogenization (multiple scales) is used to derive the effective equations governing the macro-scale wave dynamics. The constitutive coefficients are computed from the solution of micro-scale boundary-value problem for a finite number of unit cells. Eddy viscosity in a unit cell is determined by balancing the time-averaged rate of dissipation and the rate of work done by wave force on the forest at a finite number of macro stations. While the spirit is similar to RANS scheme, less computational effort is needed. Using one fitting parameter, the theory is used to simulate three existing experiments with encouraging results. Limitations of the present theory are also pointed out.     

A numerical model of wave- and current-driven nutrient uptake by coral reef communities

We developed a numerical model capable of simulating the spatial zonation of nutrient uptake in coral reef systems driven by hydrodynamic forcing (both from waves and currents). Relationships between nutrient uptake and bed stress derived from flume and field studies were added to a four-component biogeochemical model embedded within a three-dimensional (3-D) hydrodynamic ocean model coupled to a numerical wave model. The performance of the resulting coupled physical-biogeochemical model was first evaluated in an idealized one-dimensional (1-D) channel for both a pure current and a combined wave-current flow. Waves in the channel were represented by an oscillatory flow with constant amplitude and frequency. The simulated nutrient concentrations were in good agreement with the analytical solution for nutrient depletion along a uniform channel, as well as with existing observations of phosphate uptake across a real reef flat. We then applied this integrated model to investigate more complex two-dimensional (2-D) nutrient dynamics, firstly to an idealized coral reef-lagoon morphology, and secondly to a realistic section of Ningaloo Reef in Western Australia, where nutrients were advected into the domain via alongshore coastal currents. Both the idealized reef and Ningaloo Reef simulations showed similar patterns of maximum uptake rates on the shallow forereef and reef crest, and with nutrient concentration decreasing as water flowed over the reef flat. As a result of the cumulative outflow of nutrient-depleted water exiting the reef channels and then being advected down the coast by alongshore currents, both reef simulations exhibited substantial alongshore variation in nutrient concentrations. The coupled models successfully reproduced the observed spatial-variability in nitrate concentration across the Ningaloo Reef system.     

Identification of pollution sources in rivers using a hydrodynamic diffusion wave model and improved Bayesian-Markov chain Monte Carlo algorithm

● A hydrodynamic-Bayesian inference model was developed for water pollution tracking. ● Model is not stuck in local optimal solutions for high-dimensional problem. ● Model can estimate source parameters accurately with known river water levels. ● Both sudden spill incident and normal sewage inputs into the river can be tracked. ● Model is superior to the traditional approaches based on the test cases. Water quality restoration in rivers requires identification of the locations and discharges of pollution sources, and a reliable mathematical model to accomplish this identification is essential. In this paper, an innovative framework is presented to inversely estimate pollution sources for both accident preparedness and normal management of the allowable pollutant discharge. The proposed model integrates the concepts of the hydrodynamic diffusion wave equation and an improved Bayesian-Markov chain Monte Carlo method (MCMC). The methodological framework is tested using a designed case of a sudden wastewater spill incident (i.e., source location, flow rate, and starting and ending times of the discharge) and a real case of multiple sewage inputs into a river (i.e., locations and daily flows of sewage sources). The proposed modeling based on the improved Bayesian-MCMC method can effectively solve high-dimensional search and optimization problems according to known river water levels at pre-set monitoring sites. It can adequately provide accurate source estimation parameters using only one simulation through exploration of the full parameter space. In comparison, the inverse models based on the popular random walk Metropolis (RWM) algorithm and microbial genetic algorithm (MGA) do not produce reliable estimates for the two scenarios even after multiple simulation runs, and they fall into locally optimal solutions. Since much more water level data are available than water quality data, the proposed approach also provides a cost-effective solution for identifying pollution sources in rivers with the support of high-frequency water level data, especially for rivers receiving significant sewage discharges.     

Near-boundary velocity and turbulence in depth-varying stream flows

This research examined the temporal distribution of turbulent structure near a streambank toe through the progression of a flood wave in West Run (Morgantown, WV, USA). Three-dimensional velocities and water depths were measured through a 17-h flood event. Turbulence characteristics were examined: Reynolds stresses, turbulent kinetic energy, and turbulence intensities. On average, near-boundary velocity during the rising stage was less than the falling stage, likely due to the measurement location and local roughness. The velocity vectors shifted from towards bed before the flood wave to toward the streambank during progression of the flood wave. Turbulent kinetic energy increased with increasing water depth during the rising stage. Reynolds stress, τ_xz, increased with increasing water depth during the rising stage, but the majority of the stresses were negative through the storm event. Reynolds stress, τ_xy, was positive throughout the event and did not vary with depth. This work is among the first to evaluate turbulence during depth-varying flows in the field.     

Effects of irrigation water qualities on biomass and sugar contents of sugar beet and sweet sorghum cultivars

An experiment involving four qualities of irrigation water two sugar beet and three sweet sorghum cultivars was conducted in a split plot design with four replications at Rudasht Drainage and Reclamation Experiment Station in 1999. The results showed salinity of water has an adverse effect on sugar beet and sweet sorghum biomass. Sweet sorghum cultivar SSV108 had the lowest biomass under all qualities of irrigation water Sweet sorghum cultivar Rio had the maximum biomass with water qualities of 2, 5, and 8 dS m(-1). Sugar beet cultivar 7233 had the maximum biomass with 11 dS m(-1). The effect of irrigation water quality was not significant for sugar characteristics such as brix, pol and purity. However, responses of cultivars on the above parameters were significant and sugar beet cultivars had higher brix, pol and purity and lower invert sugar and starch than sweet sorghum cultivars. In conclusion, sweet sorghum cultivars are not recommended to be irrigated with saline water of more than 8 dS m(-1) for sugar production. Under such condition, they may be suitable to be grown for forage purposes.     

Dike-break induced flows: a simplified model

A simplified model for the prediction of the steady-state outflow through a breach in an inland dike is presented. It consists in the application of the mass and momentum conservation principles to a macroscopic control volume. A proper definition of the shape of the control volume enables to take the main characteristics of the flow into account and thus to compensate for the extreme simplification of the spatial representation of the model. At the breach, a relation derived from the shallow-water equations is used to determine the direction of the flow. Developments have been guided by numerical simulations and results have been compared to experimental data. Both the accuracy and the domain of validity of the simplified model are found satisfactory.     

Assessing the effect of long-term crop cultivation on distribution of Cd in the root zone

Trace elements such as cadmium (Cd) may be inadvertently added to cropland soils through application of fertilizers, irrigation water, and other amendments. These toxic trace elements pose a potentially threat to soil quality and, through the food chain transfer, to human health. A generalized soil trace element mass balance model that accounts for the interactive processes governing the reactions of trace elements in soils, and consequently removed with crop harvest and leaching out of the soil profile with irrigation water was developed in this research. The model conceptually approximates the mechanisms and kinetics of a real field cropland system. The model was used to evaluate the long-term cultivation on distribution of Cd in California croplands. Under typical California cropping practices, Cd added into the soils accumulated primarily in the plow layer while the Cd content below the plow layer was barely affected. After 100 years of continuous cultivation, the soil Cd content of the plow layer increases from the background level 0.22 mg kg⁻¹ to 0.40 mg kg⁻¹. The accumulation of Cd in the plow layer is in proportion to the external inputs and is affected by the soil and plant characteristics, and management practices. The model can be used to evaluate the environmental fates of other toxic element in soils with case specific parameters.     

Numerical study of wave–current–vegetation interaction in coastal waters

Research on interactions among wave, current, and vegetation has received increasing attention. An explicit depth-averaged hydrodynamic model coupled with a wave spectral model (CMS-wave) was proposed in this study in order to simulate the wave and wave-induced current in coastal waters. The hydrodynamic model was based on the finite volume method while the intercell flux was computed by employing the Harten–Lax–van Leer approximate Riemann solver to investigate the dry-to-wet interface, and the drag force of vegetation was modeled as the sink terms in the momentum equations. The CMS-wave model was used to investigate the non-breaking and the breaking random waves propagation in vegetation fields. Afterwards, an empirical wave energy dissipation term with plant effect was derived to represent the resistance induced by aquatic vegetation in the wave-action balance equation. The established model was calibrated and validated with both the experimental and field data. The results showed that the wave height decreased significantly along the wave propagation direction in the presence of vegetations. The sensitivity analysis for the plant density, the wave height, and the water depth were performed by comparing the numerical results for the wave height attenuation. In addition, wave and wave-induced current through a finite patch of vegetation in the surf zone were investigated as well. The strong radiation stress gradient could be produced due to the variation of the energy dissipation by vegetation effect in the nearshore zone, which impacted the direction and amplitude of the longshore current. The calculated results showed that the coupling model had good performance in predicting wave propagation and the current over vegetated water regions.