A Field Measurement Based Wind Characteristics Analysis of a Typhoon in Near-Ground Boundary Layer

: The ﬁeld measurement was conducted to observe the wind ﬁeld data of West Paciﬁc typhoon “Maria” in this research. With the application of ultrasonic anemometers installed in different heights (10 m, 80 m, 100 m) of the tower, the three dimensional wind speed data of typhoon “Maria” was acquired. In addition, vane-type anemometers were installed to validate the accuracy of the wind data from ultrasonic anemometers. Wind characteristics such as the mean wind proﬁle, turbulence intensity, integral length scale, and wind spectrum are studied in detail using the collected wind data. The relationship between the gust factor and turbulence intensity was also studied and compared with the existing literature to demonstrate the characteristics of Maria. The statistical characteristics of the turbulence intensity and gust factor are presented. The corresponding conclusion remarks are expected to provide a useful reference for designing wind-resistant buildings and structures.


Introduction
Numerous wind-sensitive structures such as high-rise buildings and large span bridges have been built in recent years, especially in China, due to the rapid economic increase. It is essential to fully study the wind characteristics for the wind resistant design of these structures. Typhoons are natural hazards that may cause substantial economic losses and pose a threat to life. The average loss in China caused by typhoons is over 200 billion yuan per year. The economic loss was estimated to be 970.7 billion yuan in 1996 and has a growing tendency since then [1]. Researchers have paid much attention to the investigation of typhoons. Specific studies have been conducted via simulation methods, wind tunnel tests, and field measurements. Computational fluid dynamics (CFD) simulations have been proven to be an efficient method compared with experimental methods. This method provides detailed wind field information and can also be well controlled under different terrain and flow conditions. However, due to technology limitations, the accuracy of the simulation results still needs to be verified by wind tunnel tests or field measurements. The wind tunnel test is the most commonly used approach to investigate the influence of wind on buildings and structures. However, it still has limitations such as the lack of consideration of detailed terrain conditions or the selection of structural parameters.
Field measurements have been regarded as the most reliable method that can directly provide real long-term wind data and structural responses even though they are costly and time consuming. Numerous researchers have conducted observations and analyses for strong typhoon processes, and thus, the wind characteristics have been well investigated [2][3][4][5][6][7][8]. Wind characteristics such as wind profile, gust factor, and turbulence intensity can be calculated 2 of 22 in different measurement conditions [9,10]. Moreover, non-stationary methods have also been applied to obtain detailed wind characteristics of typhoons [11][12][13].
Specifically, gradual increases in wind measurement research have been conducted for tall buildings in China, and significant results have been obtained. The typhoon characteristics and their influence on the wind-induced response on tall buildings have been investigated in detail in Hong Kong [14][15][16]. For instance, through field wind measurements on a 70-storey tall building in Hong Kong, the wind characteristics of four typhoons were studied, and the typhoon characteristics in the Hong Kong area were discussed [17]. Moreover, to fully study the dynamic response of large span bridges, numerous wind field measurements were performed with installed structural health monitoring systems, which have provided significant recommendations on both the design and construction of long-span bridges [18][19][20][21]. Another consideration for wind characteristic analysis is the estimation of wind power due to the increasing demand for renewable energy worldwide. The potential wind energy was studied in specific areas by the long-term record of both typhoon data and monsoon wind data [22][23][24][25][26]. Although many wind field measurements have been performed to explore wind characteristics and wind load on structures, they may show significant differences in different areas due to the features of the terrain. However, currently, few wind field experiments have been conducted in Fujian Province, even though it is typhoon-prone. In this research, to investigate the wind characteristics in southeastern Fujian Province, a wind field measurement system installed in a 100 m tower was established in Yutou Island, Pingtan County, Fujian Province. By observing the wind data during the typhoon process, the wind characteristics of the typhoon periphery such as the wind profile, turbulence intensity, gust factor, and wind power spectrum are investigated in details.

Experiment Site and Instruments
The measurement site is located on Yutou Island in Pingtan County in Fujian Province, China, with a latitude of 25 • 37 56" N and longitude of 119 • 34 45" E. Because it is surrounded by the sea on the north side, the island has suffered numerous typhoons and contains abundant wind resources. To study the wind characteristics in this area, a 100 m wind tower was built, as shown in Figure 1. Three 3D sonic anemometers as shown in Figure 2 are installed at heights of 10 m, 80 m, and 100 m. The ultrasonic systems have been widely applied in meteorological measurements for their ability to register instantaneous values of the measured characteristics of a turbulent atmosphere with high sampling rates [27,28]. In this paper, the 3D sonic anemometers applied were high-performance Wind Master Pro produced by UK Gill Company with the wind speed range: 0~65 m/s, resolution 0.01 m/s, wind direction range: 0~359 • , resolution: 0.1 • . A 10 Hz sample output frequency was used on site, indicating that one set of 3D instantaneous wind speed was recorded per 0.1 s. The ultrasonic anemometers can work at temperatures of −40~+70 • C, and a small weather station was installed at a 10 m height to monitor the working temperature of the ultrasonic anemometers with an accuracy of the thermometer of 0.01 • C. Identification codes are generated automatically to be used as the reference for data quality control. When the data are unstable or out of the range, it will be diagnosed as 0 in the recorded data. The bad or invalid data points will be eliminated and then made up by the multiple truncation variance data control method [29]. To ensure the accuracy of the 3D sonic anemometer data, six vane-type anemometers, as shown in Figure 3, were installed at heights of 10 m, 30 m, 50 m, 80 m, 90 m, and 100 m, as shown in Figure 1. The distance between the wind tower and anemometers was 1.5 m to eliminate the influence of the tower. The data acquisition system was located at a height of 6 m for data collection. To demonstrate the wind characteristics in detail, the terrain in four directions around the measurement is presented in Figure 4.

Description of the Typhoon "Maria"
The eighth typhoon in 2018, "Maria", appeared on the northwest Pacific Ocean surface approximately 1930 km east of Yilan County, Taiwan Province, at a super typhoon level. The typhoon landed on Lianjiang, Huangqi Peninsula, Fujian Province, China, at 09:10, 11 July, with a maximum wind speed of 42 m/s, a central minimum pressure of 960 hPa, and moved toward the northwest at a speed of 15 km/h. The measurement site was approximately 83.61 km away from the typhoon landing center. The landing path and measurement site are shown in Figure 5.

Description of the Typhoon "Maria"
The eighth typhoon in 2018, "Maria", appeared on the northwest Pacific Ocean surface approximately 1930 km east of Yilan County, Taiwan Province, at a super typhoon level. The typhoon landed on Lianjiang, Huangqi Peninsula, Fujian Province, China, at 09:10, 11 July, with a maximum wind speed of 42 m/s, a central minimum pressure of 960 hPa, and moved toward the northwest at a speed of 15 km/h. The measurement site was approximately 83.61 km away from the typhoon landing center. The landing path and measurement site are shown in Figure 5.

Near-Ground Wind Characteristics of the Typhoon
The record of Typhoon "Maria" was from 0:00 to 24:00 on 11 July 2018, with a duration of 24 h data. Samples were divided every 10 min [30] based on the standard time interval in China. The installed ultrasonic anemometer can automatically distinguish the invalid data caused by rainfall. Moreover, to ensure the accuracy of wind data, samples were selected according to the selected criterion, and only the samples with 98% validation data could be used for the wind characteristic analysis. The data were low pass filtered at 4 Hz to eliminate aliasing caused by high sampling frequency. Due to the high sampling frequency of the ultrasonic anemometer, the recorded wind speed data may include some bad or invalid data. The multiple truncation variance method was applied to smoothly estimate the original data as follows: First, calculate the time series difference as du(t) as: where u(t) is the tth time point of the wind speed data. The mean value of du(t) and du 2 can be expressed as: The truncation variance can be expressed as follows: The criterion to detect invalid data can be defined as: In which, c is the coefficient to determine unreasonable data. In this research, c is taken as 4, indicating that when the absolute value of difference between the mean value of the sample and the total sample was greater than four times the standard deviation, the data will be diagnosed as invalid data. When the invalid data are detected, it will be modified by the five-point interpolation method [29].

Mean Wind Data
The 10 min wind speed and wind direction measured by three sonic anemometers are indicated in Figures 6 and 7 with the distance of the typhoon center to the measurement site. The blue line is the time point when Typhoon "Maria" made landfall. As accepted by researchers, extreme winds in the eyewall regions of typhoons may cause serious damage to buildings and structures. The wind speed time history shows an 'M' shape with double peaks, which means the typhoon wind circle passed through the measurement site. The distance shows that the typhoon center first got close to the measurement site and then got far way. When the typhoon wind circle arrived, the wind speed reached a peak, and when it passed by, the wind speed reached the other peak. This "M" type indicates that the measured wind speed can represent the wind characteristic of the periphery of Typhoon "Maria". In Figure 6, the wind speed increases as the measured height increases, showing consistent varying tendencies for different heights. The maximum mean wind speed reached 26.13 m/s during the landing of a typhoon. The change in the wind direction during the typhoon process was approximately 160 • . To ensure the validity of the measured data, the 10 min average wind speed at a height of 10 m obtained from the vane-type anemometer and sonic anemometer were compared, as shown in Figure 8. The wind speed data acquired by the two types of anemometers agreed well, which indicates the reliability of the obtained wind data.

Fluctuating Wind Data
In this experiment, the output data can be time series t, and the wind speeds in three directions are u x , u y , and u z . Taking 10 min as the statistical interval, the average speed can be expressed as follows: where u x , u y are the mean wind speeds in the directions of x and y, respectively. The fluctuating wind speeds of the longitudinal, horizontal, and vertical wind speeds can be expressed as follows: where u (t), v (t), and w (t) are the fluctuating wind speeds in the longitudinal, horizontal, and vertical directions, respectively. Figure 9 indicates the fluctuating wind speed at a 10 m height during Typhoon Maria.

Mean Wind Characteristics
The wind profile is an essential parameter to describe wind characteristics in a neutrally stable atmospheric boundary layer. The power law (recommended in Code GB5009-2001) and logarithmic law are usually used in wind profile analysis. The power law can be expressed as follows: where z 1 is the reference height, which is considered to be 10 m; u(z) is the wind speed at height z; and α is the wind profile parameter. The wind profile of the measurement site was calculated based on the power law equation. Based on the measured data, the mean wind profile parameter was calculated as 0.2208, which was close to the recommended value of class C ground motion in the Chinese Code. This finding may be caused by the terrain of the measurement site. Even though the north side of the measurement site is close to sea, the south side is faced with trees, as shown in Figure 6. From Figure 9, the main direction of this typhoon is approximately the southeast direction; therefore, the roughness of the forest may slow down the wind speed. Figure 10 shows the wind profile parameters with wind speed. It is obvious that the wind profile parameter decreases as the wind speed increases. The main reason for this phenomenon may be that the wind profile coefficient is to show the influence of ground roughness in the wind speed; with the wind speed increase, the influence of ground roughness decreases, therefore, the wind profile coefficient also decreases. Figure 11 shows the wind profiles for different wind speeds during Typhoon "Maria", indicating that the wind profile significantly varies with the wind speed.

Turbulence Intensity
where σ i is the mean variance of the fluctuating wind speed and U is the mean wind speed in the statistical time interval of 10 min. Figure 12 shows the turbulence intensity at 10 m in the longitudinal, horizontal, and vertical directions. The blue circles indicate the distance between typhoon center and measurement site. Compared with the horizontal and longitudinal directions, the vertical direction showed a smaller turbulence intensity and hardly changed during the typhoon. Figure 13 shows the turbulence intensity in the three directions via the wind speed. As indicated, the turbulence intensity decreased with increasing longitudinal and horizontal directions. However, the turbulence intensity remained almost unchanged in the vertical direction.   Table 1 shows a summary of the measured turbulence intensities during Typhoon "Maria". The turbulence intensity decreases as the height increases. The ratios of the turbulence intensities I u : I v : I w at 10 m, 80 m, and 100 m are 1:0.76:0.47, 1:0.86:0.63, and 1:0.89:0.50, respectively. With increasing height, the ratio of the vertical direction increases. The measured turbulence intensity ratio I u : I v : I w at a height of 100 m is close to that recommended in the Chinese wind-resistant code I u : I v : I w = 1 : 0.88 : 0.5. However, the data showed a difference at heights of 10 m and 80 m, which indicates the difference in wind characteristics between typhoons and seasonal winds, which is the basis of recommendation in the Chinese code. The generalized extreme value (GEV) distribution has been proven to be an efficient distribution developed from extreme value theory that combines the Gumbel, Fréchet, and Weibull functions. The general probability density function of the GEV can be written as follows: where k, σ, and µ are the shape, scale, and location parameters, respectively. When k = 0, the distribution is regarded as the type I case (Gumbel distribution), with k > 0 and k < 0, and the distribution is referred to as the Type II (Frechet distribution) case and Type III (Weibull distribution) case, respectively. The turbulence intensity during this typhoon process is fitted by the GEV distribution by maximum likelihood estimates with a 95% confidence, and the corresponding probability density function is shown in Figure 14. The coefficients fitted by the GEV distribution are listed in Table 2 including the statistical characteristics such as the mean, median standard deviation, skewness, and kurtosis. To verify the accuracy of the fitting results, goodness-of-fit tests such as the Kolmogorov-Smirnov (K-S) test and the Anderson-Darling (A-D) test are conducted with a 95% confidence. It can be seen from the table that the shape values for the three heights were all above 0, indicating that the fitted distributions were Type II cases. This finding means that the turbulence intensity during Typhoon "Maria" can be described by the Freche distribution. From Table 2, it can be seen that all skewness coefficients were positive, indicating that distributions at each height were right-skewed or right-tailed. All kurtosis values were greater than 3.0, indicating that all distributions of Iu are leptokurtic and more peaked than a normal distribution with longer tails. Meanwhile, both the skewness and kurtosis coefficient have a tendency of decrease from 10 m to 80 m and decrease from 80 m to 100 m. The results show differences in the literature [31], in which the turbulence intensity agreed well with the Weibull distribution based on the observation of typhoons in Zhejiang Province, indicating the distinctive features of typhoons. Moreover, it is also due to the unique features of typhoons, the observation of typhoons, and analysis of typhoon characteristics are essential for the wind-resist design of a specific area. It has been widely accepted in the literature that the turbulence intensity decreases with increasing wind speed [31]. Tables 3-5 indicate the statistical value of the turbulence intensity at 10 m, 80 m, and 100 m, with classification by the wind speed. The mean value of the turbulence intensity decreases with the wind speed at all three heights, which agrees with the literature and previous descriptions. Moreover, the variances of the three directions at heights of 10 m, 80 m, and 90 m decrease with increasing wind speed, which agrees with [31].

Gust Factor and Peak Factor
The gust factor is a significant feature of fluctuating wind that describes the proportions of the gust wind speed to the mean wind speed. This factor is defined as the ratio of the maximum mean wind speed during gust duration t s . In the Chinese structural design code, the gust duration is also recommended as 3 s. The expression of the gust factor can be described as follows: where u(t s ) max , v(t s ) max , and w(t s ) max are the longitudinal, horizontal, and vertical maximum average wind speeds, respectively, within the time interval t s . The gust factor may be considerably affected by the measurement height, ground roughness, and time interval. In China, the time interval t s is recommended as 3 s by wind engineers to calculate the gust factor. Figure 15 shows the gust factor at heights of 10 m, 80 m, and 100 m with a time interval of 3 s. The gust factor decreases with increasing height. The same conclusion can also be seen in the statistical characteristics analysis in Table 6. Analogous to the turbulence intensity, the GEV distribution was applied to fit the gust factors at the three heights. The corresponding fitted results are shown in Figure 16. As obtained in Table 6, for the height of 10 m, the fitted shape values were all larger than 0, indicating the Type II cases for the three heights, showing agreement with the turbulence intensity. However, the shape values at 80 m and 100 m approached 0, indicating a type I case (Gumbel distribution) for the gust factors. It should be mentioned that, in Figure 14, outlier data occurred around 14:00, and the same phenomenon could also be found in the plot of turbulence intensity shown in Figure 11. From the wind speed time history, the wind speed significantly decreased around 14:00, which may be one of the reasons for the significant increase in turbulence intensity and gust factors. This also agrees with the conclusion that the turbulence intensity and gust factor decrease with the increase of wind speed ( Figure 13). Moreover, the wind direction changed around 20 • at the time segment around 14:00, which induced the fluctuation of the wind speed. Therefore, the fluctuations of wind speed and wind direction result in the outlier data in turbulence intensity and gust factor.   The time interval has proven to be an important aspect in calculating gust factors. In the Chinese Code, the gust factor is usually recommended as 3 s. In Figure 17, gust factors at different time intervals are presented. It can be seen that the gust factor decreases as the time interval increases. The calculated gust factor was far smaller at a time interval of 1 min than 3 s. Figure 17 shows the mean gust factors under time intervals of 1-60 s. The gust factor shows a linear relation with the logarithm of the time interval in both the longitudinal, horizontal, and vertical directions, as indicated in Figure 18.  The relationship between the turbulence intensity and the gust factor has attracted attention in wind engineering. However, as discussed previously, the gust factor can be affected by the time intervals. Therefore, based on the measured results, the relation between the turbulence intensity and gust factor can be expressed as follows [32]: where G u (t s ) is the longitudinal gust factor at time interval t s ; T is the basic time interval when calculating the mean wind speed, which is 10 min in this research; and k 1 and k 2 are the coefficients at t s . The values of k 1 and k 2 were suggested by Ishizaki [33] as k 1 = 0.5; k 2 = 1.0, by Choi [34] as k 1 = 0.62; k 2 = 1.27, and by Cao [35] as k 1 = 0.5; k 2 = 1.15. Based on the observed experimental data, k 1 and k 2 were fitted as k 1 = 0.3695; k 2 = 0.9224 as fitted in Figure 19, which showed significant differences from the results in the existing literature. Figure 19. Relationship between the turbulence intensity and the gust factor.
The peak factor is also a feature of the fluctuating wind that describes the ratio of the maximum wind speed fluctuation in time interval t s to the standard deviation of the fluctuating wind speed, expressed as: whereÛ(t s ) is the maximum mean wind speed within the interval t s . Figure 20 shows a relationship between the turbulence intensity and gust factor under different time intervals. Substituting Equation (14) into Equation (15), one can obtain the expression between the peak factor and the turbulence intensity:  Figure 20 indicates that the peak factors slightly changed with the wind speed. The minimum, maximum, and mean values of the peak factor were 1.40, 3.93, and 2.23, respectively. Compared with the results obtained by Wang [36], which were 2.43, 2.48 and 2.47 at 10 m, 20 m and 40 m, respectively, the present research showed a larger maximum value but a smaller average value. These changes may be due to the fluctuation in the wind speed in the presented measurement site, which can also be seen by comparing the turbulence intensities.

Turbulence Integral Length Scale
The turbulence integral length scale is defined as the average size of the turbulent eddies of flows, which describe the influence of the fluctuating wind on structures. Based on Taylor's assumption, the expression of the turbulence integral length scale can be defined as follows: (18) where L x u , L y v , and L z w are the measured turbulence integral length scales in the longitudinal, horizontal, and vertical directions, respectively. R u (τ), R v (τ), and R w (τ) are the autocorrelation functions of the fluctuating wind in the longitudinal, horizontal, and vertical directions, respectively, and τ is the time difference. The equation of the autocorrelation function can be presented as in Equation (19).
where X(t) is a random wind process. Based on Equation (19), the autocorrelation coefficients of three components at different heights are presented in Figure 23. Generally, the autocorrelation coefficients decrease with time difference τ; for the vertical components at height 100 m, it was different, which agreed with the results in the literature [1].   Table 7 shows the summary of the ratio of L x u :L y v :L z w in the literature and present research. It can be seen that the ratio between L x u :L y v in this study was close to that in [36]. However, the ratio of L z w was much larger than that in the references in the table, except in [6]. One consideration is that the terrain may significantly affect the integral length scale as the two measurement sites were both located in Pingtan County in Fujian Province, China.

Cross-Correlations and Coherence
The cross-correlations show the correlation between fluctuating wind components, which can be expressed as follows: where C R(ij) is the cross-correlation coefficient of the fluctuating wind component of I and j; R ij is the corresponding cross-correlation coefficient; and R ii and R ij are the autocorrelation of the fluctuating wind components of i and j, respectively. Figures 25 and 26 show the cross-correlation coefficients between the fluctuation wind speeds at heights of 10 and 80 m, respectively. Figures 25a-c and 26a-c show the coefficients between the u, v direction, u, w direction, and v, w direction, respectively. It can be concluded from Figures 25 and 26 that the correlations between u and v, v and w were approximately 0, which showed good agreement with the results in [1]. However, the fluctuation wind speed in the v and w directions showed a negative correlation with the coefficients −0.36 and −0.26, respectively, as seen in Figures 25b and 26b, which was comparable to the results in the neutral surface boundary layer of −0.3 [1].  The coherence model was developed by Davenport in 1961 to describe the correlation among the frequency domain, which is used as an exponential decaying model: where f is the frequency; U is the mean wind speed; C y and C z are the exponential decay coefficients in the lateral and vertical directions, respectively; and ∆y and ∆z are the distances between two measured points in the lateral and vertical directions, respectively. In this experiment, the distance in the lateral direction was 0, and then the decay model can be expressed as follows: Figures 25-27 show the coherence at different wind speeds in the three directions. It can be seen that the decay coefficient increased with increasing wind speed in all directions, which agreed with the results in the literature [1]. Comparing Figures 27-29, the coherence coefficient in the v and w directions can be better fitted than that in the u direction by the exponential decay model.

Power Spectrum Density (PSD) Function
The distribution of the fluctuating wind energy with respect to frequency, described as the power spectrum density (PDS), is significant for predicting wind loading on structures and the corresponding response. Low-frequency or larger eddies may dominate the wind energy, and high-frequency or lower eddies may be ignored in the wind load analysis. The descriptions of the PSD by the Karman model are expressed as follows [37,38]: Longitudinal direction: Horizontal and vertical directions: where f = nL x u U ,L x u is the turbulence integral scale; U is the mean wind;β is the coefficient of the friction velocity; σ 2 u = βu 2 * ; and σ u is the variance of the fluctuating wind speed components.
The Davenport spectrum, which is recommended in the Chinese Code, has also been widely used for modeling the longitudinal wind direction. The expression is described as follows: nS(z, n) where f = 1200n U(10) and U(10) is the mean wind speed at a height of 10 m. The Panofsky model has been widely used for the calculation of the power spectrum for the vertical wind speed component with the following expression [39]: Figure 30 shows the wind power spectra during Typhoon "Maria" in the longitudinal wind direction at time intervals of 10 min, 20 min, and 60 min. It can be seen that the general tendency for the wind spectra at time intervals were similar. However, for the larger time interval, the amplitudes at high frequency were larger. As the time interval of 10 min is recommended in the wind-resist design standard in China, the following analysis was based on a time interval of 10 min. larger time interval, the amplitudes at high frequency were larger. As the time interval of 10 min is recommended in the wind-resist design standard in China, the following analysis was based on a time interval of 10 min. Figure 31 shows the wind power spectra during Typhoon "Maria" in the longitudinal wind direction. The measured power spectra were compared with the spectra of the Karman and Davenport models. It can be concluded that the Karman model can be fitted better than the Davenport model in the longitudinal direction, even though the Davenport spectrum shows better fitting performance in the high-frequency area. The measured vertical spectra were fitted by the Karman model and Panofsky model as shown in Figure 32. Comparing the two models, the Panofsky model showed better results in the fitting of the measured wind spectrum at lower frequencies at a 10 m height. However, for heights of 80 m and 100 m, the measured wind spectra agreed well with the Panofsky model at both high and low frequencies.   Figure 31 shows the wind power spectra during Typhoon "Maria" in the longitudinal wind direction. The measured power spectra were compared with the spectra of the Karman and Davenport models. It can be concluded that the Karman model can be fitted better than the Davenport model in the longitudinal direction, even though the Davenport spectrum shows better fitting performance in the high-frequency area. The measured vertical spectra were fitted by the Karman model and Panofsky model as shown in Figure 32. Comparing the two models, the Panofsky model showed better results in the fitting of the measured wind spectrum at lower frequencies at a 10 m height. However, for heights of 80 m and 100 m, the measured wind spectra agreed well with the Panofsky model at both high and low frequencies.  Figure 30 shows the wind power spectra during Typhoon "Maria" in the longitudinal wind direction at time intervals of 10 min, 20 min, and 60 min. It can be seen that the general tendency for the wind spectra at time intervals were similar. However, for the larger time interval, the amplitudes at high frequency were larger. As the time interval of 10 min is recommended in the wind-resist design standard in China, the following analysis was based on a time interval of 10 min. Figure 31 shows the wind power spectra during Typhoon "Maria" in the longitudinal wind direction. The measured power spectra were compared with the spectra of the Karman and Davenport models. It can be concluded that the Karman model can be fitted better than the Davenport model in the longitudinal direction, even though the Davenport spectrum shows better fitting performance in the high-frequency area. The measured vertical spectra were fitted by the Karman model and Panofsky model as shown in Figure 32. Comparing the two models, the Panofsky model showed better results in the fitting of the measured wind spectrum at lower frequencies at a 10 m height. However, for heights of 80 m and 100 m, the measured wind spectra agreed well with the Panofsky model at both high and low frequencies.  According to the Kolmogorov turbulence theory, the fluctuating wind power spectra can be expressed by the following empirical expression: (1 + B f α ) β (27) where σ 2 u is the variance of the longitudinal wind component; A, B, α, and β are the coefficients that need to be fitted; and f is the converted frequency. The fitted power spectrum can be seen in the Figure 33. The fitted value of the empirical expression can be obtained in the following Table 8:

Conclusions
The wind characteristics of the periphery of Typhoon "Maria" were studied in this research on the basis of the observed wind speed data recorded from sonic anemometers on a 100 m tower. Based on the analysis results, the following conclusions can be summarized: (1) The field measurement was conducted with the applications of both ultrasonic and vane-type anemometers, and wind data of the periphery of Typhoon "Maria" were therefore recorded. (2) Data control methods were applied to filter the high frequency wind data to obtain more reasonable data for analysis. (3) The reliability of the obtained sonic wind speed data was guaranteed by a comparison of the mean wind speed between the data from sonic anemometers and vane anemometers. (4) The wind profile parameter calculated by the power law was 0.2208, which was larger than that recommended in the Chinese Code; the wind profile parameter decreased with increasing wind speed. (5) The turbulence intensity decreased with increase in wind speed. The gust factor showed a linear relation with the logarithm of the time interval in both the longitudinal, horizontal, and vertical directions. (6) Based on the statistical characteristics, the turbulence intensity and gust factor followed a Type II GEV distribution. The gust factor decreases with an increasing time interval and has a linear relationship with the logarithm of the time interval. The maximum and mean values of the peak factor during Typhoon "Maria" were 3.93 and 2.23, respectively. The fitted parameter of the relationship between the gust factor and turbulence intensity was comparably studied. (7) The turbulence integral length scale of the periphery of Typhoon "Maria" L x u :L y v :L z w = 1:0.72:0.59, the ratio of the vertical component was larger than that in the literature but close to those of Typhoons Naseat and Hiatang, which also occurred in Fujian Province. (8) The correlations between u and v, v and w were approximately 0, which showed good agreement with the results in the literature, and the fluctuation wind speed in the v and w directions showed a negative correlation. The decay coefficient increased with increasing wind speed in all directions, and the decay model could be better fitted better in the v and w directions.