Scientific journal
Modern problems of science and education
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 0,682

DISTRIBUTION PHYTOMASS HERBS AND TREES IN FOREST PHYTOCOENOSIS

Мазуркин П.М., Михайлова С.И.
Недостатком известных способов испытания фитоценоза является раздельное изучение древесных растений и травы. Эта раздельность анализа фитомассы отдельных компонент растительного сообщества не позволяет выявлять закономерности распределения биологического вещества по частям растительного покрова и фитоценоза в целом. Предложены метод комплексного анализа результатов измерений фитомассы, а также модели годичного прироста и опада биологического вещества в различных растениях. Ключевые слова: трава, деревья, фитоценозы, распределение фитомассы.
The disadvantage of the known methods of testing phytocoenosis is a separate study of woody plants and grasses. This separate analysis of individual components of phytomass plant communities is not possible to identify characteristics of the distribution of biological materials in portions of vegetation and phytocoenosis in general. A method for integrated analysis of the results of measurements of phytomass, and model annual increment and litterfall of biological substances in various plants. Key words: grass, trees, phytocoenoses, distribution of phytomass.

Известны способы испытания деревьев и других частей фитоценоза на содержание химических элементов, применяемые для изучения биологического круговорота химических элементов [4, с. 10-18].

Недостатком их является раздельное изучение древесных растений и травы. Эта раздельность анализа фитомассы отдельных компонент растительного сообщества не позволяет выявлять статистические закономерности распределения биологического (живого по В.И. Вернадскому) вещества по частям растительного покрова и фитоценоза в целом.

Цель статьи - повышение комплексности анализа результатов измерений фитомассы, а также моделей годичного прироста и опада биологического вещества в различных видах и частях растений.

Для сопоставления результатов измерения фитомассы растительного сообщества, по его видам растений и их структурным частям, значения фактической [4] фитомассы этих частей распределяют по кодовой шкале, ориентированной по слоям фитоценоза, размещенных «снизу вверх»: 1 - слой корней травы; 2 - мелкие корни древесных растений; 3 - крупные корни древесных растений; 4 - стволы древесных растений; 5 - крупные ветви древесных растений; 6 - мелкие ветви древесных растений; 7 - листва (хвоя) древесных растений; 8 - слой надземной части травы.

Измеренные значения фитомассы по структурным частям располагают в виде таблицы, причем по строкам располагают слои фитоценоза. По столбцам размещают виды растений или их растительные комплексы, а затем по всем строкам и столбцам вычисляют суммарные значения фитомассы, причем сумма по строке будет показывать фитомассу принятого слоя фитоценоза, а сумма по столбцу будет равна фитомассе отдельного вида или комплекса видов растений. Такие же таблицы составляют по годичному приросту и опаду биологического вещества фитоценоза.

Статистические закономерности распределения фитомассы, годичного её прироста и годичного опада, по слоям фитоценоза 1 - 8, выявляются по обобщенной формуле статистической закономерности:

f f f

f f (1)

где q - содержание фитомассы (анало-гично годичный прирост и годичный опад фитомассы) в слоях фитоценоза по структурным частям, ц/га;

q1 - закон гибели, являющийся частным случаем биотехнического закона, предложенного проф. П.М. Мазуркиным [1 - 3], и показывающий влияние почвы и корневой системы на фитомассу в частях растений, ц/га;

q2 - биотехнический закон [1 - 3], показывающий влияние древесины, то есть омертвевших клеточных структур, на фитомассу, ц/га;

q3 - биотехнический закон, показывающий влияние процессов фотосинтеза и образования в зеленых частях растений ассимилятов, ц/га;

q4 - волновая функция с амплитудой в виде закона аллометричского роста, которая показывает способность фитоценоза к адаптации к растительным условиям и при её высокой значимости (по отношению к другим составляющим модели) приводящая к коллапсу фитоценоза, ц/га;

i - номер слоя растительного сообщества (фитоценоза), то есть место или порядок распределения структурных частей у различных видов растений по принципу «снизу вверх»;

a0, a1, ...a12 - параметры статистической закономерности (1), получаемые идентификацией закономерности (1) относительно фактически измеренных значений фитомассы как конкретные числа, причем чаще всего могут быть выполнены для леса упрощающие условия: a1 =1; a6 =1; a10 =1;

a11ia13 - половина амплитуды колебательного возмущения в структуре фитоценоза по слоям (экологическим горизонтам) в виде закона показательного роста, ц/га;

a13 - половина периода у колебательного возмущения структуры фитоценоза по его слоям;

a14 - смещение (сдвиг) начала волны возмущения структуры фитоценоза относительно начала координат, находящегося в почве (при i = 0).

Сущность предлагаемого способа заключается в том, что фитоценоз рассматривается как популяция из многовидовых растений, которая приспосабливается к внешним условиям волновыми возмущениями отдельных структурных слоев и частей в общем симбиотическом росте и развитии.

При этом деревья совместно с окружающей их травой и более мелкими древесными растениями в виде кустарников и кустарничков приспосабливаются к пространству обитания через отклик на параметры внешней среды, то есть на изменения экологической (шире биологической) обстановки. При этом любой фитоценоз, как показали примеры моделирования, приспосабливается по биотехническому закону [1 - 3] и колебательным возмущением адаптируется к питательным веществам из окружающей лесной среды, в том числе симбиотическим образом отдельные растения помогают друг другу в фитоценозе.

Примеры. Все примеры для статистического моделирования изменения содержания фитомассы в различных растительных комплексах были взяты из книги [4, c.103-114, табл. 9-11, табл. 14-16, табл. 19-21].

1. Верховое осоково-сфагновое болото. Фитомасса в этом сложном растительном сообществе состоит из слоев и видов растений по табл. 1.

Для системы болотных растений была получена статистическая модель (табл. 2) фитомассы в виде формулы из двухсоставляющих:

ff (2)

Таблица 1. Содержание фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование

слоя растительного сообщества

Номер слоя

    i

Вид растения

 

Всего

Сосна

Кустарнички

Осока

Мох

Слой надземной части травы

8

-

-

0.5

27.0

27.5

Листва (хвоя) древесных растений

7

13.0

10.0

-

-

23.0

Мелкие ветви древесных растений

6

24.0

5.5

-

-

29.5

Крупные ветви древесных растений

5

-

-

-

100.0

100.0

Стволы древесных растений

4

150.0

-

-

-

150.0

Крупные корни древесных растений

3

-

-

-

-

-

Мелкие корни древесных растений

2

27.0

12.0

-

-

39.0

Слой корней травы

1

-

-

1.0

-

1.0

Итого:

214.0

27.5

1.5

127.0

370.0

Таблица 2. Содержание фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование

части растений

Слой

   i

Факт

 f

Расчетные значения (2)

Составляющие (2)

     q

f

f

         q1

      q2

надземная часть травы

листва (хвоя)

мелкие ветви

крупные ветви

стволы

крупные корни

мелкие корни

корни травы

8

7

6

5

4

3

2

1

27.5

23.0

29.5

100.0

150.0

-

39.0

1.0

27.5

20.4

39.4

91.9

153.0

135.9

37.9

0.6

-0.02

2.62

-0.01

8.09

-2.95

-

1.08

0.37

-0.07

11.39

-33.59

8.09

-1.97

-

2.77

37.00

2.72

11.00

36.36

91.10

152.79

135.85

37.92

0.63

24.80

9.38

3.06

0.81

0.16

0.02

0.00

0.00

Примечания: 1) максимальное значение относительной погрешности подчеркнуто, поэтому доверие к модели (2) будет не ниже 100 - 37.00 = 63,00 %; 2) максимальные значения составляющих набраны полужирным шрифтом.

Максимальная относительная погрешность зависит от характера распределения фитомассы по структурным частям фитоценоза. При этом количество видов растений, учитываемых во всех восьми экологических слоях фитоценоза, не ограничивается.

Годичный прирост Δq+ фитомассы характеризуется данными (табл. 3), по которым было получено статистическое уравнение (табл. 4) в виде

f
f (3)

Таблица 3. Годичный прирост фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование

слоя растительного сообщества

Номер слоя

i

Вид растения

 

Всего

Сосна

Кустарнички

Осока

Мох

Слой надземной части травы

8

-

-

0.50

27.00

27.50

Листва (хвоя) древесных растений

7

1.85

2.20

-

-

4.05

Мелкие ветви древесных растений

6

0.26

0.25

-

-

0.51

Крупные ветви древесных растений

5

-

-

-

-

-

Стволы древесных растений

4

1.60

-

-

-

1.60

Крупные корни древесных растений

3

-

-

-

-

-

Мелкие корни древесных растений

2

0.29

0.13

-

-

0.42

Слой корней травы

1

-

-

0.30

-

0.30

Итого:

4.00

2.58

0.80

27.00

34.38

Таблица 4. Годичный прирост фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование

части растений

Слой

i

Факт

f

Расчетные значения (3)

Составляющие (3)

q

f

f

q1

q2

надземная часть травы

листва (хвоя)

мелкие ветви

крупные ветви

стволы

крупные корни

мелкие корни

корни травы

8

7

6

5

4

3

2

1

27.50

4.05

0.51

-

1.60

-

0.42

0.30

27.50

4.00

0.68

0.19

1.60

6.21

0.42

0.29

-0.004

0.055

-0.172

-

0.000

-

-0.000

0.007

-0.01

1.36

-33.73

-

0.00

-

-0.05

2.33

0.00

0.00

0.00

0.04

1.56

6.19

0.39

0.00

27.50

4.00

0.68

0.15

0.04

0.02

0.03

0.29

Годичный опад Δq- фитомассы характеризуется данными (табл. 5), по которым было получено статистическое уравнение (табл. 6) в виде

f

f(4)

Сравнение фактических и расчетных значений содержания, прироста и опада фитомассы растений в верховом болоте показывает относительную простоту закономерностей (всего две составляющие). При этом в дальнейших исследованиях следует обратить внимание на точноть измерения фитомассы по слоям фитоценоза «крупные ветви» и «крупные корни».

2. Ельник сложный 83 года. Здесь можно выделить один растительный комплекс «трава + ель», то есть растительное сообщество характеризуется всего двумя видами растений. Фитомасса в этом растительном сообществе состоит из слоев и видов растений табл. 7.

Таблица 5. Годичный опад фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование слоя

растительного сообщества

Номер слоя

i

Вид растения

 

Всего

Сосна

Кустарнички

Осока

Мох

Слой надземной части травы

8

-

-

0.50

19.00

19.50

Листва (хвоя) древесных

растений

7

1.80

2.00

-

-

3.80

Мелкие ветви древесных растений

6

0.24

0.22

-

-

0.48

Крупные ветви древесных растений

5

-

-

-

-

-

Стволы древесных растений

4

0.75

-

-

-

0.75

Крупные корни древесных растений

3

-

-

-

-

-

Мелкие корни древесных растений

2

0.27

0.12

-

-

0.39

Слой корней травы

1

-

-

0.30

-

0.30

Итого:

3.06

2.34

0.80

19.00

25.22

Таблица 6. Годичный опад фитомассы в верховом сосново-сфагновом болоте, ц/га

Наименование

части растений

Слой

i

Факт

f

Расчетные значения (4)

Составляющие (4)

q

f

f

q1

q2

надземная часть травы

листва (хвоя)

мелкие ветви

крупные ветви

стволы

крупные корни

мелкие корни

корни травы

8

7

6

5

4

3

2

1

19.50

3.80

0.48

-

0.75

-

0.39

0.30

19.51

3.70

0.79

0.21

0.77

3.92

0.37

0.09

-0.007

0.097

-0.307

-

-0.018

-

0.019

0.206

-0.04

2.55

-63.96

-

-2.40

-

4.87

68.67

0.00

0.00

0.00

0.01

0.71

3.89

0.35

0.00

19.51

3.70

0.79

0.20

0.06

0.03

0.02

0.09

Таблица 7. Содержание фитомассы в ельнике сложном, ц/га

Наименование слоя

растительного сообщества

Номер слоя

i

Вид растения

 

Всего

Ель

Трава

Слой надземной части травы

8

-

6.30

6.30

Листва (хвоя) древесных растений

7

188.0

-

188.00

Мелкие ветви древесных растений

6

109.0

-

109.00

Крупные ветви древесных растений

5

229.0

-

229.00

Стволы древесных растений

4

2271.0

-

2271.00

Крупные корни древесных растений

3

732.0

-

732.00

Мелкие корни древесных растений

2

36.0

-

36.00

Слой корней травы

1

-

7.58

7.58

Итого:

3565.0

13.88

3578.88

Для ельника была получена статистическая модель (табл. 8) фитомассы в виде формулы из четырех составляющих (простая фитоценозная структура из двух видов растений приводит к их сложному их взаимодействию), то есть в виде сложной закономерности

f

f

f

f(5)

Формула (5) имеет нулевую погрешность или 100-процентную доверительную вероятность. Поэтому анализ отдельных составляющих модели, а также их сопоставление позволяют понять многое из жизни ельника.

Годичный прирост Δq+ фитомассы характеризуется данными (табл. 9), по которым было получено статистическое уравнение (табл. 10) в виде

f

f

f

f (6)

Таблица 8. Содержание фитомассы в ельнике сложном, ц/га

Слой

i

Факт

f

Расчетные значения (5)

Составляющие модели (5)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

6.30

188.00

109.00

229.00

2271.00

732.00

36.00

7.58

6.30

188.00

109.00

229.00

2271.00

732.00

36.00

7.58

-9.4е-12

6.7е-12

-2.1е-12

-8.5е-14

5.9е-12

-1.6е-12

4.3е-14

8.9е-16

-0.00

0.00

-0.00

-0.00

0.00

-0.00

0.00

0.00

82.89

56.22

38.13

25.86

17.54

11.89

8.07

5.47

0.00

0.00

0.24

89.87

2140.96

648.07

0.07

0.00

42.59

66.54

93.11

112.11

108.55

74.48

27.25

2.15

119.76

72.07

40.10

20.05

8.58

2.87

0.61

0.04

119.18

-65.24

22.48

-1.16

-3.95

2.44

-0.61

0.04

Таблица 9. Годичный прирост фитомассы в ельнике сложном, ц/га

Наименование

слоя растительного сообщества

Номер слоя i

Вид растения

Всего

Ель

Трава

Слой надземной части травы

8

-

6.30

6.30

Листва (хвоя) древесных растений

7

21.70

-

21.70

Мелкие ветви древесных растений

6

1.42

-

1.42

Крупные ветви древесных растений

5

2.88

-

2.88

Стволы древесных растений

4

29.14

-

29.14

Крупные корни древесных растений

3

9.55

-

9.55

Мелкие корни древесных растений

2

0.47

-

0.47

Слой корней травы

1

-

2.50

2.50

Итого:

65.16

8.80

73.96

Таблица 10. Годичный прирост фитомассы в ельнике сложном, ц/га

Слой

i

Факт

f

Расчетные значения (6)

Составляющие модели (6)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

6.30

21.70

1.42

2.88

29.14

9.55

0.47

2.50

6.30

21.70

1.42

2.88

29.14

9.55

0.47

2.50

5.7е-11

-2.0е-11

-1.8е-10

-5.9е-11

7.6е-12

-5.8е-12

-3.4е-11

-7.4е-12

0.00

-0.00

-0.00

-0.00

0.00

-0.00

-0.00

-0.00

4.48

4.55

4.62

4.70

4.77

4.85

4.93

5.01

0.00

0.00

0.00

1.18

23.92

6.21

0.00

0.00

0.00

0.00

0.01

0.05

0.35

1.70

4.42

2.51

33.44

17.59

8.38

3.48

1.19

0.30

0.04

0.00

1.82

17.15

-3.20

-2.95

0.79

0.19

-0.04

-0.00

В табл. 9, как и табл. 7, содержатся данные о значениях фитомассы по всем восьми экологическим слоям фитоценоза ельника. Это позволяет получить весьма точные статистические закономерности.

Годичный опад Δq- фитомассы характеризуется данными (табл. 11), по которым было получено статистическое уравнение (табл. 12) в виде

f

f

f

f(7)

Таблица 11. Годичный опад фитомассы в ельнике сложном, ц/га

Наименование

слоя растительного сообщества

Номер слоя i

Вид растения

Всего

Ель

Трава

Слой надземной части травы

8

-

6.30

6.30

Листва (хвоя) древесных растений

7

21.52

-

21.52

Мелкие ветви древесных растений

6

4.57

-

4.57

Крупные ветви древесных растений

5

3.11

-

3.11

Стволы древесных растений

4

15.25

-

15.25

Крупные корни древесных растений

3

4.83

-

4.83

Мелкие корни древесных растений

2

0.24

-

0.24

Слой корней травы

1

-

2.50

2.50

Итого:

49.52

8.80

58.32

Таблица 12. Годичный опад фитомассы в ельнике сложном, ц/га

Слой

i

Факт

f

Расчетные значения (7)

Составляющие модели (7)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

6.30

21.52

4.57

3.11

15.25

4.83

0.24

2.50

6.30

21.52

4.57

3.11

15.25

4.83

0.24

2.50

-2.5е-13

4.1е-13

1.6е-13

-6.1е-14

-3.0е-14

1.4е-14

4.7е-15

-4.4е-16

-0.00

0.00

0.00

-0.00

0.00

0.00

0.00

-0.00

4.88

4.77

4.66

4.56

4.46

4.35

4.26

4.17

0.11

1.05

6.65

22.87

32.72

11.68

0.33

0.00

31.44

26.26

21.34

16.69

12.36

8.39

4.86

1.91

59.63

41.96

27.97

17.32

9.63

4.52

1.55

0.25

32.75

41.96

14.60

-7.63

-9.57

-2.82

0.51

0.24

В ельнике, по сравнению со сосново-сфагновым болотом, имеются четкие статистические закономерности, превращаясь в функциональные зависимости с нулевой погрешностью.

3. Березняк травный. По измеренным [4] значениям фитомассы (табл. 13) получена закономерность (табл. 14) для растительного комплекса «трава + береза» в виде уравнения

f

f

f

f(8)

Таблица 13. Содержание фитомассы в березняке травном, ц/га

Наименование

слоя растительного сообщества

Номер слоя i

Вид растения

 

Всего

Береза

Трава

Слой надземной части травы

8

-

7.0

7.0

Листва (хвоя) древесных растений

7

45.0

-

45.0

Мелкие ветви древесных растений

6

-

-

-

Крупные ветви древесных растений

5

142.0

-

142.0

Стволы древесных растений

4

1500.0

-

1500.00

Крупные корни древесных растений

3

400.0

-

400.00

Мелкие корни древесных растений

2

24.0

-

24.00

Слой корней травы

1

-

15.0

15.0

Итого:

2111.0

22.0

2133.0

Таблица 14. Содержание фитомассы в березняке травном, ц/га

Слой

i

Факт

f

Расчетные значения (8)

Составляющие модели (8)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

7.0

45.0

-

142.0

1500.0

400.0

24.0

15.0

7.00

45.00

-

142.02

1500.00

400.01

23.96

15.03

0.004

0.016

-

-0.020

0.003

-0.06

0.042

-0.030

0.06

0.04

-

-0.01

0.00

-0.00

0.18

-0.20

5.07

5.85

6.76

7.80

9.01

10.40

12.01

13.87

0.00

0.00

0.34

84.37

1448.64

357.48

0.04

0.00

18.44

28.03

38.44

45.78

44.44

31.23

12.17

1.12

16.51

11.11

7.03

4.09

2.11

0.90

0.27

0.03

16.51

-11.10

7.02

-4.08

2.10

-0.89

0.27

-0.03

Формула (8) имеет максимальную относительную погрешность 0,20% или доверительную вероятность 100-0,20 = 99,80%. Поэтому анализ составляющих модели, а также их сопоставление позволяют понять многое из жизнедеятельности березняка травного.

Годичный прирост Δq+ фитомассы березняка травного характеризуется данными (табл. 15), по которым было получено статистическое уравнение (табл. 16) в виде

f

f

f

f(9)

Таблица 15. Годичный прирост фитомассы в березняке травном, ц/га

Наименование

слоя растительного сообщества

Номер слоя i

Вид растения

Всего

Береза

Трава

Слой надземной части травы

8

-

7.0

7.0

Листва (хвоя) древесных растений

7

45.0

-

45.0

Мелкие ветви древесных растений

6

-

-

-

Крупные ветви древесных растений

5

5.0

-

5.0

Стволы древесных растений

4

50.0

-

50.0

Крупные корни древесных растений

3

13.0

-

13.0

Мелкие корни древесных растений

2

1.0

-

1.0

Слой корней травы

1

-

5.0

5.0

Итого:

114.0

12.0

126.0

Таблица 16. Годичный прирост фитомассы в березняке травном, ц/га

Слой

i

Факт

f

Расчетные значения (9)

Составляющие модели (9)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

7.0

45.0

-

5.0

50.0

13.0

1.0

5.0

7.00

45.00

-

5.00

50.00

13.00

1.00

5.00

-8.0е-13

-1.4е-12

-

9.9е-14

-1.4е-14

1.8е-15

-2.7е-15

8.9е-16

-0.00

-0.00

-

0.00

-0.00

0.00

-0.00

0.00

155.60

96.98

60.44

37.67

23.48

14.63

9.12

5.68

0.83

3.93

14.51

38.01

60.52

43.84

7.55

0.04

2.53

7.07

16.94

32.98

47.51

42.53

15.97

0.72

852.44

395.96

163.39

57.35

15.92

3.05

0.30

0.01

146.89

48.84

-66.70

37.70

-13.52

2.95

-0.30

0.01

Годичный опад Δq- фитомассы характеризуется данными (табл. 17), по которым было получено статистическое уравнение (табл. 18) в виде формулы

f

f

f

f(10)

Таблица 17. Годичный опад фитомассы в березняке травном, ц/га

Наименование

слоя

растительного сообщества

Номер слоя i

 

Вид растения

 

Всего

Береза

Трава

Слой надземной части травы

8

-

7.0

7.0

Листва (хвоя) древесных растений

7

45.0

-

45.0

Мелкие ветви древесных растений

6

-

-

-

Крупные ветви древесных растений

5

1.9

-

1.9

Стволы древесных растений

4

10.5

-

10.5

Крупные корни древесных растений

3

2.8

-

2.8

Мелкие корни древесных растений

2

0.2

-

0.2

Слой корней травы

1

-

5.0

5.0

Итого:

60.4

12.0

72.4

Таблица 18. Годичный опад фитомассы в ельнике сложном, ц/га

Слой

i

Факт

f

Расчетные значения (10)

Составляющие модели (10)

q

f

f

q1

q2

q3

a

q4

8

7

6

5

4

3

2

1

7.0

45.0

-

1.9

10.5

2.8

0.2

5.0

7.0

45.0

-

1.9

10.5

2.8

0.2

5.0

-1.4е-13

-2.4е-13

-

2.6е-14

1.1е-14

2.4е-14

1.9е-15

1.8е-15

-0.00

-0.00

-

0.00

0.00

0.00

0.00

0.00

150.79

93.98

58.58

36.51

22.75

14.18

8.84

5.51

0.60

2.84

10.49

27.47

43.73

31.69

5.46

0.03

5.70

13.60

27.85

46.34

57.05

43.64

14.00

0.54

264.28

122.76

50.65

17.78

4.94

0.95

0.09

0.00

138.69

38.72

-46.95

15.74

-1.07

-0.57

0.09

-0.00

Если поделить все уравнения на общую фитомассу, то получим относительные фитомассы, а также относительные годичные приросты и опады в фитоценозах. Их сопоставление позволит сравнивать между собой структуры различных видов растительных сообществ.

Комплексные испытания древесины, листвы и корней, коры растущих деревьев, при дополнительном учете листвы и корней травы, позволят изучать экологический режим и состояние лесного и болотного ландшафта. При этом испытатель заранее знает искомую закономерность в виде обобщенной статистической модели (1), и поэтому будет прогнозировать результаты усложняющихся измерений фитомассы растительных сообществ.

СПИСОК ЛИТЕРАТУРЫ:

  • 1. Винокурова, Р.И. Закономерности физико-химических параметров древесины растущего дерева / Р.И. Винокурова, П.М. Мазуркин, Е.В. Тарасенко: Научное издание. - Йошкар-Ола: МарГТУ, 2004. - 187 с.
  • 2. Мазуркин П.М. Геоэкология: Закономерности современного естествознания: Научное изд. / П.М. Мазуркин. - Йошкар-Ола: МарГТУ, 2006. - 336 с.
  • 3. Мазуркин П.М. Рациональное природопользование: учебное пособие. В 3-х ч. Ч. 3: Экологически ответственное лесопользование / П.М. Мазуркин, С.Е. Анисимов, С.И. Михайлова; под ред. П.М. Мазуркина. - Йошкар-Ола: МарГТУ, 2006. - 359 с.
  • 4. Родин Л.Е., Ремезов Н.П., Базилевич Н.Н. Методические указания к изучению динамики и биологического круговорота в фитоценозах. - Л.: Наука, 1968. - 145 с.