Coastal defence through wave farms

The possibility of using wave farms for coastal defence warrants investigation because wave energy is poised to become a major renewable in many countries over the next decades. The fundamental question in this regard is whether a wave farm can be used to reduce beach erosion under storm conditions. If the answer to this question is positive, then a wave farm can have coastal defence as a subsidiary function, in addition to its primary role of producing carbon-free energy. The objective of this work is to address this question by comparing the response of a beach in the face of a storm in two scenarios: with and without the wave farm. For this comparison a set of ad hoc impact indicators is developed: the bed level impact (BLI), beach face eroded area (FEA), non-dimensional erosion reduction (NER), and mean cumulative eroded area (CEA); and their values are determined by means of two coupled models: a high-resolution wave propagation model (SWAN) and a coastal processes model (XBeach). The study is conducted through a case study: Perranporth Beach (UK). Backed by a well-developed dune system, Perranporth has a bar between − 5 m and − 10 m. The results show that the wave farm reduces the eroded volume by as much as 50% and thus contributes effectively to coastal protection. This synergy between marine renewable energy and coastal defence may well contribute to improving the viability of wave farms through savings in conventional coastal protection.     

A hybrid efficient method to downscale wave climate to coastal areas

Long-term time series of sea state parameters are required in different coastal engineering applications. In order to obtain wave data at shallow water and due to the scarcity of instrumental data, ocean wave reanalysis databases ought to be downscaled to increase the spatial resolution and simulate the wave transformation process. In this paper, a hybrid downscaling methodology to transfer wave climate to coastal areas has been developed combining a numerical wave model (dynamical downscaling) with mathematical tools (statistical downscaling). A maximum dissimilarity selection algorithm (MDA) is applied in order to obtain a representative subset of sea states in deep water areas. The reduced number of selected cases spans the marine climate variability, guaranteeing that all possible sea states are represented and capturing even the extreme events. These sea states are propagated using a state-of-the-art wave propagation model. The time series of the propagated sea state parameters at a particular location are reconstructed using a non-linear interpolation technique based on radial basis functions (RBFs), providing excellent results in a high dimensional space with scattered data as occurs in the cases selected with MDA. The numerical validation of the results confirms the ability of the developed methodology to reconstruct sea state time series in shallow water at a particular location and to estimate different spatial wave climate parameters with a considerable reduction in the computational effort.     

A probabilistic method for constructing wave time-series at inshore locations using model scenarios

Continuous time-series of wave characteristics (height, period, and direction) are constructed using a base set of model scenarios and simple probabilistic methods. This approach utilizes an archive of computationally intensive, highly spatially resolved numerical wave model output to develop time-series of historical or future wave conditions without performing additional, continuous numerical simulations. The archive of model output contains wave simulations from a set of model scenarios derived from an offshore wave climatology. Time-series of wave height, period, direction, and associated uncertainties are constructed at locations included in the numerical model domain. The confidence limits are derived using statistical variability of oceanographic parameters contained in the wave model scenarios. The method was applied to a region in the northern Gulf of Mexico and assessed using wave observations at 12 m and 30 m water depths. Prediction skill for significant wave height is 0.58 and 0.67 at the 12 m and 30 m locations, respectively, with similar performance for wave period and direction. The skill of this simplified, probabilistic time-series construction method is comparable to existing large-scale, high-fidelity operational wave models but provides higher spatial resolution output at low computational expense. The constructed time-series can be developed to support a variety of applications including climate studies and other situations where a comprehensive survey of wave impacts on the coastal area is of interest.     

Beach morphology and shoreline evolution: Monitoring and modelling medium-term responses (Portuguese NW coast study site)

Numerical models for shoreline evolution have been used for coastal management planning for several decades. The model calibration is a start point to project shoreline scenarios and in this aim the use of data acquired within the scope of monitoring programmes provides the opportunity to assess the models' capabilities under real condition. This work applies calibration data (retrieved from field surveys) to numerical models to predict medium-term shoreline evolution using, as a case study, a beach stretch named AC, about 3.5 km long and located downdrift of a groin on the northwest Portuguese coast. A smaller stretch AB (2.4 km long), included in the total one, which exhibits a pronounced erosive tendency usually better reproduced in shoreline evolution models, was also analysed. Based on topographic surveys, associated wave climate conditions registered between 2003 and 2008 and typical wave conditions registered over a longer wave climate time period, this work compares the calibration of two different shoreline evolution models, Long-term Configuration (LTC) and GENESIS for this period. Then, considering the 2003 topographic conditions for the models' calibration, the results of both models are discussed with respect to simulation scenarios after 10, 15 and 20 years of evolution. The 10-year evolution projections of the models are also compared to the results of a survey performed in February 2012. For the wave data calibration period (2003–2008), the average shoreline retreat of the analysed coastal stretch was reproduced with small differences (around 1% and 10% for LTC and 15% and 14% for GENESIS, considering stretches AB or AC, respectively), though local differences along the AB coastal stretch represent root mean square errors reaching up to 52% and 88% for GENESIS and LTC, respectively, and were above 118% for both models along the AC coastal stretch.     

How does multiscale modelling and inclusion of realistic palaeobathymetry affect numerical simulation of the Storegga Slide tsunami?

The ∼8.15 ka Storegga submarine slide was a large (∼3000 km³), tsunamigenic slide off the coast of Norway. The resulting tsunami had run-up heights of around 10–20 m on the Norwegian coast, over 12 m in Shetland, 3–6 m on the Scottish mainland coast and reached as far as Greenland. Accurate numerical simulations of Storegga require high spatial resolution near the coasts, particularly near tsunami run-up observations, and also in the slide region. However, as the computational domain must span the whole of the Norwegian-Greenland sea, employing uniformly high spatial resolution is computationally prohibitive. To overcome this problem, we present a multiscale numerical model of the Storegga slide-generated tsunami where spatial resolution varies from 500 m to 50 km across the entire Norwegian-Greenland sea domain to optimally resolve the slide region, important coastlines and bathymetric changes. We compare results from our multiscale model to previous results using constant-resolution models and show that accounting for changes in bathymetry since 8.15 ka, neglected in previous numerical studies of the Storegga slide-tsunami, improves the agreement between the model and inferred run-up heights in specific locations, especially in the Shetlands, where maximum run-up height increased from 8 m (modern bathymetry) to 13 m (palaeobathymetry). By tracking the Storegga tsunami as far south as the southern North sea, we also found that wave heights were high enough to inundate Doggerland, an island in the southern North Sea prior to sea level rise over the last 8 ka.     

Rockslide tsunamis in complex fjords: From an unstable rock slope at Åkerneset to tsunami risk in western Norway

An unstable rock volume of more than 50 million m³ has been detected in the Åkerneset rock slope in the narrow fjord, Storfjorden, Møre & Romsdal County, Western Norway. If large portions of the volume are released as a whole, the rockslide will generate a tsunami that may be devastating to several settlements and numerous visiting tourists along the fjord. The threat is analysed by a multidisciplinary approach spanning from rock-slope stability via rockslide and wave mechanics to hazard zoning and risk assessment.The rockslide tsunami hazard and the tsunami early-warning system related to the two unstable rock slopes at Åkerneset and Hegguraksla in the complex fjord system are managed by Åknes/Tafjord Beredskap IKS (previously the Åknes/Tafjord project). The present paper focuses on the tsunami analyses performed for this company to better understand the effects of rockslide-generated tsunamis from Åkerneset and Hegguraksla. Two- and three-dimensional site-specific laboratory experiments are conducted to study the generation, propagation, and run-up of the wave for several potential rockslide scenarios from Åkerneset. Furthermore, the two models GloBouss and DpWaves are applied for numerical simulations of the generation/propagation phase and a third model MOST is applied for numerical simulations of the near-shore propagation and inundation of the wave in selected locations. Strong emphasis is put on verification, validation, and sensitivity of the numerical models. The best match between the numerical simulations and the laboratory experiments is found for the larger scenarios with the linear dispersive solution for the propagation phase; the corresponding calculated run-up values are remarkably similar to the ones observed during the laboratory experiments.During the risk assessment it was found that the rockslide tsunami hazard (probability of impact) is higher than accepted by the Norwegian Planning and Building Act. This should at that time prevent any further development in all the exposed areas of the entire fjord system. The Act is today altered to open for specified further development in the various hazard zones. The results of the tsunami analyses are applied in risk management in terms of hazard map production and land-use planning. Two failure scenarios for each of the two unstable rock slopes are designed for the hazard zoning. The larger and less probable scenarios (1 in 5000 years) are applied for evacuation zones and routes, while the smaller and more probable scenarios (larger than 1 in 1000 years) are applied for location and design of less critical facilities accepted in the inundation zone.     

Wave climate variability in the North-East Atlantic Ocean over the last six decades

Ocean surface gravity waves play a major role in many engineering and environmental problems, both in the open ocean and in coastal zones. Therefore, it is essential to improve our knowledge on spatial and temporal variability of wave climate. This study aims at investigating this variability in the North-East Atlantic Ocean (25°W–0°W and 30°N–60° N), using a 57-year hindcast (1953–2009) obtained with a spectral wave model forced with reanalysis wind fields. The hindcast analysis reveals firstly strong seasonal fluctuations of wave climate, with winters characterized by large and long-period waves of mean direction spreading from south-west to north-west, and summers characterized by smaller and shorter-period waves originating from norther directions. From northern (55°N) to southern (35°N) latitudes, the significant wave height (Hs) decreases by roughly 40%, the mean wave direction (Mwd) rotates clockwise by about 25% while the peak period (Tp) only grows by 5%. These three parameters also exhibit a strong inter-annual variability, particularly when winter-means (from 1st of December to 1st of April) are considered. Linear trend analysis over the studied period shows spatially variable long-term trends, with a significant increase of Hs (up to 0.02 m yr^?1) and a counterclockwise shift of Mwd (up to ?0.1° yr^?1) at northern latitude, contrasting with a fairly constant trend for Hs and a clockwise shift of Mwd (up to +0.15° yr^?1) at southern latitudes. Long-term trends of Tp are less significant, with still a slight increase in the north-eastern part of the study area (up to +0.01 s yr^?1). Eventually, a comparison between the inter-annual variability of the winter-means of the three selected wave parameters and the North Atlantic Oscillation (NAO) reveals: (1) a strong positive correlation between Hs and the NAO index at northern latitudes (correlation coefficient up to R = 0.91) and a significant negative correlation at southern latitudes (up to R = ?0.6); (2) no significant correlation for Mwd north of 40°N and a clear positive correlation southward of 40°N (up to R = 0.8) and (3) a northward increasing correlation for Tp (up to R = 0.8). Long-term trends for Hs, Mwd and Tp are finally explained by a significant increase in the NAO index over the studied period.     

Application of a fuzzy inference system for the prediction of longshore sediment transport

A fuzzy inference system (FIS) and a hybrid adaptive network-based fuzzy inference system (ANFIS), which combines a fuzzy inference system and a neural network, are used to predict and model longshore sediment transport (LST). The measurement data (field and experimental data) obtained from Kamphuis [1] and Smith et al. [2] were used to develop the model. The FIS and ANFIS models employ five inputs (breaking wave height, breaking wave angle, slope at the breaking point, peak wave period and median grain size) and one output (longshore sediment transport rate). The criteria used to measure the performances of the models include the bias, the root mean square error, the scatter index and the coefficients of determination and correlation. The results indicate that the ANFIS model is superior to the FIS model for predicting LST rates. To verify the ANFIS model, the model was applied to the Karaburun coastal region, which is located along the southwestern coast of the Black Sea. The LST rates obtained from the ANFIS model were compared with the field measurements, the CERC [3] formula, the Kamphuis [1] formula and the numerical model (LITPACK). The percentages of error between the measured rates and the calculated LST rates based on the ANFIS method, the CERC formula (K_sig = 0.39), the calibrated CERC formula (K_sig = 0.08), the Kamphuis [1] formula and the numerical model (LITPACK) are 6.5%, 413.9%, 6.9%, 15.3% and 18.1%, respectively. The comparison of the results suggests that the ANFIS model is superior to the FIS model for predicting LST rates and performs significantly better than the tested empirical formulas and the numerical model.     

The influence of sea state on formation speed of alongshore variability in surf zone sand bars

The formation time of alongshore morphological variability in surf zone sand bars has long been known to differ from one beach to the other and from one post-storm period to another. Here we investigate whether the type of sea state, i.e. distant swell waves or locally generated short period wind sea, affects the formation time of the emerging alongshore topographic variability.A numerical modeling approach is used to examine the emergence of alongshore variability under different shore-normal wave forcing. A research version of Delft3D, operating on the time-scale of wave groups, is applied to a schematised bathymetry with a single bar. The model is then used to investigate several wave scenarios, examining the impact of peak period, frequency spread and directional spread on the formation time of alongshore variability.Results show that an increase in wave period has a large effect, changing the formation time up to O (250%) in case the wave period is changed from a representative value for the Dutch coast (T_p ~ 5–6 s) to an Australian South East coast value (T_p ~ 10–12 s). In contrast, modifications in the directional and frequency spread of the wave field result only in a minor change in the formation time.Examination of hydrodynamics and potential sediment transport shows that the variations in formation time are primarily related to changes in the magnitude of the time-averaged flow conditions. Variations in the magnitude of very low frequency (f < 0.004 Hz) or infragravity (0.004 < f < 0.04 Hz) surf zone flow velocities do not affect the mean sediment transport capacity. Consequently the formation speed of patterns is primarily governed by positive feedback between mean flow and morphology, and low frequency flow fluctuations are of minor importance.These findings indicate that the development of alongshore topographic variability may be faster at swell dominated open coasts, primarily due to the occurrence of longer period swell. Also, at a given site, the arrival of a long wave period swell after a storm can accelerate the emergence of variability.     

The Relative Trough Froude Number for initiation of wave breaking: Theory,experiments and numerical model confirmation

It is important to accurately locate the wave breaking region for the calculation of nearshore hydrodynamics. Energy from breaking waves drives hydrodynamic phenomena such as wave set-up, set-down, wave run-up, longshore currents, rip currents, and nearshore circulation. Numerous studies have been undertaken to describe when and where wave breakings occurs. Recent development of computer resources permits the use of phase-resolving numerical models for the study of wave propagation, transformation, and nearshore hydrodynamics. This requires new types of wave breaking criteria for the numerical model. The Relative Trough Froude Number (RTFN) is a new wave breaking criterion. This model is based on the moving hydraulic jump concept, therefore it satisfies properly posed boundary-value conditions. It has been experimentally proved that a critical RTFN at the initiation of wave breaking is consistent with and without the presence of an opposing current, but previous efforts did not investigate the theory for the critical value. This paper provides a theoretical analysis and a numerical analysis to demonstrate why the RTFN theory works as a wave breaking initiation (trigger) index. The theoretical analysis provides a universal constant for the initiation of wave breaking for all water depths assuming the Miche formula properly describes the wave breaking condition. A subroutine for wave breaking in a numerical model, FUNWAVE was modified to include the RTFN trigger. The numerical model was calibrated with data from wave tank experiments, and it was found that the critical condition is very close to the theoretical number, CTFN = 1.45. A second paper (in preparation) provides details of the theory and experiments for a second criterion for termination of wave breaking. The time scale for the establishment of the breaking region i.e., between the initiation position and termination position, depends upon the additional momentum present under turbulent condition within the breaking wave. This subject is not considered herein.     

Wave climate scatter performance of a cycloidal wave energy converter

A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power P_W = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of P_G = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate.     

Assessment of tsunami hazard for coastal areas of Shandong Province,China

Shandong province is located on the east coast of China and has a coastline of about 3100 km. There are only a few tsunami events recorded in the history of Shandong Province, but the tsunami hazard assessment is still necessary as the rapid economic development and increasing population of this area. The objective of this study was to evaluate the potential danger posed by tsunamis for Shandong Province. The numerical simulation method was adopted to assess the tsunami hazard for coastal areas of Shandong Province. The Cornell multi-grid coupled tsunami numerical model (COMCOT) was used and its efficacy was verified by comparison with three historical tsunami events. The simulated maximum tsunami wave height agreed well with the observational data. Based on previous studies and statistical analyses, multiple earthquake scenarios in eight seismic zones were designed, the magnitudes of which were set as the potential maximum values. Then, the tsunamis they induced were simulated using the COMCOT model to investigate their impact on the coastal areas of Shandong Province. The numerical results showed that the maximum tsunami wave height, which was caused by the earthquake scenario located in the sea area of the Mariana Islands, could reach up to 1.39 m off the eastern coast of Weihai city. The tsunamis from the seismic zones of the Bohai Sea, Okinawa Trough, and Manila Trench could also reach heights of >1 m in some areas, meaning that earthquakes in these zones should not be ignored. The inundation hazard was distributed primarily in some northern coastal areas near Yantai and southeastern coastal areas of Shandong Peninsula. When considering both the magnitude and arrival time of tsunamis, it is suggested that greater attention be paid to earthquakes that occur in the Bohai Sea. In conclusion, the tsunami hazard facing the coastal area of Shandong Province is not very serious; however, disasters could occur if such events coincided with spring tides or other extreme oceanic conditions. The results of this study will be useful for the design of coastal engineering projects and the establishment of a tsunami warning system for Shandong Province.     

Focused wave evolution using linear and second order wavemaker theory

In this paper, the evolution of focused waves using different paddle displacements (piston type) under laboratory conditions is presented. It is well known that in intermediate water depths, linear paddle displacements will generate spurious, free, sub and super harmonics. Thus, a second order correction to suppress these spurious free sub and super harmonics was used to generate the focused waves. The focused waves were generated in the laboratory using a linear superimposition principle, in which the wave paddle displacement is derived based on the sum of a number of sinusoidal components at discrete frequencies, whose phases are accordingly set to focus at a particular location. For this method of generation, the second order wave maker theory proposed by Schäffer [24] can be easily adopted and was used in the present study. Two different centre frequencies (f_c = 0.68 Hz and 1.08 Hz) with three different bandwidth ratios (Δf/f_c = 0.5, 0.75 and 1.0) were tested in a constant water depth, to consider both narrow and broadband spectra. These test cases correspond to wave focusing packets propagating in intermediate and deep water regions. Further, for each wave packet, two different amplitudes were considered in order to analyze non-breaking and breaking cases. Thus, by systematically generating the wave packets using the linear and second order paddle displacements, the analysis was carried out for the spectral and temporal evolution of selected long waves. The temporal evolution of the selected harmonics was analyzed using the Inverse Fast Fourier Transform (IFFT), to show the propagation of the spurious, free, long waves. Further, the variations in energy for the lower, higher and primary frequency ranges are reported for different wave paddle displacements. The analysis revealed that for the broadband spectrum the differences are more pronounced when using linear paddle displacements. We have also noticed a shift in focusing/breaking location and time (i.e. premature) due to the increase in crest height using linear displacements. The experiment data used in this paper has been provided as a supplementary, which can be used to validate the numerical models.     

Acoustic propagation analysis with a sound speed feature model in the front area of Kuroshio Extension

This paper aims to analyse acoustic-propagation character in the front area of Kuroshio Extension (KE). By analysing Argo data and the Sea surface height (SSH) data in this KEF area, a two-dimensional (2D) sound-speed feature model (SSPFM) characterising the KEF is proposed. The SSPFM has a transition zone with a width about 100 km and the sound channel changes from 1000 m south of KEF to 300 m north of KEF, resulting in a sharp gradient about 7 m/km. Along with the meandering character of the KEF axis, the sharp gradient results in a rather complicated acoustic environment in the KEF area. With reanalysis data from the hybrid coordinate ocean model, a three-dimensional (3D) sound-speed environment is established. The acoustic propagation character in the KEF area is then analysed with the 2D SSPFM and the 3D acoustic environment. Results show that the KEF affects acoustic propagation mainly by modifying the sound channel depth. Given that acoustic propagation in the KEF area is influenced mainly by the meandering KEF, with the near-real-time SSH data to locate the KEF, the 2D SSPFM is able to provide a near-real-time estimate of the underwater 3D acoustic environment.     

Offshore wind farm impacts on surface waves and circulation in Eastern Lake Ontario

A coupled wave and hydrodynamic model was applied to the Kingston Basin of eastern Lake Ontario, a region with bathymetric variability due to channels and shoals, to assess the potential impacts on surface waves and wind-driven circulation of an offshore wind farm. The model was used to simulate a series of storm events with time-varying wind forcing and validated against wave, current and water level observations. The wind farm was simulated by adding semi-permeable structures in the surface wave model to represent the turbine monopiles, and by adding an energy loss term to the fluid momentum equations in the hydrodynamic model to represent the added drag of the monopiles on the flow. The results suggest that the wind farm would have a small influence on waves and circulation throughout the wind farm area, with spatial variability due to focussing of wave energy and re-direction of the flow. Overall, the results indicate that the wave height in coastal areas will be minimally affected with changes in significant wave height predicted to be < 3%. Larger changes to the strength of circulation occur inside the wind farm region with localized changes in current magnitude of up to 8 cm s^− 1. The results of this study may help to understand the impacts of future offshore wind farms and other offshore structures in the Great Lakes.     

Numerical simulations of wave–current flow in an ocean basin

Wave–current flow is a phenomenon that is present in many practical engineering situations. Over the past several decades, this type of flow has been increasingly investigated under controlled laboratory conditions. This paper presents a numerical study of wave–current flow in the ocean basin of the LabOceano (COPPE/UFRJ). A homogeneous multiphase model based on the RANS equations and the k–ɛ turbulence model implemented in ANSYS-CFX code were used. A cross section of the ocean basin was represented. A regular wave with a height of 0.08 m and a period of 1.80 s (i.e., a wave steepness of H/L = 0.016), propagating on favourable currents, was simulated. The behaviour of the free surface elevation over time and the streamlines along the basin for wave and wave–current flows were presented. The numerical results were compared to the non-viscous theory given by the Rayleigh equation applied to the problem of wave–current interaction. Good agreement was found between the wave length estimated by the numerical results and the analytical solutions, with a deviation of less than 2%.     

Surface wave predictions in weakly nonlinear directional seas

We have employed laboratory and numerical experiments in order to investigate propagation of waves in both long and short-crested wave fields in deep water. For long-crested waves with steepness, ϵ = k_ca_c = 0.1 (a fairly extreme case), reliable prediction can be performed with the modified nonlinear Schrödinger equation up to about 40 characteristic wavelengths. For short-crested waves the accuracy of prediction is strongly reduced with increasing directional spread.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

Simulation of storm surge,wave, and coastal inundation in the Northeastern Gulf of Mexico region during Hurricane Ivan in 2004

Hurricane-induced storm surge, waves, and coastal inundation in the northeastern Gulf of Mexico region during Hurricane Ivan in 2004 are simulated using a fine grid coastal surge model CH3D (Curvilinear-grid Hydrodynamics in 3D) coupled to a coastal wave model SWAN, with open boundary conditions provided by a basin-scale surge model ADCIRC (Advanced CIRCulation) and a basin-scale wave model WW3 (WaveWatch-III). The H1wind, a reanalysis 10-m wind produced by the NOAA/AOML Hurricane Research Division (HRD), and a relatively simple analytical wind model are used, incorporating the effect of land dissipation on hurricane wind. Detailed comparison shows good agreement between the simulated and measured wind, waves, surge, and high water marks. Coastal storm surge along the coast is around 2–3 m, while peak surge on the order of 3.5 m is found near Pensacola, which is slightly to the east of the landfall location on Dauphin Island. Wind waves reach 20 m at the Mobile South station (National Data Buoy Center buoy 42040) on the shelf and 2 m inside the Pensacola/Escambia Bay. Model results show that wave-induced surge (total surge subtracted by the meteorologically-induced surge due to wind and pressure) accounts for 20–30% of the peak surge, while errors of the simulated surge and waves are generally within 10% of measured data. The extent of the simulated inundation region is increased when the effects of waves are included. Surge elevations simulated by the 3D model are generally up to 15% higher than that by the 2D model, and the effects of waves are more pronounced in the 3D results. The 3D model results inside the Pensacola/Escambia Bay show significant vertical variation in the horizontal currents. While the estuary has little impact on the surge elevation along the open coastal water, surge at the head of Escambia Bay is more than 50% higher than that at the open coast with 1.5 h delay.     

Conceptual hydrodynamic-thermal mapping modelling for coral reefs at south Singapore sea

Coral reefs are important ecosystems that not only provide shelter and breeding ground for many marine species, but can also control of carbon dioxide level in ocean and act as coastal protection mechanism. Reduction of coral reefs at Singapore coastal waters (SCW) region remains as an important study to identify the environmental impact from its busy industrial activities especially at the surrounding of Jurong Island in the south. This kind of study at SCW was often being related to issues such as turbidity, sedimentation, pollutant transport (from industry activities) effects in literatures, but seldom investigated from the thermal change aspect. In this paper, a computational model was constructed using the Delft3D hydrodynamic module to produce wave simulations on sea regions surrounding Singapore Island. The complicated semi-diurnal and diurnal tidal wave events experienced by SCW were simulated for 2 weeks duration and compared to the Admiralty measured data. To simulate the thermal mapping at the south Singapore coastal waters (SSCW) region, we first adapted a conversion of industrial to thermal discharge; then from the discharge affected area a thermal map was further computed to compare with the measured coral map. The outcomes show that the proposed novel thermal modelling approach has quite precisely simulated the coral map at SSCW, with the condition that the near-field thermal sources are considered (with the coverage area in the limit of 20 km × 20 km).