Лесоведение, лесоводство, лесные культуры, агролесомелиорация, озеленение, лесная пирология и таксация
ПРИРОДООБУСТРОЙСТВО 2' 2023
Оригинальная статья УДК 630.2*:528.94
DOI: 10.26897/1997-6011-2023-2-120-123
ПОСТРОЕНИЕ УНИВЕРСАЛЬНОЙ МОДЕЛИ ОЦЕНКИ СТЕПЕНИ ИЗМЕНЕНИЯ ПЛОЩАДИ ЛЕСНЫХ ПОКРЫТИЙ
Тахмазов Табриз, докторант
Thmzovtbriz66@gmail.com
Национальное Аэрокосмическое Агентство, г. Баку, Азербайджанская Республика
Аннотация. Статья посвящена построению универсальной модели скорости изменения площади лесных покрытий. Предложен вариант показателя скорости обезлесения, в котором площадь лесного покрытия в конце рассматриваемого интервала времени является функцией его длительности. Отмечено, что введение указанной функциональной зависимости приводит к универсальности вводимой оценки скорости процесса, так как появляется возможность оценивать не только процесс обезлесения, но и процессы лесопосадки. Проанализирован показатель усредненной по рассматриваемому интервалу времени скорости процесса обезлесения. Сформирована и решена задача вариационной оптимизации введенного усредненного показателя с учетом вновь введенного ограничительного условия применительно к вышеотмеченной функциональной зависимости. Вычислена конкретная оценка нововведенного усредненного показателя в виде смещаемой оценки по оси времени при наличии жесткой связи между конечными точками рассматриваемого временного интервала усреднения.
Ключевые слова: обезлесение, оптимизация, математическая модель, вырубка леса, биоразнообразие
Формат цитирования: Тахмазов Табриз. Построение универсальной модели оценки степени изменения площади лесных покрытий // Природообустройство. 2023. № 2. С. 120-123. Б01:
CONSTRUCTION OF A UNIVERSAL MODEL FOR ASSESSING THE DEGREE OF CHANGE IN THE AREA OF FOREST COVERINGS
Tahmazov Tabriz, doctoral student
Thmzovtbriz66@gmail.com
National Aerospace Agency, Baku, Republic of Azerbaijan
Annotation. The article is devoted to the construction of a universal model of the rate of change in the area of forest coverings. A variant of the deforestation rate indicator is proposed, in which the area of forest cover at the end of the time interval under consideration is a function of the duration of the specified time interval. It is noted that the introduction of this functional dependence leads to the universality of the introduced assessment of the speed of the process, since it becomes possible to evaluate not only the deforestation process, but also the processes of forest planting. The indicator of the speed of the deforestation process averaged over the considered time interval is analyzed. The problem of variational optimization of the introduced average indicator has been formed and solved, taking into account the newly introduced restrictive condition in relation to the above-mentioned functional dependence. A specific estimate of the new averaged indicator is calculated in the form of a shifted estimate along the time axis in the presence of a rigid connection between the endpoints of the considered averaging time interval.
Keywords: disafforestation, optimization, mathematical model, deforestation, biodiversity
Format of citation: Tahmazov Tabriz. Construction of a universal model for assessing the degree of change in the area of forest coverings // Prirodoobustrojstvo. 2023. No. 2. S. 120-123. DOI: 10.26897/1997-6011-2023-2-120-123.
10.26897/1997-6011-2023-2-120-123.
© Тахмазов Табриз, 2023
Original article
2i
&
Forest science, forestry, forest crops, agricultural afforestation, landscaping, forest pyrology and taxation
Введение. Как отмечается в исследованиях [1], обезлесение, или вырубка лесов (deforestation), является одним из важнейших вопросов в период глобального изменения климата. Обезлесение приводит к таким негативным событиям, как уменьшение биоразнообразия и увеличение эмиссии парниковых газов. При этом основные причины, приводящие к обезлесению, подразделяются на непосредственные и косвенные. К непосредственным причинам относятся в основном прямая вырубка леса и природные бедствия. К косвенным причинам относятся те факторы, которые формируются посредством социально-экономических процессов. Вместе с тем [2] обезлесение является двунаправленным процессом, то есть факторы, приводящие к обезлесению, также способны содействовать развитию лесов. Например, развитие сельскохозяйственных технологий может привести как к обезлесению, так и к развитию лесов, или, например, охрана государственного управления, введя прозрачные нормы, может улучшить сохраняемость лесов и в то же время, снижая штрафы за незаконную вырубку лесов, содействовать этому процессу.
Процесс обезлесения происходит при выполнении условия [3]:
А, + DV, > F • и, + IP, (1)
где А, - чистая стоимость сельского хозяйства; DVi - стоимость обезлесения; F i - чистая стоимость лесоводства; IPi - платежи по стимулированию; И, - коэффициент пропорциональности, И = 1,5 [3].
Как следует из выражения (1), наличие обезлесения определяется путем сравнения чистой стоимости сельского хозяйства и лесоводства. При этом в левую сторону (1) добавлена выгода ввиду обезлесения и уменьшения углеродных запасов, а в правую сторону - платежи по укреплению лесоводства.
Существует модель Чапмана-Ричардса [4], определяющая процесс изменения площади лесного покрытия в зависимости от демографических и экологических факторов:
^ = b • yb2 dp
- b3 • У
(2)
dy
где у - процент территории, не покрытой лесом;--произ-
йр
водная у по плотности населения; Ъ1,Ъ2,ЪЯ - параметры процесса.
При этом показатель у определяется как
5 - у
У = ■
S
• 100%
(3)
PRIRODOOBUSTROJSTVO 2' 2023
Интегрирование (3) дает [4]: y = a0 • (1 - a1 • exp [-a2PJ)"3
где P - плотность населения; a0, a^ a2, a4 - параметры процесса.
Как следует из выражения (4), рост плотности населения приводит к увеличению площади территории, не покрытой лесом.
Как видим из вышеприведенных модельных построений, модели (1)^(4) в основном направлены на оценку процесса обезлесения, и их возможности для процессов развития лесов детально не изучены.
С учетом вышесказанного имеет смысл разработать новую модель процесса изменения площади лесных территорий - в частности, через показатель скорости (степени) таких изменений. Далее рассмотрим возможность построения такой универсальной модели на базе известного выражения для оценки степени (скорости) процесса обезлесения.
Предлагаемый метод. Существует закон [5] совместной заинтересованности (Compoundlnterest Law), в соответствии с которым среднегодичный уровень изменения площади лесных покрытий определяется как
А2 А
где А1 и А2 - площади лесных покрытий во временные мо-
q = 1^ I -1
1_
¡2 -h
(6)
На основе (5) можно вывести формулу [6-10] для вычисления скорости обезлесения:
1 , (А
R =
• ln I ^
¡2 ¡1 V А-
(6)
Формула (6) интуитивно более понятна [11], чем формула (5). При этом оценки, полученные по формуле (6), всегда выше, чем оценки по формуле (5).
Детальные вычисления скорости обезлесения, например, приведены в работе [12]. Рассмотрим предлагаемую модификацию формулы (6) для придания ей более высокой универсальности.
Модификация показателя (6) осуществляется по следующим пунктам.
1) Выражение (6) приводится к виду
R =1 • ln А т А,
где т = t2 -1
где S - общая площадь территории.
2) Вводится на рассмотрение функция А2 = f (Т)
(7)
менты t2 и t.
Лесоведение, лесоводство, лесные культуры, агролесомелиорация, озеленение, лесная пирология и таксация
3) Принимается следующее ограничительное условие:
j f (г) dt = C; C = const
(9)
i -T1
ln-
Г„„„ - tm;„ г г A1
j 1ln dt
г
(10)
d
1ln ^^ dt + ÄA2 (г)
г A,
dA2 (г) Из (12) получаем 1
= 0
г • A2(г)
+ Л = 0
(12)
(13)
Рис. Геометрическая интерпретация условия (9) в отношении введенной функции f (г) Fig. Geometrical interpretation of condition (9) as to the introduced function f (г)
ПРИРОДООБУСТРОЙСТВО 2' 2023
Из (13) находим
1
Отметим, что условие (9) определяет универсальность введенной функции f (г), т.е. показатель А2 на интервале 0 + гтах может как расти от нуля до А2тах, так и уменьшатся в обратном порядке (рис. 1).
4) Вводится на рассмотрение среднеинте-гральная величина Я, то есть
Л -) = -ТА
С учетом (9) и (14) имеем
- 1 Т
Из (15) находим:
- j - dг = C
1 I г
1 =--ln1 max
C
Из выражений (15) и (16) находим C
5) Решается оптимизационная задача вычисления экстремала функционала (10) с учетом условия (9), при котором Дф инт достигает экстремальной величины.
Представим решение сформулированной оптимизационной задачи. С учетом выражений (9) и (10) сформулируем функционал F безусловной вариационной оптимизации:
F=-1-1 —lnd-+A Тл2-)d--C (11)
т -т ■ т a j
max mm т— 1 —mm
где X - множитель Лагранжа.
Согласно условию уравнения Эйлера [12] решение задачи должно отвечать условию:
A2 (г) = -
г1п
(15)
(16)
(17)
Таким образом, (17) представляет собой решение оптимизационной задачи (11).
Чтобы определить тип экстремума (минимум или максимум), вычислим производную (13) относительно А^т) и определим, что знак получаемого результата является минусовым, то есть при решении (17) достигается максимум функционала (11).
Вычислим значение указанного максимума функционала (10). С учетом (10) и (17) получим:
RCP.HHT =тТ ] 1ln
C
At s
где A т = -2 - — .
Примем, что
г
ln I max
i г
' ln max
dг
A
= C1;
C
CA
= C
(19)
С учетом (19) выражение (18) приведем
к виду:
R
1
Аг • C
T^ln ^ dz
J г г
(20)
В возражении (20) осуществим замену переменной интегрирования г на х = —. В этом слуг
чае (20) приобрешает следующий вид
1 Г'г
R
x ln xdx
Аг • C2
2 C2' min
С учетом известной формулы [12]
C2/
Cn/г
с............-2 x2
I x ln xdx =— ln x--
J 2 4
(21)
(22)
Тахмазов Табриз. Построение универсальной модели оценки степени изменения площади лесных покрытий
Forest science, forestry, forest crops, agricultural afforestation, landscaping, forest pyrology and taxation
PRIRODOOBUSTROJSTVO 2' 2023
С учетом (21) и (22) получим следующую формулу для вычисления #с„инт :
R„,
а
1-2ln-
T 1 t ■ mm -1 -TI mim 1-2]*-^ t 1— max -1
4-At-C2
(23)
Из выражения (23), например, при ттах = = е ■ ттП где е = 2,73 можно получить следующую оценку:
е2 (1 - 21п 1,73) + 21п1,73 +1
R„,
4z 2
или
R • z2
ср.инт max
= 0,34
(24)
Как следует из оценки (24), с ростом ттах
оценка Rс инт будет убывать по квадратичному
' т
закону учитывая исходное условие расчета тах = е,
т
тт
то есть (24) является смещаемой оценкой R
ср.инт
по оси т.
Заключение
Проведены модельные исследования известной формулы для оценки скорости обезлесения. Предложен модифицированный вариант
показателя скорости обезлесения, в котором площадь лесного покрытия в конце рассматриваемого интервала времени является функцией его длительности. Отмечено, что введение указанной функциональной зависимости приводит к универсальности вводимой оценки скорости процесса, так как появляется возможность оценивать не только процесс обезлесения, но и процессы лесопосадки.
Предложен показатель усредненной по рассматриваемому интервалу времени скорости процесса обезлесения. Сформирована и решена задача вариационной оптимизации введенного усредненного показателя с учетом вновь введенного ограничительного условия применительно к вы-шеотмеченной функциональной зависимости. Вычислена конкретная оценка нововведенного усредненного показателя в виде смещаемой оценки по оси времени при наличии жесткой связи между конечными точками рассматриваемого временного интервала усреднения.
Показано, что нововведенный усредненный показатель скорости процесса имеет максимум, который убывает по обратному квадратичному закону в зависимости от величины конечной точки рассматриваемого интервала времени.
Список использованных источников / References
1. Indarto J., Mutaqin D.J. An overview of theoretical and empirical studies on deforestation // Journal of International Development and Cooperation. Vol. 22. No 1& No 2. 2016. Pp. 107-120
2. Robinson B.E., Holland M.B., Naughton-Tre-ves L. Does secure land tenure save forests& A meta-ana-lysis of the relationship between land tenure and tropical deforestation // Global Environmental Change. 29(2014). Pp. 281-293
3. Kindermann G.E., Obersteiner M., Rametstei-
ner E., McCallum I. Carbon balance and management // Biomed Central.
4. Marzoli W.A. Appendix 3: Modeling deforestation baselines using the forest aera change model for the Calak-mul and Meseta Purepecha regions in Mexico.
5. Fearnside P.M. Deforestation in Brazilian Amazonia // Forest resources assessment 1990. Global Synthesis. FAO. Rome.
6. Puyravaud J.P. Standardizing the calculation of the annual rate of deforestation // Forest ecology and management 177. 2003. Pp. 593-596
7. Menon S., Bawa K.S. Applications of geograpyic information systems, remote sensing and a landscape ecology
Критерии авторства
Тахмазов Табриз выполнил практические и теоретические исследования, на основании которых провел обобщение и написал рукопись, имеет на статью авторское право и несет ответственность за плагиат. Статья поступила в редакцию 27.01.2023 Одобрена после рецензирования 01.03.2023 Принята к публикации 01.03.2023
approach to biodiversity conservation in the Western Ghats // Curr. Sci. vol. 73. Pp. 134-145. 1997
8. Ramesh B.R., Menon S., Bawa K.S. A vegetation based approach to biodiversity gap analysis in the Agastya-malai region, Western Ghats, India // Ambio. Vol. 26. Pp. 536-539. 1997.
9. Puyravaud J.P. Standardizing the calculation of the annual rate of deforestation // Forest ecology and Management. Vol. 177(1-3). Pp. 593-596. 2003
10. Prasad N.S. Conservation planning for the western Ghats of Kerala: assessment of habitat loss and degradation // Curr. Sci. Vol. 75. Pp. 228-235. 1998
11. Rijal S., Saleh M.B., Jaya N.S., Tiryana T. Deforestation profile of regency level in Sumatra // International Journal of sciences basic and applied research. ISSN2307-4531
12. Rijal S., Vicayana D.P., Barkey RA., Aki-muddin I., Mukhlisa A.N. Spatial patterns analysis of deforestation in palopo municipality and east luwu regency // The 4th International Conference of Indonesian Scoiety for remote Sensing. IOP Conf. Series: Earth and Environmental Science 280. 2019. Doi: 10.1088/ 1755-1315/280/1/012028
Criteria of authorship
Tahmazov Tabriz carried out theoretical and practical studies, on the basis of which he generalized and wrote the manuscript. He has a copyright on the article and is responsible for plagiarism.
The article was submitted to the editorial office 27.01.2023 Approved after reviewing 01.03.2023 Accepted for publication 01.03.2023