CC BY
29
Liter atur a
1. Skovorodin V. Ja., Purshel E.E. Issledovanie vozmojnosti formirovaniia metallokeramicheskih plenok pri finishnoi antifrikcionnoi obrabotke gilz cilindrov geomodifikatorami // Izvestie Sankt-Peterbyrgskogo gosydarstvennogo agrarnogo yniversiteta.
- 2016/-№ 42/-S. 333-340
2. Sovremenaia tribologiia: Itogi i percpektivi / Otv. red. K.V. Frolov. - M.: Izdatelstvo LKI, 2008.-480c.
3. Skovorodin V. Ja., Pankrashev A.S. Finishnaia obrabotka gilz cilindrov DVS s primeneniem antifrikcionnih materialov// Trudi vserosiiskogo - issledovatelskogo institute remonta i ekspluataciimashinno - traktornogo parka/ Redkol.: Chernoivanov V.I. [i dr.] - 2010. - №105.
- C.80-83.
4. Dunaev A.V., Zuev V.V., Vasilkov D.V. i dr. Gipotezi mehanizmov deistviia remontno -vosstanovitelnih serpenterovannih tribopreparatov// Trudi GOSNITI. - 2013. T.112. - №2. -s.134-142.
5. Kovalevskaia J.G. Hatkova A.V., Beliavskaia O.V. i dr. Vliianie nagreva na fazovie prevrasheniia v geomodifikatore treniia na osnove sloistogo gidrosilikata // Obrabotka metallov. - 2013. - №1 (58). - S 75-80.
УДК 582.79:57.033
Доктор техн. наук С.А. РАКУТЬКО (СПбГАУ, sergejl964@yandex.ru) Аспирант А.Н. ВАСЬКИН (СПбГАУ, vaskin69@mail.ru) Соискатель Е.Н. РАКУТЬКО (ИАЭП, elena.rakutko@mail.ru)
СТАТИСТИЧЕСКИЙ АНАЛИЗ ФЛУКТУИРУЮЩЕЙ АСИММЕТРИИ БИЛАТЕРАЛЬНЫХ ПРИЗНАКОВ ЛИСТЬЕВ ПЕТРУШКИ (PETROSELINUM TUBEROSUM) ПРИ ВЫГОНКЕ ПОД РАЗЛИЧНЫМ СПЕКТРОМ ИЗЛУЧЕНИЯ
Фитомониторинг, светокультура, спектр, флуктуирующая асимметрия
Энергия потока оптического излучения (ОИ) в области фотосинтетически активной радиации (ФАР) с длиной волны от 400 до 700 нм является основой обеспечения процессов фотосинтеза и получения полезной продукции в светокультуре, оказывая большое влияние на рост, развитие и физиологию растений. У растений под влиянием энергии ОИ наблюдается целый ряд физических эффектов, ведущих к регуляторным, адаптивным и другим процессам, вплоть до экспрессии генов. Единичные кванты ОИ, поглощаемые растительным организмом, запускают быстрые и хорошо заметные превращения морфофизиологического состояния растений. В целях обеспечения максимальной продуктивности выращиваемых растений при минимуме энергетических затрат применяют дополнительное облучение от искусственных источников света (ИС). При этом возможно широкое варьирование параметрами излучения: интенсивностью, продолжительностью, спектральным составом.
Стабильность развития растений проявляется во взаимодействии случайных событий в организме растения и их способности точно следовать программе, заложенной в генотипе, сопротивляясь воздействиям окружающей среды во время развития для формирования оптимального фенотипа. Недостаточное качество среды выращивания выражается в явлении нестабильности развития.
Наиболее ярким проявлением стабильности развития биообъекта на макроуровне является флуктуирующая асимметрия (ФА), заключающаяся в незначительных и случайных отклонениях параметров билатеральных (зеркальных) признаков. Уровень ФА является
минимальным лишь при оптимальных условиях среды и возрастает при любых стрессовых воздействиях [1].
Известно, что ФА тесно коррелирует с состоянием качества среды растений, растущих как в естественных [2], так и в контроллируемых [3] условиях. Феномен флуктуирующей асимметрии находит широкое применение в экологическом мониторинге природных популяций при оценке воздействия различных факторов на среду обитания живых организмов. Учет незначительных и ненаправленных отклонений от строгой билатеральной симметрии, которые выражаются тем отчетливее, чем сильнее внешние воздействия (в первую очередь антропогенное загрязнение), позволяет оценить стабильность развития организма [4].
При этом дестабилизация развития организма начинает проявляться уже на относительно низком уровне средовых нарушений, которые еще не связаны с необратимыми изменениями в организме. Это позволяет использовать ФА как индикатор стабильности развития организмов, характеризующий даже незначительные отклонения параметров среды от фонового состояния [5]. Полагая, что неоптимальность параметров ОИ, воздействующего на растения, является фактором стресса, появляется возможность принять уровень ФА за индикатор состояния растения, по которому можно оценивать эффективность и экологичность процесса выращивания растений в искусственных условиях с учетом качества потока ОИ.
Если взаимосвязь уровня ФА растений со стрессом и состоянием окружающей среды для естественных условий широко представлена в научных публикациях, то для условий светокультуры этот вопрос исследован недостаточно подробно.
Свет является важнейшим экологическим фактором, влияющим на рост растений и производство биомассы. Естественным источником излучения служит солнце. Фотоны более коротких длин волн имеют значительную энергию и могли бы оказать повреждающее действие на биологические объекты (растения), однако они задерживаются озоновым слоем стратосферы. Фотоны более длинноволнового излучения не имеют достаточной энергии, чтобы инициировать реакции фотосинтеза.
Недостаточная интенсивность света или его неудовлетворительный спектральный состав ухудшает рост и развитие растений. Как для дополнительного досвечивания растений в теплицах, так и в условиях интенсивной светокультуры необходимый уровень облученности создают с помощью натриевых, металлогалогенных или люминесцентных ламп. Эти источники света являются достаточно эффективными приборами для преобразования электрической энергии в энергию фотонного потока в целом, однако спектр их излучения далек от действительной потребности зеленого растения в энергии излучения на различных длинах волн ФАР. Появление светодиодных ИС позволят более тонко задавать спектр излучения для растений.
В результате ряда фотобиологических исследований для некоторых светокультур найдены спектральные соотношения, обеспечивающие наилучшую продуктивность выращиваемых растений. Частные реакции растений на действие излучения отдельных спектральных диапазонов могут быть предсказаны на основании хорошо известных закономерностей. Так, синий свет снижает вытягивание растения и препятствует увеличению площади листа. Зеленый свет имеет большую проникающую способность, что ведет к интенсификации фотофизиологических процессов в глубине кроны. Желтый свет приводит к сильному вытягиванию растений, удлинению как гипокотиля, так и первых междоузлий. Красный свет способствует удлинению гипокотиля и увеличению площади листьев.
Что касается общей реакции растений на действие излучения с комбинированным спектром, то ее прогнозировать труднее из-за сложного взаимодействия откликов растения на излучения отдельных диапазонов ФАР.
Цель исследования - подтверждение гипотезы о влиянии различий в спектральном составе излучения на уровень ФА листьев на примере петрушки {Petroselinum tuberosum) при выгонке зелени.
Материалы, методы и объекты исследования. Для проведения эксперимента по выгонке в условиях светокультуры заранее летом 2015 г. в открытом грунте выращивали корнеплоды корневого сорта петрушки Урожайная. Корнеплоды убрали в первой декаде октября до заморозков, обрезав ботву на уровне головки и хранили в темноте в подвальном помещении при температуре воздуха не выше + 5°С и относительной влажности воздуха 60-70%. Сохранность корнеплодов была 100%. Масса отдельного образца корнеплода варьировала от 11,46 г до 79,42 г. В наших предварительных исследованиях [6] было выявлено, что вне зависимости от спектра излучения корреляция размера полученных при выгонке листьев зелени и массы корнеплода весьма слабая, поэтому для эксперимента брали все корнеплоды в количестве 50 шт. Корнеплоды были разделены на две партии, предназначенные для выращивания под различным спектром, таким образом, что средние массы корнеплодов обеих партий были практически одинаковы (31,15±3,18 г под СД и 31,14+3,38 г под НА), т.е. в каждой партии были представлены корнеплоды всех размеров. Высадили корнеплоды 25.02.2016 г. в полипропиленовые контейнеры объемом 663 см3. Для заполнения контейнеров использовали кислый (pH 3,6) верховой сфагновый торф торфопредприятия «Пельгорское-М» Ленинградской области со степенью разложения 10% , влажностью 35% и низким содержанием основных элементов питания. Кислотность торфа нейтрализовали мелом до pH 6,2. Производили заправку торфа основными элементами питания и микроэлементами. Содержание подвижных форм элементов питания доводили до уровней, мг/л: NH4+- 20; N03+- 194,5; К+- 189,6; Са2+- 160; Mg2+- 60; Мп2+- 0,5; Си2+- 0,05; Мо6+ - 0,05; В1 - 0,05. Содержание Р5+ - 20 мг/ 100 г сухого торфа, ЕС - 1,0 мСм/см.
При посадке головку и шейку корнеплодов оставляли открытыми. Торф вокруг саженцев уплотнили, полили и поставили контейнеры в темное место с температурой воздуха +14°С для образования боковых корней. В процессе выращивания поддерживали температуру воздуха +18°С, влажность воздуха 60% [7]. Полив производили небольшими порциями воды, избегая увлажнения головки корнеплода. Срезку зелени проводили под самый корнеплод. Продолжительность облучения в сутки (фотопериод) устанавливали 8 ч.
Сравнительный эксперимент проводили в двух зонах лабораторного помещения,
разделенных светонепроницаемой шторой. Растения на рабочих столах под облучателями
располагали на поддонах размером 0,45 х 0,45 м, при этом неравномерность облучения
составляла не более 20%. Было размещено по 25 контейнеров с растениями. В процессе
эксперимента в обоих зонах поддерживали одинаковый уровень фотонной облученности - 80 2 1 „ ^ мкмоль'м" с" путем изменения высоты подвеса облучателей над верхушками растений. В
первой зоне использовали светодиодные облучатели фирмы «Ледел» марки L-fito, в
количестве 5 шт. закрепленных на штангах с шагом 0,3 м, размещенные на высоте 1,17 м над
верхушками растений. Единичная мощность облучателя 60 Вт, тип КСС - Д, состоит из 10
синих, 10 белых и 30 красных светодиодов (СД). Спектр, задаваемый соотношением
интенсивности излучения в синем к,, (400-500 нм), зеленом kG (500-600 нм) и красном kR
(600-700 нм) спектральных диапазонах ФАР, составлял: ки : к,, :kR = 31% : 9% : 60%. Во
второй зоне использовали облучатель ЖСПЗ0-400-010 «Reflux» с натриевой лампой (НА) мощностью 400 Вт, размещенный на высоте 1,07 м над верхушками растений. Спектр излучения - kB:kG:kR= 9% : 54% : 37%.
Выбор источников излучения объясняется максимальным различием их спектров. Там, где у одного источника наблюдается значительный поток в определенном спектральном диапазоне (синем и красном у СД), у другого имеется провал (НА), и наоборот. В соответствии с принятой гипотезой максимальные различия в спектре излучения должны максимальным образом повлиять на стабильность развития облучаемых растений. При одинаковой фотонной облученности по причине различного спектрального состава освещенность под СД была почти в 5 раз меньше, чем под НА (2,06 и 9,7 кЛк соответственно), при этом облученность ФАР была больше на 25% (22,1 и 17,7 Втм"2 соответственно).
В конце эксперимента (18.03.2016 г., на 22 сутки) фиксировали количество листьев в розетке и их длину, сырую массу листьев и содержание в них сухого вещества, билатеральные признаки в листьях петрушки.
Статистическая обработка результатов оценки ФА включала проверку данных на нормальность распределения величины разницы между билатеральными признаками; на присутствие направленной асимметрии и антисимметрии; на зависимость величины асимметрии признака от его размера; на коррелированность признаков; сравнение уровня ФА листьев петрушки, полученных под излучением с различным спектром.
Петрушка имеет сложные листья, состоящие из нескольких четко обособленных листовых пластинок (листочков), каждый из которых своим черешком прикреплен к общему черешку (рахису). Прикорневые и нижние стеблевые листья тройчаторассеченные с перисторассеченными или дважды-перисторассеченными сегментами, их конечные сегменты продолговато-яйцевидные или почти ромбовидные в очертании (рис. 1).
Рис. 1. Розетка листьев петрушки Рис. 2. Измеряемые билатеральные
признаки
У листьев петрушки статистическому анализу были подвергнуты четыре билатеральных признака: 1) длины первых (Л/, М/) и 2) вторых (1,2, черешков первого порядка, отходящих от рахиса, а также длины черешков второго порядка, отходящих от первого черешка первого порядка 3) наружу (£?, Кз) и 4) внутрь (1,4, листа (рис. 2). Параметры были измерены у 190 листьев петрушки, из них у 93 листьев под СД и 97 листьев - под НА.
Результаты исследования. Проверка нормальности распределения данных. Статистический анализ закона распределения значений асимметрии билатеральных признаков листьев петрушки под обоими типами ИС, выполненный с использованием критериев согласия Колмогорова-Смирнова, Шапиро-Уилка и Лиллиефорса (табл. 1), показал, что характер распределения большинства полученных выборок отличается от нормального. Это приводит к необходимости использовать непараметрические методы статистического анализа.
Таблица 1. Статистический анализ нормальности распределения значений асимметрии билатеральных признаков листьев петрушки
Признаки ИС Статистические критерии
Колмогорова-Смирнова Шапиро-Уилка Лиллиефорса
<1 Р Р Р
1 сд 0,09574 > 0.20 0,96624 0,015 <0,05
НА 0,05460 > 0,20 0,99250 0,868 > 0,20
2 сд 0,12925 <0.10 0,97848 0.127 <0,01
НА 0.11179 <0,20 0,98247 0,222 <0,01
3 сд 0,09095 > 0,20 0,96641 0,017 <0,10
НА 0.09155 > 0,20 0,95702 0,003 <0,05
4 сд 0,12224 <0,15 0,89550 0,000 <0,01
НА 0.11494 <0,20 0,97253 0,039 <0,01
Примечание. В этой и последующих таблицах жирным шрифтом выделены уровни значимости, достаточные для отклонения нулевой гипотезы.
Проверка направленности асимметрии. Отсутствие направленности асимметрии выражается в том, что распределение различий билатерального признака является статистически однородным и приблизительно симметрично расположенным отосительно нулевого значения. Проверку направленности ФА проводили с использованием критерия Уилкоксона (табл. 2).
Таблица 2. Статистический анализ направленности асимметрии билатеральных
признаков листьев петрушки
Признаки ИС Статистические критерии
Т ъ Р
1 сд 1810,500 0,258360 0,796
НА 1983,500 0,257516 0,797
2 сд 1363,500 0,324754 0,745
НА 1235,500 1,350606 0,177
3 сд 1606,500 0,968370 0,333
НА 1630,500 0,689036 0,491
4 сд 1461,500 0,760212 0,447
НА 1248,000 1,287138 0,198
Ни в одном из случаев не установлено статистически значимых (/КО,05) различий в величине признака на левой и правой стороне листа. Это означает, что имеют место флуктуации асимметрии признака вокруг нулевого среднего, что является диагностическим признаком ФА.
Тестирование на антисимметрию. Индикатором антисимметрии является отрицательный эксцесс распределения различий между сторонами (Ь-Я). При к>0 принимается гипотеза о флуктуирующем характере асимметрии.
Результаты анализа (табл. 3) показывают превышение критических значений эксцесса третьего признака под НА и четвертого под СД, однако в этих случаях эксцесс положителен, что позволяет считать, что антисимметрия у анализируемых признаков отсутствует.
Проверка размер-зависимости. Для выявления связи между величиной асимметрии признака |Ь-Я| и его средним размером (Ь+Я)/2 использовали непараметрический коэффициент ранговой корреляции Спирмена
Таблица 3 . Статистический анализ эмпирических (к) и критических (кК1,) значений эксцесса билатеральных признаков листьев петрушки
Признаки ИС Статистические к эитерии
к кКО (/>=0.05)
1 СД 0,5618 0,902
НА 0,3206 0,883
2 СД -0,2712 -0,649
НА -0,0797 -0,639
3 СД 0,6267 0,902
НА 2,4031 0,883
4 СД 7,5767 0,902
НА 0,2597 0,883
Таблица 4 . Корреляционный анализ между абсолютными значениями асимметрии и средним размером признака листьев петрушки
Признаки ИС Статистические к эитерии
Р
1 СД 0,0413 0,694
НА 0,1601 0.117
2 СД -0,0762 0,468
НА 0,0005 0,996
3 СД 0,0491 0,640
НА -0,0071 0,945
4 СД 0,0582 0,579
НА 0,2355 0,020
Анализ свидетельствует, что статистически значимая (р<0,05) корреляционная связь наблюдается только для четвертого признака под НА (табл. 4). Наличие размер-зависимого признака требует введение нормировки вида |Ь-К|/(Ь+К), которая удобна и для не проявляющих размер-зависимость признаков, т.к. позволяет получить значение асимметрии в интервале (0... 1).
Проверка коррелированности признаков и их асимметрии. Данный шаг анализа позволяет обосновать включение изучаемого признака в систему признаков, интегрально оценивающих ФА растения: чем меньше коррелируют признаки, тем более четкую и обоснованную оценку можно получить с помощью их набора. Для целей анализа использовали ранговый коэффициент Спирмена. Как и следовало ожидать, натуральные значения признаков на левой и правой стороне листа тесно скоррелированы между собой при статистической значимости ¿><0,05 (0,5927 0,8692 для СД и 0,7452 0,9097 для НА). Однако и другие натуральные значения признаков также существенно коррелированы (0,4778</^<0,8331 для СД и 0,6465</^<0,8871 для НА). Существенна корреляция и между асимметриями признаков (0,3508</^<0,7222 для СД и 0,2510 </^<0,6237 для НА). Это позволяет ограничиться для характеристики ФА листьев петрушки только одним признаком. Для удобства измерений целесообразно выбрать первый признак как имеющий наибольший размер.
Биометрические показатели. Большая продуктивность выгонки наблюдалась под НА, были получены более крупные листья с большей сырой массой. В табл. 5 приведены средние значения измеряемых биометрических показателей петрушки с указанием ошибки средней и среднеквадратичного отклонения каждого показателя.
Таблица 5 . Результаты t-теста различий средних значений показателей
Показатель Среднее Ст.откл. Р
СД НА СД НА
Кол-во листьев, шт. 7.8 7,4 1,8 2,8 0,629
Масса листьев, г 5,03 6,37 1,5 2,6 <0,05
Содержание сухого вещества, о/ /О 11,7 11,9 1,4 2,2 0,279
Длина листа, мм 145,3 182,8 28,6 43,9 <0,05
Уровень ФА, отн.ед. 0,052 0,039 0,043 0,031 <0,05
Таким образом, экспериментально выявлено влияние различного спектрального состава излучения на показатели ФА и продуктивность петрушки, выращиваемой в условиях светокультуры на зелень. Меньшие показатели ФА соответствуют большим показателям продуктивности растения [8].
Для оценки соответствия численного показателя ФА и качества среды произрастания растений существует пятибалльная шкала [9]. В соответствии с этой шкалой выделяют качества среды: условно нормальное (I балл, ФА <0,040), незначительные отклонения от нормы (II балл, ФА=0,040 - 0,044), средний уровень отклонений (III балл, 0,045 - 0,049), существенные отклонения от нормы (IV балл, ФА=0,050-0,054), критическое состояние (V балл, ФА>0,054). Сопоставление полученных данных с этой шкалой говорит о том, что если излучение НА ламп соответствует нормальному качеству световой среды при выращивании растений (I балл), то при использовании СД наблюдаются существенные отклонения от нормы (IV балл), что проявляется в снижении продуктивности светокультуры.
Выводы:
1. Выявлена существенная асимметрия листьев петрушки, выращиваемых на выгонку под излучением с различным спектральным составом. Частота встречаемости асимметрии длины наиболее крупных листочков и сегментов сложного листа петрушки составляет 53-85%.
2. Выявлена ненаправленность асимметрии билатеральных признаков и отсутствие у них антисимметрии, что позволило классифицировать наблюдаемую асимметрию как флуктуирую щую.
3. Выявлена существенная корреляция между отдельными билатеральными признаками, что дало основание выбрать для характеристики ФА листьев петрушки только один признак - длины первых черешков, отходящих от рахиса. Размер-зависимость этого признака не обнаружена.
4. Статистически достоверно у растений, выращиваемых под СД, по сравнению с использованием НА, наблюдалась меньшая масса и длина листьев. В то же время уровень ФА у листьев петрушки под ними был больший. Это подтверждает первоначальную гипотезу о том, что большие значения ФА наблюдаются в условиях, менее благоприятных для растений.
5. Показатели ФА листьев могут быть использованы для оценки качества световой среды выращивания растений петрушки. Уровень ФА может выступать диагностическим параметром приемлемости спектра источников излучения в светокультуре для выращиваемой культуры.
Литература
1. Palmer A.R., Strobeck С. Fluctuating asymmetry analysis revisited // Developmental instability (DI): causes and consequences. M. Polak, ed. Oxford University Press, New York, 2003.
2. Светлакова Т.Н., Мандрица C.A., Боронникова C.B., Суслонов А.В. Оценка изменчивости морфологических признаков Trifolium pratense L. В условиях нефтяного загрязнения почв // Университет им. В.И. Вернадского. -2019. -№1-3 (28). - С. 16-22.
3
4
5
6
7
8
9
1
2
3
4
5
6
7
8
9
Кузнецова Е.А., Челпанова О.М., Белова Е.Е., Хотулева О.В., Колонцов А.А. Оценка влияния ионов кадмия на флуктуирующую асимметрию листьев огурца посевного (Cucumis Sativus L.) // Вестник МГОУ. - 2013. - №2. - С.1-9.
Зорина А.А., Коросов А.В. Характеристика флуктуирующей асимметрии листа двух видов берез в Карелии // Экология. Экспериментальная генетика и физиология: Труды Карельс. науч. центра РАН. -Петрозаводск, 2007. -Вып. 11. -С. 28-36. Четанов Н.А., Епланова Г.В. Статистический анализ флуктуирующей асимметрии билатеральных признаков живородящей ящерицы Zootoca Vivipara // Известия Самарского научного центра РАН. -2011. -Т.13. - №1. - С. 144-152.
Ракутько Е.Н., Ракутько С.А. Рост и фотоморфогенез петрушки корневой (Petroselinum tuberosum) под оптическим излучением различного спектрального состава // Известия Санкт-Петербургского государственного аграрного университета.-2015.-№38.-С.298-304.
Потехин, Г.А., Харченко В.А., Пивоваров В.Ф. Особенности выращивания петрушки на зелень // Овощи России. - 2010. - №3 (9). - С.42-47.
Ракутько С.А., Ракутько Е.Н. Взаимосвязь флуктуирующей асимметрии листьев петрушки (Petroselinum Tuberosum) и ее продуктивности под воздействием оптического излучения различного спектрального состава // Физика - наукам о жизни. - СПб.: ФТИ им. А.Ф. Иоффе РАН, 2016. - С. 64.
Методические рекомендации по выполнению оценки качества среды по состоянию живых существ (оценка стабильности развития живых организмов по уровню асимметрии морфологических структур)/ Утверждено Распоряжением Росэкологии от 16.10.2003 № 460-р.-М., 2003.
Liter atur а
Palmer A.R., Strobeck С. Fluctuating asymmetry analysis revisited // Developmental instability (DI): causes and consequences. M. Polak, ed. Oxford University Press, New York, 2003. Svetlakova T.N., Mandrica S.A., Boronnikova S.V., Suslonov A.V. Ocenka izmenchivosti morfologicheskih priznakov Trifolium pratense L. V uslovijah neftjanogo zagrjaznenija pochv // Universitet im. V.I. Vernadskogo. -2019. -№1-3 (28). - S. 16-22.
Kuznecova E.A., Chelpanova O.M., Belova E.E., Hotuleva O.V., Koloncov A.A. Ocenka vlijanija ionov kadmija na fluktuirujushhuju asimmetriju list'ev ogurca posevnogo {Cucumis Sativus L.) II Vestnik MGOU. - 2013. - №2. - S.l-9.
Zorina A.A., Korosov A.V. Harakteristika fluktuirujushhej asimmetrii lista dvuh vidov berez v Karelii // Jekologija. Jeksperimental'naja genetika i fiziologija: //Tr. Karel's. nauch. centra RAN. Petrozavodsk, 2007. - Vyp. 11. - S. 28-36.
Chetanov N.A., Eplanova G.V. Statisticheskij analiz fluktuirujushhej asimmetrii bilateral'nyh priznakov zhivorodjashhej jashhericy Zootoca Vivipara // Izvestija Samarskogo nauchnogo centra RAN. -2011. -T.13. - №1. - S"! 144-152.
Rakutko E.N., Rakutko S.A. Rost i fotomorfogenez petrushki kornevoj (Petroselinum tuberosum) pod opticheskim izlucheniem razlichnogo spektral'nogo sostava // Izvestija Sankt-Peterburgskogo gosudarstvennogo agrarnogo universiteta.-2015.-№38.-S.298-304. Potehin, G.A., Harchenko V.A., Pivovarov V.F. Osobennosti vyrashhivanija petrushki na zelen' // Ovoshhi Rossii. - 2010. - №3 (9). - S.42-47.
Rakutko S.A., Rakutko E.N. The relationship of fluctuating asymmetry of parsley leaves (Petroselinum Tuberosum) and its productivity under the influence of optical radiation of different light quality // Physics - To Life Sciences. SPb.: Ioffe Institute. RAS, 2016. - P. 64. Metodicheskie rekomendacii po vypolneniju ocenki kachestva sredy po sostojaniju zhivyh sushhestv (ocenka stabil'nosti razvitija zhivyh organizmov po urovnju asimmetrii morfologicheskih struktur)/ Utverzhdeno Raspoijazheniem Rosjekologii ot 16.10.2003 № 460-r. M., 2003.
