Spatial analysis of provision the population of Kharkiv with civil defense facilities

Keywords: civil defense structures, cluster analysis, semivariogram, spatial autocorrelation, spatial interpolation, population density


The article analyzed the territorial aspect of the distribution and capacity characteristics of civil defense structures in the city of Kharkiv, comparing it with the population size. In the first stage, population density was calculated for specific areas within the maximum permissible radius of the nearest available shelter. The spatial characteristics of the location (level of clustering) and capacity of the defense structures were examined for these same areas. Finally, the population size was correlated with the capacity of the civil defense structures, and spatial clustering of the identified areas was conducted based on this attribute.

The geostatistical method of spatial interpolation was used to determine population density in specific areas and fill gaps in the primary data. The application of this method required the following sequential procedures: transforming the primary data according to a normal distribution, constructing a semivariogram model of the transformed variables, aggregating the model into a surface, and defining target polygons. During the calculation of the territorial provision of various types of shelters in the city of Kharkiv, the following sources were utilized: a layer of point objects from the interactive map of open data on the Kharkiv Geoportal (to determine the locations and types of protective structures), and information from the website of the non-commercial enterprise «Emergency Medical Care and Disaster Medicine Center» of the Kharkiv Regional Council (to update information on the capacity of storage facilities and anti-radiation shelters). The Global Moran’s Index and Local Moran’s Index are statistical methods used to assess spatial autocorrelation, which is the degree of clustering or spatial pattern in a variable across a defined area. In this case, they were employed to evaluate the spatial autocorrelation of the capacity of civil defense structures in specific parcels in Kharkiv. The Global Moran’s Index indicated a high level of clustering of areas based on this attribute. Using the Local Moran’s Index, parcels were classified into five object classes: High-High cluster (HH), Low-Low cluster (LL), a high-value outlier surrounded by low-values (HL), a low-value outlier surrounded by high-values (LH), and areas without cluster or outliers (non-significant).

The scientific novelty of the article, compared to related studies on a similar topic, lies in the utilization and transformation of a hexagonal grid of population density distribution in the city of Kharkiv, in accordance with the research requirements. In the conclusions, based on the results of cluster analysis and through the adaptation of the tectological principle of the weakest link to the realities of the present, a comprehensive sequential strategy for addressing the shortcomings of territorial provision of civil defense structures in ensuring the population of Kharkiv was proposed.


Download data is not yet available.

Author Biographies

Kateryna Sehida, V.N. Karazin Kharkiv National University, Svobody Sq., 4, Kharkiv, 61022, Ukraine

DSc (Geography), Professor of the Department of Human Geography and Regional Studies

Serhii Chekhov, V.N. Karazin Kharkiv National University, Svobody Sq., 4, Kharkiv, 61022, Ukraine

PhD Student of the Department of Human Geography and Regional Studies


Budynky i sporudy. Zakhysni sporudy tsyvilnoi oborony: DBN V 2.2.5-97 [Buildings and Structures. Civil Defense Protective Structures: DBN V 2.2.5-97]. [Chynnyi vid 1998-01-01]. Kyiv, 1998. 22 s. (Informatsiia ta dokumentatsiia) [in Ukrainian].

Kodeks tsyvilnoho zakhystu Ukrainy: Zakon Ukrainy vid 2 zhovtnia 2012 roku № 5403-VI [Civil Protection Code of Ukraine: Law of Ukraine dated October 2, 2012, No. 5403-VI]. Retrieved from [in Ukrainian].

Marchenko, O.M., & Tomanevych, L.M. Teoriia orhanizatsii: navch. posib [Theory of Organization: a textbook]. Lviv: LvDUVS, 2015. 359 s [in Ukrainian].

Nabir danykh pro shchilnist naselennia [Dataset on population density]. Kontur. Retrieved from [in English].

Skilky Kharkiv platyt za bezpechni zupynky [How much does Kharkiv pay for safe bus stops]. Kharkiv Today – informatsiino-analitychnyi portal Kharkova. Retrieved from [in Ukranian].

Spysok bomboskhovyshch u Kharkivskii oblasti [List of bomb shelters in Kharkiv region]. Komunalne nekomertsiine pidpryiemstvo Kharkivskoi oblasnoi rady «Tsentr ekstrenoi medychnoi dopomohy ta medytsyny katastrof». Retrieved from [in Ukrainian].

Ukryttia dlia naselennia [Shelters for the population]. Heoportal Kharkova: heoprostorovi dani infrastruktury mista. Retrieved from [in Ukranian].

Ding, Z., Dong, H., Li, Y., Yao, X., & Zhang, H. (2022). Study on the Emergency Shelter Spatial Accessibility Based on Simulation of Evacuation Process. Proceedings of the 2022 3rd International Conference on Modern Education and Information Management (ICMEIM 2022), 10-20.

Ding, Z., Dong, H., Yang, L., Xue, N., He, L., & Yao, X. (2022). Study on the Emergency Shelter Spatial Accessibility Based on the Adaptive Catchment Size 2SFCA Method. ISPRS Int. J. Geo-Inf, 11, 593.

Dou, K., & Zhan, Q. (2011). Accessibility analysis of urban emergency shelters: Comparing gravity model and space syntax. International Conference on Remote Sensing, Environment and Transportation Engineering, Nanjing, China, 5681-5684,

Dou, K., Zhan, Q., & Shiguo L. (2012). GIS-based responsibility area subdivision for metropolitan emergency shelters – Case study of Wuchang district, Wuhan city. The 6th International Association for China Planning (IACP) Conference, Wuhan, China, pp. 1-4.

How Cluster and Outlier Analysis (Anselin Local Moran’s I) works. Esri (Environmental Systems Research Institute). Retrieved from

How Kriging works. Esri (Environmental Systems Research Institute). Retrieved from

How Spatial Autocorrelation (Global Moran’s I) works. Esri (Environmental Systems Research Institute). Retrieved from

Jekl, J., & Jánský, J. (2022). Security Challenges and Economic-Geographical Metrics for Analyzing Safety to Achieve Sustainable Protection. Sustainability, 14, 15161.

Liang, Y., Xie, Zh., Chen, S., Xu, Y., Xin, Z., Yang, S., Jian, H., & Wang, Q. (2023). Spatial Accessibility of Urban Emergency Shelters Based on Ga2SFCA and Its Improved Method: A Case Study of Kunming, China. Journal of Urban Planning and Development, 149 (2).

Rizqihandari, N., & Indratmoko, S. (2017). Using Open Street Map Data for Population Distribution Model. Advances in Social Science, Education and Humanities Research, 79.

Shi, Y., Zhai, G., Xu, L., Zhu, Q., & Deng, J. (2019). Planning Emergency Shelters for Urban Disasters: A Multi-Level Location-Allocation Modeling Approach. Sustainability, 11(16), 1-19.

Su, H., Chen, W., & Wang, Z. (2020). Evaluating the Crowdedness of Urban Emergency Shelters Based on the Improved Gravity Model. IOP Conf. Ser. Earth Environ. Sci, 502, 012046.

Tang, S., Wang, J., Xu, Y., Chen, S., Zhang, J., Zhao, W., & Wang, G (2023). Evaluation of Emergency Shelter Service Functions and Optimisation Suggestions – Case Study in the Songyuan City Central Area. Sustainability, 15(9), 7283.

Transform Field (Data Management) Esri (Environmental Systems Research Institute). Retrieved from

Unal M., & Uslu C. (2016). GIS-Based accessibility analysis of urban emergency shelters: The case of Adana city. ISPRS – International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, XLII-2/W1: 95-101.

What is areal interpolation? Esri (Environmental Systems Research Institute). Retrieved from

Zhu, W., Xing, H., & Kang, W. (2023). Spatial Layout Planning of Urban Emergency Shelter Based on Sustainable Disaster Reduction. Int. J. Environ. Res. Public Health, 20, 2127.

How to Cite
Sehida, K., & Chekhov, S. (2023). Spatial analysis of provision the population of Kharkiv with civil defense facilities. Human Geography Journal, 34, 14-26.
Наукові повідомлення