Spatiotemporal Dynamics of Maximum Wind Speed Using the Wind Multiplier Downscaling Method in the Yangtze River Inland Waterway from 1980 to 2017

: Wind speed affects the navigational safety of the Yangtze River, and assessing its spatiotemporal dynamics provides support for navigation management and disaster prevention. We developed a wind multiplier downscaling method integrating the effects of land use and topography, and used meteorological station observations and European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis Interim (ERA-Interim) reanalysis data for statistical downscaling in the Yangtze River inland waterway region from 1980 to 2017. Compared with reanalysis data, the downscaling products showed improved accuracy (especially at 5–10 m/s), and are consistent with site-based interannual variability observations. Increasing maximum wind speeds in the middle–downstream area was observed from 1980 to 1990, while a decreasing trend was observed from 2010 to 2017; the opposite was observed for the upstream. Land use has signiﬁcant inﬂuence on wind speed, with a decreasing trend observed year by year for wind speed above grade 9. Although the proportion of grade 4–8 wind speed over water is small and the trend is not obvious, grade 9–10 wind speeds displayed an increasing trend from 2010 to 2017, indicating that changes in surface roughness have a signiﬁcant inﬂuence on wind speed in the Yangtze River inland waterway. decreased signiﬁcantly from 2010 to 2017. (3) V max_WM shows spatial characteristics of large nearshore wind speed (14–37 m/s), small inland wind speed (11–20 m/s), and wind speed over undulating The V max_WM value in the terrain section (Yibin and Chongqing) increases signiﬁcantly and expands westward, the center of high-value region of the Yueyang and Jiujiang section narrows and decreases, and the center of the high-value region in the and section to the in the Yangtze River may be a sensitive indicator of global climate change. It should also be noted that the Yangtze River inland waterway is highly exposed and vulnerable and is also affected by transport policy, and socio-economic and technological developments. Our study provides valuable data and knowledge support for navigational risk management, identiﬁcation of local areas at high wind speed risk, and wind emergency management in the Yangtze River. Research on the trends of extreme weather events on the Yangtze River inland waterway under the


Introduction
Large-scale human activities have led to drastic changes in the frequency of strong winds over the land surface, considerably endangering the safety of inland navigation and sustainable regional socio-economic development [1,2]. The cargo throughput of Yangtze River inland waterway transport accounts for approximately 55% of total cargo transport in China. Strong winds and associated secondary disasters (e.g., power disruption, bridge damage, ship accidents, etc.) have become the main destructive weather-and climaterelated events along the Yangtze River inland waterway transportation [3]. For example, in June 2015, a tornado caused the sinking of the "Eastern Star" passenger ship [4]. The traditional view is that wind speeds in the Yangtze River basin have generally been decreasing over the past 50 years [5], which may be related to the rising Arctic Oscillation (AO) index, urbanization, and large-scale climate change [6,7]. However, a new study found a rebounding trend in global surface wind speeds from 2010 to 2017, with a 17±2% increase in potential wind energy [8]. Climate change has significantly affected the frequency, intensity, and magnitude of variability in global wind speeds [9]. Land use/land cover changes (e.g., urbanization) can affect the surface roughness and cause small-scale wind disturbances [10]; spatiotemporal variation in wind speed at large scales is further influenced by wind speed perturbations [11]. Although the spatiotemporal variation relationships of wind speed at regional-global scales are well understood, characterizing multi-scale maximum wind speed variation due to urbanization and other surface processes in the Yangtze River inland waterway remains problematic.
Wind speed data for inland waterways are mainly derived from river-shore meteorological stations, 80-100 m meteorological observation towers, and climate data; the meteorological observations are often limited to monitoring environments, with missing data and low estimation [12]. European Center for Medium-Range Weather Forecast (ECWMF) reanalysis data can be used to characterize long time series and large-scale climate changes caused by atmospheric circulation and greenhouse gas emissions, which are also a source of data for analyzing the spatiotemporal characteristics of wind speeds at regional scales [13]. However, ECWMF reanalysis data do not sufficiently consider land surface processes during assimilation, resulting in low estimates of wind speed and low sensitivity to land use change [14].
To address these issues, scale-down methods are generally used to downscale coarsegrid reanalysis data to local values. For example, statistical downscaling methods have been used to extrapolate and interpolate riverbank observations to obtain river surface wind speeds [15], dynamic downscaling methods such as numerical simulations of weather forecasts have been adopted to obtain accurate data [16], and coupled dynamic-statistical downscaling methods (e.g., the Wind Atlas Analysis and Application program) can also be used [17,18]. These methods can effectively reduce prediction-observation errors due to missing local features in climate models [19]. Specifically, dynamic downscaling methods are numerical weather prediction models get their initial and boundary conditions from a low-resolution global forecast system and are interpolated into the model domain of interest [20], which are sensitive to the boundary conditions, often lack assimilation in the initialization of mesoscale models, require a large computational volume, and have the uncertainty of interpolation methods and large systematic spatial errors. [21,22]; traditional statistical methods are based on modeling meteorological observation factors and do not consider the influence of geographical elements such as land use and topographic environment on wind speed; coupled dynamic-statistical downscaling methods are limited by the computational capacity of computational fluid dynamics models and cannot provide large-scale and long time series products. Among them, improved statistical downscaling methods can be used to consider the effects of surface processes and effectively reduce reanalysis data errors [23], with good reliability and applicability [24], such as the wind multiplier method that integrates the effects of land use and topography [25,26]. Therefore, this method can be developed for complex geographical environments; in the case of this study, we used it to provide a high-precision wind speed dataset for the Yangtze River inland waterway.
To highlight the above key scientific issues and knowledge gaps, a regional wind multiplier downscaling method integrating high-precision land cover and topographic local effects was proposed. Using this method, reanalysis wind speed data for 1980,1990,2000,2010, and 2017 at 0.7 • × 0.7 • resolution was downscaled from a coarse grid to high spatial resolution data (1 km × 1 km: 1980-2010; 10 m × 10 m: 2017), and the spatiotemporal characteristics of the downscaled maximum wind speed we determined. We focused on assessing wind speed heterogeneity in the water-land area. The results provide scientific reference and data support for navigation management and disaster prevention and mitigation in the Yangtze River inland waterway. Specifically, (1) the land use and elevation datasets were integrated to obtain the surface roughness parameter M z and the slope parameter M h , respectively, and the wind multiplier downscaling method was used to obtain a high-resolution wind speed dataset, which is applicable to the navigable environment of the Yangtze River inland waterway; (2) the observed and downscaled wind speed data were compared and analyzed using spatial mapping methods, and the accuracy of downscaled products was assessed using three error test metrics (the Kolmogorov-Smirnov two-sample test, mean difference, and R 2 ) based on comparison to 51 meteorological stations in the study area; (3) the spatiotemporal characteristics of the downscaled maximum wind speed data under different land use types were analyzed, and the applicability and potential improvement of the wind multiplier downscaling method was considered.

Data
Based on the segment management regulations of Changjiang Waterway Bureau (http://www.cjhdj.com.cn/, accessed on 13 September 2021), the Yangtze River inland waterway (104 • 28 1 E-121 • 54 10 E, 27 • 39 30 N-32 • 42 28 N) was divided into 9 sections (Figure 1): Shanghai section, Nanjing section, Anhui section, Jiujiang section, Yueyang section, Wuhan section, Yichang section, Chongqing section, and Yibin section. In accordance with the characteristics of waterway maintenance, management, and research practices [27], the inland waterways of the Yangtze River were further grouped into three major waterway areas using ArcGIS 10.5 software (esri, Redlands, CA, USA): the upstream waterway from the Yibin section to the Chongqing section, the midstream waterway from the Yichang section to the Jiujiang section, and the downstream waterway from the Anqing section to the Shanghai section.
Atmosphere 2021, 12, x FOR PEER REVIEW 3 of 18 maximum wind speed data under different land use types were analyzed, and the applicability and potential improvement of the wind multiplier downscaling method was considered.

Data
Based on the segment management regulations of Changjiang Waterway Bureau (http://www.cjhdj.com.cn/, accessed on 13 September 2021), the Yangtze River inland waterway (104°28′1″ E-121°54′10″ E, 27°39′30″ N-32°42′28″ N) was divided into 9 sections (Figure 1): Shanghai section, Nanjing section, Anhui section, Jiujiang section, Yueyang section, Wuhan section, Yichang section, Chongqing section, and Yibin section. In accordance with the characteristics of waterway maintenance, management, and research practices [27], the inland waterways of the Yangtze River were further grouped into three major waterway areas using ArcGIS 10.5 software (esri, Redlands, CA, USA): the upstream waterway from the Yibin section to the Chongqing section, the midstream waterway from the Yichang section to the Jiujiang section, and the downstream waterway from the Anqing section to the Shanghai section.  Map of the study area ( 1 Shanghai section, 2 Nanjing section, 3 Anhui section, 4 Jiujiang section, 5 Yueyang section, 6 Wuhan section, 7 Yichang section, 8 Chongqing section, and 9 Yibin section).
Meteorological observations were obtained from the Global Surface Summary of the Day (GSOD) dataset (https://www.ncei.noaa.gov/access/search/data-search/globalsummary-of-the-day, accessed on 13 September 2021) for 1980-2017. The dataset was subjected to multi-process automated data quality control to eliminate random errors and Atmosphere 2021, 12, 1216 4 of 17 mistakes [28]. First, we obtained 51 available stations from 101 ground-based meteorological stations in the study area to satisfy the time continuity and data integrity of strong wind analysis, and then extracted the 10 m daily maximum wind speed (V max_OBS ) to validate the accuracy of the maximum wind speed downscaling product. Second, we extracted 14 representative riverine stations near the river surface from the 51 stations (Table 1) to analyze the maximum wind speed trend within the navigation segment. The reanalysis information was obtained from the ECMWF Reanalysis Interim (ERA-Interim) dataset (https://apps.ecmwf.int/datasets/, accessed on 13 September 2021). The ERA-Interim product is a uniform Gaussian grid surface grid point field (horizontal latitude and longitude grid 0.7 • × 0.7 • ); it is operated by the ECMWF Integrated Forecasting System cycle 31r2 and integrates a variety of physical processes of sea-air-land interactions. The product is credible and can effectively describe the spatiotemporal dynamics of the atmosphere, wind, and waves that affect navigation [29]. Specifically, the 10 m maximum wind speed (V max_ERA ) of the ERA-Interim dataset was extracted using the Python 3.8 software. The geographic data primarily included land use data and a digital elevation model (DEM). The 1980-2015 Landsat dataset in China (1 km × 1 km) was obtained from the Chinese Academy of Sciences Resource and Environmental Science Data Center (http: //www.resdc.cn/, accessed on 13 September 2021) [30]; this dataset is based on Landsat's Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM) images, through manual visual interpretation and multiple sources. The land use types include 6 primary types (arable land, forest land, grassland, water, residential land, and unused land) and 25 secondary types. The 2017 land use data were derived from the Finer Resolution Observation and Monitoring of Global Land Cover 10 (FROM-GLC10) dataset (http://data.ess.tsinghua.edu.cn/, accessed on 13 September 2021), which successfully produced the world's first 10-m resolution 10-classification land cover map using a finite sample stable classification method [31]. Land use data for the years 1980, 1990, 2000, and 2010 were reclassified into cropland, forest, grassland, shrubland, wetland, water, tundra, impervious surface, bare land, and snow/ice based on the classification descriptions of the 2017 land use data. The DEM was derived from the Shuttle Radar Topography Mission (SRTM) dataset (90 m × 90 m, https://cgiarcsi.community/data/srtm-90m-digital-elevation-database-v4-1/, accessed on 13 September 2021) [32] and resampled to the spatial resolution of the corresponding land use data. Then, the land use data and DEM were subjected to standard pre-processing such as projection transformation (Krasovsky_1940_Transverse_Mercator) and clipping [33] to obtain the DEM data of the study area ( Figure 1) and a land use dataset from 1980 to  Figure 2). Finally, the spatial analysis module of ArcGIS 10.5 was used to obtain the slope and land use types for the input data of the wind multiplier downscaling method.
(SRTM) dataset (90 m × 90 m, https://cgiarcsi.community/data/srtm-90m-digital-elevation-database-v4-1/, accessed on 13 September 2021) [32] and resampled to the spatial resolution of the corresponding land use data. Then, the land use data and DEM were subjected to standard pre-processing such as projection transformation (Krasov-sky_1940_Transverse_Mercator) and clipping [33] to obtain the DEM data of the study area ( Figure 1) and a land use dataset from 1980 to 2017 ( Figure 2). Finally, the spatial analysis module of ArcGIS 10.5 was used to obtain the slope and land use types for the input data of the wind multiplier downscaling method.

Wind Multiplier Method
The wind multiplier method integrates the effects of land use types and slope on wind speed by using the surface roughness coefficient (Mz) and slope coefficient (Mh) to reduce the 10 m maximum wind speed (Vmax_ERA) at the ERA-Interim grid point. Thus, downscaled maximum wind speed (Vmax_WM) datasets with 1 km and 10 m resolutions were obtained for the Yangtze River inland waterway for the period 1980-2017.
where Vmax_ERA is the maximum wind speed (m/s) at the ERA-Interim grid point. Mz was calculated by reference to the Australian wind load standard AS/NZS 1170.2 [34] for roughness (Z0); terrain category (TC) correspond to the original Mz, which is widely used and has a parametric solution of credibility. A table of the TC, elevation, and Z0 was first obtained with reference to the above corresponding relationships ( Table 2): water and grassland were classified as Class 2 features with heights of 1.5-10 m and a Z0 of 0.02 m; bareland, cropland, wetland, and shrubland were Class 3 features with heights of 0.5-10 m. Mh was used to characterize the effect of topography on wind speed, which is influenced by land steepness and the distance between a site and the nearest topographical high point [35]; there are significant local terrain effects near the tops of the mountains

Wind Multiplier Method
The wind multiplier method integrates the effects of land use types and slope on wind speed by using the surface roughness coefficient (M z ) and slope coefficient (M h ) to reduce the 10 m maximum wind speed (V max_ERA ) at the ERA-Interim grid point. Thus, downscaled maximum wind speed (V max_WM ) datasets with 1 km and 10 m resolutions were obtained for the Yangtze River inland waterway for the period 1980-2017.
where V max_ERA is the maximum wind speed (m/s) at the ERA-Interim grid point. M z was calculated by reference to the Australian wind load standard AS/NZS 1170.2 [34] for roughness (z 0 ); terrain category (TC) correspond to the original M z , which is widely used and has a parametric solution of credibility. A table of the TC, elevation, and z 0 was first obtained with reference to the above corresponding relationships ( Table 2): water and grassland were classified as Class 2 features with heights of 1.5-10 m and a z 0 of 0.02 m; bareland, cropland, wetland, and shrubland were Class 3 features with heights of 0.5-10 m. M h was used to characterize the effect of topography on wind speed, which is influenced by land steepness and the distance between a site and the nearest topographical high point [35]; there are significant local terrain effects near the tops of the mountains and along the edges of cliffs [36]. . In summary, distribution maps and parameter sets for M z and M h were obtained for the study area for different periods; for example, good spatial heterogeneity was found in both parameters in 2017 ( Figure 3). and along the edges of cliffs [36]. New Mh values were obtained from the original M slope corresponding relationship in AS/NZS 1170.2 (2011): 1.0 (< 2.9°), 1.08 (2.9°-5.7°), 1 (5.7°-11.3°), 1.32 (11.3°-16.7°), 1.48 (16.7°-24.2°) and 1.71 (≧ 24.2°). In summary, distri tion maps and parameter sets for Mz and Mh were obtained for the study area for differ periods; for example, good spatial heterogeneity was found in both parameters in 2 ( Figure 3).

Accuracy Evaluation of Downscaling Wind Speed
We compared the similarity between Vmax_OBS and Vmax_ERA or Vmax_WM by scatterp and analyzed the performance of Vmax_ERA and Vmax_WM by using the D statistic of the tw sample Kolmogorov-Smirnov (K-S) test. This method is a non-parametric test based the empirical cumulative distribution function (ECDF) and is used to assess the differe in ECDF between N data points (Y1, Y2, ..., YN) in two datasets. First, the ECDF of Vmax_  3 Anhui section, 4 Jiujiang section, 5 Yueyang section, 6 Wuhan section, 7 Yichang section, 8 Chongqing section, and 9 Yibin section).

Accuracy Evaluation of Downscaling Wind Speed
We compared the similarity between V max_OBS and V max_ERA or V max_WM by scatterplot, and analyzed the performance of V max_ERA and V max_WM by using the D statistic of the twosample Kolmogorov-Smirnov (K-S) test. This method is a non-parametric test based on the empirical cumulative distribution function (ECDF) and is used to assess the difference in ECDF between N data points (Y 1 , Y 2 , ..., Y N ) in two datasets. First, the ECDF of V max_OBS -V max_ERA and V max_OBS -V max_WM was calculated, then the ECDF was sorted in ascending order to obtain the absolute difference of the cumulative frequency series of the sample points. Finally, the D statistic of K-S was calculated for these sequences. The mean absolute error (MAE) is the average of the absolute value of the error, and the coefficient of determination (R 2 ) is the goodness fit between V max_OBS and V max_ERA or V max_WM , respectively.
where Y N is the ECDF value at point N, and n i is the number of points less than Y i . ECDF is a step function that jumps up by 1/n at each of the n data points when there are repeated wind speed data. Y i,WM , Y i,ERA , and Y i,OBS are the ECDF values of V max_ERA , V max_WM , and V max_OBS with 1/n equal probability, respectively. V i,ERA , V i,WM , and V i,OBS are the maximum wind speed values for each station in the dataset, with smaller MAE values indicating a good result as a difference between V i,WM and V i,OBS is smaller than the difference between V i,ERA and V i,OBS . The increasing R 2 value indicates an improvement when comparing a regression of V i,WM and V i,OBS to a regression of V i,ERA and V i,OBS .

Nexus Analysis of Wind Speed and Land Use
The area ratio of the downscaled maximum wind speed AV i (i = 1, ..., 10) was classified by wind speed levels in the study area (Table 3), and the area ratio of different land use types ATC j (j = 1, ..., 8) was obtained using land use data in the study area. Finally, the area ratio of wind speed A k per unit area for each land type was calculated after normalization (Equation (10)), where k is the number of wind speed-land type classification units (k = 1, ..., 80). The nexus analysis of wind speed and land use assesses the extent and impact of land use on wind speed and provided information for wind speed risk assessment.

Evaluation of Downscaling Effect of Maximum Wind Speed
The downscaled annual maximum wind speed (V max_WM ) was calculated using ERA-Interim data and Equation (1) in the study area from 1980 to 2017, and it was used to evaluate the similarity of V max_OBS in comparison with V max_ERA (Figure 4). The scatterplot of wind speed shows that the majority of the V max_OBS  Boxplots of the three statistics show comparisons of maximum wind speed between Vmax_OBS and Vmax_ERA and between Vmax_OBS and Vmax_WM ( Figure 5): the downscaled D-statistic value decreases from 0.941 to 0.294 in the upstream, from 0.471 to 0.294 in the midstream, and from 0.235 to 0.412 in the downstream, indicating that the difference between the upstream and midstream samples is small and the performance is optimal in highaltitude upstream areas where the terrain is more complex; the difference between the downstream data increases and the effect is averaged. The D-statistic value of the Yangtze River inland waterway decreased from 0.49 to 0.27, indicating that the accuracy of the downscaled product is generally better than that of Vmax_ERA. The MAE value after Boxplots of the three statistics show comparisons of maximum wind speed between V max_OBS and V max_ERA and between V max_OBS and V max_WM ( Figure 5): the downscaled D-statistic value decreases from 0.941 to 0.294 in the upstream, from 0.471 to 0.294 in the midstream, and from 0.235 to 0.412 in the downstream, indicating that the difference between the upstream and midstream samples is small and the performance is optimal in high-altitude upstream areas where the terrain is more complex; the difference between the downstream data increases and the effect is averaged. The D-statistic value of the Yangtze River inland waterway decreased from 0.49 to 0.27, indicating that the accuracy of the downscaled product is generally better than that of V max_ERA . The MAE value after downscaling decreased from 2.34 to 1.73 in the upstream, from 0.91 to 0.82 in the midstream, and from 1.34 to 1.00 in the downstream; the MAE value of the Yangtze River inland waterway changed from 2.34 to 0.745, although the change was not obvious in the downstream. The downscaled data scales tend to be concentrated, with only individual stations susceptible to extreme values. The downscaled R 2 values increased from 0.199 to 0.905 in the upper reaches, from 0.082 to 0.494 in the middle reaches, and from −0.009 to 0.898 in the lower reaches (the most significant increase); the overall R 2 value of the Yangtze River inland waterway changed from 0.026 to 0.72, indicating a significant improvement in the goodness of fit of V max_WM after downscaling. Despite the local variability and bias in the indicator analysis, the products obtained by the wind multiplier downscaling method are more suitable for characterizing the spatiotemporal variability of maximum wind speeds in the study area as compared with the ERA-Interim reanalysis data.

Spatiotemporal Characteristics of Maximum Wind Speed
The maximum wind speed values for the 14 stations in the Yangtze River channel, including Vmax_OBS, Vmax_ERA, and Vmax_WM, were extracted from downscaled wind-speed products. Temporal trend analysis shows (  The above statistical analysis shows that the upstream downscaled data have the smallest differences compared to the observed values and the best fit, but there are extreme values resulting in large differences and low concentration; the midstream downscaled data have small differences compared to the observed values and the highest concentration, but the fit is average; the downstream data have good concentration and fit, but the individual data are substantially different, resulting in lower sample quality. Despite the local variability and bias in the indicator analysis, the products obtained by the wind multiplier downscaling method are more suitable for characterizing the spatiotemporal variability of maximum wind speeds in the study area as compared with the ERA-Interim reanalysis data.
The spatial distribution of V max_WM from 1980 to 2017 reflects the spatial characteristics of large nearshore, small inland, and undulating topographical areas (Figure 7). The spatial distribution range of wind speeds on the Yangtze River inland waterway from 1980 to 1990 increased significantly, with Vmax_WM reaching over 20 m/s in the northeast direction in the Shanghai section and wind speed variation on a downscale significantly higher than the corresponding ERA-Interim data. The spatial maximum wind speed decreased by less than 16 m/s m/s downstream of the Yangtze River inland waterway from 1990 to 2000, but the spatial variation was not significant in the mid-upper stream. The spatial distribution range of wind speeds on the Yangtze River inland waterway from 1980 to 1990 increased significantly, with V max_WM reaching over 20 m/s in the northeast direction in the Shanghai section and wind speed variation on a downscale significantly higher than the corresponding ERA-Interim data. The spatial maximum wind speed decreased by less than 16 m/s m/s downstream of the Yangtze River inland waterway from 1990 to 2000, but the spatial variation was not significant in the midupper stream.

Nexus Relationship between Downscaling Wind Speed and Land Use
Changes in area by land use type were found (Figure 8): the proportion of cropland area has been decreasing year by year, accelerating from 43.39% before 2010 to 33.28% in 2017, a significant decrease of~10.11% in 8 years. Forestland increased steadily from 31.86% in 1980 to 33.28% in 2017; the proportions of grassland, shrubland, wetland, and water area remained largely unchanged. The proportion of impervious surfaces has increased annually from 1980 and accelerated after 2010, from 2.75% in 1980 to 12.18% in 2017, an increase of approximately 9.43%, indicating continued urbanization and an accelerating trend in recent years. The proportion of tundra and snow area is approximately 0, but a very small amount of glacier is present in the study area, and its effect on local wind speed variation remains to be studied.
Atmosphere 2021, 12, x FOR PEER REVIEW 13 of 18 Changes in area by land use type were found (Figure 8): the proportion of cropland area has been decreasing year by year, accelerating from 43.39% before 2010 to 33.28% in 2017, a significant decrease of ~10.11% in 8 years. Forestland increased steadily from 31.86% in 1980 to 33.28% in 2017; the proportions of grassland, shrubland, wetland, and water area remained largely unchanged. The proportion of impervious surfaces has increased annually from 1980 and accelerated after 2010, from 2.75% in 1980 to 12.18% in 2017, an increase of approximately 9.43%, indicating continued urbanization and an accelerating trend in recent years. The proportion of tundra and snow area is approximately 0, but a very small amount of glacier is present in the study area, and its effect on local wind speed variation remains to be studied. Nexus relationship analysis shows (Figure 9) that arable Ak values increase with wind speed and then decrease; the proportion of grade 4 wind speed is 0.18-0.28, the proportion of grade 5 wind speed is the largest (0.36-0.52), the proportion of grade 6 wind speed is 0.21-0.24, the proportion of grade 7 wind speed is 0.04-0.11, the proportion of grade 8 wind speed is 0.01-0.04, the proportion of grade 9 wind speed is 0-0.02. Wind speed of grade 10 or above accounts for 0.02-0.24, grade 7 for 0.04-0.11, grade 8 for 0.01-0.04, grade 9 for 0-0.02, and grade 10 or above for <0.02. The rate of change for 2010-2017 was significantly higher than that for 1980-2010; for example, the wind speed of grade 4 accelerated from 0.29 to 0.18, the wind speed of grade 5 accelerated from 0.37 to 0.52, the wind speed of grade 6 increased from 0.21 to 0.24, and the wind speed of grade 7 and above decreased year by year. The pattern of change for forestland is similar to that of cropland, with the proportion of grade 4 wind speeds ranging from 0.13 to 0.22, grade 5 wind speed ranging from 0.22 to 0.30, and grade 6 wind speed ranging from 0.18 to 0.31. From 2010 to 2017 grade 6 wind speed accelerated down from 0.20 to 0.13, grade 5 wind speed accelerated up from 0.23 to 0.30, and grade 6 wind speed accelerated down from 0.19 to 0.31; the proportion of grade 5-7 wind speed increased year by year, of which the proportion of grade 6-7 wind speed increased significantly. The proportion of wind speeds for the remaining land use types (especially above grade 8) declined significantly from 2010 to 2017, even reaching a value of 0. The proportion of grade 7 wind speed and above for grassland, shrubland, and wetland was higher than that for cropland and woodland, and the proportion of grade 5-7 wind speed increased slightly, while the proportion of wind speed in the remaining types (especially gale force 9 and above) decreased significantly, or even to 0. The proportion of grade 4-6 wind speed in the bare land decreased year by year, but the proportion of grade 7-10 wind speed increased more significantly, by approximately one order of magnitude. The proportion of wind speed above grade 11 decreased significantly, even to 0, indicating that the areas where these land types exist are prone to gale force weather. The proportion of grade 4-8 wind speed in water areas Nexus relationship analysis shows (Figure 9) that arable A k values increase with wind speed and then decrease; the proportion of grade 4 wind speed is 0.18-0.28, the proportion of grade 5 wind speed is the largest (0.36-0.52), the proportion of grade 6 wind speed is 0.21-0.24, the proportion of grade 7 wind speed is 0.04-0.11, the proportion of grade 8 wind speed is 0.01-0.04, the proportion of grade 9 wind speed is 0-0.02. Wind speed of grade 10 or above accounts for 0.02-0.24, grade 7 for 0.04-0.11, grade 8 for 0.01-0.04, grade 9 for 0-0.02, and grade 10 or above for <0.02. The rate of change for 2010-2017 was significantly higher than that for 1980-2010; for example, the wind speed of grade 4 accelerated from 0.29 to 0.18, the wind speed of grade 5 accelerated from 0.37 to 0.52, the wind speed of grade 6 increased from 0.21 to 0.24, and the wind speed of grade 7 and above decreased year by year. The pattern of change for forestland is similar to that of cropland, with the proportion of grade 4 wind speeds ranging from 0.13 to 0.22, grade 5 wind speed ranging from 0.22 to 0.30, and grade 6 wind speed ranging from 0.18 to 0.31. From 2010 to 2017 grade 6 wind speed accelerated down from 0.20 to 0.13, grade 5 wind speed accelerated up from 0.23 to 0.30, and grade 6 wind speed accelerated down from 0.19 to 0.31; the proportion of grade 5-7 wind speed increased year by year, of which the proportion of grade 6-7 wind speed increased significantly. The proportion of wind speeds for the remaining land use types (especially above grade 8) declined significantly from 2010 to 2017, even reaching a value of 0. The proportion of grade 7 wind speed and above for grassland, shrubland, and wetland was higher than that for cropland and woodland, and the proportion of grade 5-7 wind speed increased slightly, while the proportion of wind speed in the remaining types (especially gale force 9 and above) decreased significantly, or even to 0. The proportion of grade 4-6 wind speed in the bare land decreased year by year, but the proportion of grade 7-10 wind speed increased more significantly, by approximately one order of magnitude. The proportion of wind speed above grade 11 decreased significantly, even to 0, indicating that the areas where these land types exist are prone to gale force weather. The proportion of grade 4-8 wind speed in water areas was the opposite of and does not change significantly from land, with the remaining classes increasing slightly, especially the proportion of grade 9-10 wind speeds appearing from 0 to 0.01 from 2010 to 2017. The proportion of force 4-5 gales increased on construction land and decreased for the remaining types, while the proportion of force 8 gales and above declined significantly.
Atmosphere 2021, 12, x FOR PEER REVIEW 14 of 18 was the opposite of and does not change significantly from land, with the remaining classes increasing slightly, especially the proportion of grade 9-10 wind speeds appearing from 0 to 0.01 from 2010 to 2017. The proportion of force 4-5 gales increased on construction land and decreased for the remaining types, while the proportion of force 8 gales and above declined significantly.

Main Findings
Large-scale spatiotemporal variations in wind speed are perturbed by hourly and spatial-scale wind events in different underlayers, and there are pronounced regional differences in economic development and urbanization in the areas surrounding the Yangtze River inland waterway. For example, the rapid socio-economic development of the middownstream of the Yangtze River inland waterway and offshore areas can increase the uncertainty of wind speed frequency changes. We also found a significant underestimation of wind speeds in upstream areas in the reanalysis data, which may be a reflection of errors in the reanalysis data assimilation system itself. For example, the model may not

Main Findings
Large-scale spatiotemporal variations in wind speed are perturbed by hourly and spatial-scale wind events in different underlayers, and there are pronounced regional differences in economic development and urbanization in the areas surrounding the Yangtze River inland waterway. For example, the rapid socio-economic development of the middownstream of the Yangtze River inland waterway and offshore areas can increase the uncertainty of wind speed frequency changes. We also found a significant underestimation of wind speeds in upstream areas in the reanalysis data, which may be a reflection of errors in the reanalysis data assimilation system itself. For example, the model may not correctly describe local topography, climate, and surface parameters; in addition, there may be missing meteorological data and low-resolution input data. The reanalysis data represent large-scale climate changes caused by atmospheric circulation, greenhouse gases, etc. However, this is not sensitive to changes in land use type, and the results of this study confirm this limitation. However, ERA-Interim reanalysis data can be significantly improved by the wind multiplier downscaling method, and this more accurate wind speed data can be used for refined analysis of the maximum wind speed characteristics of the waterway. The wind multiplier downscaling method better integrates the combined effects of land use and topography on wind speed, and the evaluation of statistical indicators shows that the downscaling products have significantly improved accuracy and can more accurately reflect the spatial and temporal distribution characteristics of maximum wind speed in the study area.
Land use has a greater impact on overland winds than overwater winds in the study area. The maximum wind speed measured by the weather stations tends to decrease with increasing roughness of the subsurface owing to the increase in the number of buildings along with accelerated urbanization; the roughness of the water is smaller and the increase in the number of buildings around river banks and seawalls tends to form a river valley effect, resulting in a continuous decrease in the wind speed difference between water and land. Correlation analysis between land use and maximum wind speed shows that, except for bare land, the proportions of wind speed over the different land types increase and then decrease with rising wind speed, especially for grade 9 wind speed or above; the proportion was smallest in 2017. Bare land shows a trend of increasing wind speed in grade 9-10 wind speed, while the proportion of wind speed over water grade 4-8 is the opposite of that on land. The proportion is opposite to that of land, with no significant change in trend, and the rest of the classes are slightly increased, especially grade 9-10, which showed an increasing trend from 2010 to 2017. This study also demonstrates that although the risk of strong winds over water has weakened, the risk of extreme wind speeds has increased, and so the risk of damage to navigation and the region from extreme wind speeds along the Yangtze River inland waterway under climate change continues to increase.

Limitations of Study Approach
Anthropogenic activity substantially affects surface roughness, which has become a major cause of variation in the maximum wind speed in the study area. However, there are uncertainties in the analysis of the effects of land use on the maximum wind speed. For example, land use can lead to an increase in canopy resistance, which has a blocking effect on surface maximum wind speed and can change the thermodynamic properties of the surface, which in turn causes changes in surface maximum wind speed by altering local circulation and turbulent fluxes. Moreover, there is a lack of reliable data methods and assessment techniques, and the development of relevant models should be enhanced. Many stations in the Yangtze River inland waterway area are located within the transition zone between urban and suburban areas, and urban expansion not only increases surface roughness, but also leads to heat island circulation and turbulent flux perturbations, which have an important impact on the monitoring results of local wind fields.
When analyzing the impact of land use on the maximum surface wind speed based on station observations, the above thermodynamic effects should be taken into account, and existing reanalysis data should be corrected for cross-scale thermodynamic information to improve the accuracy of the observations. Moreover, the spatiotemporal resolution of wind speed products should be further strengthened, and higher resolution wind speed downscaling techniques should be developed through integration with high precision and long time series of geographical elements (e.g., land use, impermeable layer). Finally, wind speed event changes are significantly influenced by climate change, and wind speed change predictions under different climate change scenarios should be further developed to obtain more complete spatiotemporal characteristics of wind speed and its potential change trends in the study area.

Conclusions
This study proposes an integrated analysis method for wind multiplier downscaling and accuracy assessment, and a comprehensive assessment of the acquired downscaled maximum wind speed spatiotemporal characteristics. The modeling approach is generally applicable to the study area, with good accuracy and spatiotemporal representation of downscaled products, and the main conclusions obtained are as follows. (1) The trends of the site downscaled V max_WM (5.9-17.55 m/s) and observed V max_OBS (6.61-15.44 m/s) are better fitted (especially 5-10 m/s), and the downscaled accuracy is significantly improved; the upstream statistic D drops from 0.941 to 0.294 and R 2 rises to 0.905, but there are extreme values leading to large differences in maximum wind speed and low concentration (MAE decreasing from 2.34 to 1.73); midstream D drops from 0.471 to 0.294 and has the highest concentration (MAE decreasing from 0.91 to 0.82), but the fit is average (R 2 rises to 0.494); downstream data concentration (MAE decreasing to 1) and fit (R 2 from −0.009 to 0.898) are better, but individual data are too different resulting in lower sample quality (D increasing to 0.412). (2) Downscaled products and interannual rates of variability at the stations (−0.01 and −0.014 m·s −1 ·10a −1 ) are much closer, with both showing increasing trends from 1980 to 1990 (0.91 and 0.44 m·s −1 ·10a −1 ), weakening by 2010 (−0.30 and −0.13 m·s −1 ·10a −1 ), followed by significant decreases (−0.71 and −1.07 m·s −1 ·10a −1 ). Wind speeds at both midand down-stream stations increased significantly from 1980 to 1990 and decreased significantly from 2010 to 2017. (3) V max_WM shows spatial characteristics of large nearshore wind speed (14-37 m/s), small inland wind speed (11-20 m/s), and wind speed over undulating topography. The V max_WM value in the hilly terrain section (Yibin and Chongqing) increases significantly and expands westward, the center of high-value region of the Yueyang and Jiujiang section narrows and decreases, and the center of the high-value region in the Nanjing and Shanghai section narrows to the northeast. The spatial distribution range of wind speed increased from 1980 to 1990, with V max_WM values reaching >20 m/s in the Shanghai section. The maximum wind speed in the downstream region decreased below 16 m/s from 1990 to 2000, with insignificant changes in the middle and upstream regions from 2000 to 2010. The overall changes were insignificant from 2010 to 2017. The most significant differences were found in the upstream region in 2017, with V max_WM values reaching 20 m/s or more, and a significant decrease in the mid and downstream regions. This suggests that the wind multiplier downscaling method effectively integrates the effects of surface processes on wind speed, which is missing in the reanalysis data, resulting in better interannual variability characteristics.
Research on extreme weather and climate events such as maximum wind speed can improve disaster prevention and management, and the design of shipping on the Yangtze River inland waterway. For the navigable environment of the Yangtze River inland waterway, the basin environment is complex, and the Yangtze River crosses three major geomorphic terrains with large topographic drop-offs and significant differences in climatic and meteorological conditions, with an increasing trend in the frequency of extreme events [37,38]. The natural rate of climate variability in the monsoon region has a significant impact on the Yangtze River inland waterway, such as the large-scale atmospheric circulation changes caused by the East Asian monsoon fluctuations, The El Niño-Southern Oscillation (ENSO) and the Pacific Ocean interdecadal oscillation [39]. It is worth noting that the wind speed changes in the Yangtze River may be a sensitive indicator of global climate change. It should also be noted that the Yangtze River inland waterway is highly exposed and vulnerable and is also affected by transport policy, and socio-economic and technological developments. Our study provides valuable data and knowledge support for navigational risk management, identification of local areas at high wind speed risk, and wind emergency management in the Yangtze River. Research on the trends of extreme weather events on the Yangtze River inland waterway under the combined influence of climate change and human activities should be strengthened in the future, and high-precision regional model data and spatial analysis methods should be used to clarify the different influencing factors of disaster-causing risks in each section of Yangtze River inland waterway.