Decadal‐scale gravel beach evolution on a tectonically‐uplifting coast: Wellington,New Zealand

Uplift of the shoreline in tectonically‐active areas can have a profound influence on geomorphology changing the entire process dynamics of the coast as the landforms are removed from the influence of the sea. Over decadal timescales it is possible for the landforms to return to their pre‐earthquake condition and this paper examines the re‐establishment of mixed sand and gravel beaches on the coast of Wellington, New Zealand, subsequent to an uplift event in 1855. Over 60 topographic profiles were surveyed, seven sets of aerial photographs from a 67 year period were mapped and sediment size analyses conducted in order to quantify the nature of beach change following uplift, and associated relative sea level fall. These data were supported by surveys using ground penetrating radar. It is found that uplift raised the gravel beaches out of the swash zone thereby removing them from the littoral zone. Intertidal rocky reefs which occur between each embayment were also uplifted during the same event and completely interrupted the longshore transport system. Continued input of gravel material to the littoral zone allowed beaches to re‐establish sequentially along the coast as each embayment was infilled with sediment. This reconnection of the embayments with the longshore drift system is associated with the beach planform being initially drift dominated during infill but then switching to swash alignment once the embayment becomes infilled. This has resulted in shoreline accretion of over 100 m in some places, at rates of up to 4 m/yr, covering shore protection works built in the past few decades. The ability of the shore to adjust back to its pre‐uplift condition appears to be a function of the accommodation space created during uplift and the rate of sediment supply. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of beach slope breaks on nearshore groundwater dynamics

Interactions between fresh groundwater and seawater affect significantly the nearshore pore water flow, which in turn influences the fate of nutrients and contaminants in coastal aquifers prior to discharge to the marine environment. Field investigations and numerical simulations were carried out to examine the groundwater dynamics in the intertidal zone of a carbonate sandy aquifer on the tropical island of Rarotonga, Cook Islands. The study site was featured by distinct cross‐shore slope breaks on the beach surface. Measured pore water salinities revealed different distributions under the influences of different beach profiles, inland heads, and tidal oscillations. Fresh groundwater was found to discharge around a beach slope break located in the middle area of the intertidal zone. The results indicate a strong interplay between the slope break beach morphology and tidal force in controlling the nearshore groundwater flow and solute transport. The fresh groundwater discharge location was largely determined by the beach morphology in combination with the tidal force. The nearshore groundwater flow can be very sensitive to beach slope breaks, which induce local circulation and flow instabilities. As slope breaks are a common feature of beaches around the world, these results have important, general implications for future studies of nutrients transport and transformations in nearshore aquifers and associated fluxes via submarine groundwater discharge.     

Modelling the effect of waves,weathering and beach development on shore platform development

A mathematical model was used to study shore platform development. Mechanical wave erosion was dependent on such variables as tidal range, wave height and period, breaker height and depth, breaker type, surf zone width and bottom roughness, submarine gradient, rock resistance and the elevational frequency of wave action within the intertidal zone. Also included were the effects of sand and pebble accumulation, cliff height and debris mobility, and downwearing associated with tidal wetting and drying. The occurrence, location and thickness of beaches often depended on initially quite minor variations in platform morphology, but owing to their abrasive or protective effect on underlying rock surfaces, they were able to produce marked differences in platform morphology. Generalizations are difficult, but the model suggests that platform gradient increases with tidal range. Platform width also increases with tidal range with slow downwearing but it decreases with fast downwearing. Platform gradient decreases and width increases with wave energy, and decreasing rock resistance and platform roughness. With low tidal range, platform gradient is generally lower and platform width greater with beaches of fine sand than with gravel, but the relationship is more variable with a high tidal range. Platform width increases and platform gradient decreases with the rate of downwearing on bare surfaces, particularly in low tidal range environments, but the pattern is less clear on beach‐covered platforms. Platforms with large amounts of beach sediment tend to be narrower and steeper than bare platform surfaces. Platform gradient increases and platform width decreases with increasing cliff height and with decreasing cliff debris mobility. Copyright © 2005 John Wiley & Sons, Ltd.     

Wave-forced dynamics in the nearshore river mouths,and swash zones

The role of wave forcing on the main hydro-morphological dynamics evolving in the shallow waters of the nearshore and at river mouths is analyzed. Focus is mainly on the cross-shore dynamics that evolve over mildly sloping barred, dissipative sandy beaches from the storm up to the yearly timescale, at most. Local and non-local mechanisms as well as connections across three main inter-related subsystems of the nearshore – the region of generation and evolution of nearshore bars, river mouths and the swash zone – are analyzed. The beach slope is a major controlling parameter for all nearshore dynamics. A local mechanism that must be properly described for a suitable representation of wave-forced dynamics of all such three subsystems is the proper correlation between orbital velocity and sediment concentration in the bottom boundary layer; while specific dynamics are the wave–current interaction and bar generation at river mouths and the sediment presuspension at the swash zone. Fundamental non-local mechanisms are both infragravity (IG) waves and large-scale horizontal vortices (i.e. with vertical axes), both influencing the hydrodynamics, the sediment transport and the seabed morphology across the whole nearshore. Major connections across the three subsystems are the upriver propagation of IG waves generated by breaking sea waves and swash–swash interactions, the interplay between the swash zone and along-river-flank sediment transport and the evolution of nearshore sandbars. © 2019 John Wiley & Sons, Ltd.     

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated.     

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi‐steady‐state rise in the mean water‐table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed and this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water‐table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright © 2001 John Wiley & Sons, Ltd.     

A proposed mechanism for storm beach sedimentation

One of the major problems of shingle beach dynamics is the method by which coarse gravel is transferred from beach face to storm beach, the latter often lying several metres above high spring tidal levels. The mechanism which is usually proposed, cites the action of plunging breakers as being central to this problem of sediment transfer. However, the nature of net residual fluid force of plunging breakers is deemed unsatisfactory for any substantial upbeach (onshore) sediment transport during storm conditions on gravel beaches. A mechanism is proposed by which high still water levels due to high astronomical tides, onshore storm force winds and allied wave surge, promote shoaling characteristics and beach profile changes which are conducive to spilling breaker development at tidal extremities. It is the net onshore fluid force vector of the spilling breaker overtopping the beach crest which may be the cause of extreme sedimentation events on the storm beach. An example of such sediment transfer is given for a known storm beach sedimentation event at Llanrhystyd gravel beach, West Wales, during February 1974. Process variables were monitored on several days allowing the use of an inshore breaker steepness criterion, to indicate positions in the tidal regime where plunging breakers give way to spilling forms. This example serves to suggest that more attention should be given to the nature and characteristics of shoaling waves, especially in respect of breaker type, when examining problems of shingle beach dynamics and sedimentation.     

Modelling the time-varying extent of groundwater seepage on tidal beaches

The outcrop of groundwater on tidal beaches distinguishes an upper unsaturated region from a lower saturated region of the intertidal profile. Since the 1940s, it has been recognized that the extent of groundwater seepage at the beach face is one factor determining the tendency for erosive or accretionary conditions to prevail. As a primary step towards incorporating bed saturation characteristics within cross-shore sediment transport models, this paper (and accompanying program disk) details a simple model to simulate the time-varying extent of seepage face development across tidal beaches. From a comparison with field results obtained on the macrotidal Central Queensland (Australia) coast, the model appears to provide an encouraging degree of predictive capability. The model also assists in highlighting the sensitivity of seepage face development to varying beach face, tide and wave characteristics.     

Storm response of a mixed sand gravel beach ridge plain under falling relative sea levels: A stratigraphic investigation using ground penetrating radar

Beach ridge stratigraphy can provide an important record of both sustained coastal progradation and responses to events such as extreme storms, as well as evidence of earthquake induced sediment pulses. This study is a stratigraphic investigation of the late Holocene mixed sand gravel (MSG) beach ridge plain on the Canterbury coast, New Zealand. The subsurface was imaged along a 370 m shore-normal transect using 100 and 200 MHz ground penetrating radar (GPR) antennae, and cored to sample sediment textures. Results show that, seaward of a back-barrier lagoon, the Pegasus Bay beach ridge plain prograded almost uniformly, under conditions of relatively stable sea level. Nearshore sediment supply appears to have created a sustained sediment surplus, perhaps as a result of post-seismic sediment pulses, resulting in a flat, morphologically featureless beach ridge plain. Evidence of a high magnitude storm provides an exception, with an estimated event return period in excess of 100 years. Evidence from the GPR sequence combined with modern process observations from MSG beaches indicates that a palaeo-storm initially created a washover fan into the back-barrier lagoon, with a large amount of sediment simultaneously moved off the beach face into the nearshore. This erosion event resulted in a topographic depression still evident today. In the subsequent recovery period, sediment was reworked by swash onto the beach as a sequence of berm deposit laminations, creating an elevated beach ridge that also has a modern-day topographic signature. As sediment supply returned to normal, and under conditions of falling sea level, a beach ridge progradation sequence accumulated seaward of the storm feature out to the modern-day beach as a large flat, uniform progradation plain. This study highlights the importance of extreme storm events and earthquake pulses on MSG coastlines in triggering high volume beach ridge formation during the subsequent recovery period. © 2019 John Wiley & Sons, Ltd.     

Morphological and sedimentary features of sandy-muddy transitional beaches in estuaries and bays along mesotidal to macrotidal coasts

Sandy-muddy transitional beaches (SMT-Beaches), representing the transition from sandy beaches to tidal mudflats, should theoretically develop very different morphological and sedimentological characteristics in river estuaries and in semi-enclosed bays due to their contrasting dynamic sedimentary environments. Evidence, however, is rare in the scientific literature. To reveal these morphological and sedimentary differences, the sand–mud transition (SMT) boundary distribution, beach profiles, and surface and downcore sediment grain-size compositions of 27 SMT-Beaches located along mesotidal to macrotidal coasts of the western Taiwan Strait, southeastern China, were investigated. The results show that typical estuarine SMT-Beaches are mainly characterized by an ambiguous SMT, a long distance between the SMT and the coastline (31–302 m), lower SMT and inflection point altitudes (average –0.76 m and –0.04 m), and lower upper beach gradients (~0.068) with fine sand. Estuarine SMT-Beach sediments display clear interbedded mud and sand layers, implying potential SMT migrations over various timescales. By contrast, typical bay SMT-Beaches are characterized by distinct SMT, a short distance between the SMT and the coastline (11–52 m), higher SMT and inflection point altitudes (~0.24 m and ~0.35 m), and narrower upper beaches with higher gradients (~0.095) and coarse sand. Bay SMT-Beaches present relatively stable sedimentary sequences and a narrow gravel belt surrounding the inflection point and/or SMT. These morphological and sedimentary differences between the two SMT-Beach types are initially constrained by sediment supply and transport and are further affected by tide conditions and wave climate. Sediment supply and transport predominately control the sediment structures, while the tidal range strongly influences spatial variations in SMT distances. Wave climate normally drives SMT altitude variations. This study highlights the morphological and sedimentary differences in SMT-Beaches in estuaries and bays, providing important knowledge for further revealing their morphodynamic processes and potential future nourishment. © 2020 John Wiley & Sons, Ltd.     

Aeolian transport rates across raked and unraked beaches on a developed coast

Wind characteristics and aeolian transport were measured on a naturally evolving beach and dune and a nearby site where the beach is raked and sand‐trapping fences are deployed. The beaches were composed of moderately well sorted to very well sorted fine to medium sand. The backshore at the raked site was wider and the foredune was more densely vegetated and about 1 m higher than at the unraked site. Wind speeds were monitored using anemometers placed at 1 m elevation and sand transport was monitored using vertical traps during oblique onshore, alongshore and offshore winds occurring in March and April 2009. Inundation of the low backshore through isolated swash channels prevented formation of a continuously decreasing cross‐shore moisture gradient. The surface of the berm crest was dryer than the backshore, making the berm crest the greatest source of offshore losses during offshore winds. The lack of storm wrack on the raked beach reduced the potential for sediment accumulation seaward of the dune crest during onshore winds, and the higher dune crest reduced wind speeds and sediment transport from the dune to the backshore during offshore winds. Accretion at wrack seaward of the dune toe on the unraked beach resulted in a wider dune field and higher, narrower backshore. Although fresh wrack is an effective local trap for aeolian transport, wrack that becomes buried appears to have little effect as a barrier and can supply dry sand for subsequent transport. Aeolian transport rates were greater on the narrower but dryer backshore of the unraked site. Vegetation growth may be necessary to trap sand within zones of buried wrack in order to allow new incipient foredunes to evolve. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of multiple intertidal bars over semi‐diurnal and lunar tidal cycles,North Lincolnshire,England

Multiple intertidal bars are common features of wave‐dominated sandy beaches, yet their short‐term (<1 month) and small‐scale (<1 km) morphology and dynamics remain poorly understood. This study describes the morphodynamics of multiple intertidal bars in North Lincolnshire, England, during single and lunar tidal cycles under two contrasting conditions – first when significant wave height was <0·5 m and second when significant wave height frequently exceeded 1 m. The relative importance of swash, surf and shoaling processes in determining morphological change was examined using detailed field observations and a numerical model. The beach featured four intertidal bars and both cross‐shore and longshore bar morphology evolved during the field investigation, particularly under medium to high wave‐energy conditions. Numerical modelling suggests shoaling processes are most common on the seaward two bars under calm wave conditions (H_s < 0·5 m) and that surf zone processes become more common during neap tides and under more energetic (H_s < 0·5 m) conditions. Surf processes dominate the inner two bars, though swash influence increases in a landward direction. The numerical modelling results combined with low tide survey data and high‐resolution morphological measurements strongly suggest changes in the intertidal bar morphology are accomplished by surf zone processes rather than by shoaling wave or swash processes. This is because shoaling waves do not induce significant sediment transport to have any morphological effect, whereas swash action generally does not have enough scope to act as the swash zone is much narrower than the surf zone. It was found, however, that the absolute rate of morphological change under swash action and surfzone processes are of similar magnitudes and that swash action may induce a significant amount of local morphological change when the high tide mark is located on the upper bar, making this process important for bar morphodynamics. Copyright © 2007 John Wiley & Sons, Ltd.     

Downdrift coarsening of beach foreshore sediments at tidal inlets: An example from the coast of New Jersey

Analysis of grain size statistics of upper foreshore sediments on sand beaches at two tidal inlets in New Jersey, U.S.A. reveals that sediments are coarser at beaches flanking the inlets than updrift, although sediments become finer downdrift at the broad, regional scale. The local reversal of the regional trend in size grading is attributed to: (1) the offshore diversion of the finer sands along the surf zone on the ebb tidal delta, and (2) the removal of the finer sands from the inlet flank beach caused by low wave energy conditions at low stages of the tide and by deflation. Sediments thus become coarser at inlet beaches as a result of alterations in the interaction of waves with the beach and as a result of aeolian processes, not solely as a result of increased tidal current velocities as previously reported. The distance along the New Jersey barrier islands over which inlet processes are likely to affect changes in sediment size updrift averages less than 1100 m, but the impacts of inlets on the sedimentary record can be extended greater distances as a result of inlet migration.     

A cross-shore suspended sediment transport shape function parameterisation for natural beaches

A new field-based parameterisation (‘shape function’) describing the distribution of cross-shore suspended sediment transport across a beach profile is presented. Time-averaged and depth-integrated suspended sediment fluxes were measured over 39 tides at Sennen Cove, Cornwall, UK, for a range of wave conditions (offshore significant wave heights 0.1–2.5 m). The suspended sediment flux data were heuristically separated into four transport components: (1) mean flux in the surf/shoaling zone; (2) oscillatory flux in the surf/shoaling zone; (3) onshore flux in the swash/inner surf zone and (4) offshore flux in the swash/inner surf zone. Each of these transport components was related to the local water depth (h) normalised by the breakpoint depth (h_b) and the four resulting suspended transport shape functions were combined to form a total suspended load shape function. Each shape function component is scaled independently by the wave energy level through h_b. The total suspended load shape function predicts onshore sediment transport under low-energy conditions, with peaks at the breakpoint and in the swash zone, in agreement with the field observations. Under high-energy conditions the total suspended load shape function predicts onshore transport in the shoaling zone, offshore transport in the surf zone and onshore transport in the inner swash zone.     

Swash-zone morphodynamics

Hydrodynamic forcing with respect to sediment transport and morphologic change, paying particular attention to relevant swash asymmetries, is reviewed. The hydrodynamics are categorized into their individual effects: high- and low-frequency motions, bores and turbulence, in/exfiltration, shear stresses and friction coefficients. Individual effects are then related to their potential for driving or influencing sediment transport and morphological change. Additional concepts such as settling/scour lag and sediment advection that have largely been ignored are also discussed. A simple framework is presented for the morphological response of the beachface under swash zone hydrodynamic processes. The framework acknowledges that the beachface cannot be considered in isolation from the surf zone and that the two zones are strongly linked through feedback processes. The swash zone itself is also a morphodynamic system and morphological response occurs as a result of disequilibrium between the beachface gradient and asymmetries in the swash hydrodynamics. Any beachface morphological development in response to such equilibrium will have direct and indirect effects on swash hydrodynamics and sediment transport processes. It is concluded that the two issues requiring most urgent research attention with regards to swash zone sediment transport processes are the roles of sediment advection and longshore swash motion.     

Overwash processes along a gravel beach in south-east Ireland

Recent overwashing, and to a more limited extent overtopping, has occurred on a gravel-dominated beach just west of Carnsore Point in south-east Ireland. At relatively low points swash excursions have breached the duneline and extensive gravel fans have formed on landward-facing slopes. Also, at interbreach positions, swash ramps have been built-up in front of previously eroding dune cliffs, allowing overtopping at heights of up to 5 m above MHWOST, spreading a veneer of sediment over the adjacent dune ridge. In both cases events appear to be preceded by erosion and backbeach sand and gravel deposition. Stratigraphic, textural and fabric analyses assist in identifying environmental conditions that prevailed during overwash throat and fan, and overtopping ramp, sedimentation.     

Tidal influence on BTEX biodegradation in sandy coastal aquifers

A numerical study was conducted to investigate the influence of tides on the fate of terrestrially derived BTEX discharging through an unconfined aquifer to coastal waters. Previous studies have revealed that tide-induced seawater circulations create an active salt–freshwater mixing zone in the near-shore aquifer and alter the specific subsurface pathway for contaminants discharging to the coastal environment. Here the coupled density-dependent flow and multi-species reactive transport code PHWAT was used to examine the impact of these tidal effects on the aerobic biodegradation of BTEX released in a coastal aquifer and its subsequent loading to coastal waters. Simulations indicated that tides significantly enhance BTEX attenuation in the near-shore aquifer. They also reduce the rate of chemical transfer from the aquifer to the ocean and exit concentrations at the beach face. For the base case consisting of toluene transport and biodegradation, 79% of toluene initially released in the aquifer was attenuated prior to discharge with tides present, compared to only 1.8% for the non-tidal case. The magnitude of tidal forcing relative to the fresh groundwater flow rate was shown to influence significantly the extent of biodegradation as it controls the intensity of salt–freshwater mixing, period of exposure of the contaminant to the mixing zone and rate of oxygen delivery to the aquifer. The oxygen available for biodegradation also depends on the rate at which oxygen is consumed by natural processes such as organic matter decomposition. While simulations conducted with heterogeneous conductivity fields highlighted the uncertainties associated with predicting contaminant loadings, the study revealed overall that BTEX may undergo significant attenuation in tidally influenced aquifers prior to discharge.     

Tidal and meteorological forcing of sediment transport in tributary mudflat channels

Field observations of flow and sediment transport in a tributary channel through intertidal mudflats indicate that suspended sediment was closely linked to advection and dispersion of a tidal salinity front. During calm weather when tidal forcing was dominant, high concentrations of suspended sediment advected up the mudflat channel in the narrow region between salty water from San Francisco Bay and much fresher runoff from the small local watershed. Salinity and suspended sediment dispersed at similar rates through each tidal inundation, such that during receding ebbs the sediment pulse had spread spatially and maximum concentrations had decreased. Net sediment transport was moderately onshore during the calm weather, as asymmetries in stratification due to tidal straining of the salinity front enhanced deposition, particularly during weaker neap tidal forcing. Sediment transport by tidal forcing was periodically altered by winter storms. During storms, strong winds from the south generated wind waves and temporarily increased suspended sediment concentrations. Increased discharge down the tributary channels due to precipitation had more lasting impact on sediment transport, supplying both buoyancy and fine sediment to the system. Net sediment transport depended on the balance between calm weather tidal forcing and perturbations by episodic storms. Net transport in the tributary channel was generally off-shore during storms and during calm weather spring tides, and on-shore during calm weather neap tides.     

Beach erosion and recovery during consecutive storms at a steep‐sloping,meso‐tidal beach

This study analyses beach morphological change during six consecutive storms acting on the meso‐tidal Faro Beach (south Portugal) between 15 December 2009 and 7 January 2010. Morphological change of the sub‐aerial beach profile was monitored through frequent topographic surveys across 11 transects. Measurements of the surf/swash zone dimensions, nearshore bar dynamics, and wave run‐up were extracted from time averaged and timestack coastal images, and wave and tidal data were obtained from offshore stations. All the information combined suggests that during consecutive storm events, the antecedent morphological state can initially be the dominant controlling factor of beach response; while the hydrodynamic forcing, and especially the tide and surge levels, become more important during the later stages of a storm period. The dataset also reveals the dynamic nature of steep‐sloping beaches, since sub‐aerial beach volume reductions up to 30 m³/m were followed by intertidal area recovery (–2 < z < 3 m) with rates reaching ~10 m³/m. However, the observed cumulative dune erosion and profile pivoting imply that storms, even of regular intensity, can have a dramatic impact when they occur in groups. Nearshore bars seemed to respond to temporal scales more related to storm sequences than to individual events. The formation of a prominent crescentic offshore bar at ~200 m from the shoreline appeared to reverse the previous offshore migration trend of the inner bar, which was gradually shifted close to the seaward swash zone boundary. The partially understood nearshore bar processes appeared to be critical for storm wave attenuation in the surf zone; and were considered mainly responsible for the poor interpretation of the observed beach behaviour on the grounds of standard, non‐dimensional, morphological parameters. Copyright © 2011 John Wiley & Sons, Ltd.