Investigating Changes in Aeolian Sediment Transport at Coastal Dunes and Sand Trapping Fences: A Field Study on the German Coast

: For the restoration and maintenance of beach and dune systems along the coast, knowledge of aeolian sediment transport and its interaction with coastal protection measures is required. As a nature-based solution, sand trapping fences can be an integral part of coastal protection measures initiating foredune development. There are few detailed studies on aeolian sediment transport rates on coastal dunes and sand trapping fences available to date. Thus, in this work, we present the results of ﬁeld experiments conducted at the beach, coastal dune, and sand trapping fence on the East Frisian island Langeoog. The vertical sediment ﬂux proﬁle was measured by vertical mesh sand traps, and saltiphones measured the instantaneous sediment transport. A meteorological station was set up to obtain wind data. On the beach, dune toe, and dune crest, the stationary wind proﬁle can be described well by the law of the wall. Saturated aeolian sediment transport rates on the beach and dune toe were predicted by widely used empirical models. Between the sand trapping fence, these empirical transport models could not be applied, as no logarithmic wind proﬁle existed. The upwind sediment supply reduced after each brushwood line of the sand trapping fence, thereby, leading to increased deviation from the saturated conditions. recommend to extend the knowledge of the inﬂuences of sand trapping fences on aeolian sediment transport in standardized wind tunnel experiments to gain more quantitative data. To evaluate the inﬂuence of sand trapping fences on dune development for longer timescales, aerial drone investigations of the coastal dune should be conducted on the study site Langeoog. For this purpose, an analysis of a series of digital elevation models at the coastal dune and at the sand trapping fence will be performed to investigate the dune volume changes around the fence depending on each sand trapping conﬁguration. The aim of the future investigation is to evaluate the e ﬃ ciency of each sand trapping conﬁguration under di ﬀ erent prevailing wind conditions.


Introduction
Coastal protection measures on the East Frisian Islands in the German North Sea generally involve coastal dunes, coastal dikes, and associated coastal protection structures [1,2]. The main roles of the coastal dunes on the East Frisian Islands are to provide flood protection for flood-risk areas against storm surges and to secure the existence of these islands [2][3][4][5]. These islands form a natural barrier system that reduces the wave energy during storm surges on the mainland coast in Lower Saxony. The islands are, therefore, maintained to reduce flood risks for the hinterland. The maintenance of the coastal dunes at the German North Sea coast generally involves sand nourishment at the coastal dunes, beach, and foreshore [2,6]. Over the past decades, sand nourishments have been an effective method for ensuring coastal protection of the island, to bridge temporary phases of sedimentary deficit caused by sandbar relocations [2,7]. The planting of vegetation and installing sand trapping fences to stabilize the natural dynamics of the dune system represent an additional component of coastal protection on the East Frisian Islands [8].
Nature-based solutions in coastal engineering such, as installing sand trapping fences aim to use natural processes and natural resources to contribute to coastal protection. Nature-based solutions are becoming increasingly popular due to the growing environmental awareness among the population, (1) Describing the wind profile, shear velocity, and aerodynamic roughness length for the beach, dune toe, and dune crest of the foredunes as well as before and between the sand trapping fence. (2) Quantifying the vertical sediment flux profile for these measurement locations.
(3) Evaluating the performance of aeolian sediment transport models in comparison to the measured sediment transport rates under consideration of prevailing boundary conditions. (4) Evaluating the downwind evolution of instantaneous sediment transport at these measurement locations. (5) Describing the influence of sand trapping fences on aeolian sediment transport.  [41], with permission from © GeoBasis-DE/ BKG, 2020) belonging to the (b) East Frisian Islands, federal state Lower Saxony (adapted from [41], with permission from © GeoBasis-DE/ BKG, 2020) as part of the (c) island Langeoog (adapted from [42] with permission from © dl-de/by-2-0, 2020), as well as the location of the (d) study site in the Northern part of Langeoog.
The East Frisian Islands are located in the German North Sea and form a barrier island system that developed in the Holocene. Seven permanently inhabited sandy barrier islands stretch over ~90 km and are dissected by six main tidal inlets [2,6]. They are separated from the mainland (distance to mainland ~3-20 km) by a large muddy intertidal area with tidal flats. The Holocene sea level rise, the abundant sediment availability, large-scale embankments, and storm surges have contributed significantly to the development and shape of these barrier islands [43,44]. Up to now, the longshore migration of the intertidal sandbars in the coastal foreland from west to east in irregular periods of decades due to the tides, currents, wave loads, and storm surges caused a significant change in the sediment budget for the East Frisian Islands and resulted in a temporary supply or deficit of sand for these islands [1,2,43,45].
The East Frisian island Langeoog covers an area of ~20 km² and has ~20.3 km long coastal dunes stretching from the southwestern side of the island northward toward the east. No grey coastal protection structures are present. From the seaward side to inland, Langeoog's morphology consists of a natural sandy beach followed by foredunes and older dune landscapes with dunes up to 20 m. In the southwest, Langeoog has extensively large areas with salt marshes [2,6,46,47].
Langeoog can be categorized as a mesotidal barrier island. Tides in this area are semi-diurnal with a tidal range of around TR ~2.7 m at the port entrance [47][48][49]. At each tidal cycle, the tidal basins with intertidal areas and supratidal salt marshes are filled and subsequently emptied [50]. The main wave direction is from southwest to northwest [51].
The climate is temperate [52] with an annual average temperature of ~11 °C (Norderney, January 2015 to December 2019) and an annual average rainfall of ~754 mm (Norderney, January 2015 to December 2019). The wind direction are mainly from southwest to northwest characterized by wind velocities up to ~30 m/s [53]. Figure 2 shows a close-up aerial photograph of the sand trapping fence. The sand trapping fence was installed from 18 May 2020 to 19 May 2020 by the Lower Saxony Water Management, Coastal Protection and Nature Conservation Agency and was made from locally available brushwood that was put into the sand and arranged in lines parallel and orthogonal to the coastal dune [2,8]. In Figure  2a, an aerial drone photograph and (b) a close-up of the sand trapping fence at the study site of Langeoog are shown.  [41], with permission from © GeoBasis-DE/ BKG, 2020) belonging to the (b) East Frisian Islands, federal state Lower Saxony (adapted from [41], with permission from © GeoBasis-DE/ BKG, 2020) as part of the (c) island Langeoog (adapted from [42] with permission from © dl-de/by-2-0, 2020), as well as the location of the (d) study site in the Northern part of Langeoog.
The East Frisian Islands are located in the German North Sea and form a barrier island system that developed in the Holocene. Seven permanently inhabited sandy barrier islands stretch over~90 km and are dissected by six main tidal inlets [2,6]. They are separated from the mainland (distance to mainland~3-20 km) by a large muddy intertidal area with tidal flats. The Holocene sea level rise, the abundant sediment availability, large-scale embankments, and storm surges have contributed significantly to the development and shape of these barrier islands [43,44]. Up to now, the longshore migration of the intertidal sandbars in the coastal foreland from west to east in irregular periods of decades due to the tides, currents, wave loads, and storm surges caused a significant change in the sediment budget for the East Frisian Islands and resulted in a temporary supply or deficit of sand for these islands [1,2,43,45].
The East Frisian island Langeoog covers an area of~20 km 2 and has~20.3 km long coastal dunes stretching from the southwestern side of the island northward toward the east. No grey coastal protection structures are present. From the seaward side to inland, Langeoog's morphology consists of a natural sandy beach followed by foredunes and older dune landscapes with dunes up to 20 m. In the southwest, Langeoog has extensively large areas with salt marshes [2,6,46,47].
Langeoog can be categorized as a mesotidal barrier island. Tides in this area are semi-diurnal with a tidal range of around T R~2 .7 m at the port entrance [47][48][49]. At each tidal cycle, the tidal basins with intertidal areas and supratidal salt marshes are filled and subsequently emptied [50]. The main wave direction is from southwest to northwest [51].
The climate is temperate [52] with an annual average temperature of~11 • C (Norderney, January 2015 to December 2019) and an annual average rainfall of~754 mm (Norderney, January 2015 to December 2019). The wind direction are mainly from southwest to northwest characterized by wind velocities up to~30 m/s [53]. Figure 2 shows a close-up aerial photograph of the sand trapping fence. The sand trapping fence was installed from 18 May 2020 to 19 May 2020 by the Lower Saxony Water Management, Coastal Protection and Nature Conservation Agency and was made from locally available brushwood that was put into the sand and arranged in lines parallel and orthogonal to the coastal dune [2,8].
In Figure 2a, an aerial drone photograph and (b) a close-up of the sand trapping fence at the study site of Langeoog are shown. When the brushwood bushes were delivered, they were bundled by wires. As far as possible, the wires were removed from the bundles of brushwood. The investigations were conducted in front of field 1, in field 1, and in field 2 (all part of configuration 1). These fields were characterized by three bundles of brushwood per running meter orthogonal and two bundles of brushwood per running meter parallel to the coastal dune [54]. Currently, the design of the protruding brushwood height as well as the porosity of the sand trapping fences as well as their position in relation to the beach profile, and the length of the resulting fields between the sand trapping fences on the East Frisian Islands are based on empirical knowledge. Figure 3 shows the top view of the study site with all four measurement locations: the beach, the area close to the dune toe, the area on the dune crest of the foredunes, and the sand trapping fence at the dune toe. The seaward side of the foredune is in part protected by a sand trapping fence. All four measurement locations differ a lot in their prevailing boundary conditions. A closer visual impression of the locations (a) beach, (b) dune toe, (c) dune crest, and (d) sand trapping fence is given in Figure 4a. When the brushwood bushes were delivered, they were bundled by wires. As far as possible, the wires were removed from the bundles of brushwood. The investigations were conducted in front of field 1, in field 1, and in field 2 (all part of configuration 1). These fields were characterized by three bundles of brushwood per running meter orthogonal and two bundles of brushwood per running meter parallel to the coastal dune [54]. Currently, the design of the protruding brushwood height as well as the porosity of the sand trapping fences as well as their position in relation to the beach profile, and the length of the resulting fields between the sand trapping fences on the East Frisian Islands are based on empirical knowledge. Figure 3 shows the top view of the study site with all four measurement locations: the beach, the area close to the dune toe, the area on the dune crest of the foredunes, and the sand trapping fence at the dune toe. The seaward side of the foredune is in part protected by a sand trapping fence. When the brushwood bushes were delivered, they were bundled by wires. As far as possible, the wires were removed from the bundles of brushwood. The investigations were conducted in front of field 1, in field 1, and in field 2 (all part of configuration 1). These fields were characterized by three bundles of brushwood per running meter orthogonal and two bundles of brushwood per running meter parallel to the coastal dune [54]. Currently, the design of the protruding brushwood height as well as the porosity of the sand trapping fences as well as their position in relation to the beach profile, and the length of the resulting fields between the sand trapping fences on the East Frisian Islands are based on empirical knowledge. Figure 3 shows the top view of the study site with all four measurement locations: the beach, the area close to the dune toe, the area on the dune crest of the foredunes, and the sand trapping fence at the dune toe. The seaward side of the foredune is in part protected by a sand trapping fence. All four measurement locations differ a lot in their prevailing boundary conditions. A closer visual impression of the locations (a) beach, (b) dune toe, (c) dune crest, and (d) sand trapping fence is given in Figure 4a.  For the measurement location beach, the sand was relatively wet (M = 7.08-10.38 Vol%). A few shell fragments were identified on the beach, see Figure 4a. Close to the dune toe, the sand was relatively dry with a moisture content of maximum M = 1.8 Vol% and without shells or shell fragments, see Figure 4b. Many small sand ripples were recognized moving downwind. On the dune crest, dune grass of varying heights (with a medium height of about 0.5 m) and irregularities in the topography of the dune crest area of about 2 m were present, see Figure 4c. In front of the sand trapping fence, sand was relatively wet with a maximum moisture value of M = 8.73 Vol%. Sand was deposited in the form of a large, rippled sandbar, see Figure 4d, in front of the sand trapping fence. In the sand trapping fence, the moisture content of the sand varied between M = 1.80-8.90 Vol%, and no sand ripples were present.

Experimental Setup and Instrumentation
The field experiments began immediately after the completion of the installation of the sand trapping fence and were conducted from 19 May 2020 to 24 May 2020. They were focused on the aeolian transport processes and rates in correlation with the wind regime. The mean sediment transport rates were measured with three modified vertical mesh sand traps constructed according to Sherman et al. (2014) [55], covering a height of z = 0.3 m above the surface, see Figure 5a. The vertical mesh sand traps consisted of six rectangular aluminium tubes arranged one above the other. Each rectangular aluminium tube was b = 0.1 m wide, h = 0.05 m high, and l = 0.25 m long, with edges e = 2 mm thick. To enable sand to be trapped, each tube was covered with a bag made of nylon mono filament with openings of size = 50 µ m. A long nylon bag l = 0.75 m was used for the lowest tube, for all other heights the nylon back was l = 0.5 m long. A threaded rod and rubber bands fixed the nylon For the measurement location beach, the sand was relatively wet (M = 7.08-10.38 Vol%). A few shell fragments were identified on the beach, see Figure 4a. Close to the dune toe, the sand was relatively dry with a moisture content of maximum M = 1.8 Vol% and without shells or shell fragments, see Figure 4b. Many small sand ripples were recognized moving downwind. On the dune crest, dune grass of varying heights (with a medium height of about 0.5 m) and irregularities in the topography of the dune crest area of about 2 m were present, see Figure 4c. In front of the sand trapping fence, sand was relatively wet with a maximum moisture value of M = 8.73 Vol%. Sand was deposited in the form of a large, rippled sandbar, see Figure 4d, in front of the sand trapping fence. In the sand trapping fence, the moisture content of the sand varied between M = 1.80-8.90 Vol%, and no sand ripples were present.

Experimental Setup and Instrumentation
The field experiments began immediately after the completion of the installation of the sand trapping fence and were conducted from 19 May 2020 to 24 May 2020. They were focused on the aeolian transport processes and rates in correlation with the wind regime. The mean sediment transport rates were measured with three modified vertical mesh sand traps constructed according to Sherman et al. (2014) [55], covering a height of z = 0.3 m above the surface, see To enable sand to be trapped, each tube was covered with a bag made of nylon mono filament with openings of size = 50 µm. A long nylon bag l = 0.75 m was used for the lowest tube, for all other heights the nylon back was l = 0.5 m long. A threaded rod and rubber bands fixed the nylon bags to the frame. This design of the vertical mesh sand traps was used in field experiments by Campos (2018) [56]. It is a cost effective measurement tool with good efficiency results [55].
direction. The saltiphones were calibrated by checking that all six saltiphones gave nearly the same amount of counted grains hitting the microphone for a defined sediment probe. For this purpose, a weighed mass of quartz sand was poured into the housing of each saltiphone at a defined distance. Detailed information about the principle of saltiphones is given by Spaan and van der Abeele (1991) [58] and Goossens et al. (2018) [59].
In this field campaign, a total of three saltiphone stations, each station consisting of two saltiphones, were used. They were placed in such a way that the centre of the microphones was arranged at z = 0.05 m and z = 0.15 m above the bed surface. Thereby, the instantaneous sediment transport at different altitudes could be estimated. However, in weak wind conditions, it can be possible that the transport layer is so low that most particles pass underneath the lowest sensor [22].
A meteorological station, see also Figure 5b, was set up to obtain the wind velocity and wind direction at the study site. The sensors OMEGA WMS-23 consisted of a cup anemometer to register the wind velocity and a wind vane to measure the wind direction relative to north. The wind sensors were mounted on a vertical bar at predefined heights of z = 0.13, 0.17, 0.45, 0.59, 0.91, and 1.95 m and the wind vanes were at 0.24, 0.28, 0.56, 0.70, 1.02, and 2.06 m, respectively. The wind data were recorded with a measuring rate of f = 10 Hz. The measuring accuracy of the wind vane is stated with ± 3° and for the wind anemometer it is the larger value of ± 3% or 0.45 m/s [60].
An unmanned aerial vehicle (UAV, manufacturer DJI Phantom 4 with real time kinematics) was used to obtain ortho-image data, which was referenced by four ground control points. Ortho-images taken by the drone achieved an accuracy of 1 cm + 1 parts-per-million (ppm) of distance between the rover and base horizontally and 1.5 cm + 1 ppm of distance between the rover and base vertically [61]. The exact position of all measuring equipment and reference points were registered additionally using a global navigation satellite system JAVAD GNSS Receiver SigmaD with the same accuracy horizontally and vertically as the drone [62].
A sedimentological analysis according to the specifications in ISO 17892-4 [63] was conducted to determine the grain size distribution at the study site. The spatially variable moisture contents at Six saltiphones from the company Eijkelkamp, see also Figure 5a, were used to measure the instantaneous sediment transport. The saltiphones consist of a microphone, installed in a stainless-steel housing. Each sediment grain hitting the microphone is counted. The maximum cumulative number of impacts per second is generally limited to 1000. The overall measurement accuracy is given with 5% [57]. The data of the saltiphones were recorded by a data logger that was built by the IWW. Two wind vanes on the back of the tube ensured that the microphone was positioned in the wind direction. The saltiphones were calibrated by checking that all six saltiphones gave nearly the same amount of counted grains hitting the microphone for a defined sediment probe. For this purpose, a weighed mass of quartz sand was poured into the housing of each saltiphone at a defined distance. Detailed information about the principle of saltiphones is given by Spaan and van der Abeele (1991) [58] and Goossens et al. (2018) [59].
In this field campaign, a total of three saltiphone stations, each station consisting of two saltiphones, were used. They were placed in such a way that the centre of the microphones was arranged at z = 0.05 m and z = 0.15 m above the bed surface. Thereby, the instantaneous sediment transport at different altitudes could be estimated. However, in weak wind conditions, it can be possible that the transport layer is so low that most particles pass underneath the lowest sensor [22].
A meteorological station, see also Figure 5b, was set up to obtain the wind velocity and wind direction at the study site.
The sensors OMEGA WMS-23 consisted of a cup anemometer to register the wind velocity and a wind vane to measure the wind direction relative to north. The wind sensors were mounted on a vertical bar at predefined heights of z = 0.13, 0.17, 0.45, 0.59, 0.91, and 1.95 m and the wind vanes were at 0.24, 0.28, 0.56, 0.70, 1.02, and 2.06 m, respectively. The wind data were recorded with a measuring rate of f = 10 Hz. The measuring accuracy of the wind vane is stated with ± 3 • and for the wind anemometer it is the larger value of ± 3% or 0.45 m/s [60].
An unmanned aerial vehicle (UAV, manufacturer DJI Phantom 4 with real time kinematics) was used to obtain ortho-image data, which was referenced by four ground control points. Ortho-images taken by the drone achieved an accuracy of 1 cm + 1 parts-per-million (ppm) of distance between the rover and base horizontally and 1.5 cm + 1 ppm of distance between the rover and base vertically [61]. The exact position of all measuring equipment and reference points were registered additionally using a global navigation satellite system JAVAD GNSS Receiver SigmaD with the same accuracy horizontally and vertically as the drone [62].
A sedimentological analysis according to the specifications in ISO 17892-4 [63] was conducted to determine the grain size distribution at the study site. The spatially variable moisture contents at the vertical mesh sand traps, at the saltiphones, and at the meteorological station were measured by a moisture content sensor Theta-Probe ML3 to quantify the volumetric water content of the sand up to 6 cm beneath its surface [64].
In Figure 6, an aerial photograph of the measurement locations with the instrumental setup is shown. For all measurement locations, except for the dune crest, the same measuring equipment was used. On the dune crest, no sediment transport took place during the conducted experiments; thus, only wind velocity and wind direction were measured at this location. The windblown sand grains were recorded with three saltiphone measuring stations and three vertical mesh sand traps. For the beach and the dune toe, the wind tower was placed 1.0 m windward to the left, and the vertical mesh sand traps were placed 0.5 m windward to the right of the saltiphones. the vertical mesh sand traps, at the saltiphones, and at the meteorological station were measured by a moisture content sensor Theta-Probe ML3 to quantify the volumetric water content of the sand up to 6 cm beneath its surface [64]. In Figure 6, an aerial photograph of the measurement locations with the instrumental setup is shown. For all measurement locations, except for the dune crest, the same measuring equipment was used. On the dune crest, no sediment transport took place during the conducted experiments; thus, only wind velocity and wind direction were measured at this location. The windblown sand grains were recorded with three saltiphone measuring stations and three vertical mesh sand traps. For the beach and the dune toe, the wind tower was placed 1.0 m windward to the left, and the vertical mesh sand traps were placed 0.5 m windward to the right of the saltiphones.
The distance between each measuring position was 20 m. The arrangement of the instrumentation for the sand trapping fence differed with regard to Figure 6. Between the brushwood lines of the sand trapping fence the distance between wind tower and saltiphone was reduced to 0.5 m, and the distance between the individual positions was reduced to 10 m. The measuring location sand trapping fence, position 1, was located close to the dune toe in front of the sand trapping fence windward, where the wind was not yet influenced by the brushwood lines. In field 1 (position 2) and field 2 (position 3) the wind regime was influenced by the brushwood lines under the prevailing wind conditions from west during the experiments. Every 20 min, the wind tower was moved to the next measuring position. Meanwhile, the saltiphones and vertical mesh sand traps remained at their initial location. This procedure was repeated between three and five times starting at the first position again. The distance between each measuring position was 20 m. The arrangement of the instrumentation for the sand trapping fence differed with regard to Figure 6. Between the brushwood lines of the sand trapping fence the distance between wind tower and saltiphone was reduced to 0.5 m, and the distance between the individual positions was reduced to 10 m. The measuring location sand trapping fence, position 1, was located close to the dune toe in front of the sand trapping fence windward, where the wind was not yet influenced by the brushwood lines. In field 1 (position 2) and field 2 (position 3) the wind regime was influenced by the brushwood lines under the prevailing wind conditions from west during the experiments. Every 20 min, the wind tower was moved to the next measuring position. Meanwhile, the saltiphones and vertical mesh sand traps remained at their initial location. This procedure was repeated between three and five times starting at the first position again.

Wind Velocity, Shear Velocity, and Aerodynamic Roughness Length
Saturated aeolian sediment transport rates can be described by sediment transport models as a function of the shear velocity [7,22,23,25,28,65]. In most models, it is assumed that the sand is dry, the wind velocity and the sediment transport rates are steady and a logarithmic velocity profile over a flat surface is developed [37,55].
In our data set, the wind velocity profile was measured by six cup anemometers at predefined heights. An appropriate averaging time must be selected to gain reliable values for the shear velocity and aerodynamic roughness length. To account for larger boundary layer turbulences at higher altitudes, averaging over longer time scales is required [66]. However, the potential to choose shorter averaging time scales has been recognized when the highly variable aeolian sediment transport is the focus of investigation. As many researches stated [25,65,67], short averaging time scales (~10 s) allow for resolving short-term variations in surface roughness and aeolian sediment transport due to changing wind conditions. For example, Strypsteen et al. (2019) [22] and Ellis et al. [68] (2009) chose time scales for averaging of less than 1 min. In this study, 10 s as the averaging time scale was chosen.
The obtained wind data were filtered. As a criterion for valid wind measurement data, gusts defined by the German Weather Service were excluded. Therefore, the measured wind velocity must not exceed the measured 10-min average wind velocity by at least 5.0 m/s in a time interval greater than 3 s and smaller than 20 s [53]. To obtain stationary wind profiles, the following conditions were defined: The coefficient of variance of each 10-second-interval at each measurement point must not exceed the selected maximum value of ν ≤ 0.15 and the gradient of acceleration of two successive wind profiles over the selected time interval at each measurement point must not be greater than b ≤ 0.5 m/s 2 .
Under neutral atmospheric stability conditions, the wind velocity distribution over the viscous sub boundary layer for aerodynamically rough surfaces can be described [23] with the law of the wall as follows:  [69,70]. The parameters u* and z 0 were determined from experimental data through a regression analysis of the measured wind velocities at defined heights over the selected time interval of 10 s. A least-square fit of the wind velocities u [m/s] to the natural logarithm of the wind measurement height ln(z) was applied. The shear velocity u* [m/s] results from the slope of the regression line of u [m/s] on ln(z) multiplied by the Kármán constant κ [-]. The aerodynamic roughness length is equal the exponential function of the following term: the quotient of the negative y-axis interception and the slope of the regression line u [m/s] on ln(z) [25,65,71]. Only wind profiles that had a coefficient of determination greater than R 2 > 0.90 were retained for further analysis to ensure only stationary logarithmic wind profiles were analyzed in the following.
With increasing drag and lift force on the particle, there is a critical value at which grain movement is initiated. In nature, this depends on the surface moisture, surface roughness, vegetation, shells, surface slope, and fetch effect [7,35,39,72,73]. The so-called critical shear stress u* t is defined as [18]:  50 [µm] is the mean particle size (here: 210 µm). The empirical constant is assumed as A = 0.11. This value takes into account the effect of cohesion, but not the influence of the protective layer of shells or the increasing cohesion due to higher moisture contents [74][75][76]. With the given parameters, the critical shear velocity at the study site of Langeoog was calculated as u* t = 0.24 m/s.

Vertical Sediment Flux Profile and Potential Aeolian Sediment Transport Rates
The vertical mesh sand traps, see also Figure 5a, were exposed to sediment flow during a defined exposure time, and, afterward, the collected sand was weighed and divided by the inlet area and measurement duration to calculate the mean sediment flux rate. The vertical distribution of sediment transport can be described by an empirical exponential decay function [22,68,77]: where q z [kg/m 2 /s] is the sediment transport rate at a predefined height z [m], q o [kg/m 2 /s] is the extrapolated saltating sediment mass transport at the surface, and the decay rate ß [1/m] as a measure of the vertical concentration gradient. ß and q o can be obtained by a regression analysis to the data of the field experiments [78]. The regression provides generally high correlations of the coefficients, as they are only weakly related to physical aeolian parameters, such as grain size, roughness, and shear velocity. Therefore, a wide range of fitted coefficient values can be found in literature [22,78]. The total mass transport can be calculated by integrating Equation (3) as follows: where Qs [kg/m/s] is the total mass transport by saltation. It is assumed that saltating particles start moving immediately above the surface (z = 0 m) [78].
In Table 1, five representative sand transport models related to the shear velocity u* [m/s], used to calculate the saturated aeolian sand transport rates Qs [kg/m/s] for dry sand, are given. For an overview of all existing models, reference is made to Dong et al. (2003) [79] and Strypsteen (2019) [22]. With the exception of wind velocity, all parameters are considered to be constant over time. C i [-] are constants derived by wind tunnel and/or field experiments from different authors. D (here: 250 µm) [29] is the reference diameter of the sediment as a standard value for dune sand. Table 1. Summary of five chosen representative sand transport models for predicting saturated aeolian sand transport rates for dry sand (modified after [22,79]).

Refer-ence Equation Description/Assumptions Equation No
.
is chosen for naturally graded sand and the equation is only valid for grain sizes between 0.1-1.0 mm, when the threshold of motion is exceeded.
Bagnold corrected his model, introducing the term for threshold shear velocity.
Hsu (1974) Hsu related the predicted sand transport to the third power of the friction Froude number.
Van Rijn (2019) [35] The moisture content on a beach is subject to large spatial and temporal variations [35,36,39]. With increasing moisture content, the cohesion between the sand particles rises; therefore, the critical shear velocity is increased. Even low moisture contents (M~0.5 Vol%) can reduce the transport rates significantly, except when the wind velocity is intense. Saturated moisture content can reduce the instantaneous transport rate by a factor of m moisture = 2 -10 [30,36,38,76]. Therefore, to predict the sediment transport of wet sand the following equation is used: where Q s,wet [kg/m/s] is the predicted transport rate for wet sand, and Q s,dry [kg/m/s] is the saturated aeolian sand transport rates for dry sand.

Results and Discussion
According to the conducted sieve analyses, the median grain size D 50~2 18 µm was chosen as an average value for the measurement locations of the beach, dune toe, and sand trapping fence. The sieve analyses results from the measurement location of the dune crest were excluded for determining the median grain size as no sediment transport was measured on the dune crown. For more details regarding grain size distribution, d 10 and d 90 are also provided representing the grain sizes for 10% and 90% sieve retainment: D 10~1 50 µm, D 90~3 02 µm. The mean size was equal to D m~2 09 µm, and the standard deviation was σ~212 µm.
During the whole field campaign, the sand at the study site was almost continuously slightly wet with varying moisture contents around an average value of M = 6 Vol% up to a maximum moisture content of M = 14.5 Vol% (average values for all locations). The weather was sunny with a maximum temperature of T = 18 • C. The wind was strong over the selected 10-second-interval with an average wind velocity around u 2 = 16.6 m/s at a height of z = 2.0 m up to stormy conditions with a maximum of u 2 = 24.7 m/s. The wind direction came mainly from the west, which lead to coast parallel sediment transport. Aeolian sand transport occurred in the form of streamers. The average shear velocity was almost always above the threshold of movement of u* t = 0.24 m/s based on Equation (2). However, aeolian sediment transport was not measured throughout the whole field campaign.

Wind Velocity Profiles, Shear Velocities, and Aerodynamic Roughness Length
In this Section 4.1, the typical results of the measured velocity profiles and corresponding regression lines to determine the shear velocity and the roughness length for the four different measurement locations, the beach, dune toe, dune crest, and between the sand trapping fences at the study site, are presented. In Figure 7, the wind profiles to determine the shear velocity and roughness length parameters for the locations (a) beach and (b) dune toe are shown. The blue dots show the measured wind velocity at a predefined elevation, the arrows show the wind direction-in this example, the wind comes from the west, and the red line shows the regression line. The regression for all evaluated wind data fits the location beach and dune toe well with an overall coefficient of determination of R² = 0.95 (u* = 0.93 m/s) and R² = 0.88 (u* = 0.43 m/s). The lower coefficient of determination of the regression at the dune toe can be explained by the larger unevenness in the topography compared to the beach. Thereby, it is proven for the locations of the dune toe and beach that the wind profile can be displayed appropriately by the law of the wall, see Equation (1).
In Figure 8a, the wind profile for the location dune crest and, in Figure 8b, the wind profile for the location sand trapping fence in field 1 (position 2) are shown. Regression analyses for position 1 and position 3 were evaluated, but are not shown in the following Figure 8. The regression for all evaluated wind data fits the location beach and dune toe well with an overall coefficient of determination of R 2 = 0.95 (u* = 0.93 m/s) and R 2 = 0.88 (u* = 0.43 m/s). The lower coefficient of determination of the regression at the dune toe can be explained by the larger unevenness in the topography compared to the beach. Thereby, it is proven for the locations of the dune toe and beach that the wind profile can be displayed appropriately by the law of the wall, see Equation (1).
In Figure 8a, the wind profile for the location dune crest and, in Figure 8b, the wind profile for the location sand trapping fence in field 1 (position 2) are shown. Regression analyses for position 1 and position 3 were evaluated, but are not shown in the following Figure 8. The regression analysis of all wind data fit the measuring points for the location dune crest on the foredunes with an overall determination coefficient of R ² = 0.82 (u* = 0.65 m/s). Due to the strong wind, the dune grass bent on the coastal dune to varying degrees depending on the wind velocity, resulting in different roughness lengths. Even with the same vegetation cover and height, the wind flow characteristics can be different [85]. In front of the sand trapping fence at position 1, see also Figure 6, the overall coefficient of determination was R² = 0.89 (u* = 0.92 m/s). The measuring location of the sand trapping fence, position 1, was located close to the dune toe in front of the sand trapping fence windward, so that the wind was not yet influenced by the sand trapping fence.
Between the sand trapping fence at location 2 and location 3, the overall coefficient of determination was significantly lower with R² = 0.34 and R² = 0.38, respectively. On the basis of the curve fitting, between the sand trapping fence, no stationary logarithmic velocity profiles were present. In between the sand trapping fence, the wind directions and the wind velocities varied considerably, see also the wind vectors in Figure 8b. We assumed that turbulent fluctuations took place [16,86], and therefore the effects of short-time variations in the surface roughness associated The regression analysis of all wind data fit the measuring points for the location dune crest on the foredunes with an overall determination coefficient of R 2 = 0.82 (u* = 0.65 m/s). Due to the strong wind, the dune grass bent on the coastal dune to varying degrees depending on the wind velocity, resulting in different roughness lengths. Even with the same vegetation cover and height, the wind flow characteristics can be different [85]. In front of the sand trapping fence at position 1, see also Figure 6, the overall coefficient of determination was R 2 = 0.89 (u* = 0.92 m/s). The measuring location of the sand trapping fence, position 1, was located close to the dune toe in front of the sand trapping fence windward, so that the wind was not yet influenced by the sand trapping fence.
Between the sand trapping fence at location 2 and location 3, the overall coefficient of determination was significantly lower with R 2 = 0.34 and R 2 = 0.38, respectively. On the basis of the curve fitting, between the sand trapping fence, no stationary logarithmic velocity profiles were present. In between the sand trapping fence, the wind directions and the wind velocities varied considerably, see also the wind vectors in Figure 8b. We assumed that turbulent fluctuations took place [16,86], and therefore the effects of short-time variations in the surface roughness associated with large eddies cannot be displayed appropriately by a one-dimensional wind profile according to Equation (1).
In Figure 9, the relationship between the calculated shear velocity u* [m/s] and the aerodynamic roughness length z 0 [mm] at the different measurement locations are shown. The black line indicates the regression for the beach, dune toe, and sand trapping fence (position 1) with the following equation: with large eddies cannot be displayed appropriately by a one-dimensional wind profile according to Equation (1). In Figure 9, the relationship between the calculated shear velocity u* [m/s] and the aerodynamic roughness length z0 [mm] at the different measurement locations are shown. The black line indicates the regression for the beach, dune toe, and sand trapping fence (position 1) with the following equation: The regression parameters are defined as a = 0.49 and b = 4 with a determination coefficient of R² = 0.75. The dune crest data is not implemented in the regression line. The values for the aerodynamic roughness length close to the dune toe lie between z0 = 2.6 · 10 −6 -0.23 mm, on the beach z0 = 0.009-2.03 mm, and in front of the sand trapping fence between z0 = 0.003-1.79 mm. The results of the dune crest cannot be compared directly to the results of the dune toe, beach, and sand trapping fence (position 1) due to the presence of vegetation and higher bed roughness on the dune crest (see also Chapter 2). At the dune crest, the values were significantly higher (~100 times) and were subject to more scatter with aerodynamic roughness lengths of around z0 = 0.003-90.96 mm. However, they follow the same trend as the measurement locations of the dune toe, beach, and sand trapping fence (position 1). The values are in a similar range as other field experiments for estimating the aerodynamic roughness length over a coastal dune field [85].
The aerodynamic roughness length increased with the increasing shear velocity. The aerodynamic roughness length increased as a power law with the increasing shear velocity [22,25].

Vertical Sediment Flux Profile
For the measurement locations of the beach, dune toe, and sand trapping fence (position 1 and position 2), the vertical mesh sand traps over z = 0.3 m above the surface were analyzed in the following. In Figure 10, a representative example of the vertical distribution of sediment transport fluxes derived from the vertical mesh sand traps is shown for the measurement location of the beach. The trapped sediment over the height is described using an exponential decay curve as described in Equation (3). The regression parameters are defined as a = 0.49 and b = 4 with a determination coefficient of R 2 = 0.75. The dune crest data is not implemented in the regression line.
The values for the aerodynamic roughness length close to the dune toe lie between z 0 = 2.6 · 10 −6 -0.23 mm, on the beach z 0 = 0.009-2.03 mm, and in front of the sand trapping fence between z 0 = 0.003-1.79 mm. The results of the dune crest cannot be compared directly to the results of the dune toe, beach, and sand trapping fence (position 1) due to the presence of vegetation and higher bed roughness on the dune crest (see also Section 2). At the dune crest, the values were significantly higher (~100 times) and were subject to more scatter with aerodynamic roughness lengths of around z 0 = 0.003-90.96 mm. However, they follow the same trend as the measurement locations of the dune toe, beach, and sand trapping fence (position 1). The values are in a similar range as other field experiments for estimating the aerodynamic roughness length over a coastal dune field [85].
The aerodynamic roughness length increased with the increasing shear velocity. The aerodynamic roughness length increased as a power law with the increasing shear velocity [22,25].

Vertical Sediment Flux Profile
For the measurement locations of the beach, dune toe, and sand trapping fence (position 1 and position 2), the vertical mesh sand traps over z = 0.3 m above the surface were analyzed in the following. In Figure 10, a representative example of the vertical distribution of sediment transport fluxes derived from the vertical mesh sand traps is shown for the measurement location of the beach. The trapped sediment over the height is described using an exponential decay curve as described in Equation (3). The sand was caught predominantly in heights of less than z ≤ 0.05 m. When the wind was really strong with u*~1 m/s, sand was caught in heights over z > 0.15 m. This is in accordance with the literature reporting that most of the sediment transport in coastal areas occurs in a thin layer z < 0.05 m above the surface [22].
In Tables 2 and 3   The sand was caught predominantly in heights of less than z ≤ 0.05 m. When the wind was really strong with u*~1 m/s, sand was caught in heights over z > 0.15 m. This is in accordance with the literature reporting that most of the sediment transport in coastal areas occurs in a thin layer z < 0.05 m above the surface [22].
In Tables 2 and 3 Table 2. Results of the total mass flux over the measuring duration, decay rate, extrapolated saltating sediment mass transport at the surface, correlation coefficients of the regression analysis of the vertical wind profile (rounded values), shear velocity, moisture content, and aerodynamic roughness length.

Location
No.  too little sediment transport to conduct regression analysis, no vertical wind profile according to the law of the wall present * no wind data available, u * and z 0 assumed from data of the closest time interval; ** no wind data available between the sand trapping fence.
The decay rate generally increases for a constant saltating sediment mass transport at the surface with higher wind velocity or less moisture content, indicating that more sand is transported at higher elevations above the surface [56]. Considerable variation was present when comparing the decay rate and saltating sediment mass transport at the surface for different measurement locations, where the prevailing boundary conditions, like the moisture content and shear velocities, varied. Close to the dune toe, the wind regime was of medium intensity, and the sand was relatively dry, resulting in lower saltation heights, compared to the beach and sand trapping fence, when it was relatively wet; however, high wind velocities occurred. The presence of slightly wet sand increased the threshold of motion for sand, resulting in a lower total sediment flux.
Aeolian sediment transport is highly variable in both space and time. During one event, the trapped sand in two vertical mesh sand traps standing next to each other can vary a lot. This happens, if wind forces and sand transport occur in the form of streamers resulting in spatial-temporally changing threshold shear velocity [22,26,35,78].
At the measurement location of the sand trapping fence, position 3, there was too little sediment caught in the vertical mesh sand trap over the defined time interval. Thus, the vertical sediment flux profile cannot be described with the approach according to Equation (3).
The literature gives a wide range of fitting coefficients and different methods for trapping sand above the surface [22,68,78,79]. The studies at the Belgian coast at Mariakerke and Koksijde are described in the following to compare the results of this work with other authors. Strypsteen (2019) [22] measured the trapped sand with Modified Wilson And Cool Samplers (MWAC) on a dry beach, where the bottles were placed at the heights of z = 6.5, 13.5, 21, and 28.5 cm above the surface. The results were analyzed with the potential decay function as an average value out of 12 measurement locations on the beach. With an average shear velocity between u* = 0.29 -0.54 m/s and transport fluxes between Q s = 7.2 -121.8 kg/m/h, the decay rate varied between ß = 14.06 -28.57/m. The measuring time varied between minutes up to days depending on the intensity of the sediment transport [87]. Campos (2018) [56] measured the sediment transport fluxes with vertical mesh sand traps on the Belgian coast at Koksijde. With an average shear velocity between u* = 0.40-0.53 m/s and transport fluxes between Q s = 38.4-88 kg/m/h, the decay rate varied between ß = 18.13-23.72/m. In the study, each run was performed nine times with measurement durations between t = 16 min up to t = 38 min. Little to no moisture was present [56].
These two selected field campaigns measured significantly higher total mass flux but lower shear velocities and decay rates. This could be explained by the presence of slightly wet sand on the study site Langeoog, increasing the threshold of motion for the sand, resulting in a lower total sediment flux (see also Section 2) for higher shear velocities. Further reasons for these deviations may be that Campos (2018) [56] used shorter time intervals for trapping the sand in the vertical mesh traps. The trap efficiency of the measurement tool MWAC, stated as the ratio of total particle mass passing through the air per defined area to the actual value, varies between 42%-120% [59,88,89]. For the vertical mesh sand traps, there was good efficiency reported; however, the efficiency decreased for longer time periods [55,56]. The MWAC sand trap began measuring at z ≥ 6.5 cm above the surface, whereas the vertical sand trap began measuring right above the sand surface.

Observed Aeolian Sediment Transport Rates at the Measurement Location of the Beach
In Figure 11, the theoretical comparison between the field data of Langeoog beach, dune toe, sand trapping fence (position 1), and five chosen transport models derived by Bagnold (1937) [80], Kawamura (1951) [82], Bagnold (1954) [23], Kok et al. (2012) [31], and van Rijn (2019) [35] (see also Table 1) are shown. The Bagnold (1954) [23] model results in the smallest transport rates compared to the other prediction models, especially for smaller shear velocities. The orange dots show the results of the location beach, the blue dots of the location dune toe, and the brown dots of the location sand trapping fence (location 1). Campos (2018) [56] measured the sediment transport fluxes with vertical mesh sand traps on the Belgian coast at Koksijde. With an average shear velocity between u* = 0.40-0.53 m/s and transport fluxes between Qs = 38.4-88 kg/m/h, the decay rate varied between ß = 18.13-23.72/m. In the study, each run was performed nine times with measurement durations between t = 16 min up to t = 38 min. Little to no moisture was present [56].
These two selected field campaigns measured significantly higher total mass flux but lower shear velocities and decay rates. This could be explained by the presence of slightly wet sand on the study site Langeoog, increasing the threshold of motion for the sand, resulting in a lower total sediment flux (see also Chapter 2) for higher shear velocities. Further reasons for these deviations may be that Campos (2018) [56] used shorter time intervals for trapping the sand in the vertical mesh traps. The trap efficiency of the measurement tool MWAC, stated as the ratio of total particle mass passing through the air per defined area to the actual value, varies between 42%-120% [59,88,89]. For the vertical mesh sand traps, there was good efficiency reported; however, the efficiency decreased for longer time periods [55,56]. The MWAC sand trap began measuring at z ≥ 6.5 cm above the surface, whereas the vertical sand trap began measuring right above the sand surface.

Observed Aeolian Sediment Transport Rates at the Measurement Location of the Beach
In Figure 11, the theoretical comparison between the field data of Langeoog beach, dune toe, sand trapping fence (position 1), and five chosen transport models derived by Bagnold (1937) [80], Kawamura (1951) [82], Bagnold (1954) [23], Kok et al. (2012) [31], and van Rijn (2019) [35] (see also Table 1) are shown. The Bagnold (1954) [23] model results in the smallest transport rates compared to the other prediction models, especially for smaller shear velocities. The orange dots show the results of the location beach, the blue dots of the location dune toe, and the brown dots of the location sand trapping fence (location 1). Figure 11. Comparison between the measured sediment transport rates on Langeoog and the predicted sediment transport rates for five different transport models.
The measured data from Langeoog beach can be described well with the predicted transport models if the moisture content as a reduction factor is considered. For the location dune toe and sand trapping fences (position 1) the results scattered more but still showed good results. We could not clearly evaluate which equation fit the results best as only high wind velocities with shear velocities around u* = 0.70 -1.05 m/s were recorded, except for u* = 0.42 m/s, where the results scattered a lot. For further analysis, smaller shear velocities are required to indicate the best fit regression line.
An application of the transport models from Table 1 is only useful for the entrance boundary condition (upwind the sand trapping fence, position 1). Between the sand trapping fences, the Figure 11. Comparison between the measured sediment transport rates on Langeoog and the predicted sediment transport rates for five different transport models.
The measured data from Langeoog beach can be described well with the predicted transport models if the moisture content as a reduction factor is considered. For the location dune toe and sand trapping fences (position 1) the results scattered more but still showed good results. We could not clearly evaluate which equation fit the results best as only high wind velocities with shear velocities around u* = 0.70 -1.05 m/s were recorded, except for u* = 0.42 m/s, where the results scattered a lot. For further analysis, smaller shear velocities are required to indicate the best fit regression line.
An application of the transport models from Table 1 is only useful for the entrance boundary condition (upwind the sand trapping fence, position 1). Between the sand trapping fences, the transport models cannot be used as the sediment supply from above was too small and a logarithmic wind profile no longer existed (see Section 4.1).

Downwind Evolution of Saltiphone Data for the Measurement Location Beach
In Figure 12, an example of the wind and saltiphone data on the beach on 23 May 2020 for the time interval 12:35 pm to 13:10 pm is shown. This is a representative example out of four repeating measurements on the beach. The calculated shear velocities from the wind profiles and the mean wind direction are shown as an average value over 10 s. The critical shear velocity is also given, see Figure 12a. Additionally, the 95% quantiles of the saltation impacts over 10-s time intervals are shown. In Figure 12a, the black line represents the calculated shear velocity, the dotted red line indicates the critical shear velocity (u* t = 0.24 m/s, see Section 3.2.1), whereas the light orange line shows the mean wind direction. Figure 12b  transport models cannot be used as the sediment supply from above was too small and a logarithmic wind profile no longer existed (see Chapter 4.1).

Downwind Evolution of Saltiphone Data for the Measurement Location Beach
In Figure 12, an example of the wind and saltiphone data on the beach on 23 May 2020 for the time interval 12:35 pm to 13:10 pm is shown. This is a representative example out of four repeating measurements on the beach. The calculated shear velocities from the wind profiles and the mean wind direction are shown as an average value over 10 s. The critical shear velocity is also given, see Figure 12a. Additionally, the 95% quantiles of the saltation impacts over 10-s time intervals are shown. In Figure 12a   The shear velocity was always above the critical shear velocity during the whole experiment, showing continuous sediment transport for all three positions. In general, as the number of counted sediment grains on the lower saltiphone increased, so did the number of counted sediment grains on the upper saltiphone. This suggests that the height of the saltating layer increased. At position 2 (20 m downwind), lower impacts were measured compared to position 1, followed by an increase of the saltation impacts at position 3 (40 m downwind). The smaller number of impacts at the second location can likely be attributed to a higher moisture content or more shells, protecting the beach against the erosion of sand grains, resulting in higher critical shear velocities [35,90,91].
We assumed that, between location 1 and location 2, there was a sediment sink due to the high moisture contents, whereas, between position 2 and position 3, there was a sediment source with lower moisture contents. During the whole day, the moisture content was relatively high with values around M = 3.3-14.5 Vol%, showing almost saturated values. However, the moisture content was very variable in time and space. The higher moisture content resulted in higher critical shear velocities compared to the dune toe, where the sand was relatively dry.

Downwind Evolution of Saltiphone Data for the Measurement Location of the Dune Toe
In Figure 13, a representative example of the wind and saltiphone data close to the dune toe on 19 May 2020 for the time interval 12:15 pm to 13:05 pm is shown. This is an example out of three repeated measurements of the dune toe. The shear velocity was always above the critical shear velocity during the whole experiment, showing continuous sediment transport for all three positions. In general, as the number of counted sediment grains on the lower saltiphone increased, so did the number of counted sediment grains on the upper saltiphone. This suggests that the height of the saltating layer increased. At position 2 (20 m downwind), lower impacts were measured compared to position 1, followed by an increase of the saltation impacts at position 3 (40 m downwind). The smaller number of impacts at the second location can likely be attributed to a higher moisture content or more shells, protecting the beach against the erosion of sand grains, resulting in higher critical shear velocities [35,90,91].
We assumed that, between location 1 and location 2, there was a sediment sink due to the high moisture contents, whereas, between position 2 and position 3, there was a sediment source with lower moisture contents. During the whole day, the moisture content was relatively high with values around M = 3.3-14.5 Vol%, showing almost saturated values. However, the moisture content was very variable in time and space. The higher moisture content resulted in higher critical shear velocities compared to the dune toe, where the sand was relatively dry.

Downwind Evolution of Saltiphone Data for the Measurement Location of the Dune Toe
In Figure 13, a representative example of the wind and saltiphone data close to the dune toe on 19 May 2020 for the time interval 12:15 pm to 13:05 pm is shown. This is an example out of three repeated measurements of the dune toe.  The average wind velocity was u 2 = 11.6 m/s, and the average shear velocity was u* = 0.36 m/s. The wind was relatively steady coming from west (dr = 277 • ), parallel to the dune. During the time interval shown in Figure 13, the shear velocity was mostly slightly above the critical shear velocity, whereby, around 12:56 pm and 13:02 pm, the critical shear velocity was undercut, and no sediment transport was further measured. The results, therefore, indicate that the calculated critical shear velocity was a good approximation of the threshold for the initiation of sediment transport.
The lower saltiphones recorded more impacts than the upper saltiphones. As only sporadic saltation impacts were measured by the upper saltiphones, they are hard to distinguish in Figure 13; see, for example, Figure 13b around 12:17 pm. The number of impacts measured by the upper saltiphones increased due to a wind peak. This shows that most of the sediment transport occurred below z ≤ 0.15 m. This also agrees with the results of the vertical mesh sand traps close to the dune toe, where at heights of z ≥ 0.15 m above the surface little sediment were caught (see also Section 4.2).
Time intervals with increased saltation impacts corresponded temporally for all measuring positions downwind. Most grain impacts were recorded at position 1, followed by position 2 (20 m downwind) and position 3 (40 m) approximately in the same order of magnitude. The peaks in the wind velocity correspond with the peaks in the counted particles by the saltiphones.
Little spatial variations in sediment transport gradients are found downwind and reflect the spatially varying availability of sand related to the surface roughness and moisture. The moisture content of the sand close to the dune was relatively low with average values around M~1.1-1.8 Vol% and, therefore, played a minor role for supply limiting factors. However, due to the medium intense wind velocities, the bed ripples moved downwind, causing a high spatial variability of the aerodynamic roughness length (see also Figure 4). Generally, these ripples occur with wind velocities of u = 5 -10 m/s, and, with higher wind velocities, these ripples are gradually smoothed out [67].

Downwind Evolution of Saltiphone Data Influenced by Sand Trapping Fences
In Figure 14, a representative example of the time series on 24 May 2020 from 09:45 am to 10:08 am for the measurement location sand trapping fence is shown. This is a representative example out of thirteen repeating measurements. Position 1 is in front of the sand trapping fence and is treated as a boundary condition for incoming sediment that accumulates at the brushwood lines of the sand trapping fence. Position 2 (field 1) and position 3 (field 2) are in between the sand trapping fences at 10 m and 20 m, respectively, downwind (see Figure 6). All three locations belong to the sand trapping configuration 1, see Figure 2a. The average wind velocity at 2 m of height was u 2 = 20.8 m/s with a few fluctuations. The shear velocity was u* = 0.94 m/s and always above the critical shear velocity. The wind direction was relatively steady coming from the west (dr = 262 • ), which implies parallel wind to the dune. The number of counted particles decreased in the downwind direction behind the first line of brushwood of the sand trapping fence. In the first field of the sand trapping fence, less sediment was transported than in front of the sand trapping fence. In the second field, even less sediment was transported. With the measured wind conditions and the given time interval, the sand transport between the sand trapping fence could only be measured up to the second sand trap field. In field 2, the counted saltiphones were almost not measurable or zero.
In Figure 15, 95% quantiles of the counted sand grains s95 [counts/10 s] at position i relative to the boundary condition s95,position 1 [counts/10 s] in percentage over the distance downwind y [m] through the sand trapping fence for different shear velocities are shown. The 95%-quantiles of the saltation impacts were averaged over a 60-min time interval with approximately the same shear velocity. Position 1 (upwind of the sand trapping fence) was taken as the boundary condition with an input of 100% saltation impacts. In field 1 and field 2, the saltiphone data are shown relative to the boundary condition (position 1). To illustrate the decrease in the saltation impacts before and between the sand trapping fence in a simplified way, the saltation impacts are shown constant over the distance downwind y1 = 0-10 m, y2 = 10-20 m, and y3 =20-30 m-although in reality, this is not the case. The number of counted particles decreased in the downwind direction behind the first line of brushwood of the sand trapping fence. In the first field of the sand trapping fence, less sediment was transported than in front of the sand trapping fence. In the second field, even less sediment was transported. With the measured wind conditions and the given time interval, the sand transport between the sand trapping fence could only be measured up to the second sand trap field. In field 2, the counted saltiphones were almost not measurable or zero.
In Figure 15, 95% quantiles of the counted sand grains s 95 [counts/10 s] at position i relative to the boundary condition s 95,position 1 [counts/10 s] in percentage over the distance downwind y [m] through the sand trapping fence for different shear velocities are shown. The 95%-quantiles of the saltation impacts were averaged over a 60-min time interval with approximately the same shear velocity. Position 1 (upwind of the sand trapping fence) was taken as the boundary condition with an input of 100% saltation impacts. In field 1 and field 2, the saltiphone data are shown relative to the boundary condition (position 1). To illustrate the decrease in the saltation impacts before and between the sand trapping fence in a simplified way, the saltation impacts are shown constant over the distance downwind y 1 = 0-10 m, y 2 = 10-20 m, and y 3 =20-30 m-although in reality, this is not the case. For all shear velocities between u* = 0.68 -1.01 m/s, the counted saltation impacts were significantly reduced between the sand trapping fences, indicating the deposition of sand in the distance downwind. In field 1 (position 2), the sediment transport for the lower saltiphone was reduced between 5-95% compared to position 1, and, in field 2 (position 3), the sediment transport was reduced between 96-98% compared to position 1. For the upper saltiphone, the sediment transport in field 1 (position 2) was reduced between 13-44% compared to position 1, and, in field 2 (position 3), the sediment transport was reduced between 72-91% compared to position 1. At high shear velocities of up to u* = 1.01 m/s, the sand transport in field 1 was only slightly reduced for both saltiphone heights. For the upper saltiphone, the saltation impact was less reduced compared to the lower saltiphone. Due to the installation of the bundled brushwood, the sand trapping fence was relatively dense directly above the bed surface. This indicates that more sediment grains hit the brushwood and were deposited close to the bed surface.
The investigation at the sand trapping fence were conducted directly after the installation of the sand trapping fence, and, therefore, the fields caught relatively little sand. Sediment was caught predominantly by the first two brushwood lines, where sediment transport could also be measured. This implies that the first two brushwood lines of this type of sand trapping fence configuration (see also Figure 2), with the westerly longshore wind and the actual amount of trapped sand in between the sand trapping fence, were the most effective brushwood lines to trap sand. However, on Langeoog, many of the elderly installed sand trapping fences are already completely filled with sand, showing that not only the first two brushwood lines exposed to longshore wind but also the brushwood lines that are not in an exposed wind position trap sand effectively. This can be explained by the fact that strong winds from the west over north to south can also transport sediment to the sand trapping fence, where it can accumulate at the brushwood lines orthogonal and parallel to the coastal dune. In addition, with a larger amount of trapped sediment between the individual brushwood lines and longshore wind, other fields (not directly exposed to the longshore wind) can be filled with sediment as well. For all shear velocities between u* = 0.68 -1.01 m/s, the counted saltation impacts were significantly reduced between the sand trapping fences, indicating the deposition of sand in the distance downwind. In field 1 (position 2), the sediment transport for the lower saltiphone was reduced between 5-95% compared to position 1, and, in field 2 (position 3), the sediment transport was reduced between 96-98% compared to position 1. For the upper saltiphone, the sediment transport in field 1 (position 2) was reduced between 13-44% compared to position 1, and, in field 2 (position 3), the sediment transport was reduced between 72-91% compared to position 1. At high shear velocities of up to u* = 1.01 m/s, the sand transport in field 1 was only slightly reduced for both saltiphone heights. For the upper saltiphone, the saltation impact was less reduced compared to the lower saltiphone. Due to the installation of the bundled brushwood, the sand trapping fence was relatively dense directly above the bed surface. This indicates that more sediment grains hit the brushwood and were deposited close to the bed surface.
The investigation at the sand trapping fence were conducted directly after the installation of the sand trapping fence, and, therefore, the fields caught relatively little sand. Sediment was caught predominantly by the first two brushwood lines, where sediment transport could also be measured. This implies that the first two brushwood lines of this type of sand trapping fence configuration (see also Figure 2), with the westerly longshore wind and the actual amount of trapped sand in between the sand trapping fence, were the most effective brushwood lines to trap sand. However, on Langeoog, many of the elderly installed sand trapping fences are already completely filled with sand, showing that not only the first two brushwood lines exposed to longshore wind but also the brushwood lines that are not in an exposed wind position trap sand effectively. This can be explained by the fact that strong winds from the west over north to south can also transport sediment to the sand trapping fence, where it can accumulate at the brushwood lines orthogonal and parallel to the coastal dune. In addition, with a larger amount of trapped sediment between the individual brushwood lines and longshore wind, other fields (not directly exposed to the longshore wind) can be filled with sediment as well.

Conclusions
This work presents the results of field experiments conducted on the coast of the East Frisian island Langeoog, analyzing wind and sediment transport data at the beach, coastal dune, and sand trapping fence. The following conclusions can be drawn from the results presented in this work: (1) The logarithmic wind profiles for the beach, the area close to the dune toe, the dune crest of the foredunes, and upwind of the sand trapping fence (position 1) can be well described by the law of the wall. In between the sand trapping fence, no stationary logarithmic velocity profile was present. Here, the wind directions and the wind velocities varied considerably over the distance from the bottom for any given point in time. We assumed that turbulent fluctuations took place, and therefore the effects of short-time variations in surface roughness associated with large eddies cannot be displayed appropriately by a one-dimensional wind profile. The values for the aerodynamic roughness length close to the dune toe, on the beach, and in front of the sand trapping fence can be described as a function of the shear velocity. At the dune crest, where vegetation was present, the values were~100 times higher. With increasing aerodynamic roughness length, the shear velocity also increased. (2) The vertical sediment flux profile of aeolian sand transport for the measurement locations of the beach, dune toe, and sand trapping fence (position 1 and 2) can be described well with an empirical exponential decay function. The regression parameters exponential decay rate and the extrapolated saltation flux at the surface were subject to considerable variation for different measurement locations. Close to the dune toe, the wind regime was of medium intensity, and the sand was relatively dry. On the beach and between the sand trapping fence, high wind velocities occurred, and the sand was relatively wet, resulting in higher critical shear velocities and lower total sediment fluxes. (3) The sediment caught in the vertical mesh sand traps divided through the inlet area and the measurement duration could well describe the sediment transport rates on the coastal area. The results are described well by the transport models derived by Bagnold (1937) [80], Kawamura (1951) [82], Bagnold (1954) [23], Kok et al. (2012) [31], and van Rijn (2019) [35], considering the moisture content as a further reduction factor. Between the sand trapping fences, these transport models cannot be applied as no logarithmic wind profile existed and the upwind sediment supply was reduced after each brushwood line, thereby, leading to increased deviation from the saturated conditions. (4) The saltiphones are a good measurement tool to determine the moment sediment transport is initiated and, thus, for determining the value of the critical shear velocity. Small spatial variations in the sediment transport gradients were found downwind close to the dune toe. They reflected the spatially varying availability of erodible sand, which is related to the surface roughness and moisture. However, due to the medium intense wind velocities, the bed ripples were moving downwind, causing a high spatial variability of aerodynamic roughness length. At the measurement location of the beach (position 2, 20 m downwind), a smaller number of impacts were measured compared to position 1, followed by an increase of the saltation impacts at position 3 (40 m downwind). We assumed that, between location 1 and location 2, a sediment sink was present due to the high moisture contents, whereas, between position 2 and 3, we assumed a sediment source was present due to the lower moisture contents. During the whole field campaign, the moisture content was high resulting in higher critical shear velocities. (5) Sediment transport was significantly reduced by the sand trapping fence. The sand transport was highest upwind of the sand trapping fence. Downwind, in field 1 and field 2, the sediment transport for the lower and upper saltiphones was reduced. In field 2, the sediment transport was almost not measurable or zero.
This study provided the first insights into aeolian sediment transport at different locations on the coastal dune and the influence of sand trapping fences on aeolian sediment transport. We strongly recommend to extend the knowledge of the influences of sand trapping fences on aeolian sediment transport in standardized wind tunnel experiments to gain more quantitative data. To evaluate the influence of sand trapping fences on dune development for longer timescales, aerial drone investigations of the coastal dune should be conducted on the study site Langeoog. For this purpose, an analysis of a series of digital elevation models at the coastal dune and at the sand trapping fence will be performed to investigate the dune volume changes around the fence depending on each sand trapping configuration. The aim of the future investigation is to evaluate the efficiency of each sand trapping configuration under different prevailing wind conditions.