=Paper=
{{Paper
|id=Vol-3006/68_regular_paper
|storemode=property
|title=Cartographic modeling of the temperature fields of the land fund of the Novosibirsk region using satellite data of the MODIS/Terra spectroradiometer
|pdfUrl=https://ceur-ws.org/Vol-3006/68_regular_paper.pdf
|volume=Vol-3006
|authors=Polina V. Voronina,Elena A. Mamash,Igor A. Pestunov,Svetlana Ya. Kudryashova,Aleksandr S. Chumbaev
}}
==Cartographic modeling of the temperature fields of the land fund of the Novosibirsk region using satellite data of the MODIS/Terra spectroradiometer==
Cartographic modeling of the temperature fields of the land fund of the Novosibirsk region using satellite data of the MODIS/Terra spectroradiometer Polina V. Voronina1,2 , Elena A. Mamash1 , Igor A. Pestunov1 , Svetlana Ya. Kudryashova3 and Aleksandr S. Chumbaev3 1 Federal Research Center for Information and Computational Technologies, Novosibirsk, Russia 2 Novosibirsk State University, Novosibirsk, Russia 3 Institute of Soil Science and Agrochemistry of SB RAS, Novosibirsk, Russia Abstract The results of cartographic modeling of the temperature fields of soils of the land fund of the Novosibirsk region using satellite data obtained by the MODIS/Terra spectroradiometer are considered. The developed cartographic models give a clear idea of the spatial structure of the temperature fields of the soil cover of the Novosibirsk region and the qualitative changes in the temperature regime of soils in different years. According to the data of daytime and nighttime surveys, the peculiarities of the temperature distribution of the soil cover for 2001 and 2010 were established. The values of the average annual temperature of the underlying surface for 2001–2014 were calculated. It is assumed that thermal resources not reflected in the generalization of zonal zoning can be identified on cartographic models. These resources have independent ecological significance and characterize the diversity of landscape, anthropogenic and other types of climate. Keywords MODIS, satellite data processing, soil surface temperature, cartographic modeling of soil temperature fields, land fund of the Novosibirsk region. 1. Introduction Rational use of the land fund of a territory is impossible without construction of cartographic models for land surface temperature fields. The analysis of such models allows distinguishing the climatic features of the regional soil climate, revealing the uniqueness of the soil climate, and indicating their possible value for agricultural and industrial purposes [1, 2, 3, 4]. While studying the soil temperature fields the data of space imagery in the thermal infrared range can be used as a source of information [5, 6, 7]. Such images data can be obtained with the MODIS spectroradiometer installed on the Terra satellite. The regularity of the survey and the significant coverage (sometimes very inaccessible) territory is the undoubted advantage of using remote sensing (RS) data of low spatial resolution [8, 9]. For a wide scale natural-anthropogenic complex (such as the Novosibirsk region) the assessment of land fund heat supply promotes efficiency of the using of all categories of land, regardless of the intended purpose, and its SDM-2021: All-Russian conference, August 24–27, 2021, Novosibirsk, Russia " pestunov@ict.sbras.ru (I. A. Pestunov) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings http://ceur-ws.org ISSN 1613-0073 CEUR Workshop Proceedings (CEUR-WS.org) 575 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 temperature regime is a tool in the study of the ecological state of the territory and a rational use of its natural resources [10, 11, 12, 13, 14]. Cartographic modelling of the temperature fields of the land surface of the Novosibirsk region by MODIS/Terra satellite data for studying the temperature regime of soils and identifying its features is the purpose of this work. To achieve this goal, the analysis of the spatio-temporal series of temperatures obtained by processing and aggregation of satellite images for various time intervals was carried out. The average annual temperature of the land surface is a meteorological indicator. It reflects the most important physical processes occurring in atmosphere and soils. Thus, the average annual land surface temperature is an indicator of the characteristics of the soil climate and the interaction of two interpenetrating environments — soils and plants. 2. Study area The Novosibirsk Region is located in the southeast of the West Siberian Plain (55∘ N–85∘ E), occupies a relatively small and compact territory — 177.8 th. km2 , which length is 444 km from North to South and 642 km from East to West (Figure 1). According to soil and climatic zoning [10], the Novosibirsk region is characterized by a wide range of zonal climate types, the thermal resources vary from 1200 ∘ C in the southern taiga zone to 2100 ∘ C in the typical steppe zone. The region includes several zones with landscape climate types: the southern taiga on meadow-boggy; subtaiga on gray forest soils; northern forest- steppe on chernozems leached and meadow solonetzic; southern forest-steppe on chernozems Figure 1: Google Earth Pro image of the research area (Novosibirsk region). 576 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 ordinary and leached and meadow-alkaline soils; steppes on solonetzic southern and ordinary chernozems and on sod-podzolic; and steppes on gray forest soils in Salair vertical zonation system. The forest-steppe subzone of the Novosibirsk region includes Barabinsky, Priobsky and Salairsky climatic districts. The vast Barabinsk forest-steppe (117 th. km2 ) is located within a wide strip of the Ob-Irtysh interfluve and is a plain slightly inclined from northeast to southwest, slightly elevated above sea level. The relief of Baraba is hilly in the north and hryvnia or hryvnia- hollow in the south. Its characteristic geomorphological features are boggy inter-crest basins, many fresh and salt lakes, a high degree of frostiness and swampiness in the north, which sharply decreases in the south [10, 11, 12, 13]. Zonal climates of Baraba are divided into a system of local climates united in the Northern, Western, Central and Eastern regions. It’s climate- forming factors affects geographical distribution and heat supply of the of soils. Podzolic, gray forest soils and ordinary chernozems are widespread in the northern and northeastern parts of Baraba. The western and southwestern regions of the Barabinsk lowland are characterized by ordinary and solonetzic chernozems, solodized soils and malts. The center of the Baraba is occupied by meadow chernozems. The flat nature of a relief with almost ideal latitudinal bioclimatic zoning makes the Baraba forest-steppe a unique object of research on a planet. This makes the research region promising for understanding the dynamics of changes in the zoning of ecosystems in a changing climate. The Novosibirsk city (located on the Ob river) is the center of a large agglomeration, repre- senting a compact accumulation of settlements, mainly urban, locally intertwining and united in a complex multicomponent dynamic system with intensive production, transport and cultural links. The core of the agglomeration is formed by the city of Novosibirsk and adjacent towns and villages. A characteristic feature of the Novosibirsk urban agglomeration is its pronounced monocentricity, since the city occupies an area much larger than the total area of its constituent settlements. The developing urban agglomeration has a powerful impact on environmental objects including the cause of the formation of heat islands which lead to a change in the temperature regime of urbanized territories. Thus, the Novosibirsk region is a complex natural and anthropogenic complex, and its temperature regime is a tool for the study of the ecological state of the territory [14]. 3. Research materials — spatio-temporal data series One of the features of this study is the extraction and processing of satellite RS data from the archive of the Federal Research Center of ICT (FRC ICT) using the hVault technology [16]. Unlike traditional unloading of individual scenes with subsequent connection to a geographic information system (which allows visualizing data as a set of thematic layers but complicates simultaneously using a large number of images with complex processing functions), hVault technology provides virtual integration of the data from the archive into relational DBMS. This technology is based on the principle of data presentation as a set of tables containing satellite image data or derived information products, followed by analysis of spatio-temporal series by means of a DBMS. The mapping of these data to tables is as follows. Each spatial point is associated with a cortege containing geographic coordinates, a time interval of observations, and 577 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 a sequence of measured values. In this form, data extraction and transformation are performed using SQL queries. The original data storage format and their division into separate files have no matter to the user. The software module implementing this principle was developed as an extension of the free PostgreSQL DBMS [15]. In the analysis of the spatio-temporal series of the remote sensing data, the products ob- tained as a result of processing MODIS/Terra data were used. The temperature values are reconstructed from measurements of the intensity of infrared radiation recorded in channels 31 (10.78–11.28 𝜇m) and 32 (11.77–12.27 𝜇m), according to the MOD11A1 product construction algorithm [15]. The data contained in MOD11A1 are presented on a regular grid in a sinusoidal projection with a cell size of ∼1 km. Under good (cloudless) atmospheric observation conditions, the algorithm [16] provides an accuracy of surface temperature recovery within 1 K. The average values of the temperature measured at nighttime and daytime, as well as their difference at a point, were calculated by data from the FIC ICT archive for — the whole year (2002 and 2010); — for the entire time interval 2001–2014 (a period of 14 years); — the spring period from April 10 to May 30 for each calendar year from 2001 to 2014; — the spring period from April 10 to May 30 for the entire time interval 2001–2014 (a period of 14 years). The average value at a point for 14 years was calculated by average annual arrays to speed up the calculations. The study area covers 400000 points. It can be said with some convention that the calculation with the hVault technology is equivalent to simultaneous processing of 365 satellite images. 4. Results and discussion The average annual temperature of the soil surface is a meteorological element which reflects the characteristics of heat exchange throughout the year and depends on the intensity of solar radiation, the characteristics of atmospheric circulation, and the physical characteristics of the soil cover. To get an idea of the changes in the temperature regime of the Novosibirsk region we obtain the averaged surface temperature at the geographical point of the territory for a calendar year and for all the years under consideration (from 2001 to 2014) (Figure 2). According to the nighttime survey, the lowest minimum average temperature was in 2010 and 2012, and the highest maximum was in 2002 and slightly lower in 2003. According to the daytime survey, the lowest average temperature was in 2010, and the highest — in 2010 and 2012. According to the data of the regime network of the Hydrometeorological Service, the average annual soil temperature on the surface changes by 1–2 ∘ C, and at the same time, zoning is strictly traced: 0 ∘ C — in the subtaiga; 1 ∘ C — in the forest-steppe; and 2 ∘ C — in the steppe. It is noted that the surface temperatures of saline soils are 1 ∘ C lower than on non-saline automorphic soils. Extracted spatio-temporal series of the Earth surface temperature were used for calculation of the average values for nighttime and daytime surveys at a geographic point and its differences for one calendar year. Calculations were carried out for each year from 2001 to 2014. The 578 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 Figure 2: Average values and difference of soil surface temperature according to data of nighttime and daytime surveys. obtained results allow suggesting the heterogeneity of the spatial structure of Earth surface temperature and qualitative changes in the temperature regime from year to year. Figures 3 and 4 shows the distribution of the average temperature according to the data of nighttime and daytime surveys and its differences for 2002 and 2010, respectively. It can be seen that the average temperatures differ significantly in these years, the land surface temperature in 2002 is higher than in 2010. Both in the more severe and in the milder manifestation of the weather regime, the classic latitudinal zoning of the land surface temperature distribution is clearly visible. It is very typical for the southeast of Western Siberia: from north to south by zones, there is a stable increase in average negative, average long-term minimum temperatures and absolute minimum temperatures. Thus, the severity of the winter period increases from south to north, which also causes an increase in the environmental load on human and animal organisms during winter periods. Figure 5 shows the distribution of the average annual temperature of the soil surface in the system of zonal types of climate in the Novosibirsk region. First, the average value of the temperature in a pixel was determined from the data for one calendar year for 2001–2014. Then the average value for 2001–2014 was found based on the RS data for nighttime and daytime surveys. The same actions were performed for the difference between the daytime and nighttime temperatures of the underlying surface. The latitudinal zoning of the soil temperature distribution from south to north from the typical steppe zone to the southern taiga zone is clearly traced here (Figure 5, a and b). Within the zonal subdivisions, subzonal types of soil climates of the Baraba lowland can also be distinguished — podzolic and gray forest soils are widespread in the northern part, meadow chernozems soils occupy its central part, and leached and ordinary chernozems prevail in the southern subzone of the forest-steppe. Temperature fields with similar temperature values, marked on the map of the region, can characterize a variety of landscape, anthropogenic and other types of climate. Thermal resources not reflected in the generalization of zonal zoning are noted. These resources have independent ecological significance and can be identified on cartographic models. The average difference between day 579 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 a b c Figure 3: Distribution of the average temperature for 2002. Nighttime (a), daytime survey data (b), and the difference between daytime and nighttime data (c). a b c Figure 4: Distribution of the average temperature for 2010. Nighttime (a), daytime survey data (b), and the difference between daytime and nighttime data (c). a b c Figure 5: Distribution of the average annual temperature of the underlying surface of the Novosibirsk region for the period 2001–2014. Nighttime survey data (a), daytime survey data (b), and the difference between daytime and nighttime data (c). and night temperatures varies from high in the south of the region to lower in the north, again with a clearly pronounced latitudinal zoning (Figure 5, c). Soil temperature is influenced by the type of soil, the nature of its salinization, and granulo- metric composition. Thus, the average annual temperature of the underlying surface throughout the entire territory of the Novosibirsk region varies from 0 to 2 ∘ C. The annual temperatures of automorphic and semi-hydromorphic soils are positive; its variability across the territory is 1.8–2.3 ∘ C. Saline soils are, on average, colder than non-saline soils per year. The pronounced seasonal rhythm of all components of nature in the temperate zone of the Novosibirsk region is located. It determines both the seasonal variation of soil temperature and the originality of the 580 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 factors of climate formation, which is reflected in the specificity of the thermal regime of soils. The dates of the transition through zero temperature are taken as the beginning of warm and cold periods in the soil, in the first case from negative ones, in the other — from positive ones. Figure 6 shows graphs of changes in the average day and night temperatures for the region for 2002 and 2010. Note that the transition from negative daily mean values to positive ones occurs in April (in rare, anomalous years in March) and nighttime temperatures change sign in April. Although, as we can see from the charts, negative values also occur in May. A steady transition from negative to positive temperatures marks the beginning of a warm period. It is important for agricultural land use. However, it is not possible to carry out a study to establish the day when both nighttime and daytime surface temperatures become positive in all or part of the region, due to significant gaps in the analyzed spatio-temporal data series (both in space and in time). Nevertheless, we made an attempt to obtain the distribution of soil temperature at the assumed dates of the spring seasonal change in the sign of the temperature value. The great importance for the use of land in agricultural are the timing of the transition of night and day temperatures through the zero mark. Figure 7 shows the results of the analysis of the spring change in the mean temperature values near the days of the change in the temperature sign. Three monthly periods were considered: from April 10 to May 10, 2010; from April 20 to May 20, 2010; and from May 1 to May 30, 2010. For each period, the average night and day temperatures in a pixel and its differences were obtained. Average daily temperatures in these time intervals become positive almost throughout the entire territory, moreover, they significantly exceed 10 ∘ C, and we do not present them. But the average night temperatures in the pixel are still negative, but comparing the three spring periods, it can be seen that the “warming” of the soil surface occurs in the latitudinal-zonal direction. In May, almost the entire territory of the Novosibirsk Region is already in the zone of positive night temperatures (Figure 7, c). In May, the soil surface temperature fluctuates more significantly than in winter. It is clearly seen in Figure 8. a b Figure 6: Average temperature (∘ C) in the Novosibirsk region for 2002 (a) and 2010 (b). 581 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 a b c Figure 7: Distribution of average night temperature for periods: a — 10.04–10.05.2010; b — 20.04– 20.05.2010; c — 1.05-30.05.2010. Correspondence of colors to temperature values, 𝑇 ∘ C: dark blue — < 0; light blue — 0–5; green — 5–10; red — > 10. a b c Figure 8: Distribution of the average difference between day and night temperatures for the periods: a — 10.04–10.05.2010; b — 20.04–20.05.2010; c — 1.05–30.05.2010. The soil temperature in the northern regions passes to positive values faster and more intensively than in the regions located to the south. As you can see, zoning is well traced here. 5. Conclusion Soil surface temperature is an important indicator in assessing disturbances in ecological and climatic interactions. Cartographic models of thermal resources obtained from the analysis of the spatio-temporal series of Terra/MODIS satellite data allows assessing the heat supply of the land fund of the Novosibirsk region, which makes it possible to effectively use all categories of land, regardless of their intended purpose. Studying the temperature of the soil surface also makes it possible to reveal the features of the soil climate, which is important for considering ecological processes and their interaction with two interdependent environments — soil cover and vegetation. It also requires an analysis of the spatio-temporal series of the underlying surface temperature. The obtained results can serve as a basis for zoning the territory of the Novosibirsk region according to the degree of environmental safety. 582 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 Acknowledgments The study was carried out according to the state assignment of Federal Research Center for Information and Computational Technologies and Institute of Soil Science and Agrochemistry of Siberian Branch of the Russian Academy of Sciences with financial support from the Ministry of Science and Higher Education of the Russian Federation. References [1] Kudryashova S.Ya., Chumbaev A.S., Pestunov I.A., Sinyavskiy Yu.N., Chubarov D.L., Bezborodova A.N., Ermakov N.B. Cartographic modeling of soil temperature fields for Middle Siberia Transect based on conjoint analysis of automated ground-based and satellite temperature data // CEUR Workshop Proceedings. 2020. Vol. 2534. P. 124–131. ISSN:1613- 0073. Available at: http://ceur-ws.org/Vol-2534/21_short_paper.pdf. [2] Voronina P.V., Chubarov D.L., Dobretsov N.N. Analysis of time series of the Earth’s surface temperature // Proceedings of the International Conference “Computational Mathematics and Mathematical Geophysics” dedicated to the 90th anniversary of the birth of Academi- cian A.S. Alekseev. 2018. P. 108–114. (In Russ.) [3] Voronina P.V., Chubarov D.L., Dobretcov N.N. The study of the temperature regime of the Novosibirsk agglomeration according to the satellite sensing data // CEUR Workshop Proceedings. 2020. Vol. 2534. P. 327–333. ISSN:1613-0073. Available at: http://ceur-ws.org/ Vol-2534/55_short_paper.pdf. [4] Voronina P.V., Chubarov D.L., Dobretsov N.N. Temperature anomalies on the territory of Barabinskaya Lowland. The first results // Interexpo Geo-Siberia. 2019. Vol. 4. No. 1. P. 28–36. (In Russ.) [5] Voogt J.A., Oke T.R. Thermal remote sensing of urban climates // Remote Sensing of Environment. 2003. No. 86. P. 370–384. [6] Weng Q. Thermal infrared remote sensing for urban climate and environmental studies: Methods, applications, and trends // ISPRS Journal of Photogrammetry and Remote Sensing. 2009. No. 64. P. 335–344. [7] Weng Q., Quattrochi D.A. Thermal remote sensing of urban areas: An introduction to the special issue // Remote Sensing of Environment. 2006. No. 104. P. 119–122. [8] Cheval S., Dumitrescu A. The July urban heat island of Bucharest as derived from MODIS images // Theoretical and Applied Climatology. 2009. No. 96. P. 145–153. [9] Pu R., Gong P., Michishita R., Sasagawa T. Assessment of multi-resolution and multi- sensor data for urban surface temperature retrieval // Remote Sensing of Environment. 2006. No. 104. P. 211–225. [10] Soil and climatic atlas of the Novosibirsk region. Novosibirsk: Nauka, Sib. Otdelenie, 1978. 121 p. (In Russ.) [11] Voronina L.V., Gritsenko A.G. Climate and ecology of the Novosibirsk region. Novosibirsk: SGGA, 2011. 228 p. (In Russ.) [12] Kravtsov V.M., Donukalova R.P. Geography of the Novosibirsk region. Novosibirsk: Info- press, 1999.205 p. (In Russ.) 583 Polina V. Voronina et al. CEUR Workshop Proceedings 575–584 [13] West Siberian Lowland. Moscow: Geography, 1963. 262 p. (In Russ.) [14] Shokin Yu.I., Pestunov I.A., Chubarov D.L., Kikhtenko V.A., Mamash E.A., Dobretsov N.N., Voronina P.V., Sinyavsky Yu.N. Information system for satellite data acquisition, archiving and analysis as an instrument for monitoring of objects of artificial and natural origin // CEUR Workshop Proceedings. 2020. Vol. 2534. P. 229–240. ISSN:1613-0073. Available at: http://ceur-ws.org/Vol-2534/39_short_paper.pdf. [15] Shokin Yu.I., Dobretsov N.N., Mamash E.A., Kikhtenko V.A., Voronina P.V., Smirnov V.V., Chubarov D.L. An information system for acquisition, processing and access to satellite data and its applications in environmental monitoring // Computational Technologies. 2015. Vol. 20. No. 5. P. 157–174. (In Russ.) [16] Wan Z., Dozier J. A generalized split-window algorithm for retrieving land-surface tem- perature from space // IEEE Trans. Geoscience and Remote Sensing. 1996. Vol. 34. No. 4. P. 892–905. 584