Study on the Susceptibility of Drifting Snow in Ya’an–Qamdo Section of the Railway in Southwest China

: To investigate the susceptibility of drifting snow along the Ya’an–Qamdo section of the railway, which is located in a high-altitude and cold plateau in Southwest China with scarce meteorological information, the Weather Research and Forecasting Model (WRF) is used in this paper to simulate the spatio-temporal distribution of meteorological data. According to the varying terrain, the railway section from Ya’an to Qamdo is divided into two regions along 100.8 ◦ E for double-layer nested simulation. The original land use data of the WRF model are used in region 1. Due to the increased number of mountains in region 2, the original data are replaced by the MCD12Q1v006 land use data, and the vertical direction layers are densified near the ground to increase simulation accuracy. The simulated results are compared with the observation data. It is found that after densifi-cation, the results have been significantly improved. The results obtained by the WRF model can accurately simulate the change trends of temperature, rainfall, and wind speed, and the correlation coefficients are relatively high, which verifies the accuracy of WRF for simulating complex terrain regions. The simulation results further indicate that approximately 300 km of the Ya’an–Qamdo railway may experience drifting snow. Among them, no drifting snow events occur in Ya’an County, and the areas with higher probability are located at the border between Luding County and Tianquan County, followed by Kangding area. The remaining areas have a probability of less than 10%. The WRF model demonstrates its capability in the drifting snow protection of railways with limited meteorological data


Introduction
Drifting snow is a theoretical study that integrates geography, meteorology, fluid mechanics, glaciology, and other disciplines.It is common in winter and spring seasons, and mostly occurs in Xinjiang, Inner Mongolia, Tibet, and Northeast China.It can cause snow accumulation on roads and block traffic.In order to reduce the harm of drifting snow, many scholars have studied its physical properties and prevention measures.
In terms of the fundamental physical properties of drifting snow, the two-phase flow theory proposed by Bagnold [1] provides a theoretical foundation for the study of wind and snow.Owen [2] obtained the calculation formula for the height of the saltation layer through a theoretical analysis of snow particle movement and field measurements of drifting snow.Kind [3] derived equations for calculating the concentration profile of snow particles by incorporating the theory of turbulent diffusion.Budd [4] analyzed the probability distribution of the equivalent particle size of snow particles in drifting snow and found that both gamma distribution and log-normal distribution could fit the effective particle size distribution of snow particles well.Schmidt [5] proposed a critical wind speed formula based on average snow particle size using field measurement data.Kind et al. [6] discovered that using high-density particles to simulate the salting process of snow particles can better fulfill the dynamic similarity requirements compared to using low-density particles.This is because high-density particles exhibit closer resemblance to actual snow particles in terms of snow pile shape and transport rate.These observations and empirical models provide valuable references for selecting key parameters in numerical simulations and wind tunnel experiments.
In terms of drifting snow control measures, Flaga et al. [7] conducted wind tunnel experiments to study the snow particle deposition of three different forms of large sports stadium roofs and found that the shape parameters of buildings have an important influence on the design of snow loads at different locations.Therefore, shape factors cannot be ignored in building load design.Qiang et al. [8] used artificial snow particles to simulate snowfall and snow drifting on flat roofs when they occur simultaneously in a lowtemperature wind tunnel.They observed that the development of snow transport before reaching a saturated state followed a consistent pattern, regardless of whether there was concurrent snowfall or not.However, in cases where snowfall occurred concurrently with drifting, uneroded regions could be observed on the roof.Liu et al. [9] proposed equipment named "snow-wind combined experimental facility" to investigate snowdrifts around buildings and snow loads on building roofs; the experimental results prove the feasibility of this method, and the method can also be adopted on various shapes of building roofs.The wind tunnel tests using simulated material or moving snow into a low-temperature wind tunnel are not capable of accurately replicating the genuine properties and conditions of snow particles.Therefore, Li et al. [10] studied the snow accumulation in railway cuttings using a movable direct current low-temperature wind tunnel set up in the field.By adjusting the wind tunnel to simulate natural flow fields, they compared their experimental results with the detection results to verify the rationality of wind tunnel design and the accuracy of the experimental results.Due to the high cost of building low-temperature wind tunnels and the need for a large amount of complex preparation work to be carried out before experiments, some scholars have begun to use numerical simulations to study them as computing power improves.Previous researchers used computational fluid dynamics (CFDs) to study the effect of fence height, porosity, bottom clearance, and inclination angle on the snow accumulation around snow fences and their ability to prevent snow, which helps optimize their design [11][12][13].Ma et al. [14] examined the effect of different forms of windbreaks on typical roadbed snow accumulation from aspects such as flow field, snow phase concentration, and sedimentation length.The results showed that the inclination angle of the windbreak was a key parameter affecting its effectiveness.Forward-leaning guide plates were more effective in preventing drifting snow during windy and snowy weather, and the angle between guide plates and the vertical axis should be controlled at around 45 • .
Meteorological factors such as temperature, precipitation, wind speed, and wind direction have an important impact on the occurrence of drifting snow.The threshold windspeed of dry snow is generally proportional to the temperature, with lower temperatures resulting in lower starting wind speeds [15].However, when the temperature drops below −27.27 • C, the snow particles start to deteriorate, leading to an increase in the required starting wind speed [16].Temperature is also a critical parameter affecting the shape of snow particles; −5 • C is the critical temperature at which snow crystals transition from plate-like to columnar [17].Snow particles of different shapes undergo distinct development processes in a snowstorm flow.The number of suspended snow particles and the single-width snow transport rate of spherical, ellipsoidal, and cylindrical snow particles are two orders of magnitude higher than those of star-shaped and hexagonal prism-shaped snow particles [18].At a height of 10 m, a wind speed of 3-8 m/s is sufficient to cause loose snow particles to jump and move [19].When the angle between the line and the wind direction is greater than 55 • , the most snow resistance is generated, and the closer the angle is to 90 • , the greater the probability of snow resistance and the more severe the consequences [20].These meteorological factors are interconnected and have various impacts on drifting snow.
The above research mainly focuses on the occurrence and aftermath of drifting snow.However, evaluating the susceptibility before its occurrence serves as the foundation for further research.Railway lines should avoid areas with high susceptibility or set different protection measures according to different occurrence probabilities along the way to save economic costs and improve protection efficiency.Most railways pass through high-altitude and cold regions with complex terrain and scarce meteorological information, making it difficult to judge the susceptibility of drifting snow along the railway line.The Weather Research and Forecasting Model (WRF) has been widely applied in regional simulations [21][22][23] and can partially address this issue.The model can decompose the simulated region into multiple grids with an accuracy of several kilometers.By providing boundary conditions, a series of atmospheric motion equations, including Newton's second law, the first law of thermodynamics, the continuity equation, and the state equation, can be solved to achieve weather forecasting.Therefore, the WRF model can provide meteorological data with a high resolution and continuous long time series at any location in the simulation area.Yang [24] analyzed the wind speed and wind direction characteristics during a storm and snow disaster in the Mayitas traffic corridor using the WRF dataset.In a separate study, Qi [25] employed WRF to investigate route selection schemes for the Keta Railway.Additionally, Luo et al. [26] utilized WRF to analyze snowfall, wind speed, and direction along the entire Keta Railway line and derived the probability of drifting snow occurrence for the entire route.Currently, the application of WRF is limited to the study of drifting snow on the Keta Railway, and there is a need to extend this research to include the study of drifting snow on the southwest railway, which holds significant importance.
Regarding the railway in Southwest China, Ma et al. [27] studied the distribution characteristics of extreme summer precipitation and the variation mechanism of total precipitation of the railway using daily rainfall observations and ERA5 reanalysis data in summer.The results showed that the extreme summer precipitation in these areas accounted for 30% of the total summer precipitation, and the extreme precipitation in the central and western regions was lower than that in the east, but with higher precipitation frequency.Due to the scarcity of observational data, some scholars use satellite data to study the temperature and precipitation changes in the region where the railway is located [28,29].The WRF model can provide a more detailed spatiotemporal distribution of meteorological information.Gao et al. [30] evaluated the WRF model by examining the interannual variation, spatial structure, and 33-year time variation trend of surface temperature and precipitation.They concluded that the WRF model is capable of accurately analyzing significant weather phenomena on the Qinghai-Tibet Plateau.Chen et al. [31], Minder et al. [32], and Wrzesien et al. [33] have all verified the feasibility of using the WRF model to simulate precipitation in a complex topography area, and many scholars have used the WRF model to calculate southwest China [34][35][36].Therefore, the above studies have demonstrated that applying the WRF model to the southwest region has credibility and can partially compensate for the lack of meteorological data.
Existing studies on Southwest China based on the WRF model mainly focus on climate, but there are few applications in terms of engineering in Southwest China.Due to the complex terrain and wide elevation difference in Southwest China, it is very difficult to conduct an engineering field investigation.Therefore, the application of the WRF model to combine meteorology and engineering can produce a large achievement with a low cost.In order to provide an example for the study of the susceptibility of drifting snow in areas lacking measured meteorological information, this paper takes the railway section from Ya'an to Qamdo in Southwest China as an object and studies its drifting snow susceptibility based on the WRF results.This research method can provide reference for places with the same terrain and lack of meteorological data and render guidance for railway drifting snow protection in southwest China.The structure of the remainder of this present study is as follows: Section 2 describes the geographical environment of a region which southwest railway passes through, the configuration of the WRF model, and the verification of the simulation results.Section 3 discusses the calculated results of temperature, precipitation, wind speed, wind direction, and the susceptibility of drifting snow on the railway.The conclusion of this study is finally summarized in Section 4.

Study Area
The southwest railway, as shown in Figure 1, is situated on the margin of Qinghai-Tibet Plateau and crosses various geomorphic units, such as a high mountain and gorge area, western Sichuan high mountain and plain area, and the core area of the Hengduan Mountains high mountain and gorge.The section starts from Ya'an and passes through Kangding, Litang, Dege, and other places before reaching Qamdo.The length of the railway is about 700 km, and the altitude varies greatly along the way, as shown in Figure 2 (the elevation data used come from https://www.earthdata.nasa.gov/(accessed on 2 June 2023)).The complex terrain results in climatic features along the railway, and the spatial distribution of meteorological elements is significantly different.
verification of the simulation results.Section 3 discusses the calculated results of temperature, precipitation, wind speed, wind direction, and the susceptibility of drifting snow on the railway.The conclusion of this study is finally summarized in Section 4.

Study Area
The southwest railway, as shown in Figure 1, is situated on the margin of Qinghai-Tibet Plateau and crosses various geomorphic units, such as a high mountain and gorge area, western Sichuan high mountain and plain area, and the core area of the Hengduan Mountains high mountain and gorge.The section starts from Ya'an and passes through Kangding, Litang, Dege, and other places before reaching Qamdo.The length of the railway is about 700 km, and the altitude varies greatly along the way, as shown in Figure 2 (the elevation data used come from https://www.earthdata.nasa.gov/(accessed on 2 June 2023)).The complex terrain results in climatic features along the railway, and the spatial distribution of meteorological elements is significantly different.
Ya'an has a subtropical monsoon humid climate, with clear seasons, hot and rainy summers, and mild and humid winters.Tianquan area also has a subtropical monsoon climate, with warm winters and hot summers, and rain and heat occur simultaneously.Kangding to Xinduqiao area is divided into two parts by Zheduoshan Mountain.The eastern part has a warm climate with clear seasons, while the western part has cold and dry winters.Xinduqiao to Baiyu section is influenced by southwest monsoons and has cool and low-temperature summers and cold and dry winters.In Litang, Dege, and other places, the temperature is low and the annual precipitation is high and mainly in solid form.verification of the simulation results.Section 3 discusses the calculated results of temperature, precipitation, wind speed, wind direction, and the susceptibility of drifting snow on the railway.The conclusion of this study is finally summarized in Section 4.

Study Area
The southwest railway, as shown in Figure 1, is situated on the margin of Qinghai-Tibet Plateau and crosses various geomorphic units, such as a high mountain and gorge area, western Sichuan high mountain and plain area, and the core area of the Hengduan Mountains high mountain and gorge.The section starts from Ya'an and passes through Kangding, Litang, Dege, and other places before reaching Qamdo.The length of the railway is about 700 km, and the altitude varies greatly along the way, as shown in Figure 2 (the elevation data used come from https://www.earthdata.nasa.gov/(accessed on 2 June 2023)).The complex terrain results in climatic features along the railway, and the spatial distribution of meteorological elements is significantly different.
Ya'an has a subtropical monsoon humid climate, with clear seasons, hot and rainy summers, and mild and humid winters.Tianquan area also has a subtropical monsoon climate, with warm winters and hot summers, and rain and heat occur simultaneously.Kangding to Xinduqiao area is divided into two parts by Zheduoshan Mountain.The eastern part has a warm climate with clear seasons, while the western part has cold and dry winters.Xinduqiao to Baiyu section is influenced by southwest monsoons and has cool and low-temperature summers and cold and dry winters.In Litang, Dege, and other places, the temperature is low and the annual precipitation is high and mainly in solid form.Ya'an has a subtropical monsoon humid climate, with clear seasons, hot and rainy summers, and mild and humid winters.Tianquan area also has a subtropical monsoon climate, with warm winters and hot summers, and rain and heat occur simultaneously.Kangding to Xinduqiao area is divided into two parts by Zheduoshan Mountain.The eastern part has a warm climate with clear seasons, while the western part has cold and dry winters.Xinduqiao to Baiyu section is influenced by southwest monsoons and has cool and low-temperature summers and cold and dry winters.In Litang, Dege, and other places, the temperature is low and the annual precipitation is high and mainly in solid form.

Equations of WRF
In this paper, WRF (ARW) version 4.1.2is employed, and a terrain following vertical coordinates is adopted, which is expressed as where p h is the static part of the air pressure, µ = p hs − p ht , and p hs and p ht are the air pressure at the model surface and upper boundary, respectively.
Considering the horizontal pressure gradient error, the new variables are defined as the deviation from the reference state, and the atmospheric motion variables are expressed as follows p = p(z) where p is the corrected air pressure, ϕ is the geopotential height, ϕ = gz, g is the gravitational acceleration, and z is the height of the position.α is the derivative of density, which is a non-conservative variable, and is a function of the vertical profile (x, y, η), µ d = p dhs − p dht , d represents dry air, and p ′ , ϕ ′ , α ′ , µ ′ d are the original values of those variables without considering the horizontal pressure gradient error.
In the calculation process, the perturbation form of momentum equations of the WRF model can be written as: And the potential temperature equation can be expressed as: where U, V, and W are the momentum components controlling the three directions of the equation of atmospheric motion on the projection plane, respectively.Ω is the momentum of the vertical velocity field, Θ is the momentum of the potential temperature field, u, v, and w are velocity components in three directions, respectively, mx and my are magnification factors for map projections in x and y directions, respectively, m x , m y = (∆x,∆y) (∆x e ,∆y e ) , and F U , F V , F W , F Θ are forced terms caused by model physical processes, turbulent mixing, spherical projection, and Earth rotation, respectively.
The mass conservation equation and the potential height equation are: Appl.Sci.2024, 14, 475 The static equation in the perturbation form is: The diagnostic equation for complete air pressure (water vapor plus dry air) is: where p 0 is the reference pressure (usually 105 Pa), R d is the ratio constant of dry air to gas, and γ is the ratio of heat capacity c p of dry air at constant pressure to heat capacity c v of passenger.

Model Configuration
Ya'an, with a subtropical monsoon climate, is located in the western edge of Sichuan Basin.The area from Kangding to Litang belongs to the western Sichuan high mountain plain, with a wide and flat or hilly top surface, and mainly has a continental plateau-type and mountain-type climate.On the other hand, the area from Litang to Changdu belongs to the southeastern Tibetan Hengduan Mountains area, with high mountains, deep valleys, complex terrain, and mainly plateau climate.Since the terrain of Ya'an and Kangding is flatter and the climate is warmer than that of Litang Changdu, in order to adapt these differences, this paper divides the railway section from Kangding to Changdu into two regions with 100.8 • E as the boundary and performs double-layer nested simulation.The outer grid is marked as d01, and the inner grid is marked as d02.Region 1 mainly includes Ya'an to Kangding and other places, and region 2 includes Litang to Qamdo and other places.The center points of the two regions are 30.03• N, 101.97 • E and 30.6 • N, 98.62 • E, respectively.Figure 3 shows the specific simulation range.
The National Centers for Environmental Prediction (NCEP) Final Operational Global Analysis (FNL) (https://rda.ucar.edu/datasets/ds083.2/ (accessed on 13 February 2023)) with a horizontal resolution of 1 • × 1 • and a frequency of 6 h are selected as the initial conditions of the model.The simulation covers the period from 12 December 2016 to 28 February 2017, and the specific grid parameters of the two regions are shown in Table 1.In region 2, which has many mountains, the simulation accuracy is improved by using MCD12Q1v006 land use data with a horizontal resolution of 500 m instead of the original data.Moreover, to simulate the atmospheric processes within the boundary layer more precisely, the vertical direction is divided into 41 layers, with 15 layers densely distributed at a height of about 1.5 km near the ground.Table 2 shows the heights of each layer grid, and Figure 4 illustrates the elevated layers in both regions, the red line shows the position of the 15th layer.The National Centers for Environmental Prediction (NCEP) Final Operational Global Analysis (FNL) (https://rda.ucar.edu/datasets/ds083.2/ (accessed on 13 February 2023)) with a horizontal resolution of 1° × 1° and a frequency of 6h are selected as the initial conditions of the model.The simulation covers the period from 12 December 2016 to 28 February 2017, and the specific grid parameters of the two regions are shown in Table 1.In region 2, which has many mountains, the simulation accuracy is improved by using MCD12Q1v006 land use data with a horizontal resolution of 500m instead of the original data.Moreover, to simulate the atmospheric processes within the boundary layer more precisely, the vertical direction is divided into 41 layers, with 15 layers densely distributed at a height of about 1.5 km near the ground.Table 2 shows the heights of each layer grid, and Figure 4 illustrates the elevated layers in both regions, the red line shows the position of the 15th layer.Physical processes are crucial for determining the source and redistribution of momentum, heat, and water vapor in the atmospheric motion.These processes have a great impact on the overall dynamics of the atmosphere.Microphysics schemes mainly affect the precipitation simulation, among which the Purdue-Lin scheme and Thompson scheme are two relatively complex calculation schemes.These schemes involve variables such as rain, ice, snow, and graupel, making them suitable for high-resolution simulations.The planetary boundary layer scheme affects the wind simulation, among which the MRF scheme can provide ideal simulation results for the boundary layer conditions with good mixing, such as instability or convection, and the YSU scheme can simulate the turbulent motion and energy conversion in the atmosphere well and has wide applicability.A radiative parameterization scheme is used to calculate the surface heat budget.The land surface process is used to calculate the heat and water vapor fluxes at land-sea grids.The surface layer parameterization scheme is used to describe the surface characteristics and simulate the interaction between the surface and the atmosphere.Since the grid resolution in this paper is 9 km, the cumulus parameter scheme is not disabled.It is found that the simulation results of rainfall and wind speed using the Purdue-Lin scheme [37] and YSU scheme [38] are closer to the observed results, respectively.According to the characteristics of various schemes, this study combines the Purdue-Lin and Thompson microphysics schemes, as well as the YSU and MRF boundary layer schemes, which have significant impacts on rainfall and wind speed results, as shown in Table 3   Physical processes are crucial for determining the source and redistribution of momentum, heat, and water vapor in the atmospheric motion.These processes have a great impact on the overall dynamics of the atmosphere.Microphysics schemes mainly affect the precipitation simulation, among which the Purdue-Lin scheme and Thompson scheme are two relatively complex calculation schemes.These schemes involve variables such as rain, ice, snow, and graupel, making them suitable for high-resolution simulations.The planetary boundary layer scheme affects the wind simulation, among which the MRF scheme can provide ideal simulation results for the boundary layer conditions with good mixing, such as instability or convection, and the YSU scheme can simulate the turbulent motion and energy conversion in the atmosphere well and has wide applicability.A radiative parameterization scheme is used to calculate the surface heat budget.The land surface process is used to calculate the heat and water vapor fluxes at land-sea grids.The surface layer parameterization scheme is used to describe the surface characteristics and simulate the interaction between the surface and the atmosphere.Since the grid resolution in this paper is 9 km, the cumulus parameter scheme is not disabled.It is found that the simulation results of rainfall and wind speed using the Purdue-Lin scheme [37] and YSU scheme [38] are closer to the observed results, respectively.According to the characteristics of various schemes, this study combines the Purdue-Lin and Thompson microphysics schemes, as well as the YSU and MRF boundary layer schemes, which have significant impacts on rainfall and wind speed results, as shown in Table 3.

Evaluation Statistics
To verify the accuracy of the simulation results, the Pearson correlation coefficient (R) is used to measure the linear correlation degree between the simulation sequence and the observation sequence.The closer R is to 1, the higher the correlation degree.The mean bias error (MBE) is used to verify the deviation between the simulation sequence and the observation sequence.The closer the result is to 0, the smaller the deviation.The root mean square error (RMSE) is used to verify the average deviation between the simulation sequence and the observation sequence.The closer the result is to 0, the better the simulation effect.The relevant calculation formulas are as follows: Appl.Sci.2024, 14, 475 where X i and X 0 are the simulated and observed value, X i and X 0 are the simulated and observed average values, and n is the number of records.

Verify Performance for WRF Model
The source of all the meteorological data used in this paper is from the China Meteorological Administration (https://data.cma.cn/(accessed on 12 December 2022)).These kinds of data have advantages of high accuracy and stability, but the data of each station only represent the surrounding information.Taking the daily average wind speed from December 2016 to February 2017 as the verification data, the physical parameter scheme of case3 in Table 3, the original land use data of the WRF model, and the two vertical stratification methods mentioned above are used to calculate the meteorological information of Litang.Table 4 shows the evaluation index results of the two schemes.The scheme of vertical layers densified near the ground has increased the correlation coefficient by 77.72%, reduced the mean bias by 61.13%, and reduced the root mean square error by 59.04% compared with the non-density scheme.Figure 5 shows the comparison of the change trends of the simulation results and the measured results.It can be seen that the change trend of the calculation results with vertical layers densified near the ground is more consistent with the measured results.Therefore, near-surface encryption can significantly improve the simulation results of the mountain wind field.In order to validate the computed results of four cases presented in Table 3, the daily average wind speed results and precipitation results from December 2016 to February 2017 at Kangding were used for comparison.Table 5 shows the correlation coefficients of the two data.For rainfall, the Purdue-Lin scheme demonstrated a higher correlation coefficient compared to the Thompson scheme, while for wind speed, the YSU scheme outperformed the MRF scheme.Figure 6 displays the actual simulated and observed values.All four schemes accurately captured the changing trends, but the maximum and minimum wind speeds were higher than the actual results.Considering the comparison of correlation coefficients and the actual trend results, case3, shown in Table 3, appeared to In order to validate the computed results of four cases presented in Table 3, the daily average wind speed results and precipitation results from December 2016 to February 2017 at Kangding were used for comparison.Table 5 shows the correlation coefficients of the two data.For rainfall, the Purdue-Lin scheme demonstrated a higher correlation coefficient compared to the Thompson scheme, while for wind speed, the YSU scheme outperformed the MRF scheme.Figure 6 displays the actual simulated and observed values.All four schemes accurately captured the changing trends, but the maximum and minimum wind speeds were higher than the actual results.Considering the comparison of correlation coefficients and the actual trend results, case3, shown in Table 3, appeared to be more reasonable.Therefore, this study adopted the third scheme for calculations in both region 1 and region 2.  To verify the accuracy of the calculation results of case3 comprehensively, this paper further validates the results of temperature, wind speed, and rainfall calculated by the WRF model.The simulated data are the WRF calculation results under the parameter settings in Tables 1 and 2 and case3 of Table 3.The validation data used to verify the calculation results of temperature and wind speed are the ground-based meteorological observation data of Kangding and Litang stations.Because there was no precipitation in the actual situation of the two meteorological stations in the observation data of 2022, the feasibility verification of the WRF model for rainfall used the daily average observation data of Litang from December 2016 to February 2017.The comparison results are shown in Table 6, and it can be found that the correlation coefficient of temperature is above 0.80, and the correlation coefficients of wind speed and rainfall are also above 0.50.In addition, the average deviation and root mean square error values are relatively small.By comparing the correlation coefficients presented in this paper with those published in the literature for temperature [39], wind speed [40], and rainfall [41] in similar terrain, it can be observed in Table 6 that the calculated correlation coefficients for temperature and wind speed in this study are slightly higher than those reported in previous studies.The correlation coefficients for rainfall closely align with the published values.Thus, the calculation results obtained in this study can be utilized to assess the susceptibility of drifting snow, providing valuable insights.To verify the accuracy of the calculation results of case3 comprehensively, this paper further validates the results of temperature, wind speed, and rainfall calculated by the WRF model.The simulated data are the WRF calculation results under the parameter settings in Tables 1 and 2 and case3 of Table 3.The validation data used to verify the calculation results of temperature and wind speed are the ground-based meteorological observation data of Kangding and Litang stations.Because there was no precipitation in the actual situation of the two meteorological stations in the observation data of 2022, the feasibility verification of the WRF model for rainfall used the daily average observation data of Litang from December 2016 to February 2017.The comparison results are shown in Table 6, and it can be found that the correlation coefficient of temperature is above 0.80, and the correlation coefficients of wind speed and rainfall are also above 0.50.In addition, the average deviation and root mean square error values are relatively small.By comparing the correlation coefficients presented in this paper with those published in the literature for temperature [39], wind speed [40], and rainfall [41] in similar terrain, it can be observed in Table 6 that the calculated correlation coefficients for temperature and wind speed in this study are slightly higher than those reported in previous studies.The correlation coefficients for rainfall closely align with the published values.Thus, the calculation results obtained in this study can be utilized to assess the susceptibility of drifting snow, providing valuable insights.We performed linear fitting on the calculation results.The closer the points in the figure are to the diagonal line (the black dash line in figure), the closer the simulation results are to the real results.The brighter the color of the points, the more points fall in that area.Figure 7 shows the fitting results of temperature and wind speed calculation data and their changes over time.It can be found that the model can capture the temperature changes well.The simulated temperature of Kangding station is concentrated within 0 to 5 • C, and the temperature of Litang station is concentrated within −5 • C to 5 • C. The WRF simulation results of wind speed for both places are concentrated in 1-2 m/s.The simulation results of Litang station are higher than the actual recorded values; however, the WRF model accurately captured the change trend and extreme values of wind speed at both stations.Figure 8 shows the comparison between the rainfall calculation results and the actual results.The WRF model successfully captured the rainfall events.In summary, although there is some error between the data simulated by the WRF model and the actual data, the correlation of the simulation results is good, and it can accurately capture the change trend of the data, which also verifies the effectiveness of the simulation method.

Occurrence Probability Determination Method
Vionnet et al. [42] analyzed the meteorological data detected when drifting snow occurred near the Alps and found that the main meteorological feature was that snow occurred when the wind speed reached a certain range, and this model was also applied by Guyomarc'h et al. [43].Based on this discovery and previous studies, the drifting snow in this paper is defined as periods with precipitating snow with a wind speed in the range of 3-8 m/s [19].In order to ensure that the precipitation is solid, the temperature is limited to less than 0 • C. Therefore, we calculated the probability by the percentage in the number of times that met the three conditions during the simulation period.

Occurrence Probability Determination Method
Vionnet et al. [42] analyzed the meteorological data detected when drifting snow occurred near the Alps and found that the main meteorological feature was that snow occurred when the wind speed reached a certain range, and this model was also applied by Guyomarc'h et al. [43].Based on this discovery and previous studies, the drifting snow in this paper is defined as periods with precipitating snow with a wind speed in the range of 3-8 m/s [19].In order to ensure that the precipitation is solid, the temperature is limited to less than 0 °C.Therefore, we calculated the probability by the percentage in the number of times that met the three conditions during the simulation period.

Occurrence Probability Determination Method
Vionnet et al. [42] analyzed the meteorological data detected when drifting snow occurred near the Alps and found that the main meteorological feature was that snow occurred when the wind speed reached a certain range, and this model was also applied by Guyomarc'h et al. [43].Based on this discovery and previous studies, the drifting snow in

Simulation Results of Temperature Distribution
The output results of the inner grid d02 are analyzed to obtain the monthly average temperature within the simulated area, as shown in Figure 9a-f.The results show that a minimum temperature of −14 • C in the study area is calculated by the WRF model, which ensures that the snow particles do not degrade due to excessively low temperatures.The temperature distribution around the railway is relatively regular.Region 1 bounded by Kangding is divided into two parts.The eastern side of Ya'an City has warmer temperatures than the western side, and the average monthly temperature stays above 0 • C, which matches the sub-tropical humid monsoon climate of Ya'an City that does not have harsh winters.On other hand, most areas on the western side of Kangding are colder, with an average monthly temperature of around −2 • C.These temperatures change with time, and the average temperature in January is lower than that in December and February.In region 2, along the railway from Litang to Qamdo, the temperatures are frigid, always below 0 • C. Notably, the temperature in January is significantly lower than that in December and February.

Simulation Results of Precipitation Distribution
Snow is the material condition for the occurrence of drifting snow.Sufficient snow sources will increase the probability of drifting snow, and the amount of snowfall will also affect the intensity and duration of drifting snow.Figure 10 illustrates the monthly cumulative precipitation around the railway in winter.In region 1, precipitation is mainly concentrated in the Ya'an area east of Kangding and the area south of Kangding station.Precipitation is affected by temperature, and the eastern area has higher temperatures, so it has more precipitation than other areas.Notably, January exhibits the lowest average temperature and, consequently, the lowest level of precipitation.In February 2017, there was more precipitation in Kangding area, and the cumulative precipitation reached 500 mm in the southern part of Kangding.Conversely, the precipitation in region 2 is significantly lower than that in region 1.During the simulation period, only a few places have weak precipitation, and the areas with higher precipitation are far away from the railway.Moreover, the precipitation was highest in February 2017, which is similar to the result of region 1.

Simulation Results of Wind Field Distribution
Wind is the primary factor that affects the occurrence of drifting snow.The wind direction determines the movement direction of snow particles, and the wind speed determines the scale of drifting snow.Figure 11 illustrates the monthly average wind speed and direction distribution in the simulated area.The arrows in the figures represent the direction of wind speed.In region 1, the wind speed on the west side of 102 • E is higher than that on the east side.In the three winter months, the wind speed in Kangding area is always higher than that in other surrounding areas.The average wind speed in January 2017 was higher than that in December 2016 and February 2017, with a maximum wind speed of 11.41 m/s.The wind direction around the railway is variable.In the section from Ya'an to Kangding, the wind direction is mainly eastward, parallel to the railway direction, with a small angle.On the west side of Kangding, the wind direction changes more frequently.Between 101.6 • E and 102 • E, it is mainly southwest wind.Between 101 • E and 101.5 • E, it is mainly south wind and southwest wind.In region 2, the whole area has high wind speed and higher than that in region 1, with a maximum wind speed of 12.33 m/s.The average wind speed in November is higher than that in the other two months.The wind direction is mainly southwest wind, forming a large angle with the railway.
matches the sub-tropical humid monsoon climate of Ya'an City that does not have harsh winters.On the other hand, most areas on the western side of Kangding are colder, with an average monthly temperature of around −2 °C.These temperatures change with time, and the average temperature in January is lower than that in December and February.In region 2, along the railway from Litang to Qamdo, the temperatures are frigid, always below 0 °C.Notably, the temperature January is significantly lower than that in December and February.
(a) (  precipitation than other areas.Notably, January exhibits the lowest average temperature and, consequently, the lowest level of precipitation.In February 2017, there was more precipitation in Kangding area, and the cumulative precipitation reached 500 mm in the southern part of Kangding.Conversely, the precipitation in region 2 is significantly lower than that in region 1.During the simulation period, only a few places have weak precipitation, and the areas with higher precipitation are far away from the railway.Moreover, the precipitation was highest in February 2017, which is similar to the result of region 1. (a) (   12.In region 1, the maximum wind speed at P1 is 13.10 m/s, the primary dominant wind direction is south-southwest wind (SSW), accounting for 17.92% of the entire simulation period, the secondary dominant wind direction is southwest wind (SW), accounting for 16.82% of the entire simulation period, the wind speed is mainly concentrated around 2-5 m/s, and the wind speed in the probability interval of drifting snow occurrence (3-8 m/s) accounts for 61.82% of the total wind speed in the simulation period; at P2, the maximum wind speed is 17.91 m/s, the primary dominant wind direction is southwest wind (SW), accounting for 41.51% of the entire simulation period, the secondary dominant wind direction is west wind (E), accounting for 28.54% of the entire simulation period, the wind speed is mainly concentrated around 5-9 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 56.39% of the total wind speed in the simulation period; at P3, the maximum wind speed is 16.7 m/s, the primary dominant wind direction is west-northwest wind (WNW), accounting for 22.51% of the entire simulation period, the secondary dominant wind directions are west wind and southwest wind, accounting for 18.81% and 18.39% of the entire simulation period, respectively, the wind speed is mainly concentrated around 6-10 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 50.73% of the total wind speed in the simulation period.

Simulation Results of Wind Field Distribution
mines the scale of drifting snow.Figure 11 illustrates the monthly average wind speed and direction distribution in the simulated area.The arrows in the figures represent the direction of wind speed.In region 1, the wind speed on the west side of 102° E is higher than that on the east side.In the three winter months, the wind speed in Kangding area is always higher than that in other surrounding areas.The average wind speed in January 2017 was higher than that in December 2016 and February 2017, with a maximum wind speed of 11.41 m/s.The wind direction around the railway is variable.In the section from Ya'an to Kangding, the wind direction is mainly eastward, parallel to the railway direction, with a small angle.On the west side of Kangding, the wind direction changes more frequently.Between 101.6°E and 102° E, it is mainly southwest wind.Between 101° E and 101.5°E, it is mainly south wind and southwest wind.In region 2, the whole area has high wind speed and higher than that in region 1, with a maximum wind speed of 12.33 m/s.The average wind speed in November is higher than that in the other two months.The wind direction is mainly southwest wind, forming a large angle with the railway.southwest wind (SW), accounting for 26.81% of the entire simulation period, the wind speed is mainly concentrated around 4-8 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 69.06% of the total wind speed in the simulation period.Since the climate of the study area is mainly influenced by the south branch jet stream of the westerly circulation, the wind direction is mostly west or southwest wind, which is consistent with the actual situation.Moreover, the proportion of wind speed in the interval of drifting snow occurrence is relatively large, which means a higher possibility of being affected by drifting snow.

Drifting Snow Susceptibility Judgment
As we can see in Figure 13, along the railway, the places with higher probability of drifting snow occurrence are located in Kangding County and Litang County, and there In region 2, the maximum wind speed at P4 is 12 m/s, the primary dominant wind direction is south-southwest wind (SSW), accounting for 22.84% of the entire simulation period, the secondary dominant wind direction is west wind (W), accounting for 20.17% of the entire simulation period, the wind speed is mainly concentrated around 2-4 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 55.59% of the total wind speed in the simulation period.At P5, the maximum wind speed is 14.90 m/s, the primary dominant wind direction is southwest wind (SW), accounting for 38.79% of the entire simulation period, the secondary dominant wind direction is southwest wind (SW), accounting for 26.81% of the entire simulation period, the wind speed is mainly concentrated around 4-8 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 69.06% of the total wind speed in the simulation period.Since the climate of the study area is mainly influenced by the south branch jet stream of the westerly circulation, the wind direction is mostly west or southwest wind, which is consistent with the actual situation.Moreover, the proportion of wind speed in the interval of drifting snow occurrence is relatively large, which means a higher possibility of being affected by drifting snow.

Drifting Snow Susceptibility Judgment
As we can see in Figure 13, along the railway, the places with higher probability of drifting snow occurrence are located in Kangding County and Litang County, and there is a long railway line passing through the drifting snow susceptible areas in Litang, Batang, and Baiyu counties, with a total length of about 321 km that may experience drifting snow along the line.The entire area of Ya'an City has almost zero probability of drifting snow occurrence; about 12 km of railway in the adjacent area of Luding County and Tianquan County may experience drifting snow, with a probability of 15%; in Kangding County, there are two sections of railway that may experience drifting snow phenomenon, located at 101.9  E has a probability of 6%.The railway located in Kangding County has a slightly higher probability of drifting snow phenomenon than other areas, and the length of the railway affected is larger, so the railway in this section needs to strengthen its protection against drifting snow.

Discussion of Occurrence Probability
As the WRF model relies on input boundary conditions, it can simulate diverse climatic scenarios by altering the input meteorological data to accommodate climate change.Climate variations, such as alterations in temperature and precipitation patterns, can exert influences on snow conditions.For instance, elevated temperatures can expedite snow melting rates, while shifts in precipitation can impact the formation and distribution of snow, ultimately affecting the occurrence of drifting snow.Given the WRF model's capacity to adjust to climate change, it can similarly accommodate shifts in drifting snow conditions due to climate variations.As a result, the derived probability of drifting snow occurrence based on the WRF model remains dependable.
This study utilizes temperature, snowfall, and wind speed information derived from WRF simulations to assess the probability of drifting snow around the railway.The calculated results have been validated against measured data and data from other published papers, which enhances the credibility of the findings in the paper.However, there are two errors in this paper: one is the calculation error of the WRF model itself, and the other is the error caused by the judgment model of drifting snow.If these two aspects are improved, the calculation results will be more accurate.
of about 4%; about 50 km of railway in Litang County at 99.8° E-100.2°E has a probability of 8%; about 87 km of railway passing through Batang County and Baiyu County at 99.15° E-99.7°E has a probability of 4% for drifting snow phenomenon; about 24 km of railway in Jomda County at 98.25° E-98.45°E has a probability of 6%; about 71 km of railway passing through Qamdo County and Chaya County at 97.50° E-97.90°E has a probability of 5%; about 12 km of railway in Chaya County and Basu County at 97.35° E-97.40°E has a probability of 6%.The railway located in Kangding County has a slightly higher probability of drifting snow phenomenon than other areas, and the length of the railway affected is larger, so the railway in this section needs to strengthen its protection against drifting snow.The accuracy of the WRF calculation results is mainly affected by boundary conditions, parameter scheme combinations, and model parameters.The initial conditions are reanalysis data, that is, the historical weather information obtained by integrating radar satellite, aircraft, ship, and station observation data and numerical forecasting products using data assimilation technology.Therefore, the accuracy of the data used for assimilation will affect the accuracy of the reanalysis data and then affect the calculation results of the WRF model.The combination of different parameter schemes of the WRF model will produce different calculation results, and one model parameter of WRF may affect multiple meteorological variables, so the selection of model parameters needs to consider multiple meteorological variables.However, WRF has many parameter schemes and model parameters.With the current computing power, it is difficult to obtain the optimal parameter scheme combination and model parameters through thousands of long-time calculations, which will also affect the calculation results of WRF.
Meteorological information such as wind and snow are the main conditions that affect drifting snow occurrence, but the terrain around the railway, vegetation types, and the form of the roadbed of the railway also have a certain impact on it.Since the purpose of this paper is mainly to calculate the probability of drifting snow from the perspective of meteorological information, snow and wind data are mainly taken into account in the judgment model.However, terrain and other factors are not considered, so the occurrence probability calculated based on the drifting snow judgment model in this paper may be higher than the actual result, and the calculation result in this paper is more secure.
To reduce errors of possibility caused by various factors, improvements can be made in those aspects, and the data of temperature, precipitation, and wind speed in FNL can be replaced with observation data of weather stations to improve the accuracy of boundary conditions and further improve the accuracy of WRF simulation.Sensitivity analysis of the parameter scheme and model parameters can be carried out to find the optimal scheme and parameter under the condition of saving computing resources as much as possible.In addition to the main factors such as wind and snow, secondary factors such as terrain and vegetation are taken into account to establish a more comprehensive evaluation model of drifting snow susceptibility.

Conclusions
The railway from Kangding to Qamdo is divided into two regions for simulation in this paper.In region 2 with complex terrain, the vertical direction near the ground is densified to 15 layers at a height of about 1.5 km, and the MCD12Q1v006 land use data are used to replace the original data.The susceptibility of drifting snow along the railway is judged according to the simulation results.The conclusions are as follows: (1) Among the simulations of temperature, rainfall, and wind speed, the simulation result of temperature is the best, followed by wind speed and rainfall.The WRF model can simulate the change trend of the three meteorological elements very accurately during the calculation period.(2) For mountainous areas with complex terrain, using the method of increasing the number of layers in the vertical direction and densifying near the ground can significantly improve the simulation result of wind speed by the WRF model.(3) The regional climate of the Ya'an to Qamdo railway is influenced by the south branch jet stream of the westerly circulation.The temperature is high in the east and low in the west, the rainfall is more in the east and less in the west, and the wind speed is high in the west and low in the east.The average wind speed is concentrated around 5-8 m/s, and the dominant wind direction in winter is mainly westerly and southwesterly.(4) The WRF model can be applied to determine the susceptibility of drifting snow in areas lacking meteorological information: the total length of Ya'an-Qamdo section is about 700 km, of which about 321 km of railway may experience drifting snow.The railway is more likely to be affected by drifting snow disaster when passing through Kangding County, Luding County, and Tianquan County where they border each other; although there are places with higher probability of drifting snow occurrence in the simulation area, these places are far away from the railway and will not affect it.
In the term of railway protection with high susceptibility, the meteorological information calculated by WRF can be extracted and taken as the boundary condition of computational fluid dynamics (CFDs) so as to more accurately study the snow condition of various roadbeds and the protective effect of various forms of snow fence to achieve the best protection of the railway.
The WRF model is easy to access, has a fast calculation rate, and provides high reliability of results.It can be applied to any location lacking but requiring long time series and high-resolution meteorological information, creating significant economic value at a relatively low cost.Therefore, this computational method has enormous potential for applications in engineering, industry, and other fields.

Figure 3 .
Figure 3. Simulation range for both regions.(a) The simulation range of region 1; (b) the simulation range of region 2.

Figure 3 .
Figure 3. Simulation range for both regions.(a) The simulation range of region 1; (b) the simulation range of region 2.

Figure 5 .
Figure 5. Wind speed variation trends at different ETA levels.

Figure 5 .
Figure 5. Wind speed variation trends at different ETA levels.

Figure 6 .
Figure 6.Comparison of simulation results in Kangding.(a) Comparison of precipitation result; (b) comparison of wind speed result.

Figure 6 .
Figure 6.Comparison of simulation results in Kangding.(a) Comparison of precipitation result; (b) comparison of wind speed result.

Figure 7 .
Figure 7.The results of the WRF fitting process and a comparison between the simulated and observed values.(a) The fitted line of temperature in Kangding; (b) the variation in temperature in Kangding; (c) the fitted line of temperature in Litang; (d) the variation in temperature in Litang; (e) the fitted line of wind speed in Kangding; (f) the variation in wind speed in Kangding; (g) the fitted line of wind speed in Litang; (h) the variation in wind speed in Litang.

Figure 7 .Figure 7 .
Figure 7.The results of the WRF fitting process and a comparison between the simulated and observed values.(a) The fitted line of temperature in Kangding; (b) the variation in temperature in Kangding; (c) the fitted line of temperature in Litang; (d) the variation in temperature in Litang; (e) the fitted line of wind speed in Kangding; (f) the variation in wind speed in Kangding; (g) the fitted line of wind speed in Litang; (h) the variation in wind speed in Litang.

Figure 9 .
Figure 9.The distribution of average temperature in the simulation period.(a) Temperature of region 1 in December; (b) temperature of region 2 in December; (c) temperature of region 1 in January; (d) temperature of region 2 in January; (e) temperature of region 1 in February; (f) temperature of region 2 in February.

Figure 9 .
Figure 9.The distribution of average temperature in the simulation period.(a) Temperature of region 1 in December; (b) temperature of region 2 in December; (c) temperature of region 1 in January; (d) temperature of region 2 in January; (e) temperature of region 1 in February; (f) temperature of region 2 in February.

Figure 10 .
Figure 10.The distribution of cumulative precipitation in the simulation period.(a) Precipitation of region 1 in December; (b) precipitation of region 2 in December; (c) precipitation of region 1 in January; (d) precipitation of region 2 in January; (e) precipitation of region 1 in February; (f) precipitation of region 2 in February.

Figure 10 .
Figure 10.The distribution of cumulative precipitation in the simulation period.(a) Precipitation of region 1 in December; (b) precipitation of region 2 in December; (c) precipitation of region 1 in January; (d) precipitation of region 2 in January; (e) precipitation of region 1 in February; (f) precipitation of region 2 in February.According to the different wind speeds and directions along the railway, five points, P1 (102 • E, 30.13 • N), P2 (101.6 • E, 29.98 • N), P3 (101.1 • E, 30.03 • N), P4 (99.50 • E, 30.37 • N), and P5 (97.7 • E, 31.15 • N), are selected as representatives to study the wind speed and direction distribution in their surrounding areas.The wind rose diagram and the wind speed probability density diagram are shown in Figure12.In region 1, the maximum wind speed at P1 is 13.10 m/s, the primary dominant wind direction is south-southwest wind (SSW), accounting for 17.92% of the entire simulation period, the secondary dominant wind direction is southwest wind (SW), accounting for 16.82% of the entire simulation period, the wind speed is mainly concentrated around 2-5 m/s, and the wind speed in the probability interval of drifting snow occurrence (3-8 m/s) accounts for 61.82% of the total wind speed in the simulation period; at P2, the maximum wind speed is 17.91 m/s, the primary dominant wind direction is southwest wind (SW), accounting for 41.51% of the entire simulation period, the secondary dominant wind direction is west wind (E), accounting for 28.54% of the entire simulation period, the wind speed is mainly concentrated around 5-9 m/s, and the wind speed in the probability interval of drifting snow occurrence accounts for 56.39% of

Figure 11 .
Figure 11.The distribution of average wind field in the simulation period.(a) Wind field of region 1 in December; (b) wind field of region 2 in December; (c) wind field of region 1 in January; (d) wind field of region 2 in January; (e) wind field of region 1 in February; (f) wind field of region 2 in February.

Figure 11 .
Figure 11.The distribution of average wind field in the simulation period.(a) Wind field of region 1 in December; (b) wind field of region 2 in December; (c) wind field of region 1 in January; (d) wind field of region 2 in January; (e) wind field of region 1 in February; (f) wind field of region 2 in February.

Figure 12 .
Figure 12.Wind rose diagram and the wind speed probability density graph.(a) the wind rose diagram of P1; (b) the wind speed probability density graph of P1; (c) the wind rose diagram of P2; (d) the wind speed probability density graph of P2; (e) the wind rose diagram of P3; (f) the wind speed probability density graph of P3; (g) the wind rose diagram of P4; (h) the wind speed probability density graph of P4; (i) the wind rose diagram of P5; (j) the wind speed probability density graph of P5.

Figure 12 .
Figure 12.Wind rose diagram and the wind speed probability density graph.(a) the wind rose diagram of P1; (b) the wind speed probability density graph of P1; (c) the wind rose diagram of P2; (d) the wind speed probability density graph of P2; (e) the wind rose diagram of P3; (f) the wind speed probability density graph of P3; (g) the wind rose diagram of P4; (h) the wind speed probability density graph of P4; (i) the wind rose diagram of P5; (j) the wind speed probability density graph of P5.

Figure 13 .
Figure 13.Probability of drifting snow along the railway.(a) The susceptibility of drifting snow in region 1; (b) the susceptibility of drifting snow in region 2.Figure 13.Probability of drifting snow along the railway.(a) The susceptibility of drifting snow in region 1; (b) the susceptibility of drifting snow in region 2.

Figure 13 .
Figure 13.Probability of drifting snow along the railway.(a) The susceptibility of drifting snow in region 1; (b) the susceptibility of drifting snow in region 2.Figure 13.Probability of drifting snow along the railway.(a) The susceptibility of drifting snow in region 1; (b) the susceptibility of drifting snow in region 2.

Table 1 .
The grid parameters of two regions.

Table 2 .
Vertical height of each mesh layer within 1.5 km.

Table 1 .
The grid parameters of two regions. .

Table 2 .
Vertical height of each mesh layer within 1.5 km.

Table 3 .
Schemes of physical parameters for WRF simulation.

Table 4 .
Evaluation results of different ETA levels.

Table 4 .
Evaluation results of different ETA levels.

Table 5 .
Comparison of results of different schemes.

Table 6 .
Evaluation results of temperature, wind speed, and rainfall.

Table 6 .
Evaluation results of temperature, wind speed, and rainfall.
• E-102 • E and 101.6 • E-101.8 • E, respectively, with lengths of 20 km and 23 km, respectively, and probabilities of 12% and 8%, respectively; about 22 km of railway in Yajiang County at 101.2 • E-101.4 • E may experience drifting snow, with a probability of about 4%; about 50 km of railway in Litang County at 99.8 • E-100.2 • E has a probability of 8%; about 87 km of railway passing through Batang County and Baiyu County at 99.15 • E-99.7 • E has a probability of 4% for drifting snow phenomenon; about 24 km of railway in Jomda County at 98.25 • E-98.45 • E has a probability of 6%; about 71 km of railway passing through Qamdo County and Chaya County at 97.50 • E-97.90 • E has a probability of 5%; about 12 km of railway in Chaya County and Basu County at 97.35 • E-97.40