The Comparison of Soil Agrochemical and Biological Properties in the Multi-Cropping Farming Systems

Multi-cropping systems play an important role in improving the quality of soil properties. A field experiment was carried at the Experimental Station of Vytautas Magnus University Agriculture Academy (Lithuania) in 2017 to 2019. The aim of the study was to compare agrophysical and biological properties of the soil in the multi-cropping systems of sole (spring barley, spring wheat, pea, caraway), binary (spring barley–caraway, spring wheat–caraway, pea–caraway) and trinary (spring barley–caraway–white clover, spring wheat–caraway–white clover, pea–caraway–white clover) crops. In the second and the third years of caraway cultivation, when solely caraway was grown, the total nitrogen content was significantly lower than in binary and trinary crops (8.5% and 17.4%, respectively). The results indicated that the highest organic carbon content was in the third year of caraway cultivation in trinary crop when caraway was grown with peas and white clover. In the third year, the highest saccharase and urease activity was found in trinary crop where caraway was grown with spring barley and white clover. A strong positive correlation was observed between the content of saccharase and urease and the total nitrogen, organic carbon, and potassium available in the soil. The results of the study suggest that multi-cropping is important for soil conservation and the sustainability of agro-ecosystems.


Introduction
Over the last five decades, advances in agriculture have made it possible to meet the world's highest demands for food, feed, and fibre [1]. However, sustaining the needs of an ever-increasing population is challenging due to widespread urbanisation, severe land degradation, and climate change [2]. Successful adaptation to and mitigation of climate change through agricultural management requires the development of simple, costeffective, and largely scalable approaches. Therefore, an agricultural management strategy is essential [3]. The solution to these challenges is to achieve the long-term sustainable use of resources, taking into account the ecological and economic aspects of sustainability. In essence, eco-efficiency is about achieving more agricultural production with fewer resources [4]. Multi-cropping systems could increase crop diversity and avoid vulnerability to biotic stress [5], as well as ensure the sustainability of agriculture and increase crop and food production while reducing land use [6,7].
Multi-cropping is the cultivation of two or more crops in the same field based on ecological principles [8][9][10]. The main differences between mono-and multi-cropping are that the cultivation of agricultural crops differs in the length of the growing season as well as in biological and agronomic characteristics [11,12]. According to Lizarazo et al. [12], multi-cropping reduces nutrient leaching into deeper layers of the soil, as well the study was to determine and compare the agrochemical (total nitrogen, organic carbon, available phosphorus, available potassium) and biological (root biomass, saccharase, and urease activity) soil properties in sole, binary, and trinary crops.

Total Nitrogen Content in the Soil
First year of caraway vegetative season. At the end of the growing season in 2017, the total nitrogen content was significantly higher in the soil under binary crops of spring barley with caraway and pea with caraway and a trinary crop of pea with caraway and white clover by 15.8, 17.4, and 19.3%, respectively, compared to the sole crops of the aforementioned plants ( Table 1). The soil with the sole caraway crop was not significantly different from the binary and trinary crops in terms of total nitrogen content. Second year of caraway vegetative season. After caraway harvesting in 2018, the total nitrogen content in the soil with caraway after spring barley (binary crop) and caraway with white clover after spring barley and spring wheat (trinary crop) was significantly higher-13.8, 12.1, and 13.0%, respectively-compared to the sole crops of the same plants ( Table 1). The total nitrogen content in the soil with the sole caraway crop was found to be significantly (8.5 to 13.7%) lower than in binary and trinary crops.
Third year of caraway vegetative season. After caraway harvesting in 2019, the content of total nitrogen in the soil with caraway after spring barley and spring wheat with and without white clover was significantly higher (19.4, 14.5, 21.3, and 18.2%, respectively) compared to the soil with the sole crops ( Table 1). The total nitrogen content in the bare fallow left after caraway harvesting was found to be significantly (9.2 to 17.4%) lower than in the soil, where caraway was grown with and without white clover.

Organic Carbon Content in the Soil
First year of caraway vegetative season. The organic carbon content in the soil under binary crops of spring barley with caraway, spring wheat with caraway, and pea with caraway and trinary crops of spring wheat with caraway and white clover and of pea with caraway and white clover was found to be significantly higher (by 16.3, 16.7, 9.0, 17.8, and 14.6%, respectively) than in the soil with the sole crops (Table 2). the organic carbon content of the soil with the sole caraway crop was found to be significantly higher than that of the binary crop of caraway with peas and the trinary crop of caraway with spring barley and white clover, 8.2 and 9.4%, respectively. Second year of caraway vegetative season. The organic carbon content in the soils where caraway grew after spring barley with spring wheat (binary crop) and after spring wheat with white clover (trinary crop) was found to be significantly higher, by 16.0, 15.3, and 13.3%, respectively, than in the soils with sole crops (Table 2). the organic carbon content in the soil with sole caraway crop was not significantly different from that of the binary and trinary crops.
Third year of caraway vegetative season. After caraway harvesting in 2019, the organic carbon content in the soils where caraway grew after spring barley and spring wheat without and with white clover was found to be significantly, 12.9, 16.4, 13.6, and 15.7%, higher, by 12.9, 16.4, 13.6, and 15.7%, respectively, compared to sole crops ( Table 2). The organic carbon content in the soil of bare fallow after caraway harvesting was found to be significantly lower, by 15.6 to 18.8%, respectively, than in the soil covered with binary and trinary crops.

Available Phosphorus Content in the Soil
First year of caraway vegetative season. At the end of the growing season in 2017, the available phosphorus content in the soil under binary crop of spring barley with caraway and trinary crop of spring barley with caraway and white clover was found to be significantly higher, by 34.9 and 27.1%, respectively, than those of the sole crops of the aforementioned plants ( Table 3). The available phosphorus content in the soil with sole caraway crop was found to be significantly lower compared to its binary crop with spring barley and trinary crop with spring barley and white clover by 15.9 and 10.7%, respectively. Second year of caraway vegetative season. After caraway harvesting, the available phosphorus content in the soil where caraway was grown after spring barley without white clover and together with white clover was found to be significantly higher, by 61.1 and 44.6%, respectively, than in the soil where spring barley was grown alone ( Table 3). The available phosphorus content in the soil with sole caraway crop was not significantly different from that of the binary and trinary crops.
Third year of caraway vegetative season. After caraway harvesting in 2019, the available phosphorus content in the soil where caraway was grown after spring barley and spring wheat without white clover and with it was found to be significantly higher, by 33.1, 30.1, 19.7, and 24.8%, respectively, than in the soil under sole cropping ( Table 3). The available phosphorus content in bare fallow soil after caraway harvesting was found to be significantly lower, by 8.5 to 23.7%, respectively, than in binary and trinary crops, except when caraway was grown after pea without and with white clover.

Available Potassium Content in the Soil
First year of caraway vegetative season. At the end of the growing season in 2017, the available potassium content in the soil with trinary crop of spring barley with caraway and white clover was significantly higher, by 14.9%, compared to the soil with a sole crop ( Table 4). The soil with a sole caraway crop did not show any significant difference in the available potassium content compared to binary and trinary crops.
Second year of caraway vegetative season. After caraway harvesting, there was no significant difference in the available potassium content in the soil of sole, binary, and trinary crops (Table 4).
Third year of caraway vegetative season. The available potassium content in the soil where caraway grew after spring barley, spring wheat, and pea without white clover (binary crops) and with clover (trinary crops) was found to be significantly higher than in the soil of sole crops, 46.2, 32.9, and 30.2 and 62.8, 40.7, and 36.2%, respectively (Table 4). Available potassium content in the soil of bare fallow after caraway harvesting was not significantly different compared to binary and trinary crops.

Root Dry Biomass
First year of caraway vegetative season. When spring barley was grown in binary crop with caraway intercrop and when growing spring barley, spring wheat, and pea in trinary crop with caraway and white clover intercrop, significantly higher, by 78.5, 85.8, 60.4 and 53.0%, plant root mass was determined compared to sole crop of the aforementioned plants ( Figure 1). Second year of caraway vegetative season. In 2018, when caraway was grown after spring wheat harvesting without white clover (binary crop) and with white clover after spring barley, spring wheat, and peas (trinary crop), the plant root mass was significantly higher compared to that of sole crop-2.1, 1.9, 2.5, 1.7, and 2.5 times, respectively ( Figure 1).
Third year of caraway vegetative season. In 2019, when caraway was grown after spring wheat harvesting (binary crop) and with white clover (trinary crop), the plant root mass was found to be significantly higher, by 74.4 and 55.0%, respectively, compared to sole crop ( Figure 1). The plant root mass in bare fallow left after caraway harvesting was found to be significantly lower, by 5.1 to 9.0 times, respectively, compared to that of binary and trinary crops.

Saccharase Activity in the Soil
First year of caraway vegetative season. The highest saccharase activity was found in the soil of binary crop of spring barley with caraway ( Figure 2). This could be influenced by the formation of a more abundant plant root system. Compared to the sole crop, except caraway, and binary crop of spring wheat and pea with caraway intercrop, the saccharase activity was found to be significantly higher, by 30.2, 34.7, 32.9%, 27.6, and 65.8%, respectively. The soil saccharase activity was found to be significantly higher, by 42.7%, in the trinary crop of pea with caraway and white clover compared to that in the binary crop of pea with caraway.
Second year of caraway vegetative season. The activity of enzyme saccharase in the soil in the second year of caraway growth was found to be higher than in the first year ( Figure 2). This could be influenced by the higher content of organic matter in the soil. The activity of the enzyme saccharase in the soil where caraway was grown together with white clover after spring barley and pea was found to be significantly higher compared to sole crop-by 34.8 and 39.0%, respectively. The saccharase activity in the soil where caraway was grown with white clover after peas was found to be significantly higher, by 21.6%, than that in the soil where caraway was grown without clover. The activity of enzyme saccharase was found to be significantly lower, from 34.0 to 68.7%, in the soil with a sole crop of caraway than in the soil where caraway was grown after pea (binary crop) and after spring barley, spring wheat, and peas in combination with white clover (trinary crop).
Third year of caraway vegetative season. The activity of enzyme saccharase in the soil in the third year of caraway cultivation was found to be higher than in the first and second years of cultivation ( Figure 2). In the soil where caraway was grown after peas without white clover, and together with clover after spring barley, spring wheat, and peas, the saccharase activity was found to be significantly higher than in the soil where caraway was grown as the sole crop-64.0, 79.0, 43.0, and 56.2%, respectively. The saccharase activity in the soil where caraway was grown with white clover after barley was found to be significantly higher, by 42.4%, compared to the soil where it was grown without clover. The activity of enzyme saccharase in the soil with binary and trinary crops was found to be significantly higher, by 2.0 to 3.2 times, than in bare fallow left after caraway harvesting.

Urease Activity in the Soil
First year of caraway vegetative season. The activity of enzyme urease in the soil with binary crop of spring wheat with caraway and trinary crops of spring wheat and pea with caraway or white clover intercrops was found to be significantly higher compared to sole crop-16.7, 16.7 and 2.5 times, respectively (Figure 3). The activity of enzyme urease in the soil with sole crop of caraway was found to be significantly lower, by 7.5 to 12.5 times, than that of caraway grown in binary and trinary crops.
Second year of caraway vegetative season. The activity of enzyme urease in the soil in the second year of caraway cultivation was also found to be higher than in the first year ( Figure 3). Urease activity in the soil with a trinary crop of caraway with white clover after spring wheat and pea was found to be significantly higher compared to sole crop by 3.7 and 2.0 times, respectively. The activity of enzyme urease activity in the soil where caraway grew with white clover after spring wheat and pea was found to be significantly higher, by 83.3 and 71.4%, than that in the soil of binary crop. Urease activity in the soil with sole crop of caraway was significantly lower, by 2.2 to 2.4 times, than in the soil of the trinary crop after spring wheat or pea.
Third year of caraway vegetative season. The activity of enzyme urease in the soil was found to be significantly higher than in the first and second years of caraway cultivation (Figure 3). The highest urease activity was found in the soil of trinary crop of spring barley with caraway and white clover. The urease activity in the soil in which caraway was grown without white clover (binary crop) and together with it (trinary crop) after spring barley, spring wheat, and pea was found to be significantly higher compared to sole crop-2.3, 4.3, and 3.8 and 4.7, 6.7, and 3.3 times, respectively. The urease activity in the soil where caraway was grown with white clover after spring barley was found to be significantly higher, by 2.0 times, than that in the soil where caraway was grown without clover. The urease activity in bare fallow left after caraway harvesting was found to be significantly lower, by 6.5 to 14.0 times, compared to the soil where caraway was grown without white clover (binary crop) and in combination with white clover (trinary crop).

Principal Component Analysis (PCA)
PCA showed different correlations between individual parameters depending on the multi-cropping system used. Two groups of correlated parameters were identified in the first year (2017) of caraway cultivation (Figure 4). The first group of correlated parameters included total nitrogen, organic carbon, root biomass, and the activity of soil enzyme urease. The second group of correlated parameters included available phosphorus, available potassium, and the activity of soil enzyme saccharase. In the second year (2018) of caraway cultivation, the activity of the soil enzyme urease was more dependent on the agrochemical soil properties (total nitrogen, available potassium, and organic carbon) and root biomass than on the activity of enzyme saccharase.
In the third year (2019) of caraway cultivation, the activity of the soil enzymes saccharase and urease correlated with the levels of total nitrogen, organic carbon, and available potassium.

N, OC, P, K Contents in the Multi-Cropping System
Studies have shown that the soil nutrient balance was better maintained in soils with binary crops with caraway intercrop and trinary crops with caraway and white clover intercrop compared to sole crops. The total nitrogen, organic carbon, available phosphorus, and potassium in the soil increased the most in binary and trinary crops in the third year of caraway cultivation compared to the initial contents.
The results of other authors showed that the fixed nitrogen content in leguminous plants in multiple crops depended on several factors such as species, plant morphology, crop density, and total nitrogen in the soil [33]. The results of Andersen et al. [34] show that when forming multiple crops of legumes and other crops, the benefits of N 2 fixation are partially lost if the selected plants have a stronger competing capacity than leguminous plants. Hauggaard-Nielsen et al. [35] stress that multi-cropping can have a positive impact on nutrient retention. Despite competition between plants, significantly higher total nitrogen levels were found in the soil of the second and third years of caraway cultivation compared to the sole crop ( Table 1).
The soil organic carbon content increased significantly by 3.8-17.8% in 2017-2019 in the multi-cropping system ( Table 2). The study confirmed that the availability of nutrients in the soil was significantly higher in binary and trinary crops compared to sole crops. Data from Adamu and Yusuf [36] showed the same trend of increasing organic carbon in multiple crops compared to sole crops. This can be explained by the fact that root exudates create favourable conditions for the nutrition of microorganisms and thus for plant nutrition. The soil is also positively affected by the multiple crop root system, which promotes the activity of micro-organisms. Studies have shown that carbon availability increased in the carbon cycle [37]. More nutrients important for plant mineral nutrition were released, while a denser soil cover reduced the rate of mineralisation and inhibited nutrient leaching.
The available phosphorus also depends on the organic matter content in the soil [38]. In the presence of available phosphorus deficiency and drought conditions, it has been observed that crops grown in mixtures absorb phosphorus released in the soil from leguminous plants directly, thus allowing for plant resilience [39]. In addition, the availability of organic phosphorus is also increased by a variety of leguminous crops. Muofhe, Dakora [40] found that growing legumes and cereals together increased the uptake of phosphorus due to the release of lactic anions that chelated Fe 3+ and subsequently released phosphorus from FePO 4 . Studies have shown that leguminous plants released higher levels of phosphorusmobilising exudates [41], which could also support phosphorus nutrition of other plant species in multiple crops with lower phosphorus uptake capacity [42]. Our experimental data confirmed that the phosphorus content in the soil was significantly increased with the cultivation of binary and trinary crops compared to sole crops (Table 3). This apparent effect may have been due to the higher accumulation of residues in the topsoil layer that were more resistant to decomposition. The intertwining of caraway roots with pea and white clover roots had a positive effect on the release of phosphorus from the soil, while at the same time, the remaining decomposing root parts contributed to the organic matter content in the soil.
In the third year of caraway cultivation, the available potassium was significantly higher in the soil with binary and trinary crops compared to sole crops ( Table 4). The results of Nasar et al. [43] show that the cereal-legume multi-cropping system significantly increased the available potassium in the soil compared to the sole crop. The reason for this increase in available potassium is not fully understood; however, it may be due to the effect of cation antagonism between potassium and calcium when roots interact in multiple-crop conditions [37,44,45]. In addition, most of the potassium is accumulated in plant by-products (straw) [46]. Therefore, it can be argued that a thicker layer of plant residues in the third year of caraway cultivation may also affect the higher amount of potassium in the soil.

Root Biomass in the Multi-Cropping System
In the case of multi-cropping, resource uptake and competition between plant roots during the growing season is more evenly distributed, with peak nutrient uptake occurring at different stages. The results of the study showed that in the first year of caraway cultivation, the plant root mass was significantly higher in binary and trinary crops compared to sole crops-by 53.0 to 85.8%, respectively. A similar trend was observed in both the second and third years of caraway cultivation, with a significantly higher-from 1.7 to 2.5 times-plant root mass formed in binary and trinary crops compared to sole crops ( Figure 4). Papa et al. [47] found that roots of different plant species took up moisture and nutrients at different rhythms and intensities. The interactions between plants and interspecific competition promote plant rooting [17]. Competition is also largely avoided during vegetation [43,48]. However, Bellostas et al. [49] received contradictory results: with a two-species multi-crop, the intertwining of plant roots in the early stages of growth led to negative competitive effects. According to the authors, two weeks after sowing the binary crop, plant dry matter production was reduced by 15-20% compared to the sole crop [49]. It can be argued that a more abundant and denser root system has a greater suction power, resulting in the better uptake of nutrients from the soil and supplying them to plants, making them better prepared for adverse environmental conditions.

Soil Enzyme Activity in the Multi-Cropping System
Changes in the saccharase activity in the soil depended mainly on the root mass and the amount of crop residues. The highest levels of enzyme saccharase were found in trinary crop in the third year of caraway cultivation (Figure 2). This means that multiple crops can stimulate the growth and development of plant roots, improving the biological and chemical properties of the soil. The highest saccharase activity was in the arable soil layer. It is mainly related to organic carbon, available phosphorus and potassium, CO 2 release, and plant residues. Cui et al. [50] found that the activity of the soil enzyme saccharase in multiple crops was significantly correlated with the soil ammonium and nitrate nitrogen. This means that pH and total nitrogen were the main factors influencing soil enzyme activity [50]. The increase in saccharase activity indicates an acceleration of hydrolytic carbohydrate degradation and intensification of organic matter mineralisation processes in the soil [51]. This implies that multiple crops can improve soil nutrient cycle by improving the activity of soil saccharase enzymes.
The studies on urease activity showed that the reduction in the tillage intensity in the soil in the third year of caraway cultivation, where caraway was grown with white clover after barley (trinary crop), may have created a better biochemical soil environment, leading to a higher stable enzyme concentration (Figure 3). In fact, soil enzymes are responsible for nutrient cycles and availability [20,52]. In addition, the enzyme urease acts by hydrolysing the C-N bond of amides and releasing ammonia during nitrogen mineralisation [53]. Thus, higher activity of this enzyme in trinary combination of spring barley, caraway, and white clover suggests that crop residues had a positive effect on enzyme urease activity and subsequently improved the nitrogen availability and soil biological properties. Moreover, trinary crops combinations with a predominance of leguminous plants increased the levels of this enzyme in the third year of caraway cultivation. No mineral fertiliser was applied and the root system of the leguminous plants transformed the rhizosphere of the plants grown together. Biological nitrogen fixation by leguminous plants acidified the rhizosphere and released higher levels of soil enzymes (urease and other enzymes) as well [54]. Our data showed that the content of saccharase and urease was the most correlated with the total nitrogen, organic carbon, and available potassium in the soil (Figure 4).
The texture of the topsoil is sandy loam, and the agrochemical properties are the following: pH KCl -6.70; organic carbon (OC)-0.91-1.08%; plant-available phosphorus (P 2 O 5 )-213-318 mg kg −1 ; and potassium (K 2 O)-103-125 mg kg −1 . Thus, the above-listed soil characteristics show that the topsoil richness in organic carbon (OC) is low, but the available phosphorus (P 2 O 5 ) corresponds to group V (very high content) and potassium (K 2 O) to group III (average concentration) according to the evaluation of agrochemical properties of Lithuanian soil [56]. A 1-factor field experiment with 10 treatments was set up in 2017 ( Figure 6). The experiment was performed in four replications, and a randomized complete block design (RCBD) was used. The size of each experimental plot was 60 m 2 (5 m × 12 m), one replication block was 600 m 2 ; 2 m buffer rows were left between the individual blocks. There were 40 plots in total in the experimental field. Figure 6. Experimental treatments: Four treatments were formed by sole crops, three treatments were combinations of caraway with the other crops as binary crops, and three treatments were trinary crops with white clover added to binary crops.
After the main crops-spring barley, spring wheat, and peas in 2017 and 2018-and after sole caraway crop in 2018, the harvested plots were disked and deeply ploughed. In the spring of 2018 and 2019, spring barley "Orphelia KWS" (180 kg ha −1 ) was sown (Figure 7). In the first year of caraway vegetation (2017), sole spring barley and spring wheat crops and binary crops with undersown caraway were fertilised with ammonium nitrate at the rate of 180 kg ha −1 , and trinary crops with undersown with caraway and white clover at the rate of 150 kg ha −1 .
In the second (2018) and the third (2019) years of caraway vegetation, caraway was not fertilised, and no plant protection products were used.

Meteorological Conditions
April was cold and humid ( Figure 8). May was very dry (hydrothermal coefficient (HTC) was 0.29). The amount of precipitation in June was close to the long-term average precipitation. July was cool with HTC of 1.53 (optimal humidity). In August, HTC was 1.00 (optimal humidity). September was warm with HTC 2.26 (excess humidity). In October, precipitation exceeded the long-term average of precipitation several times. In January 2018, the temperature was 2.2 • C higher than the long-term average. February and March were colder than usual. The temperature in April was 3.3 • C higher compared to the long-term average. The temperature in May was 4.0 • C higher than the long-term average, and the precipitation was 44.1 mm lower than usual. The monthly HTC was 0.33 (very dry), and the plants lacked moisture for growth. The temperature in June was 1.4 • C lower than the long-term average, and HTC was 1.10 (optimal humidity). The temperature in July was 1.4 • C higher than the long-term average, and the precipitation was 40.9 mm higher than usual, and HTC was 2.20 (excess humidity). August was warm with an HTC of 1.12 (optimal humidity). Temperatures were higher than usual in the autumn and winter months. April's temperatures were 2.2 • C higher compared to the long-term average, and HTC was 0.03 (very dry). In May, HTC was 0.92 (optimal humidity). The temperature in June was 4.4 • C higher compared to the long-term average, and HTC was 0.80 (insufficient humidity). July was cool with an HTC of 1.12 (optimal humidity). In August, the HTC was 1.21 (optimal humidity).
The Lithuanian climate zone is at humid continental climate (Köppen-Geiger code: Dfb) [57]. Dfb: Warm summer, temperate climate. Wet, rainfall is evenly distributed throughout the year. The average temperature does not exceed 22 • C in any month. Frost is usually absent for 3-5 months per year. In summer, heat (30 • C and above) rarely lasts longer than a week. This subtype of the climate is widespread throughout the territories of Lithuania and the rest of the Baltic States, only, in this area, the seasonal variation of temperatures is mitigated by the effects of the sea.

Analytical Methods
The agrochemical properties of the soil were determined before the installation of the experiment and after harvesting. Soil samples were collected in 15 plots, about 300 g from the 0-25 cm ploughing layer of each experimental plot using a soil auger. Soil pH was determined potentiometrically in 1 n KCl extract, soil total nitrogen (N) was determined by the Kjeldahl method, available phosphorus P 2 O 5 (P) and available potassium K 2 O (K) (mg kg −1 soil) were calculated using the Egner-Rim-Domingo (A-L) method, and organic carbon (%) (OC) was determined by incineration of samples at 900 • C using a Heraeus incinerator.
The plant root dry biomass was determined by the method of small monoliths (10 × 10 × 10 cm) [58] after the main crop harvest in the first year of caraway vegetative season and after the caraway harvest in the second and third years of caraway vegetative season. Samples were taken from two soil layers: 0-10 cm and 10-20 cm. The roots were washed through sieves and dried in a drying oven at 105 • C. Plant root mass was calculated in absolute dry matter content t ha −1 .
The activity of soil hydrolases (urease and saccharase) was determined: ureaseaccording to the method of Hofmann and Schmidt (1953); saccharase-according to the method of Hofmann and Seegerer (1950), modified by A. I. Chunderova (1973) [59] after the main crop harvest in the first year of caraway vegetative season and after the caraway harvest in the second and third year of caraway vegetative season. Soil samples were taken from 15 plots in each field with a soil drill at a depth of 0 to 25 cm. Samples of natural moisture were dried in boxes at laboratory temperature.

Statistical Analysis
The research data was statistically evaluated by one-way analysis of variance (ANOVA) using ANOVA software from the statistical analysis package SELEKCIJA [60]. The Duncan criterion was used to assess the significance of the differences. Differences between the averages of treatments marked with different letters are significant at the 95% confidence level (p < 0.05). There are no significant differences when p > 0.05.
Principal component analysis (PCA) was used to determine the correlation between different parameters in the multi-cropping system. PCA creates new artificial variables (principal components) based on the variables (features) that were analysed. Its main assumption is the possibility of visualizing the relationships of individual variables on a two-dimensional graph showing the coordinate system of the first two principal components. Based on the position of the vectors in space, it can be determined which features are correlated with each other. The smaller the angle between the vectors, the stronger the positive correlation. When the vectors are aligned on the same line but in opposite directions, there is a strong negative correlation between the variables. However, when the vectors are at an angle close to 90 degrees, no correlation occurs. Statistical analysis was performed using Statistica 10 software package (TIBCO Software Inc., Palo Alto, CA, USA).

Conclusions
The soil of binary and trinary crops had a better nutrient balance compared to the soil of sole crops. The total nitrogen, organic carbon, available phosphorus and potassium in the soil increased the most in binary and trinary crops during the third year of caraway cultivation compared to the initial level.
In the second and third years of caraway cultivation, significantly higher levels of total nitrogen and available phosphorus were found in the soil of multi-cropping system compared to sole crop. The soil organic carbon content increased significantly-from 3.8 to 17.8%-in the multi-cropping system in 2017 to 2019. In the third year of caraway cultivation, the available potassium content was found to be significantly higher in the soil of binary and trinary crops compared to the soil of sole crops.
In the first, second, and third years of caraway cultivation, plant root biomass in binary and trinary crops was significantly higher than that in sole crops.
The highest levels of saccharase and urease enzymes were found in trinary crop in the third year of caraway cultivation. The levels of saccharase and urease enzymes were most strongly correlated with the levels of total nitrogen, organic carbon, and available potassium in the soil.
The results of these studies suggest that multi-cropping is important for soil conservation and sustainable agro-ecosystems.

Conflicts of Interest:
The authors declare that they have no conflict of interest.