A new multi-level model for calculating the carbon footprint of agroecosystem products is proposed. The concept of “final carbon footprint” is introduced, which includes both direct CO₂ emissions from the operation of tractors, combines, oxidation of soil humus, CO₂-eq. during the transformation of nitrogen fertilizers in the soil, and indirect CO₂ emissions – carbon dioxide release into the atmosphere during the production of tractors, combines, tillage equipment, mineral fertilizers etc.
Based on the results of field experiments on gray forest soils in the Southern Moscow region, it is shown that when applying average doses of mineral fertilizers to field crops, the indirect CO₂ emissions are comparable to the CO₂ input from organic fuel oxidation when machinery is operating in the field. At higher doses of fertilizers, the indirect emissions are significantly greater than the CO₂ emissions from machinery operation. In order of increasing CO₂ emissions per 1 ha of sowing, crops on gray forest soils are arranged as follows: corn for silage > barley > winter wheat > clover.
Clover is a carbon-negative crop (−1.7 t/ha CO₂), i.e., CO₂ sequestration in the soil exceeds all CO₂ emissions from hay crop production. The final carbon footprint for grain crops, calculated using the standard method, was as follows: for winter wheat (with a fertilizer dose of N40P40K40) – 116 kg CO₂ per 1 centner of grain, for barley (with a dose of N60P40K40) – 104 kg CO₂ per 1 centner of grain. The final carbon footprint, taking into account the aftereffects of predecessors, was: for winter wheat (predecessor: two-year clover) – 48 kg CO₂ per 1 centner of grain; for barley (predecessor: silage corn) – 113 kg CO₂ per 1 centner of grain.

The results of our phenological observations of flowering rhythms of 37 herbaceous and 21 woody plants from the collections of the Botanical Garden of Moscow State University in 2019 are presented for three phases, namely, the beginning, peak and the end of flowering. The average timing of the “peak flowering” phase onset for the period of 2009–2019 was calculated, and it has turned out that the vast majority of the studied plants bloomed in 2019 earlier than usually. To assess the effect of summary active temperatures (SAT) on plant development, their sums were calculated by months (according to the data of MSU Meteorological Observatory) with thresholds of 0°С, 5°С, and 10°С in comparison with their average values. Our comparison of the data obtained has shown a direct effect of SAT on the timing of the onset of flowering phenophases, and also confirmed the importance of not only the sum of positive temperatures, but also the sum of temperatures exceeding the threshold of 10°C for flowering of the primary flowers. For a significant advance in the average timing of flowering, such plants need to receive an excess of heat about a month before flowering. An increased amount of active temperatures even immediately before and during flowering has a beneficial effect on the timing of flowering. At the same time, SAT values above 10°C are important.

A scheme of formation of the sedimentary blanket of the South Kara basin is considered, which can further be used for numerical reconstruction of its thermal history. The scheme is based on our analysis of the literature information on the structure and geological history of the Barents-Kara region. This information included an interpreted seismic profile crossing the studied area, drilling data from four wells located along the profile (the University, Rusanov, Leningrad and Kharasavey ones), measurements of the heat flow and deep temperatures in the basin. The proposed scheme considers the formation of the basin as a series of sedimentation stages with various combinations of clay shales, siltstones and sandstones and sedimentation in the Cretaceous and Paleogene with their subsequent erosion in the Miocene. The erosion amplitude is estimated by the observed change in the porosity of sedimentary rocks with depth. The initial heat flow in such a model should correspond to the flow of modern axial zones of continental rifting or be lower for the areas remote from the corresponding segments of the Late Permian-Early Triassic continental rifting system.

The thermal evolution of the permafrost in the sedimentary section of the Tyumen superdeep SG-6 well has been numerically reconstructed using the ICE2020 software package, which is part of the GALO flat basin modeling system. The thermal evolution of the sedimentary strata in the last 3.5 My is considered as the final stage of the basin modeling, whose formation began with continental rifting in the Late Permian. Abrupt climate changes in the late Pliocene–Holocene led to a decrease in the rock temperature by 15–20°C in the upper 1–1.5 km of the SG-6 sedimentary section. The maximum thickness of the permafrost in the study area was about 711 m, reached 2.6 Mya. The maximum thickness of the permafrost for the last ice age (23–18 thous and years ago) was 412 m, reached about 14.5 thousand years ago. According to our modeling, the modern base of the permafrost is at the depth of 311 m and is degrading with the rate about 13 m/1000 y. The results of our calculations with a database of climatic data limited to the last 50 and 100 thousand years differ markedly from the modeling results with the complete database for the last 3.5 My.