A Mathematical Model of Regional Socio-Economic Development of the Russian Arctic Zone

This work touches upon the tasks of describing regional socio-economic development. The nature of the considered problem indicates the fact that the only relevant tool here is mathematical modeling. In this paper, the application of mathematical modeling is considered for the problem of managing regional development. The results of calculations based on the regional dynamic model that passes through the hierarchy of instabilities (the correspondence of the same stationary points of the model to different parameters) are presented. These instabilities lead to increasingly complex structures.


Introduction
This paper deals with the tasks of describing the socio-economic development of regions on the basis of constructed descriptions.For the development of the regional economy and tourism, the interested reader can refer to the input-output analysis with the resource of the tourism industry described in [1].Another example of the development of a region is the modeling of a "green" economy: a dynamic model of the "green" economy, a general extended map, and a reduced-parametric cognitive map of the development of the "green" economy in the economic space of the region [2].The nature of the analyzed problem indicates the fact that the only applicable tool here is mathematical modeling.Indeed, while studying various kinds of natural-science problems and developing technical systems, machines and technologies, one can rely both on the physical experiment and on the centuries-old experience of constructing models in physics, mechanics, and other branches of natural science.The same approach is barely applicable in the sphere of economic problems, not to mention the problems of effective management in tasks of that sort, where this "experimental" approach is impossible.The conceptual view of applications of mathematical modeling in the considered class of problems [3] is very instructive, since such modeling incorporates all mathematical modeling technologies mentioned by A.A. Samarskii almost 40 years ago.To calculate the parameters of investment attractiveness of the Arctic region, it is also possible to use the mathematical modeling considered by the authors in [4] for poorly researched hydrocarbon objects.Moreover, studying the problem on the level of modern advanced approaches shows that, in the class of considered problems, the implementation of mathematical modeling should be based on computer technologies for the simple reason that the object, or the problem itself, is a unique subject for research, not only as systems of "big data" but also as multi-scale systems containing elements of uncertainties [5][6][7][8].
The application of mathematical modeling approaches is considered with the example of a model problem of regional development management.The core of the model is the state function describing regional development and represented by the vector function s(t) with the components [x(t), y(t), z(t)], where x(t) is the population of the region, y(t) is the number of jobs in the real sector of the regional economy, and z(t) is the indicator of energy supply in the region.The vector u(t) is the controlling vector that determines the required values for the state vector.The controlling vector has m-components [u 1 (t), u 2 (t), . . . ,u m (t)], each of which characterize one or another controlling factor.For example, u 1 (t) is the coefficient of demographic activity, u 2 (t) is the coefficient of people's anti-motivation to childbearing, u 3 (t) is the energy-supply coefficient of the region, etc.

The Deterministic Model of Management of Socio-Economic Development of the Region
The models of nonlinear dynamics represent a mathematical toolset that enables one to describe and analyze the development of various dynamic systems of different complexity levels-from mechanical systems up to social ones-based on objective laws [6,[9][10][11][12].In [6], the classification of the considered systems and models according to the complexity degree and to the possible expected results of modeling is presented (see Table 1).The concept of managing the socio-economic development of the region is based on a synthesized model of the management process (the inverse problem of the management process).The model introduces three main indicators of the region: s 1 (t) is the regional population, s 2 (t) is number of jobs in the real sector of the regional economy, and s 3 (t) is the indicator of the regional energy supply [13-17].
The system-forming mathematical base of the dynamic model is a system of differential equations describing the changes of the three key indicators introduced.Through the coefficients (there are nine of them all in all), the mechanisms for implementing regional policies directed to achieve the goal of public administration are formed [18].
The development of the system of differential equations can be built on the natural-scientific and the socio-economic interpretation of the meaning of regional performance indicators s i (i = 1, 2, 3).
We will present the construction of the regional development model in a concise form.We will start with the fact that the derivative of the population of the region ds 1 /dt is the rate of change in the population and it is naturally connected both with the population and with the number of jobs the real economy of the region s 2 (the indicator of economic development) and with the indicator of energy supply in the region s 3 .The rate of change for s 1 is proportional to the regional population itself, i.e. the higher the population is, the greater the growth rate is: where u 1 is the demographic activity coefficient.We will determine the influence of indicators s 2 and s 3 on the population growth rate.The number of jobs in the real sector of the regional economy is determined with the minimum number of workers in the production sector required to produce a certain assortment of goods and services [19][20][21][22][23][24][25][26][27][28].
For a given value of s 1 the number of jobs the real economy s 2 will reduce the population growth rate by u 2 s 1 s 2 , where u 2 is the coefficient of people's anti-motivation to childbearing.Additionally, for a given value of s 1 , the energy supply s 3 will increase the regional population growth rate by the amount u 3 s 1 s 3 , where u 3 is the ratio of energy supply of the region.In other words, the more energy that enters the region there is, the higher the regional population growth rate is [29,30].
Thus, the first differential equation can be written in the following form: Turning now to the change in the number of jobs in the real sector of the region s 2 , we note that the derivative ds 2 /dt is the rate of change of the economic development indicator.At the same time, the rate of change of the economic development indicator is inversely proportional to the number of jobs in the real sector; that is, the more jobs created in the real sector there are, the more difficult it is to increase the number of job opportunities: where u 4 is the coefficient of people's interest in economic development.For a given value of s 2 , if people are interested in the economic development, an increase in the regional population s 1 will increase the growth rate of economic development in the real sector by the amount u 5 s 1 s 2 , where u 5 is the coefficient of people's interest in economic development.The manifestation of such a property is natural, since it makes up the basis of the society's self-preservation.The reverse manifestation of this property will lead to the self-destruction of the society.For a given value of s 2 , the energy supply of the region s 3 will increase the economic growth rate s 2 by amount u 6 s 2 s 3 , where u 6 is the coefficient of energy supply of the workplaces.That is, the more energy there is for the development of the real sector of the economy, the higher the economic growth is [29].Thus, the differential equation for the indicator s 2 takes the following form: Turning to the indicator of energy supply of the region s 3 , we point out that its derivative ds 3 /dt is the rate of change in the region's energy supply.It is quite obvious that the rate of change in energy consumption is proportional to the amount of energy consumed; that is, the more energy consumed in the society there is, the higher the growth rate of the energy supply is: where u 7 is the development factor of the regional energy supply.
For a given value of s 3 , the growth of regional population s 1 will reduce the rate of change in the energy consumption by amount u 8 s 1 s 3 , where u 8 is the conformity ratio of the population with the energy supply.As the population increases, the rate of change in energy supply decreases.For a given value s 3 , an increase in the number of jobs in the real sector s 2 will reduce the rate of change in energy supply by amount u 9 s 2 s 3 , where u 9 is the conformity ratio of the economic development and the energy supply.In other words, the higher the growth rate of the real sector of the regional economy is, the lower the energy supply per workplace is [29].
Hence, the differential equation for the s 3 indicator can be written in the following form: Thus, the management model of sustainable socio-economic development of the region has the following form (all the indicators are presented in the form of relative, dimensionless values): where the demographic indicator is the value , where s * 1 (t) is the population at time t, and s * 1 (t 0 ) is the population at the initial time point t 0 .The economic development indicator is , where s * 2 (t) is the number of jobs of the real sector of economy at time t, and s * 2 (t 0 ) is the number of jobs of the real sector of economy at time t 0 .Eventually, the energy supply indicator is , where s * 3 (t) is the current energy supply indicator at time t, and s * 3 (t 0 ) is the current energy supply indicator at time t 0 .The coefficients u i (i = 1, 2, . . ., 9) in the system of Equation ( 7) are the following factors: u 1 is the demographic activity coefficient; u 2 is the coefficient of people's anti-motivation to childbearing; u 3 is the energy supply coefficient; u 4 is the coefficient of people's interest in economic development; u 5 is the coefficient of the real sector economic development; u 6 is the coefficient of energy supply per workplace; u 7 is the energy supply coefficient of the region; u 8 is the conformity ratio of the population with the energy supply; u 9 is the conformity ratio of the economic development with the energy supply.

The Stochastic Model of Managing the Socio-Economic Development of the Region
The nine coefficients involved in the system of differential equations are the components of the vector of public administration: u = [u 0 1 , u 0 2 , u 0 3 , u 0 4 , u 0 5 , u 0 6 , u 0 7 , u 0 8 , u 0 9 ].The state of the region vector for the model is Such an approach to describing the regional management makes it possible to base the management concept on the formulation of the following optimal management problem [30][31][32].
Let us formulate the management problem under the assumption that the controlling vector u is sought on the number set of its components.In this case we are dealing with the problem of parametric optimization.Thus, it is required to determine the vector of the control parameters u, under which the state of the region vector s(t) = [s 1 (t), s 2 (t), s 3 (t)] can be described by the system of Equation (7)  under the following conditions: • s 1 (t) is the population of the region, which is not defined (the indicator is free); • s 2 (t) = y * is fixed (one gives the fixed number of jobs in the real sector of the economy); • the minimum of the variable s 3 is given: The system of Equation ( 7) [33][34][35][36][37] can only partially describe the development dynamics of the region, because is deterministic and therefore does not take into account fluctuations that are constantly arising in the system.These fluctuations are formed by different types of markets (goods, capital, and labor) against the background of competition and limited resources.For these reasons, at the next stage of the research, one moves on to implementation of statistical methods, which also enables more accurate forecasting and the creation of possible future scenarios.However, statistical methods generally require a significant amount of statistical data and some additional studies to identify a set of influencing factors, with a subsequent search for dependencies [38].From our point of view, the main problem is the determination of the fullest possible number of factors that are used for the further analysis.This will lead to a significant change in the type of model used and will increase the number of parameters used.In addition, the problem of the availability of statistical data on all the possible parameters still remains [39].
Another important point is the evaluation of the values of integral economic indicators, for example, the generalized index of industrial growth.These indicators are especially nonlinearly dependent on the solution of the system of Equation (7).Determining the values of such indicators requires a large-scale simulation of the trajectories space, whereas the modeling scape depends on the desired accuracy.
In order to simplify the analysis, we propose limiting the number of parameters and ignoring the entire set of parameter space.The model will be constructed only on the basis of the factors selected by the experts.It is proposed that a factor of the stochastic nature of the described processes is introduced into the system of Equation ( 7) for the purpose of taking the stochastic effect into consideration.In this case, due to the transition to a system of stochastic differential equations with a small number of parameters, it is possible to model the behavior of the system not only for a certain set of statistically distinguished values of the parameters.This will also allow us to calculate the integral parameters mentioned above.For solving the problem, the so-called weak methods are used for short time periods, and strong methods (for example, Monte Carlo method) are used for larger time periods.The latter require considerable computational resources, since, in spite of the relative simplicity of the computations performed for calculating of a separate trajectory, a considerable number of independent experiments must be carried out.The required number of simulated trajectories depends not only on the time frame that serves the basis for the simulation but also on the desired accuracy and the type of functionality, in the case of calculating the integrated indicators.
The authors proposed a modified system of Equation ( 7), which takes into account possible random fluctuations in the calculated parameters.
where in the second equation an additional term is introduced.The term shows that the higher the population is, the higher the competition on the labor market will be; moreover, the more incentives to stick to the job there are, the lower the personnel turnover will be.In the system of Equation ( 9), a, b, c, d, e, f are calibration coefficients that can be obtained from statistical data processing.The traditional approach to solving the system of Equation ( 9) is the use of powerful computational resources to determine a certain sample from the entire space of trajectories.It is necessary to emphasize that the solution of the above-stated problems at the current level involves a high amount of computational experiments.Preliminary evaluation of the needed amount of calculation clearly shows that the capabilities of even powerful workstations are not sufficient.
In contrast to the behavior of the previously introduced deterministic model, the stochastic character of real systems plays an important role in practice.As it is one of the factors leading the system to non-deterministic behavior, one can consider, for example, an economic crisis that leads to the loss of stability, bifurcation, and the emergence of new stable states of the system.
The presence of fluctuations, as a rule, leads to errors in estimating the real parameters of the system and, eventually, to forecasting errors for the observed phenomena.In order to demonstrate the behavioral difference of deterministic and stochastic systems, we added only one stochastic term to the second equation of the original system (Equation ( 9)).Having introduced the term, we tried to relate the magnitude of fluctuations in the number of jobs to the available working-age population.It was also assumed that the stochastic differential is an Ito differential [7,8,40].
In this case, as was stated above, one looks for the minimum of the s 3 parameter: We present here the results of calculations for the deterministic and the stochastic models.The solution of the deterministic model is shown in Figure 1.where in the second equation an additional term is introduced.The term shows that the higher the population is, the higher the competition on the labor market will be; moreover, the more incentives to stick to the job there are, the lower the personnel turnover will be.In the system of Equation ( 9), a, b, c, d, e, f are calibration coefficients that can be obtained from statistical data processing.The traditional approach to solving the system of Equation ( 9) is the use of powerful computational resources to determine a certain sample from the entire space of trajectories.It is necessary to emphasize that the solution of the above-stated problems at the current level involves a high amount of computational experiments.Preliminary evaluation of the needed amount of calculation clearly shows that the capabilities of even powerful workstations are not sufficient.
In contrast to the behavior of the previously introduced deterministic model, the stochastic character of real systems plays an important role in practice.As it is one of the factors leading the system to non-deterministic behavior, one can consider, for example, an economic crisis that leads to the loss of stability, bifurcation, and the emergence of new stable states of the system.
The presence of fluctuations, as a rule, leads to errors in estimating the real parameters of the system and, eventually, to forecasting errors for the observed phenomena.In order to demonstrate the behavioral difference of deterministic and stochastic systems, we added only one stochastic term to the second equation of the original system (Equation ( 9)).Having introduced the term, we tried to relate the magnitude of fluctuations in the number of jobs to the available working-age population.It was also assumed that the stochastic differential is an Ito differential [7,8,40].
The state vector of the region s(t) = [x(t), y(t), z(t)] was defined as the solution of Equations ( 7) and ( 9    It should be noted that, due to the properties of the Ito stochastic integral [7,8,40], the averaged trajectory of the solution for the system of Equation ( 9) will coincide with the trajectory of the deterministic system (see Figure 1).However, it is necessary to pay attention to Figure 2, where the behavior of one fixed trajectory of the solution of the proposed stochastic system is shown.The comparison of the trajectories in Figures 1 and 2 clearly shows that, during significant time intervals, the trajectories differ substantially.Thus, there is a need to evaluate additional parameters of the system in order to describe its possible behavior more accurately.For example, it would be useful to evaluate decision parameters such as the evolution of variance, confidence intervals, and the spatial distribution of the trajectories.It should be noted that, due to the properties of the Ito stochastic integral [7,8,40], the averaged trajectory of the solution for the system of Equation ( 9) will coincide with the trajectory of the deterministic system (see Figure 1).However, it is necessary to pay attention to Figure 2, where the behavior of one fixed trajectory of the solution of the proposed stochastic system is shown.The comparison of the trajectories in Figures 1 and 2 clearly shows that, during significant time intervals, the trajectories differ substantially.Thus, there is a need to evaluate additional parameters of the system in order to describe its possible behavior more accurately.For example, it would be useful to evaluate decision parameters such as the evolution of variance, confidence intervals, and the spatial distribution of the trajectories.

Figure 1 .
Figure 1.Solution for the deterministic model.

Figure 1 .
Figure 1.Solution for the deterministic model.

Figure 2 .
Figure 2. Solution for the stochastic model.

Table 1 .
Systems and models according to the complexity degree.