This paper describes a study of the daily and yearly variability of one of the main characteristics of sporadic E layers (Es), the critical frequency (foEs). Our analysis is based on ionograms recorded by GIRO network ionosondes. The study estimates the spatial and temporal variability of the Es layer parameters and generates their statistics at different seasons and phases of the solar cycle. In turn, the statistics of the Es layer parameters can provide an assessment of the capability of the HF propagation forecasting. Further, maps of the distribution of Es layer critical frequency are produced. This paper implements the “cloud” model of Es layer in HF propagation model (Northern Ionosphere Model & Ray Tracing, NIM-RT) that can accurately reproduce many features observed in experimental measurements. Within this framework, a number of vertical sounding ionograms with the presence of Es layer are simulated based on the NIM-RT software.