1 SAND BAR BEHAVIOR: OBSERVATION S AND MODELING By LUCIANO ABSALONSEN A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILO SOPHY UNIVERSITY OF FLORIDA 2012
2 2012 Luciano Absalonsen
3 To my m om d ad grandmother and s ister
4 ACKNOWLEDGMENTS I would like to thank Capes and Fulbright for providing the scholarship that allowed me to come to UF to pursue my PhD. Also the committee members that are taking their time to read, criticize and give suggestions to improve the quality of this dissertation. All the friends (Andy, Berkay, Chloe, Cihan, Leslie, Miao, Shih Feng,Tianyi, Uriah and other coastal fellows) that share d their culture experience and knowledge and made these four years in Gainesville more interesting. To Amy, for her help, incentive and good friendship during this difficult journey pursuing the PhD. Special acknowledgement to my advis e r, Dr. Robert Dean, who was always available to help me with all my questions regarding the dissertation (and there were many) and valuable suggestions and guidance He ha s also been a great inspiration as a person for his great character and passion for coastal engineer ing I hope that I was able to assimilate at least 1 % of his qualities during this period that I was his student Um outro agradecimento especial para minha fam lia (Me, Pai, V e Lela) que sempre me apoia em todos os momentos (a maior parte de alegria) e i ncentiva nesta minha vida de estudante profissional. Sem vocs, com certeza, isto no seria possvel.
5 TABLE OF CONTENTS page ACKNOWLEDGMENTS ................................ ................................ ................................ ...... 4 LIST OF TAB LES ................................ ................................ ................................ ................ 7 LIST OF FIGURES ................................ ................................ ................................ .............. 8 LIST OF ABBREVIATIONS ................................ ................................ .............................. 11 ABSTRACT ................................ ................................ ................................ ........................ 12 CHAPTER 1 INTRODUCTION ................................ ................................ ................................ ........ 14 2 LITERATURE REVIEW ................................ ................................ .............................. 16 2.1 Beach Erosion ................................ ................................ ................................ ...... 16 2.2 Equilibrium Beach Profile ................................ ................................ ..................... 17 2.3 Differences Between Winter and Summer Profiles ................................ ............. 19 2.4 Conceptual Models ................................ ................................ ............................... 22 2.4.1 Conceptual Model for Systems with Multiple Sand Bars ........................... 25 2.4.2 Field Obs ervations of the Conceptual Model with Multiple Bars ............... 26 2.4.3 Conceptual Models Application for Coastal Management ........................ 28 2.5 Breaki ng Characteristics ................................ ................................ ...................... 29 2.5.1 Wave Breaking Types ................................ ................................ ................. 30 2.5.2 Beach Slope and Breaking Point ................................ ............................... 31 2.5.3 Cross Shore Transport Associated with Different Break Types ............... 33 2.6 Influence of Beach Profile on Wave Setup ................................ .......................... 35 2.7 Sand Bar Formation Theories ................................ ................................ .............. 36 2.8 Wave Height Decay inside the Surf Zone ................................ ............................ 38 3 MOTIVATION ................................ ................................ ................................ .............. 40 4 FIELD STUDIES ................................ ................................ ................................ ......... 42 4.1 Beach Profile Data ................................ ................................ ................................ 42 4.2 Definition of Sand Bars ................................ ................................ ......................... 42 4.3 Wave Data ................................ ................................ ................................ ............ 45 4.4 Description of Walton and Volusia Counties ................................ ....................... 46 4.4.1 Walton County ................................ ................................ ............................ 47 4.4.2 Volusia County ................................ ................................ ............................ 49 5 QUALITATIVE MODEL: WAVE HEIGHT AND NUMBER OF SAND BARS RELATIONSHIP ................................ ................................ ................................ .......... 58
6 5.1 Wave Analysis ................................ ................................ ................................ ...... 58 5.2 Average Number of Bars and Average Profile ................................ .................... 59 5 .3 Bar Response to Average Wave Energy ................................ ............................. 60 5.3.1 Walton County ................................ ................................ ............................ 60 5.3.2 Volusia County ................................ ................................ ............................ 61 5.4 Summary of Findings from the Qualitative Examination ................................ ..... 62 6 NUMERICAL MODEL TO PREDICT THE EVOLUTION OF SAND BARS IN THE NEARSHORE ZONE ................................ ................................ ......................... 67 6.1 Model Description ................................ ................................ ................................ 67 6.2 Comparison with Florida Data ................................ ................................ .............. 71 6.3 Walton and Volusia County Groups ................................ ................................ ..... 72 6.4 Walton County ................................ ................................ ................................ ...... 73 6.4.1 Comparison between the Groups ................................ .............................. 73 6.4.2 Sand Bars versus Wave Height ................................ ................................ 74 6.5 Volusia County ................................ ................................ ................................ ...... 76 6.6 Hurricane Season Effects ................................ ................................ ..................... 77 7 SHORT TERM PREDICTION OF THE NUMBER OF SAND BARS ........................ 89 7.1 Dally, Dean and Dalrymple Model and Modifications ................................ ......... 90 7.2 Role of Sand Bars in the Breaking Process ................................ ........................ 91 7.3 Prediction of the Number of Sand Bars ................................ ............................... 93 8 CONCLUSIONS ................................ ................................ ................................ ........ 102 LIST OF REFERENCES ................................ ................................ ................................ 104 BIOGRAPHICAL SKETCH ................................ ................................ .............................. 110
7 LIST OF TABLES Table page 4 1 Number of profiles with zero, one, two and three bars and the average number of bars in Walton and Volusia Counties in all surveys from 1972 to 2007. ................................ ................................ ................................ ....................... 46 4 2 Results of the LI applied to all surveys in Walton County, evaluating the inner and outer bars linearity ................................ ................................ .......................... 49 4 3 Results of the LI applied to all surveys in Volusia County, evaluating the inner and outer bars linearity ................................ ................................ .......................... 50 6 1 Dates that surveys were conducted at Walton County ................................ ......... 73 6 2 Statist ics and coefficient values for Walton County. ................................ ............. 75 6 3 Statistics and coefficient values for Volusia County. ................................ ............ 77 7 1 Results obtain ed with the model simulations for four surveys in Walton County. ................................ ................................ ................................ .................. 94 7 2 Results obtained with the model simulations for four surveys in Volusia County. The values indicate the simulations tha t correctly predicted, over predicted and under predicted the number of bars ................................ .............. 95
8 LIST OF FIGURES Figure page 2 1 Schematic representation of wave b reakpoint in dissipative and Reflective Beaches. ................................ ................................ ................................ ................. 39 4 1 Location of the two study areas in the State of Florida. Walton County in the west coast and Volusia County in the east coast. ................................ ................. 51 4 2 Criteria utilized to define the sand bars inside the three zones in the beach profile.. ................................ ................................ ................................ .................... 52 4 3 The 13 profiles of Group 1 at Walto n County during the 1997 survey (black) and the average profile (red). ................................ ................................ ................ 53 4 4 Location of the buoys 42039 and 41009 used to extract the wave data for Walton and Volusia Counties, respectively. ................................ .......................... 54 4 5 Beach survey before Hurricane Ivan. Profiles with one sand bar (blue crosses) and profiles with two sand bars (diamonds). The profiles with two bars have the inner bars represented wit h blue diamond and the outer bar with red diamond.. ................................ ................................ ................................ .. 55 4 6 Beach survey after Hurricane Ivan. Profiles with one sand bar (blue crosses) and profiles with two sand bars (diamonds). The profiles with two bars have the inner bars represented with blue diamond and the outer bar with red diamond. ................................ ................................ ................................ ................. 56 4 7 S and bar depths and widths from all profiles in Walton (red circle) and Volusia (blu e cross) Counties. Calculated linear trend for both Counties (straight lines). ................................ ................................ ................................ ........ 57 5 1 Average beach profiles at Walton and Volusia Counties representing the differences in beach profile slope. ................................ ................................ ......... 63 5 2 Measured profile examples in Volusia County. A) Profile with one bar. B) Profile with two bars. C) Profile with three bars. ................................ ................... 64 5 3 Results from the qualitative model for Walton County. A) Comparison between average wave energy and average number of bars. B) Wave energy and indication of Hurricanes .. ................................ ................................ ............... 65 5 4 Res ults from the qualitative model for Volusia County. A) Comparison between average wave energy and average number of bars. B) Wave energy and indication of Hurricanes.. ................................ ................................ ................ 66
9 6 1 Theoretical examp le demonstrating the model assumptions. A) Wave height. B) The model results (red dashed line) and N eq (solid blue line). The periods where the model is adjusting with the N eq are highlighted in gray. ....................... 79 6 2 Expected number of sand bars. A) U sing N eq as a function of Dean Number. B) U sing N eq as a function of Profile Parameter, for Volusia (blue) and Walton (Red) Counties. The coefficients used in the equations are also presented. ...... 80 6 3 Division of the groups with same dates in the seven surveys conducted in Walton County. ................................ ................................ ................................ ....... 81 6 4 Comparison between average number of ba rs from measured data in the five groups at Walton County. ................................ ................................ ...................... 82 6 5 Comparison between different N eq functions (dashed lines) used to calculate the number of bars in the nearshore zone for Walto n County and measured data (red solid line). ................................ ................................ ................................ 83 6 6 Quantitative model results for Walton County. ................................ ...................... 84 6 7 Comparison between the number of bars during the initial nine months of simulation for the five groups in Walton County. ................................ ................... 85 6 8 Comparison between different N eq functions (dashed lines) used to calculate the number of bars in the nearshore zone for Volusia County and measured data (red solid line). The first value is the same for all functions, since it was used as an initial condition. ................................ ................................ .................... 86 6 9 Quantitati ve model results for Volusia County. ................................ ..................... 87 6 10 Seasonal results obtained from the quantitative model. A) Walton County. B) Volusia County. One year periods (black line) and average season behavi or (red line).. ................................ ................................ ................................ ................ 88 7 1 Representation of the breaking and stability criteria (dashed lines) used in the Dally, Dean and Dalrymple model.. ................................ ................................ 97 7 2 Three steps of the simulation showing the sand bar importance on the breaking processes and wave height decay inside the nearshore zone. ............. 98 7 3 Model application to predict th e number of sand bars in the nearshore zone. The wave height (red line), breaking criterion (upper black dashed line), stability criterion (lower black dashed line) and beach profile (black solid line).. ................................ ................................ ................................ ....................... 99
10 7 4 Simulation where the second breakpoint is associated with an area with abrupt change in the profile slope and not to a sand bar. The wave height (red line), breaking criterion (upper black dashed line), stability criterion (lower black dashed line) and beach profile (black solid line). ........................... 100 7 5 Example of the breakpoint associated with a developed bar/trough system in Volusia County. The wave height (red line), breaking crit erion (black dashed line), stability criterion (lower black dashed line) and beach profile .. ................ 101
11 LIST OF ABBREVIATION S FDEP Florida Department of Environmental Protection LIDAR Light D etection and R anging MHW Mean High Water NDBC National Data Buoy Center NOAA National Oceanic and Atmospheric Adm inistration
12 Abstract of Dissertation Presented to the Graduate School of the University of Florida in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy SAND BAR BEHAVIOR : OBSERVATION S AND MODELING By Luciano Absalonsen August 2012 Chair: Robert G. Dean Major: Coastal and Oceanographic Engineering The relation between the number of sand bars in the nearshore zone and the forcing wave conditions w as evaluated using observations and numerical model ing at Walton and Volusia Counties, in the State of Florida. The decadal behavior of the number of bars was well represented using a simple quantitative model that correlates the wave for cing and number of bars T he Dean and Profile non dimensional parameters, used to determine the number of bars in equilibrium with the wave forcing, presented the best correlation between model and observations, indicating that wave and sediment characteri stics are important to predict the number of bars in the nearshore zone. The model succeeded in representing the rapid increase in the number of bars during major storms or hurricanes and the slower decrease in number of bars during calm periods. The break point theory used in the short term model captured the relationship between the inner and outer bar s where the outer bar tend s to be more linear, since only high waves are able to effectively transport sediments in to great er depths. The inner bar morpholo gy is more three dimensional where alongshore differences are more evident. The offshore migration or increase in the outer bar depth was observed to be
13 the trigger for the inner bar formation. The short term model represented the number of sand bars when the surveys were conducted close to periods of high energy when the sand bars were formed or migrated to the position and depth according to the wave forcing that generated the system with multiple bars. It was also observed that the inner and outer bars a re interconnected systems with multiple bars and the study of this system has to be conducted including all bars at the same time.
14 CHAPTER 1 INTRODUCTION Sediment transport inside the nearshore zone is highly dynamic where morphological features and hydr odynamics are intimately associated. Prominent morphological features encountered on sand y beaches are the shoreline, berm and bars. The shoreline is usually the main feature in which coastal managers, stakeholders and the general public are interested. Th e principal reason for th is special interest is that several states and countries have a significant part of their economy directly or indirectly relat ed to the use of this area and shoreline changes, especially associated with beach erosion, can cause sig nificant financial losses Travel and tourism i n the United States is the single n ational largest industry, employer and foreign revenue earner and U.S. beaches are the leading tourist destination (Houston, 1996) Gable (1997) argued that tourism is the m ost important source of external revenue in the Caribbean and that all tourism development s w ere located in the coastal zone. This is not just observed in the USA and Caribbean Islands, but it is a worldwide phenomenon and several other examples for differ ent countries could be show n here. Miami Beach, in the State of Florida, is a good example wh ere as a result of the waves were breaking on coastal structures. Consequently, visitor numbers dec lined and the facilities had deteriorated Following beach nourishment and infrastructure improvements in the from foreign tourists alone is $2.4 billion, about 50 times the $52 million cost of the 20 year project. Miami Beach has more than twice (21 million) the combined number of
15 tourist visits to Yellowstone (2.6 million), the Grand Canyon (4 million) and Yosemite (3.3 million) N ational P arks combined (Houston, 2002). C learly beaches are vital to tourism based economies and careful attention has to be directed to these areas to reduce or avoid problems associated with beach erosion. Nearshore zones are also important for marine and coastal organisms, e.g. turtles have th eir nests at beaches and migratory birds use these areas to rest and to obtain their food. The preservation of this habitat is fundamental for the sustainability of the entire ecosystem. Comprehension of beach dynamics requires monitoring of its natural mo rphological changes through time, according to different forcing time scales Certain knowledge of the natural behavior of every beach type can be obtained by means of a field monitoring program long enough for the beach to attain different energetic state s and in which their efficiency in protection can be comprehended ( Benavente et al., 2002) Understanding the processes governing this system is fundamental to maintaining it attractive to visitors, to diminish the costs associated with coastal management programs and construction of structures or beach nourishment projects to protect landward sites and restore recreational areas. Shoreline and sand bar behaviors are interconnected responding to the changes in waves, tides and currents. S and bars strongly i nfluence the cross shore hydrodynamics and consequently the sediment transport in the nearshore zone. The sand bar behavior inside the nearshore zone their characteristics the interaction between inner and outer bars and their response to wave conditions in different time scales are the main focus of this dissertation.
16 CHAPTER 2 LITERATURE REVIEW It is well known that sediments are transported offshore during condition s of high energy, through wave driven circulation and this sediment migrates onshore du ring periods of mild energy, forming the winter/storm and summer/normal profiles. Torrey Pines Beach, California, is one of several example s of this behavior. The beach undergoes typical seasonal changes in configuration owing to changes in wave climate. During summer wave conditions the beach has a 30 to 60 m wide backshore, a relatively steep foreshore and a pronounced berm. Winter storm waves overtop the summer berm and erode the backshore, thus reducing the exposed beach width. The winter beach is typi fied by a gently sloping beach face that in places exten d shorewa rd to the toe of the sea cliff (Winant et al., 1975) Another example is Duck, North Carolina, one of the most studied beaches in the world, where this cycle was observed to repeat through ma ny years. Longshore uniform sand bars approximately parallel to the shoreline with little bathymetric changes were observed at Duck after an energetic period between August 25 th and October 12 th 1994 (Gallagher et al., 1998) Dynamically equilibrated coa sts with a balance between onshore and offshore sediment migration and predominately cross shore transport have minim al sediment lost and this cycle can repeat over long periods of time. Problems arise when the sediment transport is not balanced usually wi th gradients of alongshore or offshore sediment transport and resulting in beach erosion. 2.1 Beach Erosion The m ain agents responsible for sediment imbalance are commonly considered as natural or anthropogenic causes. Natural causes are frequently associ ated with sea
17 level rise, influencing worldwide beaches o n a global scale and o n the local scale the natural migration of inlets can cause erosion o n adjacent beaches. Anthropogenic causes include the construction of structures perpendicular to the coast ( e.g. groins and jetties), parallel to the coast (e.g. seawalls) on top of the dune system s or during periods where the shoreline was advanced. Dam construction reduce s the amount of sediment available from rivers to the coast. T he High Aswan Dam, built in 1964 in Egypt, diminished the Nile River sediment discharge to near zero and initiated severe beach erosion near the Damietta and Rosetta outlets (Frihy and Dewidar, 2003) Meckel (1975) reported that dams at the Colorado River trapped most of the sedimen ts that otherwise would have been available for delta construction and sediment transport at the coast. Anthropogenic influences usually reduce the amount of sediment available to the coastline and diminish the natural seasonal cycle, creating areas with s carce amount s of sediment and consequently erosion. Natural causes related to inlet natural migration are of reduc ed importance present ly since many inlets are fixed for navigational purposes. An example is observed in the State of Florida, where many in lets are fixed for navigation However, this action result s in erosion downdrift and deposition updrift of the inlets, since part of the longshore sediment transport is blocked by the jettie s Davis and Barnard ( 2003) demonstrated that proximity to inlets can affect shoreline trend rates for several k ilo m eters downdrift of inlets 2.2 Equilibrium Beach Profile One of the first attempts to understand the cross shore behavior of coastal zones was to predict the general beach cross shore configuration, sin ce a wide range of
18 morphologic features were observed during beach surveys, and it appears that different agents are acting on the beach to create them Bruun (1954) defined an equilibrium beach profile as the statistical average profile which maintains i ts form apart from smaller fluctuations including seasonal fluctuations. The theory used to develop this concept is based on constructive and destructive forces that can balance on beaches exposed to same, or more realistically, similar hydrodynamic condit ions for long periods. After the balance, different sediment sizes are moved to different portions of the beach, where they are in equilibrium according to the balance between constructive and destructive forces. The application can be extended when analyz ing long term trends of sea level. Future scenarios with different sea levels can be predicted using this method, since it is based on the average equilibrium forces at the beach and these forces can be estimated by numerical models. Coarser sediments are generally found in shallower areas, since they are exposed to higher energy and there is a decrease in sediment size in deeper zones, where the wave energy is reduced ( Charles et al., 1994 ; Wang et al., 1998). Different sediment sizes and wave conditions were evaluated to determine their importance in shaping the sea bed. The equilibrium beach profile was created combining these conditions and developing a general equation to represent the cross shore beach profile The equilibrium beach profile that rela tes sediment size with distance from the shore using the 2/3 power, was verified with good agreement in beach profiles at Denmark and USA (Bruun, 1954 and Dean, 1977). Application of the equilibrium beach profile to ten beaches in New Zealand with arbitra ry cross shore sediment size suggests that the theory is effective to some degree
19 in reflecting profile characteristics due to sediment size variations (Dean et al., 1993). T his function predicts the general beach profile shape, required as a requisite usi ng the statistical definition, however the main restriction of its application when compar ed to field data is that it is not capable of predict ing sand bars or troughs (Dean and Dalrymple, 2002) To include the presence of sand bars in the equilibr ium b each profile, Inman et al. (1993) divided the beach profile in two segments using the breakpoint to distinguish the onshore and offshore areas. For each segment they proposed a similar function described above and the best fit was found using a 0.4 power, instead of 2/3. However this method can predict only profiles with one bar and a total of seven parameters is predicted using this approach and applicability f o r engineer i ng or scientific purposes is limited due to the large amount of parameters. 2.3 Differences B etween W inter and S ummer P rofiles Winter and summer profiles present remarkable differences in their configuration, morphology and time response to wave forcing. T he time scale s for the creation of winter and summer profiles are considerably different. In the first case, offshore sediment transport and sand bar creation occurs i n a time scale of hours. During hurricanes, or events with large wave heights, beach eros ion, offshore sand migration and sand bar formation occurs in few hours. During extreme events, say a 1 00 or 500 years storm, the sediments can migrate offshore to location s deeper than the depth of closure to region s where the net cross shore sediment t ransport is weak and these sediments may not return to the more active beach Kennedy et al. (2010) verified infilling of seventeen borrow pits in the east
20 and west coast of Florida, during the 2004 2005 hurricane season, only one of which was considered to be landward of the depth of closure. Borrow pit depths are greater than adjacent areas and the forces acting to remove sediments from them are very weak, creating permanent sediment deposits if the sand enter s the borrow pits, and consequently sediment deficit at the beach. T herefore, most of the sediment infilling was the result of sediment transport from outside the depth of closure. Careful analysis is important when selecting borrow area locations for beach nourishment, otherwise beach nourishment lo ngevity can be shorter than predicted since part of the placed material can return to the borrow area during high energy events through cross shore sediment transport. O nshore sediment transport occurs over a considerably lo ng er time scale than offshore t ransport and long periods without significant increase s in wave height are required for sediments to migrate shoreward Sunamura and Takeda (1984) presented one of the first attempts to quantify onshore migration of sand bars using two dimensionless quanti ties based on the offshore wave height, wave period and sediment size. Sand bars onshore migration rate of 3.5 m/day were observed in Tairua Beach, New Zealand, during four days of observations. The rate was slightly higher since it occurred mainly during periods with low tide, when the water level above the bar were smaller than 0.25 m M ost of the material transported was by bedload transport, since the waves were low and the sediment median diameter is 0.6 mm (van Maanen et al., 2008). Other than the differences in response time between winter and summer profiles, there is a remarkable difference in the nearshore morphology I n winter profiles
21 s ediments migrate to a specific depth in the beach profile, associated with the breakpoint and create linear sand bars, with approximately two dimensional morphologies. In the summer profile, o ffshore bar characteristics are more variable and alongshore differences become more evident where bars and beaches show more three dimensional features. Longshore disco ntinuities of two to three km lengths in bar position exist in the Dutch North Sea Coast, which are caused by the crescentic bar morphology (Stive et al., 1996) Bar splitting, which involves a longshore bar developing a forked or bifurcated appearance, wa s observed at Wanganui, New Zealand. After bifurcat ing the seaward segment that is considered the main bar segment migrates farther offshore while the inner segment moves into the landward trough and completely detaches from the original bar (Shand, 2007) Other three dimensional morphologies are transverse bars, formed when horn areas associated with rhythmic morphology weld onto the beach with a capability to transfer sediments onshore when they detach from the inner bar and realign in the longshore curr ent direction (Konichi and Holman, 2000). These features have an extensive variety of length scales, in the west coast of Florida (Franklin County) with average breaker height of 6 cm, their spacing may vary between 64 to 218 m and cross shore length s betw een 107 and 640 m (Niedoroda and Tanner, 1970). A geomorphologic review of intertidal and subtidal bars distinguishing seven nearshore types providing horizontal and vertical scales is presented by Wijnberg and Kroon ( 2002). These different signatures ar e associated with differences in wave energy, nearshore slope and tidal range.
22 C ross shore migration rates can vary in the longshore direction with waves focusing or diverging creating complex morphologies and erosional hot spots can be generated at the shoreline (Ribas et al., 2003). Rip currents and beach cusps are morphodynamic beach characteristics often associated with mild wave conditions. Rip currents are narrow seaward directed jets oriented approximately shore normal, that originate in the surf z one and broaden outside the breaking region, forming the rip head. A review of rip currents is presented by MacMahan et al. (2006), where theories of formation, comparisons between model and observed values, scaling and flow structures are described. Two f eatures that are interesting to emphasize is that the spacing between rip currents is the order of 100 m and they can achieve velocities greater than 2 m/s, being a constant hazard at the beach. I n the State of Florida more people fall victim to rip curren ts than hurricanes and tornadoes (Lascody, 1998). Beach cusps are crescentic morphological features with regular longshore shapes characterized by seaward protruding features (horns) and gently sloped landward embayments located between two horns I n Melbo urne Beach, Florida east coast, the cusps wavelengths varied between 25 m and 45 m (van Gaalen et al., 2011). Th e cycle of rapid offshore bar migration and slow onshore movement are observed in places continuously monitored over periods usually greater than one year. If the system is dynamically stable, without cross shore or longshore sediment transport losses, the sand bars migrate inside the nearshore zone and no shoreline erosion is evident over long periods of time. 2. 4 Conceptual Models M orphodyn amic stud ies of the nearshore zone started with conceptual models classifying different beaches or the same beaches during different periods of the year.
23 Their main aim is to group beaches with similar morphological features (based on the nearshore plan vi ew) and hydrodynamic forcing to facilitate understanding of nearshore morphodynamics. The se classifications are useful for coastal planners, engineers and scientists to quickly identify systems under investigations and initial plans can be developed based on them. However, coastal zones are highly dynamic and it is important to emphasize that these classifications are generalizations and site specific modifications are expected for different coasts. Conceptual models mostly used in the literature to classif y beaches are described below. Wright and Short (1984) developed a model based on observations at 26 beaches in Australia that evolved forming similar features. Their model used the non dimensional parameter developed by Dean (1973) where the breaking wav e height (H b ), wave s ) were correlated according to E quation 2 1. The model presents six stages changing from the most reflective stage (D < 1) through 4 intermediate stages (1 < D < 6) to the most dissipative stage (D > 6). A general decrease in the mean sediment diameter and beach slope occur s going from reflective to dissipative stages. The relation between an increase in sediment size and an increase in beach face slope was previously observed by Wiegel (1950) and Bascom ( 1951) with data from different beaches in USA. Other than the grain size, the degree of beach exposure to incident waves affect the beach slope i n which less exposed beaches have steeper slope for a given median grain size, since more onshore sediment tra nsport occurs in this condition.
24 Nearshore zones with more three dimensional shapes are expected during intermediate stages. These features occur when alongshore differences rhythmic or irregular become more evident. Rip currents are often observed during intermediate stages cutting the sand bars and increasing alongshore heterogeneities. Two dimensional features are evident in more dissipat iv e stages that are characterized by bars which are linear, continuous and relatively parallel to the coast. In the reflective stage nearshore zones without sand bars are expected, where a common morphologic feature is the presence of beach cups. This conceptual model was developed using traditional beach survey methods and visual observations. The traditional survey methods consist of the use of level and rod and later substituted by global positioning systems (GPS) surveys. Both allow nearly instantaneous beach characterization, usually not during high ly energetic conditions, since large wave heights may not allow be ach survey s to be conducted The use of video cameras to monitor nearshore zones introduced a continuous remo t e tool to monitor the nearshore zone. Lip p mann and Holman (1990) analyzed data derived from video cameras installed at Duck and observed similar behavior as Wright and Short They developed another conceptual model with eight stages based on the earlier model proposed by Wright and Short (1984) and subdivid ed two intermediate stages, creating one classification more appropriate to Duck. Both m orphodynamic models described above were developed on microtidal beaches, where waves were the main hydrodynamic forcing, however tide effects are evident in a large portion of coastlines in all latitudes and continents. Davies (1964)
25 proposed a classifica tion according to the tide amplitude, where tidal ranges smaller than 2 m are considered microtides, between 2 and 4 m are mesotides and greater than 4 m are macrotides. Short (1991), reviewed the literature on these environments and pointed out the diffe rences between microtidal and meso macrotidal regions. Masselink and Short (1993) proposed another conceptual model including the tide hydrodynamic effects, using the Dean number and the relative tidal range (that consists of the tidal range divided by H b ) This m odel is not discussed in detail since the areas studied in this dissertation are on the Florida coastline and influenced by microtides, however it is important to keep in mind that the major influence of the tide is to change the breakpoint pos itio n and intensity throughout a cycle, extending breaking processes effects to wider areas inside the nearshore zone. These conceptual models were based only on beaches with one offshore bar in the nearshore zones and multiple sand bars were not considered i n the models. Other conceptual models including the presence of multiple bars are discussed in the next section. 2. 4 .1 Conceptual M odel for S ystems with M ultiple S and B ars Nearshore zones with more than one sand bar are common in many places in the world however few studies are dedicated to these zones, and understanding their behavior is essential for coastal management and application of correct solutions to mitigate coastal problems. A conceptual study at the long barrier island of Terschelling in the Netherlands observed nearshore bar behavior over a time scale of years and concluded that typically a bar goes through three different stages during it s life cycle (Ruessink and Terwin d t, 2000). Generation in the inner nearshore zone in shallow
26 waters, of fshore migration through the surf zone to intermediate depths and finishing the cycle decaying in the outer nearshore zone. The offshore bar migration allows higher waves to reach the inner bar, since the same wave heights can break farther inshore. During this process a new bar is created near the shore with relatively fast offshore migration rate, separating from the shoreline, indicating that the trigger of a new cycle is the outer bar offshore migration and decay. The outer bar may act as an energy fil ter to the inner nearshore zone, determining the local wave conditions at the shoreward located bar. Waves with large heights dissipate part of their energy by breaking on the outer sand bars and smaller waves break after passing the bar. The sheltering pr ocess shows that the morphodynamic behavior inside the nearshore zone in systems with multiple bars is interconnected. Segregating the studies focusing only on one bar and not on the entire system can result in inappropriate interpretations. 2. 4 .2 Field O bservations of the Conceptual Model with Multiple Bars Some locations where the cycle of sand bar formation, offshore migration and decay w as observed using frequent and consecutive surveys are described in the following paragraphs. Pape et al. (2010) ob served this cycle at the Gold Coast (Australia) and Egmond (Netherlands) with periods of 5 and 15 years, respectively. Similarly, at Wanganui (New Zealand) a system with usually two subtidal bars experienced net offshore migration with the mean life cycle of approximately three years (Shand et al., 1999). Kuriyama (2002) described the same cycle at Horns (Japan), however the cycle is even faster at this location, being completed in only one year.
27 Nearshore slope differences is one explanation to changes in the b ar cycle period. L onger cycles are expected o n mild ly sloping beaches, since the bar travels greater distances in the nearshore zone until the breaking waves stop influenc ing their mobility. Other than the trajectory length, Shand et al. (1999) sugge sted that mild ly slop ing beaches tend to be more exposed to high waves, which would generate higher bars, with large volumes, that require long time for seaward migration. The comparison between Hasaki (nearshore slope of 0.008) and results from the Nethe rlands (nearshore slope between 0.004 and 0.008) confirm the mild slope hypothesis, since the cycle in Hasaki is shorter than in the Netherlands. However, according to the nearshore slope, the duration time at Hasaki should be longer than those at Duck (0. 0097) and Wanganui (0.0088) and the opposite is happening. This result indicates that although the nearshore slope influence s the duration time, it is not enough to fully explain these differences. Longshore currents were suggested as another agent contrib uting to changes in the duration time (Kuriyama, 2002) however uncertainties still remain about all agents influencing the cycle duration Inner and outer bars can present different behaviors in time and space in systems with multiple bars. Ruessink and Terwin d t (2000) proposed that the inner bar is more stable than the outer bar, since the most energetic waves are breaking and dissipating their energy o n the outer bar, and only part of the energy is reaching the inner bars. However the opposite behavior was observed at Gold Coast, where the inner bar is more dynamic than the outer bar. It was also observed that the cycle duration can change in the same location. At Duck, it was observed that one bar stayed in the outer
28 bar phase for seven years and the ot her bars spent only two years in the same outer bar phase (Plant et al., 1999). Those conclusions illustrate the complexity of nearshore zones, especially when composed of multiple bars. Conclusive answers for the explanation of the cycle duration and inn er and outer bar behavior are still under investigation, but it is clear that several variables are acting on the system and that these variables often have their effects superimposed and/or interconnected, creating a challenge to full comprehension of thi s zone. Small changes in morphology can quickly grow when currents, waves and sediment strongly interact, like in the nearshore zone of sandy beaches. 2. 4 .3 Conceptual Models Application for Coastal Management Morphodynamic conceptual models provide a desc ription of the qualitative beach morphology and hydrodynamic forcing associated with generat ing different plan form s. However their main focus is not in quantify ing the processes associated with the bars formation or shoreline changes, being restricted to qualitative analysis or the ir combination with quantitative methodologies. Benedet et al. (2004) applied the Wright and Short (1984) model to the east coast of Florida and classified it in terms of four susceptibility ratings (extreme, high, moderate and low). However some modifications were necessary to account for the perched nature of Floridian morphotypes on reefs the impacts of morphodynamic feedbacks from sediment starvation, hardground exposure and bimodal nature of beach sediment grain size distri bution. Benavente et al. (2002 ) used average seasonal values of the Dean Number, Surf Similarity and Surf Scaling Parameters (the last two are described later on ( section 2.5.1 ) to construct a susceptibility matrix for the SW coast of Spain and indicated t hat it is a useful initial tool in coastal hazard stud ies related to storms. They
29 also concluded that seasonal beaches, characterized by an alternation between two or more stages through the year are more resistant to storms than unif orm unchanging beaches where b each seasonality appears to be a key variable in the severity of coastal damage by storms. Physical processes acting inside the nearshore zone must be qualitatively and quantitatively evaluated to fulfill the limitations described by the conceptu al models and achieve better results in the prediction of nearshore processes. 2. 5 Breaking Characteristics The most important processes modulating beach bathymetry is the sediment transport caused by the breaking waves. Morphologic features with differe nt spatial scales can be formed during this process. The wave height inside the nearshore zone, maximum wave height and position where the waves start to break are often important parameters to be determined for coastal and shore protection structures or coastal management plans for specific locations. These values are not only related to the wave characteristics (wave height and period) but also to the profile bathymetry. Therefore, changes in beach slope, sediment size morphologic features and sand bars presence (including quantity in the cross shore profile ) or absence alter the breaking point location and the wave height decay inside the nearshore zone. In some situation s the waves might not break and reflect offshore, due to a combination of differen t factors. The main reason is when seawalls are constructed and the beach is eroded to them, so the waves approaching the shore reach the seawalls directly and are mostly reflected offshore. Additional scour during storms decreasing the sediment amount is observed in front of seawalls (Dean, 1986) On beaches without
30 structure s, more wave reflection occurs o n reflective beaches, where the beach slope is steep due to the coarse sediments. Sediment concentrations and sediment availability to be transported b y suspended or bedload inside the nearshore zone are directly influenced by breaking processes. Understanding this process is fundamental to predict ing the locations where sediments tend to accumulate/spread inside the nearshore zone and associated positio n where sand bar s /trough s are expected to be formed Breaking characteristics also affect the shoreline position, since both are directly related to the cross shore sediment transport and differences in breaking characteristics cause onshore or offshore tr ansport building or eroding the beach 2. 5 .1 Wave Break ing Types Distinct wave breaking types are observed in different beaches or in the same beach but during different periods in the year. It is noticed that some places have waves breaking in wide zone s and at other beaches the waves break in a narrow zone with their impact being restricted to a specific location. It is also evident that these waves change their sinusoidal shape approaching the breakpoint and differences in shape and amount of turbulenc e produced is variable during breaking and broken waves In order to differentiate these breaking processes and classify the waves according to their main characteristics, breakers have been organized as spilling, plunging and surging breakers. A brief des cription of these types provided by Wiegel (1964) is presented below: and a gradual break across the nearshore zone. Plunging breakers curl over the top of
31 the crest and t his mass of water plung es down when the crest front becomes steep and then concave. Surging breakers peak up as if to break in the manner of plunging breaker, but then the wave base surges up the beach face with the resultant disappearance of the collapsin g wave crest (Wiegel, 1964). C ollapsing breaker s are sometimes distinguished as a fourth group and their characteristics are between plunging and surging. Important criteria to differentiate breaker groups are beach slope and wave steepness. Iversen (1952 ) analyzing waves with different offshore steepness concluded that at large values of deep water wave steepness, the breakers on all beach slopes were spilling. At smaller values of deep water steepness the waves tend to plunge with greater tendencies on b eaches with steeper slopes. At the extreme lower values, the breakers tended to surge. Battjes (1974) combined the beach slope and offshore wave steepness effects to create the Surf Similarity Parameter This parameter quantitatively distinguishes the fou r breaker types qualitatively described previously. The Surf Scaling Parameter (Guza and Inman, 1975) is related to the Surf Similarity Parameter, where the same variables were used when it was develop ed 2. 5 .2 Beach Slope and Breaking Point Offshore wav es propagating towards the shore start to shoal in the nearshore zone due to the bathymetry increas ing their height until the point that they become u nstable At this point there is a transition from irrotational to r ota tional motion and the waves break. Beach profile relevance to the breaking process was first observed by McCowa n (1 891 ), using solitary wave theory who found that waves would break when the ratio
32 between wave height and water depth (H/h) exceeds the critical value of 0.78. This value was o bserved later to var y according to the beach slope Iversen (1952) used a wave channel changing the beach slope from 1:10 to 1:50 and a rigid bottom and observed that steep beach profiles have larger critical values than mild profiles. Weggel (1972) analyz ing data from wave tanks with regular, monochromatic and progressive waves also observed that this ratio is greater on steep beaches and the breaking characteristics changes to more plunging breakers. It is expected, by this observation, waves breaking clo ser to the shoreline in reflective beaches, while waves with similar characteristics break farther offshore in more dissipative stages, only considering plane beaches. A schematic representation of the breaking point position in relation to changes in the beach slope is presented on Figure 2 1. Other than the location where waves start to break, the beach slope influences the ratio between offshore wave height and the wave height at the breaking point. Iversen (1952) using a laboratory wave channel with di fferent plane sloping bottoms observed that breakers on the 1:10 slope are approximately 40 % higher than those in the 1:50 slope. He reported that these relatively steep slopes may have produced excessive reflection, particularly for longer waves and cons equently measured erroneous wave heights. Wave s teepness also changes the break point position, where an increase in wave steepness increase s the critical value (0.78) and waves break more onshore. Stive (198 4 ) conducted an e xperiment with two different wav e steepness es in a wave flume using periodic waves with minimum second harmonic on a plane concrete slope of 1:40 and observed that increasing the offshore wave steepness (from 0.01 to 0.032) the
33 breakpoint moved approximately three meters onshore (wave fl ume total length is 55 m). Komar and Gaughan (1972) used l aboratory and field data to develop a relationship between wave steepness and the ratio between the breaking wave height and offshore wave height based on Airy wave theory. An effort to categorize beach profiles between winter and summer, according to the wave steepness using laboratory measurements ha s indicated that the wave steepness transition between summer and winter profiles is approximately 0.025, whereas smaller value results from large sca le tests. For sustained wave steepness conditions greater than this value, a storm beach profile occur s and for sustained wave steepnesses below this value, a normal profile occur s (Dean, 1973). The relationships developed between waves and beach slope con tributed to distinguish beach profiles using only offshore wave conditions. 2. 5 .3 Cross S hore Transport Associated with Different Break Types Dyhr Nielsen and Sorensen (1970) described qualitatively the breakpoint bar model as a combination of undertow, d ue to radiation stress arguments, transporting mobilized sand offshore within the surf zone and converging with sand transported onshore due to the skewness of nonlinear waves just outside the surf zone to form a bar near the wave breaking point. The shape stabilize s when the slope gets so steep that the bottom current cannot transport the sand up over it, or when there is equilibrium between the amount of sand carried up by the bottom current and the sand returned in suspension at higher levels. Large cha nges in concentration occur over times shorter than a wave period and at spatial scales shorter than a wavelength (Gallagher et al., 1998), and different types of breaking waves have the capacity to suspend variable amount s of sediments.
34 Therefore, s edimen t suspension observed in the surf zone is spatially and temporally intermittent and breaking characteristics are important for further sediment transport processes. The heuristic model developed by Dean (1973) assumes that breaking waves suspend the sedime nts in the water column. If the time required for the sediments to settle to the bottom is greater than half of the wave period, the sediments are deposited offshore since the average water velocity is directed offshore. If it is smaller than half of the w ave period the net transport is onshore since the particles are predominantly influenced by onshore velocities. The breaking waves associated with those profiles govern the sediment transport inside the surf zone. Spilling waves, usually associated with m ild beach slopes and wide surf zones, start to dissipate their energy farther offshore than the plunging breakers. The turbulence generated by these different breaking processes causes differences in the suspended sediment concentration, since it is propor tional to the turbulence intensity in the surf zone. Ting and Kirby (1994) observed that during spilling waves the sediments suspended by the breaking waves have large turbulence decay time and only a small portion of the energy is dissipated over one wave period. When the next wave approaches the same position, the sediments are still suspended and have an offshore orbital velocity component, moving offshore and creating erosional profiles with one or more sand bars. Plunging waves, on the other hand, are associated with steep beach profiles and the waves do not dissipate their energy throughout the surf zone. They approach the shore and break in a specific location and create highly turbulent zones, however, turbulence decay occurs over a relatively short time. The sediments
35 stay in suspension for a short period and when the following wave approaches the same location the sediments are already deposited onshore. In summary, plunging waves are breaking processes that contribute to build back the profiles dur ing periods with mild energy and spilling breakers act to transport sediment offshore. Differences in the wave shape influence the current characteristics inside the nearshore zone. The waves approaching the shore change from sinusoidal to more peaked wit h shallow trough shapes, becoming more asymmetric inside the nearshore zone. The effects of varying wave asymmetries is that for small wave asymmetries, the longshore current profile is displaced substantially landward from the current generated without ro llers and the longshore current on the bar trough increases (Lippmann et al., 1995) However the effects of longshore currents are not the main focus of the current research. 2. 6 Influence of B each P rofile on W ave S etup Cross shore bathymetric effects on wave setup are a focus of investigation since it is desired to predict the water level, mainly during high ly energetic events, to determine areas that might be affected during flood s and quantify the risks during these events. Low areas, e.g. Florida, are e specially concern ed about changes in the water level, since small variations can inundate extensive areas and cause severe financial damages. Different slopes and presence of absence of bars cause changes in the wave setup at the coast. Stephen s et al. (2 011) observed through numerical simulations of artificial barred beaches with different shapes that when the bar is deep, substantial wave energy can propagate inshore and dissipate on the beachface, and vice versa when the bar is shallow. Furthermore, a m ean water level generated over barred profiles with steep slope offshore of the bar is greater than over that with gentle slope in
36 the same zone because the wave energy is dissipated more rapidly, producing a faster rate of change in radiation stress. A co mparison between the offshore bar slope and the beach face slope is also relevant. Beachface slope becomes more influential at places where the offshore bar slope has lower gradients than beachface slope. It occurs because the water level relatively increa ses as the bar deepens. Raubenheimer et al. (2001) obtained similar conclusions using numerical model ing and observations at Duck, suggesting that setup near the shoreline depends on the entire surf zone bathymetry and increases with decreasing surf zone b each slope. The sand bar depth can influence the wave setup, where shallow bar crests favor higher setups since there is a decrease in the radiation stress and the setup gradient increases. 2. 7 Sand Bar Formation Theories Two major theories of bar formatio n are often discussed in the literature. The first theory associates bar formation with infragravity waves and the second assumes wave breaking as the responsible mechanism for bar formation. The common point between these two theories is that incident wav es alter the flow field and consequently the bed morphology. However, there are several differences between those theories and it was observed that infragravity waves play a second ary role in the formation of sand bars. In th e infragravity wave theory, it is suggested that the sediments migrate from the nodes towards the anti nodes, since the orbital velocity below the nodes scour the sediments and they deposit below anti nodes (Short, 1975 ; Holman and Bowen, 1984 ). Secondary importance of infragravity wav es to cross shore sediment transport were noted by Aagaard et al. (2002) where they excluded their influence assuming only the effects of incident waves and mean currents to the sediment fluxes.
37 The second theory is that bar formation is a result of sedime nt convergence at the wave breakpoint. Seaward of the breakpoint, sediment moves onshore due to orbital asymmetries induced by the incident waves; shoreward of the breakpoint, sediment is transported offshore due to offshore directed mean currents (underto w) generated by breaking. Undertow velocit ies depend upon local cross shore radiation stress gradients and scale inversely with local water depth (Svendsen, 1984). Undertow velocities are expected to be larger over steep beaches than on gently sloping bea ches (Longuet Higgins, 1983). The wave breaking intensification on bars during storms causes a maximum in the undertow just sh oreward of the bar crest (Galla gher et al., 1998). Similarly, Thornton et al. (1996) a nalyzing data at Duck during the DELILAH exp eriment in 1990, observed that the undertow i s a maximum over the bar crest and a minimum in the trough. Aagaard et al. (2002) proposed a model that predicts onshore directed sediment transport for large bed shear stresses in relatively deep water occurr ing on gently sloping beaches. With increased breaking intensity in shallow water and for relatively steep nearshore slopes, undertows increase and the sediment transport becomes offshore directed. The lower waves during summer do not break over the bar an d the accumulating effect under the bars is replaced by an onshore bottom current that erode s the bar configuration and change s the winter profile to summer profile. Hanes and Huntley (1986) observed that in oscillatory flows the flow field acceleration may play an important role in suspending sediments where boundary layer
38 separation is likely to occur during the flow decelerating phase. In this situation, suspension appears to be stronger during the onshore directed phase of the wave motion than during the offshore motion and may be determined by fluid acceleration more than velocity. Field observations and numerical model results conducted by Elgar et al. (2001) suggested that the onshore sediment transport and sand bar migration are related to cross s hore gradients in skewed fluid accelerations associated with pitched forward non linear waves. Additionally, the region of strongly skewed accelerations moves shoreward with the bar, suggesting that feedback between waves and evolving morphology can result in continuing onshore bar migration when mean flows are weak. Hoefel and Elgar (2003) using a numerical model including gradients in acceleration skewness were able to predict the offshore and onshore sediment migration over a 45 day period at Duck. 2. 8 Wave Height Decay i nside the Surf Zone The waves approaching the shore shoal until they start to break in shallow waters according to the parameters described in the previous section. After the breaking point the waves transfer their energy generating tu rbulence into heat and the wave height decrease s in the nearshore zone. N umerical models were developed to describe th e wave height decay in plan e sloping and barred beaches and their results compared with laboratory and field data Thornton and Guza (198 3) developed a model similar in concept to Battjes and Janssen (1978) where the wave decay is calculated using the change in energy flux. The energy dissipation is considered primarily due to the conversion of organized potential and turbulent kinetic ene rgy to turbulent kinetic energy, using the bore dissipation function. This model calculates the wave height decay in monotonically
39 depth decreasing beaches, however the stabilization and following shoaling process are not described using this model. Dally, Dean and Dalrymple (1985) develop ed a similar wave height decay model however a stability crite rion was used to allow the wave breaking stabilization and the initiation of a second shoaling process in the inner surf zone. Th e inclusion of this term is es pecially important on barred and dissipative beaches where multiple breakpoints are commonly observed (more details about this model are presented on section 7.1) Figure 2 1. Schematic representation of wave breakpoint in dissipative and Reflective Be aches
40 CHAPTER 3 MOTIVATION M orphodynamics in the nearshore are extremely complex, sometimes the outer and inner bars present distinct behavior s making the predictability of these systems extremely difficult. Understanding the relationship be tween the wave energy and number of sand bars inside the nearshore zone at different time scales is important to better characterize this environment. It is expected that if the same wave energy acts on a beach for a long period of time, the beach achieve s an equilibrium profile and the number of bars stay s constant. The present research uses this assumption to develop a qualitative and a quantitative model to predict the average number of sand bars in a large scale (decadal) in the nearshore zone based on incident wave energy and beach characteristics. As reflected by the literature review most of the researches are concern ed with the nearshore morphologic response to individual storms or short periods of bar recovery, or in microscale solutions, therefor e the development of mesoscale models, require knowledge of more generic properties where specific conditions are simplified Current mesoscale models mostly focus on shoreline changes and neglect subaqueous changes in the beach morphology. The beach mo rphology, mainly in multiple bar systems, is a result of the beach long term response to the forces acting on it. Therefore the integrated effects of several storm recovery cycles have to be understood and long term series of forcing conditions have to be analyzed to explain the variability of the number of bars in the nearshore zone. Such knowledge is probably best obtained from the analysis of long time series of the relationship between
41 morphology and forcing conditions where new insights of the sand ba r behavior and quantitative estimates for coastal changes can be obtained On the other hand some of the short scale influence s of systems with multiple bars are essential The knowledge of the interaction between inner and outer bars in systems with mult iple bars is still poor ly understood and better insights between their relation s should be addressed to better understand this complex system. Using the assumption that sand bars can focus waves with different heights to break in a narrow zone, the short s cale relation between bars and breakpoints can be used to predict the number of bars and their interaction inside the nearshore zone.
42 CHAPTER 4 FIELD STUDIES 4.1 Beach Profile Data Beach profiles in Walton and Volusia Counties, on the east an d west coasts of Florida respectively, (Figure 4 1 ), documented by the Florida Department of Environmental Protection (FDEP) were used in the present research. Profiles are usually surveyed every 300 m in the alongshore direction based on physical monument s placed landward of anticipated beach fluctuations The cross shore distance is variable, however ideally the surveys extend to the depth of closure. Some of these monuments have been eroded or damaged through time, since most of them were placed approxim ately 40 years ago, and virtual monuments that are based only on the coordinate position have replaced them. Some of the Light detection and ranging (Lidar) data available on the http:/ /www.csc.noaa.gov/digitalcoast/data/coastallidar website were used mainly to observe alongshore differences in Walton and Volusia Counties, since this technique provide s high density data, with good precision (horizontal and vertical error in the order of deci metres) and reduction in survey time when compared with traditional methods. Few surveys using Lidar were conducted for this area and when Lidar data is available the FDEP extract s the profiles at the monuments from this source. 4.2 Definition of San d Bars C riteria to define sand bars were established for the profile analysis, thus allowing a large number of profiles to b e evaluated and quantitative comparisons to b e made during the analysis. The measured cross shore profile length s are variable, as m entioned previously, therefore only profiles deeper than six m were evaluated to
43 extract the sand bars and six m depth were excluded With this threshold depth it is expected that the sand bars in the profiles are correctly measured even in profiles with more than one bar, and potential errors offshore bars were minimized. This limit is greater than the depth of closure calculated as 5 and 5.5 m by Dean (2002) for Walton and Volusia Counties, respectively. A total of 1,10 5 and 5 40 profiles deeper than six m were analyzed in Walton and Volusia Counties respectively. More information and concerns associated with the depth of closure w ere presented by Nicholls et al. (1996). Beach profile zonation was considered the best approach for the sand bar height definition and three different zones were established according to their depths. Bar height was considered assuming the vertical dista nce from the bar crest to the landward bar trough. For depths smaller than 4.5 m the morphological features with heights greater than 0.15 m were considered as sand bars. Features with heights greater than 0.9 m were considered as bars in depths between 4. 5 and 7.5 m. For depths greater than 7.5 m, only heights greater than 1.5 m were considered. In these three zones the bar width had to be greater than 25 m. This division was based on the theory that bar formation occurs under breaking waves and features n ot capable of alter ing the hydrodynamic forces were not considered A theoretical profile showing the division used for defining sand bars characteristics is present ed on Figure 4 2 These limits are somewhat subjective a nd other depth limits could be prop osed, which as more appropriate for other application s or location s Kuriyama (2002) considered morphological features greater than 0.5 m, from bar crest to landward trough for the
44 entire profile, however no explanation for this threshold value was given. Usually researches on sand bar characteristics do not specify these criteria, with sand bars defined visually Even if uncertain the limits used here were consistent and were held constant during the analysis, allowing a quantitative way to extract the p rofile morphology. Another method to define sand bars is to calculate the equilibrium beach profile and extract bar properties with respect to the equilibrium beach profile. This method was discarded since the area evaluated in both Counties is large and the surveys were not conducted always on the same monuments, which would require apply ing equilibrium beach profiles that were not always obtained at the same location. Profiles surveyed at different monuments, but o n the same day, were averaged in order t o demonstrate general behavior for an area instead of just one representative profile. The issue with averaging beach profiles from different locations is where they proce dure was necessary, the mean high water (MHW) level was used to combine and average the profiles from different monuments, shifting the profiles in the cross shore direction and making the MHW position coincide on all them. One problem associated with thi s approach is that sometimes different features in a profile, for example bar and trough, from different profiles coincide in the same location and the averag ing damps the results masking these features. The average of 13 profiles from Walton County illust rate s dampening effects when bars and troughs are located in the same cross shore position (Figure 4 3). Even if this situation was
45 observed in the data, the average profiles were judged representative of the general trends occurring in the group and the m ain morphological features preserved. 4.3 Wave D ata Offshore wave heights were obtained from NOAA/NDBC website http://www.ndbc.noaa.gov using buoys 42039 and 41009, for Walton and Volusia Counties, respectively. Buoy 42039 located southeast of Pensacola in 307 m water depth has data available from December 1995 to April 2012. During this period a small quantity of data were missing (3.4%). The buoy used for Volusia County is 41009, located east of Cape Canaveral in 44 m water depth and data from this station is available from August 1988 to April 2012 (Figure 4 4 ). Waves with periods greater than 7.5 s can start to be influenced by the bottom in the buoy 41009 (considering the deep water limit as half of the wave length ). The first beach survey at Volusia County where wave data was available was in May 2001, and only data starting in 2001 were included in the analysis. During this period 0.6 % of the data are missing due wave gauge inoperability. Significant wave heigh ts from these two buoys were converted to wave energy using E quation s 4 1 and 4 2 Where is the water density g is the gravitational acceleration and H rm s and H s are the rout mean square and significant wave height respectively (4 1) (4 2)
46 4.4 Description of Walton and Volusia Counties An overview of all the available surveys in Walton and Volusia with emphasis on the number of sand bars in a profile, bar crest depth crest position in relation to the shore and the interaction between inner and outer bars is provided to characterize the nearshore zone in these two Counties These parame ters partially describe sand bar response s to different wave energy conditions and assist by introduc ing the areas that are analyzed in the following chapters. Table 4 1. Number of profiles with zero, one, two and three bars and the average number of bars in Walton and Volusia Counties in all surveys from 1972 to 2007. Walton County Volusia County Year 0 bar 1 bar 2 bars 3 bars Ave rage Year 0 bar 1 bar 2 bars 3 bars Ave rage 1973 0 40 3 0 1.1 1972 3 25 33 4 1.6 1975 15 23 2 0 0.7 1988 30 130 63 8 1.2 1981 1 69 15 0 1.2 2001 0 10 17 9 2 1995 0 114 12 0 1.1 2002 0 41 4 0 1.1 1996 5 78 42 1 1.3 02/2006 0 4 35 2 1.9 1997 5 110 12 0 1.1 09/2006 1 9 31 0 1.7 1998 0 32 19 3 1.5 2007 0 1 72 8 2.1 04/2004 4 112 7 0 1 11/2004 10 72 45 0 1.3 2 005 4 44 75 4 1.6 2007 10 81 35 1 1.2 Total(%) 4.9 70.1 24.2 0.8 Total(%) 6.3 40.7 47.3 5.7 The bar crest position and depth were obtained for all profiles with sand bars and the number of bars varied from profiles with no bar to three bars. The first surveys in 20 years, increasing their frequency after 1995. The majority of the profiles at Walton County presented only one bar (70.1 %) and profile s with three bars were observed on a few occasions In Volusia County, profiles with two bars (47.3 %) were more frequent
47 followed by profiles with one bar (40.7 %) and more profiles with three bars were observed when compared with Walton County. Table 4 1 shows the distribution of profiles from zero to three bars throughout the years. In order to investigate the sand bar linearity or if they have two or three dimensional characteristics, a linearity index (LI) was defined and applied ( E quation 4 3 ) V alues approaching unity indicate linear or tendency to two dimensional morphologies and the de crease in th is value is associated with an increase in the bar three dimensionality where alongshore differences in the morphology are more evident. (4 3) ( 4 4) In E quation 4 3 and 4 4 b is the bar crest depth; the same equation was used replacing b with the distance from the bar crest to the shore also referred to as bar position. Sand bars in profiles with two bars were divided according to their positio n and depth to define differences in the inner and outer bars stabilities since different wave conditions are affecting the inner and outer bar s and their response to the same offshore wave conditions can vary The third bar ( the offshore bar) in profiles with three bars was removed from the analysis since their occurrence was less frequent and 4.4.1 Walton County ion of which bar was more uniform and two dimensional. The LI results presented a large variability, however
48 after major storms, the outer bar at Walton County was more linear than the inner bar, since the LI were larger after these events. This was observ ed in the 1998 and November 2004 surveys, after hurricanes Georges and Ivan, respectively. One example where the outer bar became more linear after an event with high wave heights (Hurricane Ivan) is presented on Figures 4 5 and 4 6. Before the Hurricane o nly seven of the 127 profiles had two sand bars (diamonds on Figure 4 5) and the inner bar, when observed, was located in very shallow waters (blue diamonds on Figure 4 5) After Hurricane Ivan, most of the profiles presented two sand bars (diamonds on F igure 4 6), being the majority of profiles in this survey. It is evident that the inner bar, only observed in few profiles before Ivan, was created during this event. The outer bar changed to deeper water as a response to the high waves, however the positi on stayed almost constant with slightly bar offshore migration. The LI of the outer bar depth increased after Hurricane Ivan (from 0.981 to 0.99 6 ) and it was another indication that sand bar formation and their depths are directly related to the wave energ y acting prior to the survey. The LI for the outer bar position also increased after Ivan (from 0.968 to 0.977). However the inner bar became more three dimensional after Ivan (with a decrease in the LI) indicating that the beach morphology can still prese nt three dimensional features even after Hurricanes or major storms. Similarly, Gallagher et al. ( 1998 ) observed significant bathymetric inhomogeneities developed during a storm at Duck This fact was associated with the presence of megaripples with wave l engths of few meters that are present during both high and low wave energy periods. The results of Equation 4 3 and 4 4 applied to all surveys in Walton County are presented on Table 4 2.
49 Table 4 2. Results of the LI applied to all surveys in Walton Coun ty, evaluating the inner and outer bars linearity Walton County Survey Year Inner Depth Outer Depth Inner Distance Outer Distance 1973 0.9 88 0.9 5 5 0.9 90 0.9 75 1975 0. 962 0.9 63 0.9 61 0.99 2 1981 0.9 85 0.9 89 0.9 71 0.9 77 1995 0. 964 0. 961 0.9 94 0.9 76 19 96 0.9 72 0.998 0.9 73 0.990 1997 0. 935 0.9 47 0.9 94 0.9 43 1998 0. 997 0.99 0 0.9 84 0.9 77 04/2004 0.9 99 0.9 81 0.9 99 0.9 68 11/2004 0.9 94 0.99 6 0.9 77 0.9 77 2005 0. 956 0.9 88 0.9 77 0.9 89 2007 0.9 97 0.9 84 0.9 91 0.9 7 9 4.4.2 Volusia County The inner and outer bar depths presented opposite behaviors, an increase in two dimensional morphologies in the outer bar occurred at the same time that an increase in three dimensional morphologies in the inner bar, observing two consecutive surveys. Also, a decrease in the inner bar LI (more three dimensional) was observed when the outer bar became more two dimensional (higher LI). The LI for the bar distance presented a similar behavior, the exceptions occurred in the last two surveys that the LI increase for the inner and outer bars at the same time. Shand et al. (1999) found that sand bars are higher on more exposed beaches with mild beach slopes. However this was not observed when comparing the profiles with the milder slope in Volusia and steeper slopes in Walton County The average bar height was greater in Walton County (1.01 m) than in Volusia County (0.71 m). All the
50 bar heights and widths in both Counties are presented on Figure 4 7, where it is clear that the sand bars are smaller and narrower in Volusia. Table 4 3. Results of the LI applied to all surveys in Volusia County, evaluating the inner and outer bars linearity Volusia County Survey Year Inner Depth Outer Depth Inner Distance Outer Distance 1972 0.9 91 0.9 63 0.9 86 0.9 78 1988 0.9 83 0.9 76 0.9 93 0.9 47 2 001 0.9 06 0.998 0.990 0.995 2002 0.9 99 0.9 78 0.9 89 0.9 68 02/2006 0.9 69 0.99 5 0.9 81 0.9 89 09/2006 0.9 90 0.99 2 0.9 86 0.9 91 2007 0.9 96 0.99 7 0.9 78 0.98 2
51 Figure 4 1. Location of the two study areas in the State of Florida Walton County in the we st coast and Volusia County in the east coast.
52 Figure 4 2. Criteria utilized to define the sand bars inside the three zones in the beach profile Only bars with heights greater than 0.15 m were considered from 0 to 4.5 m; heights greater than 0.9 m from 4.5 to 7.5 m and heights greater than 1.5 m on depths greater than 7.5 m.
53 Figure 4 3. The 13 profiles of Group 1 at Walton County during the 1997 survey (black) and the average profile (red).
54 Figure 4 4 Loca tion of the buoys 42039 and 41009 used to extract the wave data for Walton and Volusia Counties, respectively.
55 Figure 4 5 Beach survey before Hurricane Ivan. Profiles with one sand bar (blue crosses) and profiles with two sand bars (diamonds). The profiles with two bars have the inner bars represented with blue diamond and the outer bar with red diamond. A) Bar crest depth. B) Bar crest distance to the shore.
56 Figure 4 6 Beach survey after Hurricane Ivan. Pr ofiles with one sand bar (blue crosses) and profiles with two sand bars (diamonds). The profiles with two bars have the inner bars represented with blue diamond and the outer bar with red diamond. A) Bar crest depth. B) Bar crest distance to the shore. D ep th and distance limits distinguishing the inner and outer bars to calculate the Linearity Index (dashed horizontal line).
57 Figure 4 7 Sand bar depth s and wid t hs from all profiles in Walton (red circle) and Volusia (blue cross) Counties. Calculated linear trend for both Counties (straight lines).
58 CHAPTER 5 QUALITATIVE MODEL: W AVE HEIGHT AND NUMBE R OF SAND BARS RELATIONSHIP A qualitative approach was used to determine the correlation between the average number and behavior o f sand bars influenced by different wave conditions This approach is based on the theory of winter and summer conditions developed by Bagnold (1947), on which breaking waves are the main force for bar formation and the sediments move offshore during energ etic conditions and migrate onshore during longer periods with mild energy. Bar formation and consequently an increase in the average number of bars is expected with the increase in wave energy. The opposite behavior occurs during calm periods, when the en ergy decreases and the average number of bars decreases. The effects of different hurricanes that impacte d Florida and increased the wave energy in Walton and Volusia Counties were qualitatively observed in this approach. 5.1 Wave Analysis The average wave energy between two consecutive surveys was used as the characteristic wave energy during that period for the qualitative model, showing the wave energy influencing the beach during those periods. These intervals are not constant and their durations change d according to when beach surveys were obtained in each county. The comparison between the wave energy in consecutive periods indicates if an increase or decrease in the average number of bars would be expected. Wave energy periods for Walton and Volusia Counties are differentiated with W and V followed by numbers according to the beach survey. The wave energies for the first and last time spans were calculated, however it is not possible to interpret them since no beach profiles were conducted before and after them and t he sand bar
59 morphodynamic responses to the wave conditions are unknown Also the first period in Walton County, prior to W1, has less than three months of data, and its interpretation as a period with high average wave energy can be biased by the lack of additional data. 5.2 Average Number of Bars and Average Profile The number of sand bars extracted from the profiles was average d resulting in one number of bars for each survey ( T able 4 1) The average was chosen to minimize potential errors associated with beach three dimensional morphologies, since the profiles ar e obtained along a single cross shore line based on the monuments at the beach and not always the same monuments were surveyed for all years This potential source of error is disc ussed later. All the profiles from Wa lton (1,105 ) and Volusia (5 40 ) Counties were averaged to determine the general nearshore characteristics in both c ounties especially emphasizing the beach slope (Figure 5 1 ). The average beach profile at Walton is ste eper (1/65) than at Volusia (1/90) County. The beach slope differences are consistent with the coarser mean sediment size of 0.3 mm in Walton County compared to 0.18 mm in Volusia (Charles et al., 1994) It is also consistent with the greater average wave energy at Volusia, which tends to flatten the beach profile, transport sediments offshore and create sand bars. This large number of profiles includes surveys not used in the model since they were conducted prior to availability of wave data therefore com parisons between hydrodynamics and morphologies were not possible. The greater number of bars in the Volusia profiles observed in the data was expected as the beach slope i s milder in compar ison with Walton. The waves are expected to break farther offsho re and as the distance to the shore is still large they can
60 stabilize, shoal inside the nearshore zone and break for a second time, as proposed by Dolan and Dean (1985) for systems with multiple bars This process can repeat several times and profiles with more than one bar are the consequence of multiple breakpoints. Beach profiles with one, two and three sand bars characterizing the situation described above are presented on Figure 5 2 all from Volusia County. 5.3 Bar Response to Average Wave Energy 5.3. 1 Walton County The average number or bars in Walton County correlated with the wave forcing (R 2 = 0.41 ) Periods with low (W1, W3 and W6) and high wave energy (W4 and W5) were found allowing a qualitative evaluation between wave forcing and sand bar resp onse. When the average wave energy decreased between two consecutive periods, the average number of bars also decreased. On the other hand, the periods with an increase in the wave energy were associated with an increase in the average number of bars. Wav e energy increase was mainly associated with the passage of six Hurricanes (Earl, Georges, Ivan, Dennis, Katrina and Ike) in Walton County (Figure 5 3 ). The effects of the first two Hurricanes (Earl September 3 rd and Georges September 29 th 1998) observed during this period were captured by the beach survey in November 1998. The sand bars responded to the increase in the wave energy from 1,723 to 2,335 J/m 2 with an increase in the average number of bars from 1.1 to 1.5. Hurricanes Ivan, Dennis and Katrina happened over a short period of time between September 2004 and August 2005 (time span W4 and W5) and wave energies greater than 1 .3*10 5 J/m 2 were observed during the first two hurricanes and resulted in significant beach profile changes. This large increa se in average wave energy, from 1,909 to 2,313 J/m 2 and
61 later to 2,811 J/m 2 (45 % more than during calm periods) increased the average number of bars from 1.0 (W3) to 1. 3 (W4) remaining high (1. 6 ) in the W5 period. The Hurricane ermined, since n o survey was conducted subsequently, however an other increase in the average number of bars would be e xpected after this hurricane. A reduction in the average number of sand bars occurred after calm periods W1, W3 and W6. This decrease is a ssociated with onshore sediment transport, where the bars migrate back to the beach. The most significant was during W3, when the average number of bars decreased from 1.5 to 1. 0 as a result of more than five years without significant increases in wave ene rgy. 5.3.2 Volusia County Only five surveys were conducted in Volusia County, result ing in a smaller quantity of data to be compared with the model (Figure 5 4 ). However the same methodology that was applied for Walton was used for Volusia County. Hurric anes Frances, Jeanne and Wilma occurred in time span V2, however a small decrease in the average wave energy occurred from V1 (2,630 J/m 2 ) to V2 (2,470 J/m 2 ). The small wave energy decrease, even with three hurricanes, can be explained by the long period, almost four years, where the data were averaged, decreasing the wave energy peak effects during time span, where a decrease in the average number of bars was expect ed in the 2006 survey, rather a significant increase from 1.1 to 1.9. The coefficient of determination (R 2 ) between the average wave energy and average number of bars was extremely low (0.001) at Volusia. This situation highlighted the necessity of frequent bea ch surveys in order to apply the qualitative approach with better accuracy. W here beach surveys are not available more frequently data can be averaged over long periods where wave
62 conditions can change considerably and the averages become non representati ve of the entire situation. Also as the hurricanes, especially hurricane Wilma, occurred close the wave forcing that happened during the closest period. The morphologic indicates that the beach morphology is not simply a result of the hydrodynamic conditions acting in the present moment but also during past events. Kriebel and Dean (1993) applied a convolution method for time dependent beach cross shore respon se to comparison between predicted and measured time scale for beach response. This also indicates the importance of when beach surveys are conducted, since it is expe cted that the beach morphology represent s the conditions occurr ing dur in g periods closer to the survey, with this period being associated with the energy at the beach. If energetic conditions happen, this period is shorter, otherwise it is longer. The res ponse of the average number of sand bars during the next two consecutive time spans V3 and V4 occurred as expected, decreasing from 1.9 to 1.7 at the calm time span V3 and increasing from 1.7 to 2.1 during the more energetic time span V4. 5.4 Summary of Fi ndings from the Qualitative Examination The qualitative approach, based on a simple methodology, showed reasonably good results when comparing the average wave energy and the average number of sand bars. It also supported a long term correlation between th e wave energy and the number of bars observed on sandy beach profiles with different grain sizes and exposed to different wave energy conditions. Hurricane influences were well represented at Walton ; however the same results were not achieved at Volusia C ounty. The time when the survey s w ere conducted and the period between two consecutive surveys were
63 found to be key to the achievement of good results with this methodology. The positive feedback between waves and sand bars obtained using the qualitative a pproach motivated the development of a simple numerical model to quantitatively predict the number of bars in the nearshore zone. Similar numerical model s have not been presented before in scientific reports, and details about the model ar e described i n ch apter six Figure 5 1. Average beach profiles at Walton and Volusia Count ies representing the differences in beach profile slope.
64 Figure 5 2. Measured profile e xamples in Volusia County A) P rofile with one bar. B) Profile with two bars. C) Prof ile with three bars.
65 Figure 5 3. Results from the qualitative model for Walton County. A) Comparison between average wave energy and average number of bars. B) Wave energy and indication of H urricanes The vertical dashed lines shows when be ach surveys were conducted and the highlighted periods indicate when the wave energy decreased between consecutive periods
66 Figure 5 4. Results from the qualitative model for Volusia County. A) Comparison between average wave energy and averag e number of bars. B) Wave energy and indication of H urricanes The vertical dashed lines shows when beach surveys were conducted and the highlighted periods indicate when the wave energy decreased between consecutive periods
67 CHAPTER 6 NUMERICAL MODEL TO PREDICT THE EVOLUTION OF SAND BA RS IN THE NEARSHORE ZONE Most of the available qualitative and quantitative models describing the nearshore characteristics are developed for beaches with only one sand bar. The increase in the number of sand bars, resulting in nearshore zones with two or multiple sand bars, increases the complexity in this zone since morphodynamic interactions occur inside the nearshore zone The basis forming th e model developed here, simplifies these interactions, establishing a relation between the wave forcing and the quantity of sand bars in the nearshore zone, where this relation is dependent on beach and wave characteristics. 6. 1 Model Description The model to predict sand bar evolution was evaluated using the same profile database from the FDEP described previously, in Walton and Volusia Counties. The model considers that the number of bars in equilibrium, N eq ( i.e., the number of bars expected in the nearshore zone) is a function of offshore wave and beach characteristics To decide which variables are important, two criteria suggested by Wijnberg (2002), were used. Firstly, the location of change in the variable value should match the location of change in the behavior and, secondly, the change in variable value should be equally abrupt as the change in the behavior. This implies that it is implicitly assumed that a gradual change in the extrinsic conditions results in a gradual change in morphologic behavior. Consequently, although the relation between input conditions an d morphologic response is probably non linear, it is assumed that the non linearity does not cause abrupt changes in the morphologic response.
68 The model assumes alongshore uniformity, although it is known that this is not representative during calm condit ions, it is a reasonable approximation when the model objective is the prediction of the number of sand bars over a large temporal scale of decades where influences from smaller perturbations can be omitted The physical variables in this situation are a function of the cross shore distance and time. Calculation for long time periods encompasses greater variations in the forcing as observed i n Figures 5 3 and 5 4, and is impractical for predictions using transport models based on local quantities. In cont rast, transport models for macro scale require formulations based only on simple descriptors of the forcing, but must include information on the principles underlying key features of the larger scale morphology (Larson and Kraus, 199 2 ). In their paper the authors propose the use of quasi equilibrium models that describe the transformation between different beach states representing key morphologic features where the rate of change is determined by simple or bulk parameters that characterize the forcing in m oving between states. The parameters frequently identified during the literature review that control the sand bar behavior inside the nearshore zone are the waves and the morphologic characteristics. Some parameters including one or more agents influencing the sand bars were applied in the data analysis. These characteristics are: the offshore wave height (H 0 ), wave period (T) and sediment size which is represented in terms of settling velocity and beach slope. N eq was tested using four different functions which are parameters frequently used in nearshore process studies First, the function was based only on the offshore wave height. Wave steepness (H 0 /L 0 ) was also tested to include the wave period influence, where L 0 is the offshore wave length.
69 Two non dimensional parameters were tested. The first being the Dean number (D), based on the settling time of suspen ded sediments by breaking waves proposed by Dean (1973). Assuming that for short settling times the sediment is influenced by onshore orbital velo cities and deposited onshore and the sediment is influenced by offshore velocities and deposited offshore during long settling times. These two processes create the normal and storm profiles, respectively. N eq as a function of the Profile (P) number (gH 0 2 / 3 T) proposed by Dalrymple (199 2 ) was the last function tested. This is also a non dimensional number and it was developed based on results reported by Larson and Kraus (1992) obtained in a large wave tank, resulting in a threshold value (~10,400) able to distinguish normal and storm profiles. In contrast with the Dean number which was based on small scale laboratory data, the Profile number was obtained using large wave tank data. In summary, N eq was tested as a function of H 0 H 0 /L 0 D and P, using E quati on 6 1 (6 1) w here The coefficient k 1 is included to represent site specific characteristics. It is dependent on the beach characteristics, as surf zone width, sediment size and beach slope, where those three parameters are known to be inter rela ted. The rate of change of the number of bars is based on the difference between the predicted and equilibrium number of bars as follows (Equation 6 2 ) (6 2) The model includes that bar formation occurs more rapidly than the decay. The coefficient ( k 2 ) during periods of bar decay (N b > N eq ) is smaller than the coefficient ( k 3 )
70 during periods of bar formation (N b < N eq ). As discussed previously, this account s for observations that during energetic waves, bars form faster and during calm periods the bars migrate slowly onshore and the number of bars decrease s Calibration methods were applied to determine the minimum error for the three coefficients used in the model ( k 1 k 2 and k 3 ). The same values of k 2 and k 3 were used in Walton and Volusia Counties, o nly K 1 was changed, since the latter defines the expected number of bars which depends on morphological characteristics and the first two are related to sediment transport time response s Figure 6 1 shows a conceptual example to demonstrate the model fun ctionality. It is idealized showing an increase in the wave height (from one to two meters) that remains constant and after a period the wave height decrease s back to the original value (Figure 6 1 A ) A profile without sand bars is used as an initial condi tion. During the first stage (P1) the number of bars is not in equilibrium with the incident wave conditions and the number of bars increase s until it stays in equilibrium with the expected number of bars associated with the current forcing (Figure 6 1 B ) When the wave height increases, there is a second period that the model is not in equilibrium (P2) A new equilibrium is achieved with a second increase in the number of bars. T he wave height decrease s and the model achieve s the third equilibrium with a d ecrease in the number of bars (P3) The rate of increase and decrease in the number of bars is not the same. Faster rates are observed during periods of bar formation and slower rates during the periods of bar decay (P1 and P2 < P3) similarly to the natu r al behavior of sand bars
71 6 2 C omparison with Florida Data The numerical model results were compared with data obtained from Walton and Volusia Counties. These two Counties were selected because they are influenced by different wave conditions and also dis tinct morphologic features with coarser sediments at Walton The average beach profile at Walton is steeper than at Volusia County (Figure 5 1 ) indicating that forcing and morphology in a mesoscale observation are distinct. To quantitatively represent tho se differences in the numerical model, an empirical correlation, based on the dataset described on Section 4 .1 between the expected number of bars and wave energy was used (Figures 6 2 A and 6 2 B ). These relations were obtained using different values of k 1 for Walton and Volusia and using N eq as a function of Dean and Profile numbers. The number of bars in equilibrium is larger in Volusia County due to the wider surf zone, finer sediment and higher waves. The correlations using the other two functions, H 0 and H 0 /L 0 are not shown because the correlation with the model results was lower which are described later (s ections 6 4 and 6 5 ) The division between a profile with more or less than one bar using the Profile parameter was P = 3.09 *10 5 and P = 1.89 *1 0 5 for Walton and Volusia Count ies respectively. These values are considerably higher than the threshold value of 10,400 for the division between summer and winter profiles proposed by Dalrymple (199 2 ), however that value was obtained using large wave tan k data versus the values obtained here based on field data at Florida. The division between storm and summer profiles using D are between four and five, and D = 6.4 and 3.8, for Walton and Volusia Counties were obtained using the empirical relationship.
72 6 3 Walton and Volusia County Groups Beach surveys were not conducted on the same day in the entire C ounty because the shoreline in these two Counties, e specially Volusia, is extensive and s ometimes the surveys were stopped due to high wave conditions. For t hese reasons, beach segments were surveyed in different days and waves with considerable different heights occurred during one survey period. To avoid data bias due to different wave forcing during the same survey, five areas where separated and analyzed i ndivid ually in Walton County (Figure 6 3 ). The five areas had all their profiles surveyed on the same day. Another concern was avoid ing alongshore differences in the profiles, mainly during calm periods whe n beach features become more three dimensional and alongshore differences in the sand bar morphology are likely to occur Only shoreline lengths greater than three km were used as groups and it is believe d that the average number of bars over the three km shoreline minimizes the problems associated with a longshore variation s in the profile w ithin a group making the averages representative of the group for the time scale used in the model One continuous beach segment, between R 165 and R 195, with approximately nine km was surveyed four times at Volusia, instead of five surveys in the qualitative This segment is located 4.8 km south of Ponce de Leon Inlet, where the shoreline change fluctuations decrease to the ambient County fluctuations (Absalonsen and Dean, 2011), diminishing the inlet influence.
73 6 4 Walton County 6 4 .1 Comparison b etween the Groups The five groups in Walton County showed a similar behavior of the measured average number of bars (Figure 6 4 ). Group 1 during the 2005 survey, was the only group surveyed before Hurricane Katrina and the only one with a decrease in the average number of bars between the second survey in 2004 and the 2005 survey, since ber of bars associated with Hurricane Katrina. All the other groups were surveyed after Hurricane Katrina and presented an increase in the average number of sand bars during that survey. The average number of bars varied considerably between the groups in the 1996 survey used as an initial condition This difference was the result of an increase in wave height ( to 2.5 m) that occurred between February 23 rd and March 4 th The dates when the surveys were conducted are presented on T able 6 1. Table 6 1. Dat es that surveys were conducted at Walton County Survey Dates on Different Groups Walton County S1 S2 S3 S4 S5 S6 S7 Group 1 3/5/96 11/19/97 11/18/98 4/15/04 8/15/04 8/24/05 7/21/07 Group 2 3/4/96 11/21/97 11/18/98 4/15/04 8/16/04 9/12/05 7 /21/07 Group 3 2/23/96 12/7/97 11/17/98 4/15/04 8/17/04 9/17/05 7/21/07 Group 4 2/22/96 12/11/97 11/17/98 4/15/04 8/18/04 9/19/05 7/21/07 Group 5 2/22/96 1/24/98 11/17/98 4/15/04 8/19/04 9/20/05 7/21/07 C hanges in decadal coastal behavior on the Dutch Coast showed that l arge man made structures appear to have effects that far exceed the spatial extent of known local effects of these hard structures ( Wijnberg, 2002) However this is not the reason of the changes between groups in Walton County, where no large hard structures are present along the coastline.
74 6 4 .2 Sand Bars v ersus Wave Height The results for group three using different functions for N eq are presented i n Figure 6 5 T he entire results from the model are not shown i n Figure 6 5 but only t he comparison between model results for the same time that beach surveys were conducted. Comparing the results using the four different functions, the best results were obtained using N eq as a function of the Profile number (R 2 = 0.75 and Var = 0.039). The first survey available was not used to calculate R 2 and variance, since it was the initial condition and therefore the model values were set equal to measured values. R 2 and variance results using different functions are summarized i n T able 6 2. The shape of all curves, using N eq as a function of the four different variables is similar, increasing and decreasing at the same moment when two consecutive survey dates are compared. The only exception is the 2005 survey, where N eq as a function of P decreased w hile the model results using the other functions predicted an increase. Even in this case, the prediction was the closest to the measured value. The results in Figure 6 6 compare the time evolution of the number of bars predicted by the model versus the av erage number of bars observed in Walton County. Measured offshore wave heights are also included (Figure 6 6 A ). Only the model results using N eq as a function of the Profile number are showed, since this parameter presented the best agreement with the obse rvations. Model results are not sensitive to the initial conditions applied. The results for the five groups differ only during the first seven months, after which high waves occur and the model results are the same, demonstrating the insensitivity to ini tial conditions of the number of bars and greater dependence on wave characteristics. The model
75 evolution during the first seven months is expanded i n Figure 6 7 where it become s clearer that the model results are more dependent on wave conditions and th at high waves tend to cause the model results to converge to the same value. Plant et al. (2004) observed in their hindcast prediction of the beach profile evolution at Duck, that the large number of bathymetric calibration data and the relatively small nu mber of free model parameters led to efficient optimization, which did not depend strongly on the initial parameter values used, albeit their model was used in a different time scale, maximum of 17 days of profile evolution hindcast. Wijnberg and Kroon (20 02) reported that the relaxation time of the nearshore system is small compared to the forcing event duration. In this case, each forcing sequence is uniquely related to a particular morphologic response sequence; i.e. independent of the initial morphologi cal condition, a given forcing sequence produces one and the same morphological sequence. In reality, the duration of events is often not long enough to let the bar system reach an equilibrium with the conditions (e. g., Bauer and Greenwood, 1990; Holman a nd Sallenger, 1993). When the relaxation time is large compared to the duration of the event, the chronology in the forcing events becomes important for understanding the evolution of the bar system. Table 6 2. Statistics and coefficient values for Walto n Count y Walton County Group 3 H 0 H 0 /L 0 Dean Profile Variance 0.092 0.094 0.062 0.039 R 2 0.581 0.656 0.548 0.747 k 1 0.680 0.597 0.395 0.002 k 2 0.004 0.003 0.008 0.002 k 3 0.320 1.880 2.400 1.900
7 6 The simple macroscale model represented the rapid increase in the expected number of bars after events with high energy. During the three energetic events associated with hurricanes Earl and Georges (1999), Ivan (2004) and Dennis and Katrina (2005) the model predicted a maximum average of 2.1 bars after h urricanes Dennis and Katrina. It also shows the slow decrease in the average number of bars after periods with high energy. This process is evident comparing the mild curve slope during periods with smaller energy and the abrupt increase during periods wit h high represented in the qualitative approach was preserved in the results using the numerical model. 6 5 Volusia County The four functions for N eq were compared using the data a t Volusia. Larger coefficients of determination were found using all functions ( T able 6 3 ). The small quantity of beach surveys is one factor contributing to these larger values. N eq as a function of Dean Number showed the best correlation and smallest var iance. Figure 6 8 presents the comparison between survey data and the four different functions. N eq as a function of the wave steepness was the only situation with a different behavior from the others decreasing the predicted number of bars from the Septe mber 2006 to the 200 7 survey The previous problem in the qualitative approach, with the misinterpretation of the long period when H urricanes Frances and Jeanne (2004) and later on Wilma (2005) altered the wave conditions, was represented well using the n umerical model with the large number of bar variations during this period. The m odelled value increased, with a peak of 2.4 during Frances and Jeanne. The changes in the beach morphology, moving
77 the sediment offshore and increasing the number of bars were predicted by the model and confirmed in the next survey (Figure 6 9) Table 6 3. Statistics and coefficient values for Volusia County Volusia County H 0 H 0 /L 0 Dean Profile Variance 0.021 0.056 0.010 0.013 R 2 0.860 0.840 0.990 0.620 k 1 1.210 9.410 0.510 0.002 k 2 0.004 0.003 0.008 0.002 k 3 0.320 1.880 2.400 1.900 The values of k 2 and k 3 were fixed for both Counties, however they changed according to the function used (Tables 6 2 and 6 3). At Volusia, the coefficients were 0.008 day 1 and 2.4 day 1 respectively, using N eq as a function of D. If the difference between the N b and N eq is one, a period of bar decay, it takes 125 days to the bar decay; on the other hand if this difference is 1, a period of bar formation, it takes less than half a day (0.41 days) to create one bar. Allowing the use of variable coefficients ( k 2 and k 3 ) for each County increased the correlation between modelled and measured data (these results are not presented). However, it is judge d that th ese coefficients are related to the bar mobility and keeping the same value for both c ounties is more representative of conditions happening in nature. 6 6 Hurricane Season Effects Seasonal effects can be determined by the model representing the changes associated with different wave heights during the year. An increase in the predicted
78 number of bars during the hurricane season and a decrease during the following months resulted in the well know storm and normal profiles. The results from the model appear to have a seasonal signal (Fi gures 7 5 and 7 8 ). To investigate this effect in more detail the model results were divided in one year periods and those periods were superimposed. A total of 13 periods were created in Walton and seven periods in Volusia County (Figure 6 10 ). The hurri cane season in the Atlantic Coast extends officially from June 1 st to November 30 th however 96 % of the major hurricane days occurs from August to October (Landsea, 1993). T he hurricane influences are evident evaluating these results, where a major increa se in the number of bars modelled occurs during this period. In Walton County t he peak number of bars predicted is 2.1 which occurred after hurricane Ike in 2004. After September there is a tendency of a slow decrease in the average number of bars. The se asonal average characteristics for the 13 periods is also plotted on Figure 6 10 and even with dampened values, expected when averages are applied, the same behavior is observed with an increase from 0.9 to 1.2 bars during the hurricane season. At Volusia, the s mall values in May 2002, distinct for the rest of the data, are related to the initial condition applied by the model, using the average number of bars in the first survey. The average number of bars in this survey was considerably low er in compariso n with other surveys in the same region in the County. Seasonal effects were even clearer in Volusia County. A m inimum average value of 1.55 bars was calculated in the beginning of August with an increase to 1.95 bars one month later, with a maximum ave rage number of bars predicted occurring in the
79 beginning of November (1.98). The value stays high until February, t hen starts to decay during a long period of six months, until the beginning of the hurricane season. Two peaks in the model results occurred close to each other, associated with Hurricane Frances (2.36) and Hurricane Jeanne (2.4). Figure 6 1. Theoretical example demonstrating the model assumptions. A ) Wave height B ) T he model results ( red dashed line ) and N eq (solid blue line) The period s where the model i s adjusting with the N eq are highlighted in gray P1 and P2 are periods of bar formation and P3 a period of bar decay
80 Figure 6 2 Expected number of sand bars A) U sing N eq as a function of Dean Number. B) U sing N eq as a function of Profile Parameter, for Volusia (blue) and Walton (Red) Counties The c oefficients used in the equations are also presented.
81 Figure 6 3 Division of the groups with same dates in the seven surveys conducted in Walton County.
82 Figure 6 4 Co mparison between average number of bars from measured data in the five groups at Walton County.
83 Figure 6 5 Comparison between different N eq functions (dashed lines) used to calculate the number of bars in the nearshore zone for Walton County and measured data (red solid line) The first value is the same for all functions, since it was used as an initial condition.
84 Figure 6 6 Quantitative model results for Walton County. A) Wave height measurements from buoy 42039. B) Predicted numb er of bars using N eq as function of Profile Number and the average number of bars measured in the five groups at Walton County
85 Figure 6 7 Comparison between the number of bars during the initial nine months of simulation for the five groups in Walt on County.
86 Figure 6 8 Comparison between different N eq functions (dashed lines) used to calculate the number of bars in the nearshore zone for Volusia County and measured data (red solid line) The first value is the same for all functions, since it was used as an initial condition.
87 Figure 6 9 Quantitative model results for Volusia County. A) Wave height measurements from buoy 41009. B) Predicted number of bars using N eq as function of Dean Number and the average number of bars meas ured in Volusia County.
88 Figure 6 10 Seasonal results obtained from the quantitative model A) Walton County. B) Volusia Count y One year periods (black line) and average season behavior (red line). The period when 96 % of the major hurricanes occur i n Florida (gray highlight)
89 CHAPTER 7 SHORT TERM PREDICTIO N OF THE NUMBER OF S AND BARS It was observed that the sand bars at Walton and Volusia Counties are relatively linear and the ir depths are located in a relatively narrow area of the nearshore zone especially after major storms Therefore another method that correlates the number of sand bars with the number and location of the break ing point s in the nearshore zone was applied for the same data set described in the previous chapter This method is focused o n the short time scale, or on the sand bar response to wave energy conditions that created the bars A different focus from the method applied in the previous chapter that described the evolution of the number of bars over a longer time scale (yea rs or decade) is examined here The Dally, Dean and Dalrymple (1985) wave breaking model was used to calculate the wave height decay in the nearshore zone. This model was applied since it was used in previous research ( Dally, 1992; Rakha and Kamphuis, 1995 ) with good agreement with large wave tank and field data. The model also requires few input conditions, which are the offshore wave conditions and the beach morphology An advantage of this model when compared with other models that calculate the wave dec ay in the nearshore zone is that it is able to predict multiple breakpoints and these breakpoints can b e associated with sand bar position s Battjes and Janssen (1978) proposed another method to calculate the wave height decay in the nearshore zone, howeve r th at method is not able to predict wave reformation and multiple break points, and it is not possible to correlate the number of breakpoints with the number of sand bars in systems with multiple bars.
90 7 .1 Dally, Dean and Dalrymple Model and Modifications The Dally, Dean and Dalrymple (1985) model is a two dimensional model that calculates changes in wave height in the nearshore zone due to wave shoaling and breaking, induced wave setup and bottom stress. The model considers that waves propagating from sl oping to horizontal bottom will not instantaneously stop breaking just because the bottom become horizontal (as dictated by the 0.78 criterion), but breaking continues until a stable wave height is attained. The energy dissipation rate per unit plan area i s assumed to be proportional to the difference between the local energy flux and the stable energy flux, according to E quation 7 1. Analytical solutions can be obtained using this method for planar and regular sloping beaches and a numerical solution of E q uation 7 1 can be used for more complex profiles (7 1) Where EC g and E s C gs are the local and stable energy fluxes, K is a dimensionless decay coefficient and h is the water depth. The two conditions on Equation 7 1 apply after the first break and before the first break Equation 7 1 is always zero. Th e model utilizes two threshold values to predict the initial breaking position and the position where the waves stabilize and start to shoal again (stability criteri on ). Both criteria depend on the water dept h. In contrast with the original model, these two threshold values were not constant in the applications herein The breaking criteri on the threshold value that determines where the waves start to break, is associated with the beach slope, therefore a lar ger breaking criterion was used in the simulation at Walton, since the profiles at Walton are steeper than at Volusia.
91 The stability criterion varied also according to the beach slope. The stability criterion for profiles with steep slope s are smaller, since the waves suddenly break in an narrow zone with plunging breaker characteristics and dissipate larger amount s of energy in a specific zone ( Ting and Kirby, 1994 ) Wave stabiliz ation and a second shoal ing and breaking process happen s less frequently i n beaches with steep slope. Figure 7 1 shows the relation of the beach profile with the breaking and stability criteria where the importance of the bar troughs in stabilizing and commencing the second shoaling process is shown 7 .2 Role of Sand Bars in the Breaking Process S imulation s w ere conducted to investigate the role of the sand bar as a trigger to the breaking waves and to evaluate focusing effects caused by the bar crests on the breaking waves where waves with different heights break in a re latively narrow zone on the bar crest This phenomenon was observed by studies using video cameras to monitor the beach (Lippmann and Holman, 1990). An example of this simulation is showed in Figure 7 2 for a profile with three bars (this was the maximum n umber of bars observed in a single profile in the database ) This simulation consisted of three separate steps where the sand bars, from the outer to the inner bar, were subsequently included in the simulations after each step and the breakpoint position c ompared with the bar crest position Actual beach profiles with three bars from Volusia County were used in these simulations. To initiate the first step, the offshore wave height was determined to shoal in the nearshore zone and to break over the offshore bar crest. An offshore wave height of 1.5 m was used in this simulation and the wave shoal ed until the offshore bar crest (835 m from the MHW ), where it achieve d the breaking criteri on with a wave height of 2.1 m and commenced
92 break ing This wave dissipat es it s energy until the stability criterion was reached i n the bar trough (at 863 m). After this point the wave re shoal ed and br o k e again with 1.5 m height at 986 m and the stability criterion attained and the wave d id not shoal a third time. Howev er, the location of the second break in the model is 35 m onshore of the second bar crest (Figure 7 2A) The second step was initiated including the second and third bars (the first or inner bar was not present in this simulation) and it is show n i n Figure 7 2 B The wave height follows the same distribution until the first breaking point, however after that, the second bar changes the beach profile and forces the wave to break offshore (934 m) when compared with the first step (986 m), reinforcing the conce ntration of breaking waves on the top of sand bars. The bathymetr ic change, including the second bar, triggered the third breaking processes that occurred at the original third bar location (1 024 m) T he third breaking is only possible because the beach s lope is gentle (1:48) and there is enough distance for the wave to shoal if this profile was steeper the wave would reach the shoreline before shoal ing to a breaking height The last simulation, the third step, used the original profile, with three ba rs, and the results were very similar to the simulation without the first bar The only difference is that the third break was 10 m farther offshore than with the second step T he fourth file has only three bars (Figure 7 2C) The same procedure, a simulation with three steps, was applied on other profiles with three bars at Volusia and similar results were obtained d emonstrat ing the sand
93 bars potential to focus the wave energy during t he simulations and the good correlation between the number of breaking points and number of sand bars 7 3 Prediction of the Number of Sand Bars The simulation presented i n section 7 .2 confirmed that the sand bars are able to focus the waves breaking on th eir crest and predict correct ly the number of breaking points on the profile can be used as an indicati on of the number of bars in the nearshore zone. This is the assumption used by the video camera systems, like Argus and CamEra (developed in USA and New Zealand, respectively), to monitor the nearshore zone and sand bar behavior (Lippmann and Holman, 1990 ; Aarninkhof et al., 2003) The second application of the Dally model was to predict the number of sand bars in the cross shore profile. In this applicat ion only the surveys conducted o n the same day (described i n section 6.3 ) were used for the analysis. The simulations were conducted similarly as described in the previous section, forcing the waves to break at the offshore bar, however the bars were not r emoved and the original profile (with one, two or three bars) were used. B reakpoints in depths shallower than 0.6 m were not considered in the analysis, since these shallow breakpoints are not directly associated with sand bar formation and they are more correlated with swash processes. Figure 7 3 shows one example where the model was applied in a profile with two sand bars in Walton County and the model predicted two breaking points. It is important to reinforce that the third breakpoint, at 0.2 m water d epth in this case, was removed from the calculation and not correlated with a sand bar location. The same process was repeated using profiles obtained during f our surveys at Walton County to evaluate the model capability to predict the number of bars in co nditions with high and low waves. The model predict ed well the number of sand bars
94 that occurred in the surveys conducted after Hurricanes Ivan and Katrina, with 71.9% and 70.7% of agreement between the number of bars and the number of breakpoints consider ing all five groups Table 7 1 summarizes the results obtained for each of the five groups and also for all the groups together. These values are consider ed favorabl e since there are differences in the morphology between the profiles. An over prediction of the number of breakpoints instead of under prediction, was observed in the majority of the situations where the number of breakpoints and the average number of bars were different This can be explained because the beach profiles used here are a snapshot representing the energy balance between constructive and destructive forces acting at the beach in th e moment prior to the survey and morphologic features that are evolving to became a bar or that are decaying after being a bar are not considered, using j ust the profile snapshot. This difference is also associated with the abrupt limits used to defin e sand bars i n this dissertation since no visual interpretation s or adjustments w ere applied to compensate transitory stages in the data. Table 7 1. Results obtained with the model simulations for four surveys in Walton County. The results are divided according to the beach profiles that were surveyed in the same date. The values indicate the simulations that correct ly predicted, over predicted and under pred icted the number of bars Surveys 1 Apr Nov 04 Aug 05 Jul 07 Groups 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 Correct 7 6 9 8 5 7 9 9 8 8 5 10 8 10 8 5 5 2 3 4 Over 0 2 0 1 3 0 0 0 1 0 0 0 1 0 1 5 1 0 0 0 Under 6 3 2 2 2 6 2 2 2 3 9 1 2 1 2 3 5 9 8 7 Total (%) Correct 62.5 71.9 70.7 33.3 Over 10.7 1.8 3.4 10.5 Under 26.8 26.3 25.9 56.1 The model is sensitive to changes in the beach slope and the over prediction of sand bars was usually associated with locations with abrupt change s in the beach
95 profile slope where sand bars were not present. One profile was predicted to have two bars by the model (two breakpoints), but is ac tually a profile with only one bar ( Figure 7 4 exemplifies this situation ) The wave shoaled until the offshore bar at 353 m and started to break, the wave shoaled for a second time after reaching the stability criterion and broke the second time at 549 m, in a location without a sand bar but with an abrupt change in the beach slope. Table 7 2. Results obtained with the model simulations for four surveys in Volusia County. The values indicate the simulations that correctly predicted, over predicted and un der predicted the number of bars Survey Jan 02 Feb 06 Sep 06 Oct 07 Correct 6 20 26 18 Over Predict 1 1 0 4 Under Predict 24 10 5 9 Total (%) Correct 19.4 64.5 83.9 58.1 Over Predict 3.2 3.2 0.0 12.9 Under Predict 77.4 32.3 16.1 29.0 The other two surveys to which the model was applied occurred after calm periods and the sand bar depth and position w ere not necessarily associated with the offshore wave con ditions (as in the two surveys described previously), but they were influenced by different wave conditions. Smaller wave heights pass directly without breaking on the outer bar causing minimum effects on them. However these waves influence the behavior of the inner bar, changing their depths and positions. The model results for pre Ivan still presented a good correlation between the number of bars and the number of breakpoints, where 62.5 % of the situations were correctly predicted. However the results in the August 2007 survey were only correct to 33.3 % of the profiles, indicating a methodology limitation for conditions associated with bar formation or migration to deeper areas. When there are large periods between the bar formation and the
96 application o f the model (when surveys are available), the conditions that occurred between them act to modify the profile to the new forces. At Volusia County the surveys were conducted for long period s after major energetic events. The results for all surveys at V olusia County are presented i n T able 7 2. The 2002 survey ha s most of the profiles with only one bar and the model correctly predicted only 19.4 % of the profiles in this situation. Most of the model results, again, over predicted the number of breakpoints during this survey. The presence of a well developed bar trough and steep slope right after the trough in most of the se profiles created the conditions to the waves stabilize at the trough and break for the second time in the steep slope onshore of the tr ough (Figure 7 5 ).
97 Figure 7 1. Representation of the breaking and stability criteria (dashed lines) used in the Dally, Dean and Dalrymple model. The offshore waves start to break after attain ing the breaking criterion and stop breaking and start to sh oal after attain ing the stability criterion. Both criteria are linearly dependent on depth
98 Figure 7 2. Three steps of the s imulation showing the sand bar importance on the break ing processes and wave height decay inside the nearshore zone. A) The f irst step only with the offshore bar B) T he second step including the second bar C) T he third step including the inner bar and applying the simulation in the profile with all bars.
99 Figure 7 3. Model application to predict the number of sand bars in the nearshore zone. The wave height (red line), breaking criteri on ( upper black dashed line) stability criterion (lower black dashed line) and beach profile (black solid line). The third break ( close to the shore) is in shallow waters (smaller than 0.6 m ) and was not considered.
100 Figure 7 4. Simulation where the second b reakpoint is associated with an area with abrupt change in the profile slope and not to a sand bar. The wave height (red line), breaking criteri on ( upper black dashed line) stability criterion (lower black dashed line) and beach profile (black solid line).
101 Figure 7 5. Example of the breakpoint associated with a developed bar / trough system in Volusia County. The wave height (red line), breaking criteri on (black dashed line) stabil ity criterion (lower black dashed line) and beach profile (black solid line). The third break (close to the shore) at 0.15 m water depth was not considered since it is shallower than the 0.6 m limit
102 CHAPTER 8 CONCLUSIONS Sand bar dynamic s we re discussed in terms of long and short time scales, using observation and modeling, to describe their response to different wave conditions and the interaction between bars in systems with multiple bars. In the long term scale, t he simple quantitative mac roscale model addressed previous limitations and showed a considerable improvement compared with the qualitative approach. The period between surveys is irrelevant using th e quantitative method, since the model is based on wave and beach characteristics an d accounts for Neq the number of bars in equilibrium for the current conditions, as a function of the Profile Parameter and Dean Number presented the best results for Walton and Volusia Count ies respectively. The se parameters e mphasiz e that the sediment size plays an important role in determining the average number of bars, since the other two The hurricane season influences are evident for t he bar formation and to the increase in linear morphologies In both counties an increase in the average number of bars occurred mainly during a short period of three months. During the re minder of the year a dec rease in the number of bars was observed, a ssociated with onshore movement of the sand bars when the wave energy is mild. The comparison between the quantitative model at Volusia and Walton Counties showed maximum and mean number of bars greater at Volusia, which was expected and previously observ ed in the measured data since Volusia County has a smaller sediment size, milder beach slope and is exposed to higher wave energy.
103 The short te r m model (Dally, Dean and Dalrymple, 198 5 ) based on the assumption that bars are formed under the breakpoints p redicted the number of bars in the nearshore zone when the beach profiles were conducted close to the wave forcing that generated and influenced the sediment transport. The predictability decreased when surveys were conducted long periods after major storm s and different wave conditions influenced the profiles. The inner and outer bars behavior s are close ly associated in nearshore zones with multiple sand bars. It was observed that t he outer bar acts as a filter to the high waves and it is mainly influenced by major storms. T he inner bar is influenced by a different set of smaller wave heights, since i t is protected by the outer bar The study of the nearshore zone with multiple sand bars should be conducted including all the bars, otherwise misinterpretatio n can cause erroneous solutions for coastal problems or coastal management plans.
104 LIST OF REFERENCES Aagaard, T., Black, K.P., Greenwood, B., 2002. Cross shore suspended sediment transport in the surf zone: a field based parameterization. Marin e Geology 185, 283 302. Aarninkhof, S.G.J., Turner, I.L., Dronkers, T.D.T., Caljouw, M., Nipius, L., 2003. A video based technique for mapping intertidal beach bathymetry. Coastal Engineering 49, 275 289. Absalonsen, L., Dean, R.G., 2011. Characteristics o f the shoreline change along Florida sandy beaches with an example for Palm Beach Coun ty. Journal of Coastal Research 27, 16 26. Bagnold, R.A., 1947. Sand movement by waves: some small scale experiments with sand of very low density. Journal of the Inter na tional Coastal Engineering 27, 447 469. Bascom, W.N., 1951. The relationship between sand size and beach face slope. Transact ions American Geophysical Union 32 (6), 866 874 Batt j es, J. A. 1974. Surf similarity. Proceedings of the 14th International Confer ence on Coastal Engineering, 1978. ASCE, Copenhagen, Denmark, pp. 466 479. Battjes, J. A. Janssen, J. P.F.M. 1978. Energy loss and set up due to breaking of random waves Proceedings of the 16 th International Co nference on Coastal Engineering,1978. ASCE, H amburg, Germany, pp. 569 587. Bauer, B.O., Greenwood, B., 1990. Modification of a linear bar trough system by a sta nding edge wave. Marine Geology 92, 177 204. Benavente, J., Del Ro, L., Anfuso, G., Gracia, F., Reyes, J., 2002. Utility of morphodynamic ch aracterisation in the prediction of beach damage by storms. Journal of Coastal Research 36, 56 64. Benedet, L., Finkl, C., Klein, A., 2004. Morphodynamic classification of beaches on the Atlantic coast of Florida: geographical variability of beach types, b each safety and coastal hazards. Journal of Coastal Research 39, 360 365. Bruun, P., 1954. Coast erosion and the development of beach profiles US Army Corps of Engineers, Beach Erosion Board, Tech. Memo. 44. Charles, L., Dean, R.G., Malakar, S.B., 1994. Sediment data for Florida's east coast University of Florida Department of Coastal and Oceanographic Engineering. UF Institutional Repository. Wetlands General Research Collections. Dally, W.R., 1992. Random breaking waves: field verification of a wave b y wave algorithm for engineering application. Coastal Engineering 16, 369 397.
105 Dally, W.R., Dean, R.G., Dalrymple, R.A., 1985. Wave height variation across beaches of arbitrary profile. Journal of Geophysical Research 90, 11917 11927. Dalry mple, R.A., 1992 Prediction of storm/normal beach profiles. Journal of Waterway, Port Coastal and Ocean Engineering 118, 193 200. Davies, J., 1964. A morphogenic approach to world shorelines. Zeitschrift fr Geomorphologie 8, 127 142. Davis Jr., R.A., Barnard, P., 2003. Morphodynamics of the barrier inlet system, west central Florida. Marine Geology 200, 77 101. beach profile concepts with New Zealand data. Marine Geology 109, 253 266. Dean, R.G., 1973. Heuristic models of sand transport in the surf zone. Proceedings of the Conference in Engineering Dynamics in the Surf Zone 1973. Sydney, Australia, p p. 2 08 2 1 4 Dean, R.G., 1977. Equilibrium beach profiles: US Atlantic and Gulf coasts. Department of Civil Engineering and College of Marine Studies, Ocean Engineering Report 12, University of Delaware Newark, Delaware Dean, R.G., 1986. Coastal Armoring: Effects, Principles and Mitigation. P roceedings of the 20th Coastal Engineering Confer ence, 1986. ASCE, Taipei, Taiwan, pp. 1843 1857. Dean, R.G., 2002. Beach nourishment: theory and practice. Advanced Series on Ocean Engineering, Vol. 18. World Scientific P ress, River Edge, NJ Dean, R.G., Dalrymple, R.A., 2002. Coastal processes with eng ineering applications. Cambridge University Press, New York NY Dolan, T., Dean, R. G. 1985. Multiple longshore sand bars in the upper Chesapeake Bay. Estua rine, Coastal and Shelf Science 21, 727 743. Dyhr Nielsen, M., Sorensen, T., 1970. Sand transport p henomena on coasts with bars. Coastal Engineering 855 865. Elgar, S., Gallagher, E.L., Guza, R., 2001. Nearshore sandbar migration. Journal of Geophysical Research 106, 11623 11627. Frihy, O., Dewidar, K., 2003. Patterns of erosion/sedimentation, heavy min eral concentration and grain size to interpret boundaries of littoral sub cells of the Nile Delta, Egypt. Marine Geology 199, 27 43. Gable, F., 1997. Climate change impacts on Caribbean coastal areas and tourism. Journal of Coastal Research 49 69.
106 Gallaghe r, E.L., Elgar, S., Guza, R., 1998. Observations of sand bar evolution on a natural beach. Journal of Geophys ical Research 103, 3203 3215. Guza, R.T., Inman, D.L., 1975. Edge waves and beach cusps. Jou rnal of Geophysical Research 80 2997 3012. Hanes, D.M., Huntley, D.A., 1986. Continuous measurements of suspended sand concentration in a wave dominated nearshore environm ent. Continental Shelf Research 6, 585 596. Hoefel, F., Elgar, S., 2003. Wave induced sediment transport and sandbar migration. Science 299, 1885 1887. Holman, R., Bowen, A., 1984. Longshore structure of infragravity wave motions. Journal of Geophysical Research 89 (C4), 6446 6452 Holman, R., Sallenger Jr, A., 1993. Sand bar generation: a discussion of the Duck experiment seri es. Journal of C oastal Research 76 92. Houston, J.R., 1996. International tourism and US beaches. Shore and Beach 96 (2) 1 3 Houston, J.R., 2002. The economic value of beaches a 2002 update. Shore and Beach 70, 9 12. Inman, D.L., Elwany, M.H.S., Jenkins, S.A., 1993. Shor erise and bar berm profiles on ocean beaches. Journal of Geophysical Research 98, 18181 18199. Iversen, H., 1952. Waves and breakers in shoaling water Proceedings of the 3 rd International Conference on Coastal Engineering 1952. ASCE, Massachusetts, USA pp. 1 12. Kennedy, A., Slatton, K., Starek, M., Kampa, K., Cho, H., 2010. Hurricane Response of Nearshore Borrow Pits from Airborne Bathymetric Lidar. Journal of Waterway, Port, Coastal, and Ocean Engineering 136, 46 58. Komar, P.D., Gaughan, M.K., 1972. A iry wave theory and breaker height prediction Proceedings of the 13th International Conference on Coastal Engineering. 1972. ASCE, Vancouver, Canada, pp. 405 418. Konicki, K.M., Holman, R., 2000. The statistics and kinematics of transverse sand bars o n an open coast. Marine Geology 169, 69 101. Kriebel, D.L., Dean, R.G., 1993. Convolution method for time dependent beach profile response. Journal of Waterway, Port, Coastal and Ocean Engineering 119, 204 226.
107 Kuriyama, Y., 2002. Medium term bar behavior and associated sediment transport at Hasaki, Japan. Journal of Geophysical Research 107 (C9) 3132. Landsea, C.W., 1993. A climatology of intense (or major) Atlantic hurricanes. Monthly Weath er Review 121, 1703 1713. Larson, M., Kraus, N.C., 1992. Dynamics of longshore bars. Proceedings of the 23 rd International Conference on Coastal Engineering, 1992. ASCE, Venice, Italy, pp. 2219 2219. Lascody, R.L., 1998. East central Florida rip current program. National Weather Digest 22, 25 30. Lippmann, T., Holman, R., 1990. The spatial and temporal variability of sand bar morphology. Journal of Geophysical Research 95, 11575 11590. Lippmann, T., Thornton, E., Reniers, A., 199 5 Wave stress and longshore current on barred profiles. Coastal Dynamics Gdansk, Poland, pp. 4 01 412. Longuet Higgins, M.S., 1983. Wave set up, percolation and undertow in the surf zone. Proceedings of the Royal Society of London. A. Mat hematical and Physical Sciences 390, 283 291. MacMahan, J.H., Thornton, E.B., Reniers, A.J.H.M., 2006. Rip curren t review. Coastal Engineering 53, 191 208. Masselink, G., Short, A.D., 1993. The effect of tide range on beach morphodynamics and morphology: a conceptual beach mod el. Journal of Coastal Research 9, 785 800. McCowan, J., 189 0 On solitary waves. Philosophi cal Magazine, Journal of Science, London, Edinburgh, Dublin. 32, 5. Meckel, L., 1975. Holocene sand bodies in the Colorado delta area, northern Gulf of California. Deltas Models for Exploration. Houston: Geological Society. 239 265. Nicholls, R.J., Birkeme ier, W.A., Hallermeier, R.J., 1996. Application of the depth of closure concept. Proceedings of the 25 th International Conference on Coastal Engineering 1996. ASCE, Orlando, USA pp. 3874 3887. Niedoroda, A.W., Tanner, W., 1970. Preliminary study of trans verse b ars. Marine Geology 9, 41 62. Pape, L., Plant, N., Ruessink, B., 2010. On cross shore migration and equilibrium states of nearshore sandbars. Journal of Geophysical Research 115 ( F03008 )
108 Plant, N.G., Holland, K.T., Puleo, J.A., Gallagher, E.L., 200 4. Prediction skill of nearshore profile evolution models. Journal of Geophysical Research 109 ( C01006 ) Plant, N.G., Holman, R.A., Freilich, M.H., Birkemeier, W., 1999. A simple model for interannual sandbar behavior. Journal of Geophysical Research 104, 15. Rakha, K., Kamphuis, J., 1995. Numerical and hydraulic model study of wave decay on a shelf beach. Coastal E ngineering 24, 275 296. Raubenheimer, B., Elgar, S., Guza, R., 2001. Field observations of wave driven setdown and setup. Journal of Geophysical Research 106, 4629 4638. Ribas, F., Falqus, A., Montoto, A., 2003. Nearshore oblique sand bars. Journal of Geophysical Research 108, 3119. Ruessink, B., Terwindt, J., 2000. The behaviour of nearshore bars on the time scale of years: a conceptual model M arine Geology 163, 289 302. Shand, R., 2007. Bar splitting: system attributes and sediment budget implications for a net offshore migrating bar system. Journal of Coa stal Research 50, 721 730. Shand, R.D., Bailey, D.G., Shepherd, M.J., 1999. An inter site comparison of net offshore bar migration characteristics and environmental conditions. Journal of Coastal Research 750 765. Short, A., 1975. Multiple offshore bars and standing waves. Journal of Geophysical Research 80, 3838 3840. Short, A., 1991. Macro me so tidal beach morphodynamics: An overvi ew. Journal of Coastal Research 7, 417 436. Stephens, S.A., Coco, G., Bryan, K.R., 2011. Numerical Simulation of Wave Setup over Barred Beach Profiles: Implications for Predictability. Journal of Waterway, Port, Coas tal and Ocean Engineering 137, 7. Stive, M. J.F 1984. Energy dissipation in waves breaking on gentle slopes. Coastal Engineering 8, 99 127. Stive, M.J.F., Guillen, J., Capobianco, M., 1996. Bar migration and duneface oscillat ion on decadal scales. Proceedi ngs of the 25th International Co nference on Coastal Engineering, 1996. ASCE, Orlando, USA, pp. 2884 2896. Sunamura, T., Takeda, I., 1984. Landward migration of inner bars. Marine Geology 60, 63 78. Svendsen, I.A., 1984. Mass flux and undertow in a surf zon e. Coastal Engineering 8, 347 365.
109 Thornton, E.B., Guza, R., 1983. Transformation of wave height distribution. Journal of Geophysical Research 88, 5925 5938. Thornton, E.B., Humiston, R.T., Birkemeier, W., 1996. Bar/trough generation on a natural beach. Jo urnal of Geophysical Research 101, 14. Ting, F.C.K., Kirby, J.T., 1994. Observation of undertow and turbulence in a laboratory surf zone. Coastal Engineering 24, 51 80. van Gaalen, J.F., Kruse, S.E., Coco, G., Collins, L., Doering, T., 2011. Observations o f beach cusp evolution at Melbourne Beach, Florida, USA. Geomorphology 129, 131 140. v an Maanen, B., de Ruiter, P., Coco, G., Bryan, K., Ruessink, B., 2008. Onshore sandbar migration at Tairua Beach (New Zealand): Numerical simulations and fie ld measuremen ts. Marine Geology 253, 99 106. Wang, P., Davis, R.A., Kraus, N.C., 1998. Cross shore distribution of sediment texture under breaking waves along low wave energy coasts. Journal of Sedimentary Research 68, 497 506. Weggel, R.J., 1972. Maximum breaker heigh t. Journal of the Waterways, Harbors a nd Coastal Engineering 98, 529 548. Wiegel, R., 1950. Trafficability Studies of Selected Beaches on the Pacific Coast of the United States and Oahu, TH: By RL Wiegel. University of California, Department of Engineering Wiegel, R.L., 1964. Oceanographical engineering. Englewood Cliffs, NJ. Wijnberg, K.M., 2002. Environmental controls on decadal morphologic behaviour of th e Holland coast. Marine Geology 189, 227 247. Wijnberg, K.M., Kroon, A., 2002 Barred beaches. Geomo rphology 48, 103 120. Winant, C.D., Inman, D.L., Nordstrom, C.E., 1975. Description of seasonal beach changes using empirical eigenfunctions. Jou rnal of Geophysical Research 80 (15), 1979 1986 Wright, L., Short, A., 1984. Morphodynamic variability of surf zones and beach es: a synthesis. Marine Geology 56, 93 118.
110 BIOGRAPHICAL SKETCH Luciano Absalonsen received his B.Sc. in oceanography at FURG, Brazil, where he dis covered oceanography and his especial interest in coastal processes arouse. He Elrio Toldo Jr., focusing his research on the study of shoreline evolution and understand ing processes responsible for beach erosion or advancement. He met his Ph.D. adviser, Dr. Robert G. Dean, during a beach nourishment course offered in Brazil and came to the University of Florida, as a Fulbright scholar in August 2008, developing his resea rch on the behavior of sand bars. Luciano is honored to receive his Ph.D. in coastal engineering and hopes to be working in this fascinating world for the rest of his career.