**Aim.** The widespread use of global navigation satellite systems (GNSS) has led to the development of new methods designed to determine and accumulate the index of ionosphere ionization (VTEC). Using these data it is possible to significantly improve the accuracy and reliability of determining the coordinates of the observation point. Therefore, the task of constructing a model of ionization index is relevant. **Method. **To construct a spatial model, we used the spherical Legendre functions of the first kind of real order with integer degree as a basic system of functions. We found the magnitude of order using the Sturm-Liouville theory since it depended on the size of the investigated region. Such system of functions form two orthogonal systems of functions in the region under study (the sphere segment), but does not have recurrence relations between functions, therefore, it is necessary to use function expansion in a hypergeometric series to find them. Also in order to find unknown coefficients of the model it is necessary to use the Tikhonov regularization parameter, since the matrix of normal equations will not be stable. For calculating the time model of the ionosphere the coefficients of different spatial models were expanded in series of power polynomials. **Results.** Based on the data of the ionization parameter values obtained of 19 permanent stations of the ZAKPOS network using the Trimble Pivot Platform software, the spatial-temporal model of this parameter was constructed using the Legendre spherical functions up to the 3rd order as well as with power polynomials up to 3rd order. The standard deviation between the measured and model values of the VTEC parameter does not exceed 1*TECU*. **The scientific novelty and practical significance.** We developed algorithm for construction of the space-time model of the ionosphere parameter. A ionosphere model of high resolution is obtained, which can be used to solve geodetic tasks in order to provide the necessary accuracy in determining the coordinates of the point, as well as to study and forecast the space weather.

1. Abdelazeem M., Celik R., Rabbany A. EI. On the development of a regional ionospheric correction model for low-cost single frequency GNSS users. The 10th International conference on mobile mapping technology, Cairo, Egypt, 2017.

2. Dzhuman B. B. Aproksymatsiya anomaliy syly vahy metodom ASHA na terytoriyu Arktyky [Approxi¬mation of gravity anomalies by method of ASHA on Arctic area]. Geodesy, cartography and aerial photography, 2014, no. 80, pp. 62–68.

3. Dzhuman B. B. Pro pobudovu modeli lokal'noho hravitatsiynoho polya [On the constraction of local gravitational field model]. Geodynamics, no. 1(14), 2013, pp. 29–33.

4. Gao Y., Liu Z. Precise Ionosphere Modeling Using Regional GPS Network Data. Journal of Global Positioning Systems. 1, 2002, pp. 18–24.

https://doi.org/10.5081/jgps.1.1.18

5. Haines G. V. Computer programs for spherical cap harmonic analysis of potential and general felds. Haines. Comput. Geosci. 14, 1988, pp. 413–447.

https://doi.org/10.1016/0098-3004(88)90027-1

6. Haines G. V. Spherical cap harmonic analysis. J. Geophys. Res. 90, 1985, pp. 2583–2591.

https://doi.org/10.1029/JB090iB03p02583

7. Hobson E. W. The Theory of Spherical and Ellipsoidal Harmonics. New York: Cambridge Univ. Press, 1931, 476 p.

8. Kelvin L, Tait P. Treatise on natural philosophy. New York: Cambridge Univ. Press., 1896, 852 pp.

9. Liu J., R. Chen, J. An, Z. Wang and J. Hyyppa. Spherical cap harmonic analysis of the Arctic ionospheric TEC for one solar cycle. Journal of Geophysical Research, vol. 119, 2014, pp. 601–619.

https://doi.org/10.1002/2013JA019501

10. Marchenko A., Dzhuman B. Regional quasigeoid determination: an application to arctic gravity project. Geodynamics, no. 1(18), 2015, pp. 7–17.

11. Ohashi M., Y. Sato, A. Yamada, Y. Kubo and S. Sugi¬moto. Studies on Spherical Cap Harmonic Analysis for Japanese Regional Ionospheric Delays and its Prediction. Proceedings of the 47th ISCIE Inter¬national Symposium on Stochastic Systems Theory and Its Applications Honolulu, Dec. 5–8, 2015

12. Schaer S. Mapping and predicting the earth's ionosphere using the global positioning system. PhD thesis, Astronomical Institute, University of Berne, Switzerland, 1999, 205 p.

13. Schmidt M., M. Fengler, T. Mayer-Gurr, A. Eicker, J. Kusche, L. Sanchez, S. C. Han. Regional gravity modeling in terms of spherical base functions. J. Geod. V. 81, 2007, pp. 17–38.

https://doi.org/10.1007/s00190-006-0101-5

14. Tykhonov A. N., Arsenyn V.Ya. Metody reshenyya ne¬korrekt¬nykh zadach [Methods for solving incorrect tasks]: 2-nd edition. Nauka. Hlavnaya redaktsyya fyzyko-matematycheskoy lyteratury, 1979.

15. Yankiv-Vitkovska L. M., Savchuk S., Pauchok V. The determination and procedure transformation of the ionosphere parameters with GNSS-observations. Geodesy, cartography and aerial photography, 2015, no. 82, pp. 5–12.

16. Yankiv-Vitkovska L. M., Savchuk S. H., Pauchok V. K., Matviichuk Ya. M. and Bodnar D. I. Recovery of the Spatial State of the Ionosphere Using Regular Definitions of the TEC Identifier at the Network of Continuously Operating GNSS Stations of Ukraine. Journal of Geodesy and Geomatics Engineering, (2016) 1–9. З.37–48.

17. Yankiv-Vitkovska L. M. Metodyka vyznachennya parametriv ionosfery u merezhi suputnykovykh stantsiy zakhidnoyi Ukrayiny [A procedure for the determination of ionosphere parameters on the basis of the gnss network in western Ukraine]. Kosmichna nauka i tekhnolohiya. 2013, 19, no. 6, pp. 4–-52.

18. Yankiv-Vitkovska L. M. Metodyka userednennya danykh dlya pobudovy rehional'noyi modeli ionosfery [The technique of averaging data for construction of the regional ionosphere model]. Geodesy, cartography and aerial photography, 2014, no. 79, pp. 35–41.