Near-surface wind field characteristics of the desert-oasis transition zone in Dunhuang,China

The desert-oasis transition zone(DOTZ)serves as a buffer area between the desert and oasis.Understanding its wind field characteristics is of great significance for the prevention and control of aeolian disasters in the oasis.In this study,we used meteorological data during 2013–2019 from the portable meteorological stations at five sites(site A on the edge of the oasis,sites B,C,and D in the DOTZ,and site O in the desert)in Dunhuang,China to analyze the near-surface wind field characteristics and their causes,as well as to reveal the key role of the DOTZ in oasis protection.The results showed that the mean wind speed,frequency of sand-driving wind,and directional variability of wind decreased from west to east within the DOTZ,and wind speed was significantly affected by air temperature.The terrain influenced the prevailing winds in the region,mainly from northeast and southwest.Only some areas adjacent to the oasis were controlled by southeasterly wind.This indicated that the near-surface wind field characteristics of the DOTZ were caused by the combined effects of local terrain and surface hydrothermal difference.At site D,the annual drift potential(DP)was 24.95 vector units(VU),indicating a low wind energy environment,and the resultant drift direction(RDD)showed obvious seasonal differences.Additionally,the DOTZ played an important buffering role between the desert and oasis.Compared with the desert,the mean wind speed in the oasis decreased by 64.98%,and the prevailing wind direction was more concentrated.The results of this study will be useful in interpreting the aeolian activity of the DOTZ in Dunhuang.

VALIDATION OF NEAR-SURFACE WINDS OBTAINED BY A HYBRID WRF/CALMET MODELING SYSTEM OVER A COASTAL ISLAND WITH COMPLEX TERRAIN

The results from a hybrid approach that combines a mesoscale meteorological model with a diagnostic model to produce high-resolution wind fields in complex coastal topography are evaluated.The diagnostic wind model(California Meteorological Model,CALMET) with 100-m horizontal spacing was driven with outputs from the Weather Research and Forecasting(WRF) model to obtain near-surface winds for the 1-year period from 12 September 2003 to 11 September 2004.Results were compared with wind observations at four sites.Traditional statistical scores,including correlation coefficients,standard deviations(SDs) and mean absolute errors(MAEs),indicate that the wind estimates from the WRF/CALMET modeling system are produced reasonably well.The correlation coefficients are relatively large,ranging from 0.5 to 0.7 for the zonal wind component and from 0.75 to 0.85 for the meridional wind component.MAEs for wind speed range from 1.5 to 2.0 m s-1 at 10 meters above ground level(AGL) and from 2.0 to 2.5 m s-1 at 60 m AGL.MAEs for wind direction range from 30 to 40 degrees at both levels.A spectral decomposition of the time series of wind speed shows positive impacts of CALMET in improving the mesoscale winds.Moreover,combining the CALMET model with WRF significantly improves the spatial variability of the simulated wind fields.It can be concluded that the WRF/CALMET modeling system is capable of providing a detailed near-surface wind field,but the physics in the diagnostic CALMET model needs to be further improved.

Temporal-spatial cross-correlation analysis of non-stationary near-surface wind speed time series

Temporal-spatial cross-correlation analysis of non-stationary wind speed time series plays a crucial role in wind field reconstruction as well as in wind pattern recognition.Firstly,the near-surface wind speed time series recorded at different locations are studied using the detrended fluctuation analysis(DFA),and the corresponding scaling exponents are larger than 1.This indicates that all these wind speed time series have non-stationary characteristics.Secondly,concerning this special feature( i.e.,non-stationarity)of wind signals,a cross-correlation analysis method,namely detrended cross-correlation analysis(DCCA) coefficient,is employed to evaluate the temporal-spatial cross-correlations between non-stationary time series of different anemometer pairs.Finally,experiments on ten wind speed data synchronously collected by the ten anemometers with equidistant arrangement illustrate that the method of DCCA cross-correlation coefficient can accurately analyze full-scale temporal-spatial cross-correlation between non-stationary time series and also can easily identify the seasonal component,while three traditional cross-correlation techniques(i.e.,Pearson coefficient,cross-correlation function,and DCCA method) cannot give us these information directly.

Linking moisture and near-surface wind with winter temperature to reveal the Holocene climate evolution in arid Xinjiang region of China

An increasing number of palaeo-climatic records have been reported to identify the Holocene climate history in the arid Xinjiang region of northwest China. However, few studies have fully considered the internal linkages within the regional climate system, which may limit our understanding of the forcing mechanisms of Holocene climate change in this region. Here, we systematically consider three major issues of the moisture/precipitation, temperature and near-surface wind relevant to the Holocene climate history of Xinjiang. First, despite there still has debated for the Holocene moisture evolution in this region, more climatic reconstructions from lake sediments, loess, sand-dunes and peats support a long-term regional wetting trend. Second, temperature records from ice cores, peats and stalagmites demonstrate a long-term winter warming trend during the Holocene in middle-to high-latitudes of Asia. Third, recent studies of aeolian sedimentary sequences reveal that the near-surface winds in winter gradually weakened during the Holocene, whereas the winter mid-latitude Westerlies strengthened in the Tienshan Mountains. Based on this evidence, in the arid Xinjiang region we propose an early to middle Holocene relatively cold and dry interval, with strong near-surface winds;and a warmer, wetter interval with weaker near-surface winds in the middle to late Holocene during winter. Additionally,we develop a conceptual model to explain the pattern of Holocene climate changes in this region.From the early to the late Holocene, the increasing atmospheric COcontent and winter insolation,and the shrinking of high-latitude continental ice-sheets, resulted in increasing winter temperatures in middle to high latitudes in the Northern Hemisphere. Subsequently, the increased winter temperature strengthened the winter mid-latitude Westerlies and weakened the Siberian high-pressure system,which caused an increase in winter precipitation and a decrease in near-surface wind strength. This scenario is strongly supported by evidence from geological records, climate simulation results, and modern reanalysis data. Our hypothesis highlights the important contribution of winter temperature in driving the Holocene climatic evolution of the arid Xinjiang region, and it implies that the socio-economic development and water resources security of this region will face serious challenges presented by the increasing winter temperature in the future.

Research on entropy weight variation evaluation method for wind power clusters based on dynamic layered sorting

This paper presents an evaluation method for the entropy-weighting of wind power clusters that comprehensively evaluates the allocation problems of wind power clusters by considering the correlation between indicators and the dynamic performance of weight changes.A dynamic layered sorting allocation method is also proposed.The proposed evaluation method considers the power-limiting degree of the last cycle,the adjustment margin,and volatility.It uses the theory of weight variation to update the entropy weight coefficients of each indicator in real time,and then performs a fuzzy evaluation based on the membership function to obtain intuitive comprehensive evaluation results.A case study of a large-scale wind power base in Northwest China was conducted.The proposed evaluation method is compared with fixed-weight entropy and principal component analysis methods.The results show that the three scoring trends are the same,and that the proposed evaluation method is closer to the average level of the latter two,demonstrating higher accuracy.The proposed allocation method can reduce the number of adjustments made to wind farms,which is significant for the allocation and evaluation of wind power clusters.

The effects of wind-driven waves and ocean spray on the drag coefficient and near-surface wind profiles over the ocean

By introducing a wave-induced component and a spray-induced component to the total stress, a mathematical model based on the Ekman theory is proposed to detail the influence of wind-driven waves and ocean spray on the momentum transport in a marine atmosphere boundary layer（MABL）. An analytic solution of the modified Ekman model can be obtained. The effect of the wave-induced stress is evaluated by a wind wave spectrum and a wave growth rate. It is found that the wave-induced stress and spray stress have a small impact compared with the turbulent stress on the drag coefficient and the wind profiles for low-to-medium wind speed. The spray contribution to the surface stress should be much more taken into account than the winddriven waves when the wind speed reaches above 25 m/s through the action of a ＂spray stress＂. As a result, the drag coefficient starts to decrease with increasing wind speed for high wind speed. The effects of the winddriven waves and spray droplets on the near-surface wind profiles are illustrated for different wave ages, which indicates that the production of the spray droplets leads the wind velocity to increase in the MABL. The solutions are also compared with the existed field observational data. Illustrative examples and the comparisons between field observations and the theoretical solutions demonstrate that the spray stress has more significant effect on the marine atmosphere boundary layer in the condition of the high wind speed compared with wave-induced stress.

Field Observations of Near-Surface Wind Flow Across Expressway Embankment on the Qinghai–Tibet Plateau

Crushed rock layers(CRLs),ventilation ducts(VDs)and thermosyphons are air-cooling structures(ACSs)widely used for maintaining the long-term stability of engineered infrastructures in permafrost environments.These ACSs can effectively cool and maintain the permafrost subgrade’s frozen state under climate warming by facilitating heat exchange with ambient air in cold seasons.As convection is a crucial working mechanism of these ACSs,it is imperative to understand the near-surface wind flow(NSWF)across a constructed infrastructure,such as an embankment.This article describes a yearlong field observation of the NSWF across an experimental expressway embankment,the first of its kind on the Qinghai–Tibet Plateau(QTP).The wind speed and direction along a transect perpendicular to the embankment on both the windward and leeward sides and at four different heights above the ground surface were collected and analyzed.The results showed that the embankment has a considerable impact on the NSWF speed within a distance of up to ten times its height,and in the direction on the leeward side.A power law can well describe the speed profiles of NSWF across the embankment,with the power-law indices(PLIs)varying from 0.14 to 0.40.On an annual basis,the fitted NSWF PLI far away from the embankment was 0.19,which differs substantially from the values widely used in previous thermal performance evaluations of ACSs on the QTP.Finally,the significance of the NSWF to the thermal performance of the ACSs,particularly the CRLs and VDs,in linear transportation infrastructure is discussed.It is concluded that underestimating the PLI and neglecting wind direction variations may lead to unconservative designs of the ACSs.The results reported in this study can provide valuable guidance for infrastructure engineering on the QTP and other similar permafrost regions.

Aeolian dust transportation and deposition by near-surface winds in arid and semiarid China

Dustfall collections were carried out in April and May 2001 and in March 2002 at six sites in northern China.Our results showed that the total deposition of dust fractions 【250 μm in diameter and the deposition of Fe both decreased exponentially with increasing distances from the source areas,and that the half-attenuation distance (HAD) for dust deposition was about 229 km in this re-gion.The HAD was closely related to the grain-size distribution of the dust,and the 15 to 20 μm fractions had the longest HAD.However,the fractions 【15 μm in diameter can be easily adsorbed to coarse particles and deposited after only short distances,and the HAD for the fractions 15 to 100 μm in diameter showed a power relationship with the grain-size distribution.The HAD for Fe deposition was 233 km,which was a little longer than that of total dust deposition,which suggests that the Fe content is higher in fine particles than in coarse particles,as previous studies have suggested.In addition,our analysis showed that under the control of current climatic conditions,the coarse fractions in dust derived from northwestern China cannot be transported over long distances,instead,it is transported primarily by near-surface winds (【3 km above the ground).The Fe in aeolian dust generated from arid and semiarid regions of China and deposited in the North Pacific region is usually transported by the upper westerlies.

Near-Surface Wind-Speed Stilling in Alaska during 1984-2016 and Its Impact on the Sustainability of Wind Power

Based on wind-speed records of Alaska’s 19 first-order weather stations, we analyzed the near-surface wind-speed stilling for January 1, 1984 to December 31, 2016. With exception of Big Delta that indicates an increase of 0.0157 m·s–1·a–1, on average, all other first-order weather stations show declining trends in the near-surface wind speeds. In most cases, the average trends are less then?–0.0300?m·s–1·a–1. The strongest average trend of?–0.0500?m·s–1·a–1 occurred at Homer, followed by?–0.0492?m·s–1·a–1 at Bettles, and?–0.0453?m·s–1·a–1 at Yakutat, while the declining trend at Barrow is marginal. The impact of the near-surface wind-speed stilling on the wind-power potential expressed by the wind-power density was predicted and compared with the wind-power classification of the National Renewable Energy Laboratory and the Alaska Energy Authority. This wind-power potential is, however, of subordinate importance because wind turbines only extract a fraction of the kinetic energy from the wind field characterized by the power efficiency. Since wind turbine technology has notably improved during the past 35 years, we hypothetically used seven currently available wind turbines of different rated power and three different shear exponents to assess the wind-power sustainability under changing wind regimes. The shear exponents 1/10, 1/7, and 1/5 served to examine the range of wind power for various conditions of thermal stratification. Based on our analysis for January 1, 1984 to December 31, 2016, Cold Bay, St. Paul Island, Kotzebue, and Bethel would be very good candidates for wind farms. To quantify the impact of a changing wind regime on wind-power sustainability, we predicted wind power for the periods January 1, 1984 to December 31, 1994 and January 1, 2006 to December 31, 2016 as well. Besides Big Delta that suggests an increase in wind power of up to 12% for 1/7, predicted wind power decreased at all sites with the highest decline at Annette (≈38%), Kodiak (≈30%), King Salmon (≈26%), and Kotzebue (≈24%), where the effect of the shear exponents was marginal. Bethel (up to 20%) and Cold Bay (up to 14%) also show remarkable decreases in predicted wind power.

Analysis of Near-Surface Layer Meteorological Elements for a Summer Fog in Ninglang Airport

Based on hourly observation data from 2016 to 2019 and hourly automatic meteorological observation data from 2018 to 2019, we analyzed the annual, monthly and daily variation characteristics of fog days in Ninglang Airport, and diagnosed a case of summer fog. The results show that the average annual fog days in Ninglang Airport are 128 days and the annual characteristic of fog days is bimodal. Analyzed on fog cases in 2019, it was a summer fog process lasting more than 20 hours and visibility less than 500 m. The data from the automatic observational stations around Ninglang airport can describe the near-surface layer features of fog and illustrate the progress of fog. The meteorological elements of automatic observational stations changed suddenly 2 hours before fog, so the data of the automatic observational stations can be used to forscast fog.

Analysis and Simulation of the Stratospheric Quasi-zero Wind Layer over Korla, Xinjiang Province, China

The stratospheric quasi-zero wind layer (QZWL) is a transition region with low zonal wind speeds in the lower stratosphere at an altitude of ~20 km. The zonal wind direction above the QZWL layer is opposite to that below the QZWL layer and the north –south wind component is small. The atmospheric wind field near the stratospheric QZWL is an important factor affecting the flight altitude and dynamic control of stratospheric airships. It is therefore necessary to study the stratospheric QZWL to provide better environmental information for these aircraft. High-resolution radiosonde data were used to analyze the characteristics of the stratospheric QZWL over Korla, Xinjiang Province, China. A weak wind layer in which the wind direction suddenly reversed from westerly to easterly was observed at ~20 km in the lower stratosphere, characteristic of the stratospheric QZWL. The Weather Research and Forecasting model was used to simulate the profiles of the horizontal wind speed and direction over Korla. The forcing effect of each diagnostic term in the equation on the zonal wind speed was analyzed. The results showed that the advection term was the dominant factor forcing the zonal wind speed. The wave term had a secondary forcing role, although the forcing effect of the wave term on the zonal wind speed was significant in some regions.

Impacts of a wind stress and a buoyancy flux on the seasonal variation of mixing layer depth in the South China Sea

The seasonal variation of mixing layer depth （MLD） in the ocean is determined by a wind stress and a buoy- ance flux. A South China Sea （SCS） ocean data assimilation system is used to analyze the seasonal cycle of its MLD. It is found that the variability of MLD in the SCS is shallow in summer and deep in winter, as is the case in general. Owing to local atmosphere forcing and ocean dynamics, the seasonal variability shows a regional characteristic in the SCS. In the northern SCS, the MLD is shallow in summer and deep in winter, affected coherently by the wind stress and the buoyance flux. The variation of MLD in the west is close to that in the central SCS, influenced by the advection of strong western boundary currents. The eastern SCS presents an annual cycle, which is deep in summer and shallow in winter, primarily impacted by a heat flux on the air-sea interface. So regional characteristic needs to be cared in the analysis about the MLD of SCS.

The research on boundary layer evolution characteristics of Typhoon Usagi based on observations by wind profilers

Vertically exploring the characteristics of the typhoon boundary layer（TBL） plays an important role in recognizing typhoon structure. The boundary layer radial direction and tangential wind characteristics of Typhoon Usagi based on the observational data of three boundary layer wind profiler stations along the route of Typhoon Usagi（No. 1319） and by combining with sounding data. The results show that：（1） maximum tangential wind appears in the vicinity of the eye area of Usagi, and it basically maintains a height of around 1 800 m when Usagi keeps a strong typhoon level, with the rapidly decreasing strength of Usagi after it lands, the speed of the maximum tangential wind and its vertical range both decrease;（2） the height of the maximum tangential wind is close to that of the inflow layer top of the typhoon, and is greater than that of the boundary layer estimated on the basis of Richardson number or potential temperature gradient, while the height of mixed layer judged on the basis of the signal-to-noise ratio（SNR） or its gradient is usually low;（3） the the boundary layer height can reach higher than2 100 m before Usagi lands. When the typhoon level or above is achieved, the boundary layer height observed by various stations does not change much, basically staying at between 1 200 and 1 600 m. With the decreasing strength of Usagi after its landfall, the boundary layer height rapidly drops.

Structures and Characteristics of the Windy Atmospheric Boundary Layer in the South China Sea Region during Cold Surges

An observational analysis of the structures and characteristics of a windy atmospheric boundary layer during a cold air outbreak in the South China Sea region is reported in this paper. It is found that the main structures and characteristics are the same as during strong wind episodes with cold air outbreaks on land. The high frequency turbulent fluctuations （period 〈 1 min） are nearly random and isotropic with weak coherency, but the gusty wind disturbances （1 rain〈period 〈 10 min） are anisotropic with rather strong coherency. However, in the windy atmospheric boundary layer at sea, compared with that over land, there are some pronounced differences： （1） the average horizontal speed is almost independent of height, and the vertical velocity is positive in the lower marine atmospheric boundary layer; （2） the vertical flux of horizontal momentum is nearly independent of height in the low layer indicating the existence of a constant flux layer, unlike during strong wind over the land surface; （3） the kinetic energy and friction velocity of turbulent fluctuations are larger than those of gusty disturbances; （4） due to the independence of horizontal speed to height, the horizontal speed itself （not its vertical gradient used over the land surface） can be used as the key parameter to parameterize the turbulent and gusty characteristics with high accuracy.

Gustiness and coherent structure under weak wind period in atmospheric boundary layer

Statistical analysis of turbulent and gusty characteristics in the atmospheric boundary layer under weak wind period has been carried out.The data used in the analysis were from the multilevel ultrasonic anemometer-thermometers at 47 m,120 m,and 280 m levels on Beijing 325 m meteorological tower.The time series of 3D atmospheric velocity were analyzed by using conventional Fourier spectral analysis and decompose into three parts:basic mean flow(period > 10 min),gusty disturbances(1 min < period < 10 min)and turbulence fluctuations(period < 1 min).The results show that under weak mean wind condition:1)the gusty disturbances are the most strong fluctuations,contribute about 60% kinetic energy of eddy kinetic energy and 80% downward flux of momentum,although both the eddy kinetic energy and momentum transport are small in comparison with those in strong mean wind condition;2)the gusty wind disturbances are anisotropic;3)the gusty wind disturbances have obviously coherent structure,and their horizontal and vertical component are negatively correlated and make downward transport of momentum more effectively;4)the friction velocities related to turbulence and gusty wind are approximately constant with height in the surface layer.

Non-Gaussian Lagrangian Stochastic Model for Wind Field Simulation in the Surface Layer

Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.

Robust Solution for Boundary Layer Height Detections with Coherent Doppler Wind Lidar

Although coherent Doppler wind lidar(CDWL)is promising in detecting boundary layer height(BLH),differences between BLH results are observed when different CDWL measurements are used as tracers.Here,a robust solution for BLH detections with CDWL is proposed and demonstrated:mixed layer height(MLH)is retrieved best from turbulent kinetic energy dissipation rate(TKEDR),while stable boundary layer height(SBLH)and residual layer height(RLH)can be retrieved from carrier-to-noise ratio(CNR).To study the cause of the BLH differences,an intercomparison experiment is designed with two identical CDWLs,where only one is equipped with a stability control subsystem.During the experiment,it is found that the CNR could be distorted by instrument instability because the coupling efficiency from free-space to the polarization-maintaining fiber of the telescope is sensitive to the surrounding environment.In the ML,a bias up to 2.13 km of the MLH from CNR is found,which is caused by the CNR deviation.In contrast,the MLH from TKEDR is robust as long as the accuracy of wind is guaranteed.In the SBL(RL),the CNR is found capable to retrieve SBLH and RLH simultaneously and robustly.This solution is tested during an observation period over one month.Statistical analysis shows that the root-mean-square errors(RMSE)in the MLH,SBLH,and RLH are 0.28 km,0.23 km,and 0.24 km,respectively.

Boundary-Layer Wind Structure in a Landfalling Tropical Cyclone

In this study, a slab boundary layer model with a constant depth is used to analyze the boundary-layer wind structure in a landfalling tropical cyclone. Asymmetry is found in both the tangential and radial components of horizontal wind in the tropical cyclone boundary layer at landfall. For a steady tropical cyclone on a straight coastline at landfall, the magnitude of the radial component is greater in the offshoreflow side and the tangential component is greater over the sea, slightly offshore, therefore the greater total wind speed occurs in the offshore-flow side over the sea. The budget analysis suggests that： （1） a greater surface friction over land produces a greater inflow and the nonlinear effect advects the maximum inflow downstream, and （2） a smaller surface friction over the sea makes the decrease of the tangential wind component less than that over land. Moreover, the boundary layer wind structures in a tropical cyclone are related to the locations of the tropical cyclone relative to the coastline due to the different surface frictions. During tropical cyclone landfall, the impact of rough terrain on the cyclone increases, so the magnitude of the radial component of wind speed increases in the offshore-flow side and the tangential component outside the radius of maximum wind speed decreases gradually.

Wind Structure in an Intermediate Boundary Layer Model Based on Ekman Momentum Approximation

A quasi three–dimensional, intermediate planetary boundary layer (PBL) model is developed by including inertial acceleration with the Ekman momentum approximation, but a nonlinear eddy viscosity based on Blackadar’s scheme was included to improve the theoretical model proposed by Tan and Wu (1993). The model could keep the same complexity as the classical Ekman model in numerical, but extends the conventional Ekman model to include the horizontal accelerated flow with the Ekman momentum approximation. A comparison between this modified Ekman model and other simplified accelerating PBL models is made. Results show that the Ekman model overestimates (underestimates) the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow due to the neglect of the acceleration flow, however, the semi–geostrophic Ekman model overestimates the acceleration effects resulting from the underestimating (overestimating) of the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow. The Ekman momentum approximation boundary layer model could be applied to the baroclinic atmosphere. The baroclinic Ekman momentum approximation boundary layer solution has both features of classical baroclinic Ekman layer and the Ekman momentum approximate boundary lager.

Horizontal structure of convergent wind shear associated with sporadic E layers over East Asia

At present,the main detection instruments for observing sporadic E(Es)layers are ground-based radars,dense networks of ground-based global navigation satellite system(GNSS)receivers,and GNSS radio occultation,but they cannot capture the whole picture of the horizontal structure of Es layers.This study employs the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension model(WACCM-X 2.1)to derive the horizontal structure of the ion convergence region(HSICR)to explore the shapes of the large-scale Es layers over East Asia for the period from June 1 to August 31,2008.The simulation produced the various shapes of the HSICRs elongated in the northwest-southeast,northeast-southwest,or composed of individual small patches.The close connection between Es layer critical frequency(foEs)and vertical ion convergence indicates that the HSICR is a good candidate for revealing and explaining the horizontal structure of the large-scale Es layers.