New Method for Estimation of Aeolian Sand Transport Rate Using Ceramic Sand Flux Sensor (UD-101)

In this study, a new method for the estimation of aeolian sand transport rate was developed; the method employs a ceramic sand flux sensor (UD-101). UD-101 detects wind-blown sand impacting on its surface. The method was devised by considering the results of wind tunnel experiments that were performed using a vertical sediment trap and the UD-101. Field measurements to evaluate the estimation accuracy during the prevalence of unsteady winds were performed on a flat backshore. The results showed that aeolian sand transport rates estimated using the developed method were of the same order as those estimated using the existing method for high transport rates, i.e., for transport rates greater than 0.01 kg m−1 s−1.

Aeolian sand transport mechanisms observed in the field under unsteady wind conditions differ significantly from those observed under steady wind conditions. Since the 1990s, high-frequency sampling systems such as piezoelectric mass flux sensors [17][18][19][20][21][22], high-temporal-resolution sediment traps [23][24][25][26][27], and microphonic flux sensors [28] have been employed for the detection of aeolian sand OPEN ACCESS flux at frequencies of several hundreds of hertz and higher. Wind tunnel and field experiments performed using these systems have improved our understanding of sand transport mechanisms under unsteady wind conditions. However, further investigations are required, especially under complex field conditions.
The ceramic sand flux sensor UD-101 (Chuo Kosoku), which is a piezoelectric sand flux sensor that can measure the flux of aeolian sand impacting on its surface, was developed by Kubota et al. [29,30] (see Figure 5c). The sensor detects piezoelectric signals generated by collisions between aeolian sand grains and the sensor surface using the same principle as that employed in Sensit [17] and Safires [19]. A significant difference between UD-101 and the other piezoelectric sensors is that UD-101 is unidirectional and does not self-orient toward the wind direction, while both Sensit and Safires are omnidirectional. However, UD-101 has the following advantages over the other piezoelectric sensors: (i) The sensor is calibrated by performing both wind tunnel and field experiments.
(ii) The sensor records an impact value of 0 or 1 at 10 kHz; it samples once per second, and theoretically, it can detect a maximum of 10,000 impact counts in the sampling interval, i.e., the sensor has high temporal resolution. (iii) The sensor is relatively inexpensive. Udo et al. [22] carried out short-term field observations at the Hasaki Beach in Japan using UD-101 during the period 12-16 January 2005. The following features were inferred from the UD-101 data: (i) During periods without rainfall and in the presence of longshore winds (conditions similar to those employed in wind tunnel experiments), the aeolian sand flux at a certain height above a flat ground surface increased significantly with the wind velocity and approximately equaled the flux estimated in wind tunnel experiments using Kawamura's [3] equation. (ii) The flux decreased significantly with an increase in precipitation, i.e., with an increase in the moisture content of the sand surface; however, even during periods with rainfall, flux was detected during strong wind conditions. (iii) The flux increased with a decrease in the angle between the sensor direction and wind direction. Further, Udo [31] carried out relatively long-term observations during the period from April 2005 to January 2006 at the Hasaki Beach and showed that a large amount of sand was transported during typhoons. This result implies that UD-101 can be employed for performing long-term flux measurements even during storm events. As the next step toward making long-term flux measurements in fields possible, it is necessary to develop a measurement method for not only the aeolian sand flux at a certain height above the ground surface but also the aeolian sand transport rate (e.g., the total aeolian sand flux between the ground surface and the maximum height of the aeolian sand layer). In particular, the accurate estimation of the rate of seasonal, yearly, or decadal sediment transport due to both littoral drift and aeolian transport at beaches would be invaluable for coastal management. However, the aeolian transport rates estimated by most researchers have not been validated because of the lack of practical estimation/measurement methods; examples of such estimations are estimations made by using conventional equations along with the mean wind data [32] or by integrating areas of beach morphological changes from a boundary under the assumption that the sand transport rate at the boundary is zero [33].
The objective of this study was to develop a new method employing UD-101 for the estimation of the aeolian sand transport rate [29,30]. We applied the method to field data obtained during the period from 22 September to 6 December 2005 [31] and assessed the estimation accuracy to examine the suitability of the method for the direct estimation of long-term transport rates in the field. Figure 1 shows a time series of the instantaneous horizontal wind velocity u and blown-sand impact count n measured at the Hasaki Beach; the median sand grain size d 50 at the beach was 0.2 × 10 -3 m [22,31]. n shows a positive relationship with u during periods without any rainfall, but decreases significantly after a rainfall event. The sensor cannot detect the impact of blown sand if wet sand adheres to its surface after a rainfall event; however, it can detect the impact if this sand is subsequently removed by heavy rainfall. The sensor does not detect the impact of raindrops.  [22,31]; the data were recorded for 5 min at a frequency of 1 Hz (a) during intervals without rainfall on 12 January 2005 and (b) after a rainfall event on 25 September 2005. The wind velocity was measured using an ultrasonic anemometer at a height of 0.9 m above the flat ground surface (z u = 0.9 m). The impact count was measured using UD-101 at z s = 0.04 m. The horizontal distance between the wind and saltation sensors was 1 m. The median sand grain size d 50 of the ground surface at the beach was 0.2 × 10 -3 m. Many studies have shown that the vertical wind velocity profile in a wind tunnel is logarithmic [34,35]. In previous studies, the threshold wind velocity u t has been estimated from both a wind profile law and Bagnold's [1] equation. The logarithmic wind profile law is expressed as:

Characteristics of UD-101 Data
Here, u * is the wind shear velocity, κ is the von Kármán constant (= 0.4), z u is the height at which u is measured, and z 0s is the saltation roughness length. Saltation is the hopping motion of sand grains and the primary mode of aeolian sand transport in beach intertidal and supratidal zones. Aeolian sand transport leads to an increase in the aerodynamic roughness length. u t is obtained by substituting u * = u *t in Equation (1). The threshold wind shear velocity u *t is estimated from where a is the numerical transport rate coefficient (= 0.1), ρ is the air density, σ is the sand density (= 2,650 kg m -3 ), g is the gravitational acceleration, and d is the sand grain size. z 0s is estimated from the relation proposed by Charnock [36] and Chamberlain [37]: where c 0 is a constant (= 0.16 for saltation on beaches [38]).

Relationship between Impact Counts and Aeolian Sand Flux at a Measurement Height
Kubota et al. [30] performed experiments in a 20-m-long wind tunnel using the experimental setup that was employed by Hotta et al. [39]; the sand used in the experiments had a d 50 value of 0.25 × 10 -3 m. The authors showed that the impact counts approached a constant value asymptotically at frequencies greater than 8 kHz. They also compared the sand mass flux q s measured in the wind tunnel using an 8-kHz sensor [30] with q measured using a vertical-distribution-type sand trap [39]. q s at the measurement height z s was estimated from n using the following equation [22]: where t 0 is the sampling period of n (= 1 s) and d s is the diameter of the sensor (= 0.012 m). The vertical-distribution-type sand trap was developed by Hotta and Horikawa [40], and it can be used to measure the flux at 37 different heights in the range 0-0.505 m. q was corrected for trap efficiency. The wind tunnel experiments performed by Kubota et al. [30] and Hotta et al. [39] showed that q s (z s ) for d 50 = 0.25 × 10 -3 m was related to q(z s ) measured using the vertical-distribution-type sand trap [22], which was closely approximated as:

Relationship between Aeolian Sand Flux at a Measurement Height and Total Sand Flux
The relationship between the sand flux at z s , q(z s ) [kg m -2 s -1 ], and the total sand flux Q [kg m -1 s -1 ] is shown in Figure 2. Hotta et al. [39] employed the same wind tunnel as that used by Kubota et al. [30], with 0.17 m s -1 < u * < 1.77 m s -1 and 0.14 × 10 -3 m < d 50 < 0.68 × 10 -3 m. Here, only the data for d 50 = 0.25 × 10 -3 m are shown. Q was found to be approximately proportional to q(z s ), as expressed by: where the regression coefficient R 2 of the approximation was 0.96, 0.99, 1.00, 1.00, 0.99, 0.98, 0.97, and 0.97 for z s = 0.005, 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, and 0.14 m, respectively. As shown in Figure 3, the constant A shows a positive relationship with z s , and it can be expressed as: The relationships between z s and R 2 and between z s and A (shown in Figures 2 and 3, respectively) indicate that z s should be set in the range 0.02-0.08.

Summary of Field Measurements
Field measurements were performed at the Hasaki Beach using UD-101. The beach faces the Pacific Ocean and is located in Ibaraki Prefecture [35°50′25″N, 140°45′42″E; Figure 5(a)] [22,31]. The beach is of the dissipative type and has a foreshore slope ranging from 0.02 to 0.05 and a width of 100 m [ Figure 5 deployed, and they were directed toward the north, east, south, and west. The wind direction in the horizontal plane is expressed clockwise from north, from 0° to 360°. The data of u, θ, and n were recorded on a data logger system (FieldPoint cFP-2000, National Instruments) at a frequency of 1 Hz and were automatically downloaded on a weekly basis from the logger. The precipitation Pr was measured at 10-min intervals at the Choshi Meteorological Observatory (35°44′18″N, 140°51′24″E) by the Japan Meteorological Agency. In addition, in this study, z s was measured manually almost every day in order to calculate Q using Equation (10). The instantaneous z s was calculated by linear interpolation. The z s value of the sensor oriented toward the wind direction tended to be large and that of the sensor directed in the opposite direction tended to be small; for example, for north winds, the z s value of the sensor directed northward tended to be large and that of the sensor directed southward tended to be small.
The measurement results are shown in Figure 6. Data could not be obtained from 28 October to 13 November 2005 because of network problems. Large amounts of sand were transported during the transport events on 24-26 September due to the impact of typhoon 0517, which resulted in a total precipitation of 39.5 mm; similar events involving the transport of a large amount of sand occurred on 11-12 October (period without rainfall); on 18-20 October (due to the impact of typhoon 0520), when the total precipitation was 6 mm; and on 14-15 November (period without rainfall) [31].

Evaluation of the Accuracy of the New Method by Comparing Q with Q B , Q K , Q O , and Q L
The Q value estimated from n was compared with that obtained using the conventional equations for a constant wind velocity in order to examine the accuracy of Q estimated from the UD-101 data. Here, the commonly employed equations of Bagnold [1], Kawamura [3], Owen [5], and Lettau and Lettau [6] (for Q B , Q K , Q O , and Q L , respectively) were used [13,16,[41][42][43]: Here, D is the reference sand grain diameter (= 0.25 × 10 -3 m) and c B , c K , c O , and c L are constants with values of 1.8, 2.78, 0.25 + w 0 /3u * (where w 0 is the fall velocity), and 4.2, respectively. The sand transport rate was estimated to be zero when u * < u *t . u * was estimated from Equation (1) using the measured value of u; u *t was obtained using Equation (2) and w 0 was determined from Rubey's [44] equation: where s is the relative density of the sand (= σ/ρ) and ν is the kinematic viscosity of air. These equations were derived for a constant wind velocity, and hence, Q estimated using the UD-101 data was averaged over 10 min. The 10-min mean sand flux Q mean during periods without rainfall was compared to Q B , Q K , Q O , and Q L . Figure 7 shows a comparison of Q mean values obtained from n with Q B , Q K , Q O , and Q L values derived from the equations given above, for periods without rainfall. Data for n = 0 and u * < u *t (i.e., u < 4.05 m s -1 ) were excluded. Q mean showed a clear positive relationship with Q B , Q K , Q O , and Q L , and the estimated values followed the order Q O < Q B < Q L < Q K . The best agreement was observed between Q mean and Q O when 0.02 < Q mean < 0.10 and between Q mean and Q B when Q mean > 0.10. Figure 8 shows double logarithmic plots of Q mean versus Q B , Q K , Q O , and Q L for east-northeast (15° < θ mean < 105°), south-southeast (105° < θ mean < 195°), west-southwest (195° < θ mean < 285°), and north-northwest (285° < θ mean < 375°) winds. For small Q mean values, i.e., less than 0.01 kg m -1 s -1 , the Q mean value was overestimated by Q B , Q K , Q O , and Q L . The probable causes for the overestimation are the low measurement accuracy of n for small Q mean , the low estimation accuracy of z 0s for small u * , variations in the wind velocity, and the high moisture content at the beach surface. Raupach [38] showed that c 0 in Equation (3) depended on u * , especially for small u * . A decrease in the instantaneous wind velocity below the threshold velocity, even when the mean wind velocity is greater than the threshold velocity, can lead to an overestimation. The moisture content of the backshore is nonzero even during periods without rainfall because of the proximity of the backshore to the sea. Udo et al. [22] demonstrated that the sand transport flux at a height, q mean , decreased by more than one order of magnitude when the moisture content increased during rainfall periods. Table 1 lists the linear regression relationships between Q mean estimated from n and Q B , Q K , Q O , and Q L for the east-northeast, south-southeast, west-southwest, and north-northwest winds (see Figure 8). While a strong correlation was observed for the north-northwest wind (in the longshore direction), a clear correlation was not observed in the case of the other winds because only a small number of Q mean values were greater than 0.01 kg m -1 s -1 . Udo [31] showed the dependency of q mean on θ mean and that q mean was high when θ mean was in the longshore direction (especially in winter); however, q mean was low when θ mean was seaward (in late summer and winter) and was even lower when θ mean was landward (in spring and autumn). The differences between Q mean and Q B , Q K , Q O , and Q L were found to be independent of θ mean .  east-northeast (15° < θ mean < 105°), south-southeast (105° < θ mean < 195°), west-southwest (195° < θ mean < 285°), and north-northwest (285° < θ mean < 375°) winds. Q B is either zero or greater than 1.76 × 10 -3 kg m -1 s -1 because Equation (11) does not include the term (u * -u *t ). The dashed line denotes perfect agreement. Table 1. Linear regression expressions for the relationships between Q mean estimated from n, and Q B , Q K , Q O , and Q L for east-northeast, south-southeast, west-southwest, and north-northwest winds. R 2 is the regression coefficient and N is the number of data points. See Figure 8 for related plots.  Hotta and Horikawa [40] performed field experiments at beaches by using a trench-type trap and a vertical-distribution-type trap to estimate the aeolian sand transport rate and showed that the possible ranges of c B and c K are 0.5-2.0 and 1.0-3.0, respectively. This result suggests that the actual values of Q B and Q K are possibly one-third of those shown in Figures 7 and 8. All the field results showed that the sand transport rate estimated using the new method together with the UD-101 data was of the same order as that estimated using the conventional equations, although several differences were observed between the aeolian sand transport rate in the wind tunnel and that in field experiments [38,40,45,46]. The new method is useful for the direct estimation of the sand transport rate in a field analysis. Further analysis using the UD-101 data can provide additional insights into aeolian sand transport mechanisms in fields.

Conclusions
In this study, we developed a new method for the estimation of the aeolian sand transport rate and demonstrated the usefulness of UD-101 for field measurements of the aeolian sand transport rate. The transport rates estimated from the UD-101 data were reasonable and showed a clear positive relationship with the rates estimated using the conventional equations. Further, for Q mean > 0.01 kg m -1 s -1 , the transport rates estimated using the UD-101 data were of the same order as those estimated using the conventional equations. The best agreement was observed between Q mean and Q O when 0.02 < Q mean < 0.10 and between Q mean and Q B when Q mean > 0.10. It is necessary to develop a method for accurately determining the constant values in the conventional equations. The cause for the disagreement between Q mean and Q B , Q K , Q O , and Q L for small Q mean should be investigated in a future research. In field analyses, the estimation of large transport rates under strong winds is important and the proposed estimation method is expected to be useful in obtaining the transport rate.