On the influence of boundary conditions on the instability of geostrophic currents

全文:

1. Введение

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

С целью такого исследования мы будем анализировать неустойчивость океанского геострофического течения с параболическим вертикальным профилем скорости общего вида [Скороходов и Кузьмина, 2022; Кузьмина и др., 2023]. Это позволит при определенных параметрах модели рассмотреть и частный случай течения с линейным вертикальным профилем скорости, которое является типичным для изучения неустойчивости, возникающей из-за вертикального сдвига (см., например, [Eady, 1949; Miles, 1965; Stern, 1975; Педлоски, 1984; Калашник, 2001; Шакина, 2013]).

Для описания неустойчивости зонального течения конечного поперечного масштаба в ограниченном по вертикали слое будем опираться на модель, основанную на уравнении потенциального вихря в квазигеострофическом приближении с учетом вертикальной диффузии импульса и массы. Модельное уравнение, вывод которого подробно представлен в работах [Kuzmina, 2016; Кузьмина и др., 2018, 2019], сводится к спектральному уравнению четвертого порядка (аналог уравнения Орра–Зоммерфельда). Принципиально важным является то, что четвертый порядок уравнения при классической параметризации диссипации обусловлен учетом в модели вертикальной диффузии массы, а не вертикальной диффузии импульса. Именно это является аргументом для выбора условий на горизонтальных границах слоя: кроме равенства нулю вертикальной скорости возмущений, дополнительно полагается равенство нулю потоков плавучести. Можно предположить, что эти, физически обоснованные граничные условия, занижают влияние трения на неустойчивые возмущения. Поэтому представляет интерес рассмотреть данную задачу с условиями равенства нулю возмущений горизонтальной скорости (условия прилипания) вместо условия равенства нулю потоков плавучести. Так, в работах [Скороходов и Кузьмина, 2018; Кузьмина и др., 2019], в которых анализировалась неустойчивость течения с параболическим профилем скорости (максимум скорости течения был в центре слоя) подробно обсуждались два вида граничных условий. В статьях [Скороходов и Кузьмина, 2018, 2019] представлены и некоторые результаты расчетов неустойчивых возмущений для течения большого поперечного масштаба с учетом условий прилипания, показывающие, что условия прилипания стабилизируют возмущения с малыми волновыми числами. Учитывая, что в данных работах изучались только конкретные виды течений, целесообразно провести более общий анализ влияния условий прилипания на неустойчивость геострофических течений. В связи с этим в настоящей работе будет проведено сравнение влияния на растущие возмущения геострофического течения более общего вида двух типов граничных условий. Будет также представлено обсуждение параметризации трения в задачах описания квазигеострофических движений в океане.

2. Модельные задачи

Для удобства изложения статьи кратко повторим постановку задачи [Кузьмина и др., 2023]. Уравнение потенциального вихря в безразмерном виде для зонального течения с параболическим вертикальным профилем с максимумом в центре слоя было выведено и подробно анализировалось в [Кузьмина и др., 2018, 2019]. Для скорости основного потока с вертикальным параболическим профилем общего вида оно выводится аналогично, и без учета бета-эффекта имеет вид:

$\left(\frac{\partial }{\partial t}+U\left(\text{z}\right)\frac{\partial }{\partial \text{x}}\right)\left(\frac{{\partial }^2\text{p}}{\partial {\text{z}}^2}+\text{Bu}\cdot \Delta \text{p}\right){-}_{\text{zz}}^{}\frac{\partial p}{\partial \text{x}}=$

$=\frac{1}{\text{R}}\left(\frac{{\partial }^4\text{p}}{\partial {\text{z}}^4}+\text{Pr}\cdot \text{Bu}\cdot \frac{{\partial }^2\Delta p}{\partial {\text{z}}^2}\right)$, (1)

где p – возмущение давления, $U\left(z\right)=1-a{z}^2+bz$, $a=\frac{s_1{H}_0^2}{U_0}$, $b=\frac{s_2{H}_0}{U_0}$, s₁, s₂ – положительные константы, характеризующие величины вертикальных линейного и постоянного сдвигов скорости, U₀ – скорость течения в середине слоя, 2H₀ – вертикальный масштаб слоя, $\text{R}=U_0{H}_0^2/\left(KL\right)=\text{Pe}{H}_0/L$, $\text{Pe}=U_0{H}_0/K$ – число Пекле (аналог числа Рейнольдса), K – коэффициент вертикальной диффузии плавучести, $\text{Bu}=H_0^2{N}_0^2/\left(f^2{L}^2\right)$ – число Бургера (см., например, [Cushman–Roisin, 1994]), N₀ – частота плавучести в фоновом состоянии, L – горизонтальный поперечный масштаб зонального геострофического течения, f – параметр Кориолиса, Pr – число Прандтля, $\Delta p=\frac{{\partial }^{2}p}{\partial x^2}+\frac{{\partial }^{2}p}{\partial y^2}$. Ось x направлена на восток (зональная координата); ось y – на север (меридиональная координата), ось z – вверх. Скорость зонального течения направлена вдоль координаты x. Диапазон изменения переменной y в рассматриваемой области: $y\in \left[\mathrm{0,1}\right]$. Диапазон изменения переменной z: $z\in \left[-\mathrm{1,1}\right]$. Для определенности положим, что скорость течения на нижней границе слоя, $z=-1$, равна нулю, то есть $U_0=s_1{H}_0^2+s_2{H}_0$, и, следовательно, a + b = 1.

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

Безразмерные переменные и искомая функция вводились по формулам: x* = x / L, y* = y / L, z* = z / H₀, t* = tU₀ / L, p* = p / p₀, p₀ – размерная константа, $\left[{\text{p}}_{\text{0}}\right]={ì}^2{ñ}^{-2}$ (давление нормировано на отсчетную плотность). Звездочки, которыми отмечены безразмерные переменные и искомая функция, в (1) опущены.

Условия равенства нулю вертикальной скорости на горизонтальных границах слоя и условия равенства нулю возмущений скоростей на боковых границах течения сводятся к следующим уравнениям:

$\left(\frac{\partial }{\partial t}+\left(1-a{z}^2+bz\right)\frac{\partial }{\partial x}\right)\frac{\partial p}{\partial z}-$

$-\left(-2az+b\right)\frac{\partial p}{\partial x}=\frac{1}{R}\frac{{\partial }^{3}p}{\partial z^3}$, $z=\pm 1$, (2)

$\frac{\partial p}{\partial x}=0$, $y=0, y=1$. (3)

Дополнительное условие равенства нулю потоков плавучести на границах слоя имеет вид:

$\frac{{\partial }^{2}p}{\partial z^2}=0$, z = –1, z = 1. (4)

И, наконец, альтернативное к условию (4), граничное условие прилипания имеет вид:

$\frac{\partial p}{\partial x}=\frac{\partial p}{\partial y}=0$, z = –1, z = 1. (5)

Мы будем рассматривать две следующие задачи: уравнение (1) с граничными условиями (2)–(4) (задача I) и уравнение (1) с граничными условиями (2), (3), (5) (задача II).

Следуя стандартному подходу исследования неустойчивости течений (см., например, [Eady, 1949; Lin, 1955]) подставим в (1)–(5) выражение для возмущения давления в виде: $p=F\left(z\right)e^{ik\left(x-ct\right)}\text{sin}\left(\pi ny\right)$. Здесь мы ограничимся случаем n = 1: масштаб возмущения в меридиональном направлении равен масштабу фронта или течения (максимально растущая мода). После подстановки получим соответствующее уравнение и граничные условия в спектральном виде:

 $\left(1-a{z}^2+bz-c\right)\left(\frac{d^{2}F}{d{z}^2}-\text{Bu}\cdot \left(k^2+{\pi }^2\right)F\right)+2aF=$

$=\frac{1}{ik\text{R}}\left(\frac{d^{4}F}{d{z}^4}-\text{Bu}\cdot \text{Pr}\cdot \left(k^2+{\pi }^2\right)\frac{d^{2}F}{d{z}^2}\right)$ , (6)

$\left(1-a+b-c\right)\frac{dF}{dz}-F\cdot \left(-2a+b\right)=\frac{1}{ik\text{R}}\frac{d^{3}F}{d{z}^3}$, z = 1, (7)

$-c\frac{dF}{dz}-F\cdot \left(2a+b\right)=\frac{1}{ik\text{R}}\frac{d^{3}F}{d{z}^3}$, z = –1, (8)

$\frac{d^{2}F}{d{z}^2}=0$, z = –1, z = 1, (9)

F = 0, z = –1, z = 1, (10)

где k – безразмерное волновое число вдоль координаты x, c = c₁ + i · c₂ – комплексная фазовая скорость, i – мнимая единица. Возмущения растут со временем, то есть являются неустойчивыми, если c₂ > 0. Задача I: уравнение (6) и граничные условия (7)–(9). Задача II: уравнение (6) и граничные условия (7), (8), (10).

При a = 1, b = 0 задача I сводится к задаче, рассмотренной в работах [Скороходов и Кузьмина, 2018; Кузьмина и др., 2019], при a = 0, b = 1 она сводится к задаче [Кузьмина и др., 2020], которая подробно рассматривалась для слоя от z = 0 до z = 1.

3. Результаты расчетов модельных задач I и II

Для поиска собственных функций (СФ) и собственных значений (СЗ) модельных задач I и II использовался новый высокоточный аналитико-численный метод решения уравнения 4-го порядка с малым параметром при старшей производной [Skorokhodov and Kuzmina, 2024], более эффективный, чем разработанные ранее методы для решения аналогичных задач (см., например, [Скороходов, Кузьмина, 2018, 2021]). Основное отличие данного метода от предыдущих заключается в том, что на отрезке интегрирования $z\in \left[-\mathrm{1,} 1\right]$ [–1, 1] вводится система трех, а не двух, окрестностей, в каждой из которых решение строится в виде степенных разложений, причем гладкая сшивка разложений позволяет вычислять СФ и СЗ с высокой точностью. Для эффективности решения настоящих задач в широком диапазоне изменения волнового числа при проведении расчетов использовались построенные асимптотические разложения (аналитические решения) при $k\to 0$. Это позволило уменьшать время численного счета и получать собственные значения при произвольно малом изменении волнового числа.

Учитывая, что нас интересует, прежде всего, неустойчивость течения, мы будем представлять в системе координат c₁ и c₂ только такие распределения СЗ при непрерывном изменении волнового числа k (так называемые, траектории СЗ), на которых наблюдаются собственные значения с положительной мнимой частью, c₂ > 0 (подробнее о траекториях СЗ см., например, [Orszag, 1971; Reddy et al., 1993; Trefethen, 1997; Shkalikov, 2004; Скороходов, 2007 а, б]).

На первом этапе для построения траекторий оценивались на основе анализа асимптотических разложений СФ и СЗ при $k\to 0$. Отметим, что спектры СЗ задач, подобных задачам I и II, являются счетными и дискретными для любого k и значений параметров задачи (см., например, [Shkalikov, 2004]). На втором этапе, с целью анализа, в первую очередь, неустойчивости течения, из множества СЗ выбирались СЗ с максимальными величинами c₂. Прежде всего, выбирались СЗ с  $c_2\ge 0$ и те СЗ с c₂ < 0, модули которых имели минимальные значения. Как правило, выбирались десять СЗ. Эти СЗ упорядочивались в соответствии с уменьшением величины c₂, а именно: первое СЗ имело максимальное значение c₂, второе СЗ имело следующее по величине значение c₂, и так далее. Наконец, численно рассчитывались с большой точностью траектории выбранных СЗ при почти непрерывном увеличении волнового числа до значения k = 2500. Обозначение траекторий принималось в соответствие с упорядоченным набором СЗ: первая траектория – это траектория первого СЗ, вторая траектория – это траектория второго СЗ, и так далее. Для любого случая анализа задач I и II рассчитывались не менее десяти траекторий. Методы построения траекторий и спектров собственных значений для задач I и II подробно изложены в работах [Скороходов и Кузьмина, 2018, 2021; Skorokhodov and Kuzmina, 2024].

На рис. 1 представлены две первые траектории СЗ задачи I при Bu = 0.0001 (широкое течение), R = 10, Pr = 1 и различных значениях параметра а. Результаты свидетельствуют о том, что существуют неустойчивые возмущения с большими длинами волн (малыми волновыми числами, k << 1).

 

Рис. 1. Первые (сплошные линии) и вторые (штриховые линии) траектории СЗ для задачи I при Bu = 0.0001, R = 10, Pr = 1 и различных значениях параметра а: (а) a = 1 ; (б) а = 0.5; (в) а = 0.2. Точки на траекториях показывают значения СЗ при k = 0.1 (4), k = 1 (3), k = 100 (2), k = 1000 (1)

 

При a = 1 (рис. 1а), когда максимум скорости среднего течения расположен в центре слоя, фазовая скорость неустойчивых возмущений на двух траекториях может превышать максимальную скорость течения (U_max = 1). Такую неустойчивость необходимо отличать от бароклинной неустойчивости и неустойчивости критического слоя (подробнее см. [Кузьмина и др., 2019]). При а = 0.5 (рис. 1б) неустойчивые возмущения с фазовой скоростью, превышающей максимальную скорость потока (U_max = 1.125), наблюдаются также на двух траекториях, причем на второй траектории такие возмущения наблюдаются на больших волновых числах от k = 9.5 до k = 115. При a = 0.2 безразмерная максимальная скорость потока равна U_max = 1.6, и согласно рис. 1в неустойчивых возмущений с фазовой скоростью, превышающей максимальную скорость потока, нет. Подробное описание указанных возмущений в зависимости от отклонения максимума скорости потока от центра слоя дано в работе [Кузьмина и др., 2023].

Расчеты первых двух траекторий СЗ для задачи II при тех же значениях параметров, что и ранее, показаны на рис. 2. На первой траектории (рис. 2а) нет неустойчивых собственных значений с фазовой скоростью, превышающей максимальную скорость потока, причем неустойчивые возмущения на этой траектории наблюдаются только при k > 2.8. На второй траектории неустойчивые возмущения наблюдаются при k > 7.4, причем в диапазоне $7.5\le k\le 33$ неустойчивые возмущения имеют фазовую скорость больше максимальной скорости потока. При a = 0.5 (рис. 2б) на первой траектории неустойчивые СЗ возникают при k > 1.5. На второй траектории неустойчивых СЗ не наблюдается. При a = 0.2 неустойчивые возмущения есть только на первой траектории при k > 0.9. Таким образом, полученные результаты свидетельствуют о существенном сужении спектра неустойчивых возмущений, прежде всего за счет стабилизации длинноволновых возмущений (малые k).

 

Рис. 2. То же, что и на рис. 1, но для задачи II

 

Особый интерес представляет сравнение результатов расчетов задач I и II для течения с линейным вертикальным профилем скорости (a = 0): рис. 3 и 4.

 

Рис. 3. Первая (сплошная линия) и вторая (штриховая линия) траектории СЗ для задачи I при Bu = 0.0001, R = 10, Pr = 1, а = 0. Стрелками показаны направления изменения СЗ при увеличении волнового числа k

 

Рис. 4. Траектории СЗ для задачи II, на которых есть области СЗ с положительной мнимой частью (Bu = 0.0001, R = 10, Pr = 1, а = 0). Стрелками показаны направления изменения СЗ при увеличении волнового числа k

 

Рассмотрим подробнее на рис. 3 изменение значения c₂ на первой траектории. При волновых числах k, лежащих в диапазоне изменения $0<k\le 4$, c₂ $\approx$ 0.557. При k > 4 значение c₂ меняется медленно с ростом k: так при k = 10 величина c₂ = 0.574, и только при k ≈ 120 уменьшается до значения 0.001. Далее наблюдается слабое увеличение и последующее уменьшение c₂ с ростом k. Максимальное значение c₂ = 0.011 при больших волновых числах (при k > 120) достигается при k = 740.

Описание неустойчивости течения с линейным вертикальным профилем скорости с учетом диффузии массы и импульса было подробно представлено в работе [Кузьмина и др., 2020], но для слоя от z = 0 до z = 1. Согласно этой работе область на первой траектории рис. 3 с постоянной фазовой скоростью (c₁ = 1) следует отнести к неустойчивым бароклинным возмущениям. Области на первой и второй траектории с положительными c₂ и фазовыми скоростями, отличными от единицы, следует отнести к неустойчивости критического слоя.

Расчеты задачи II при а = 0 показали, что неустойчивые СЗ могут наблюдаться на 4-х траекториях. Для удобства первая траектория и дополнительная к ней траектория представлены на двух первых фрагментах рис. 4. На последнем фрагменте рис. 4 представлены еще две траектории СЗ, на которых наблюдаются области СЗ с положительными значениями c₂, свидетельствующие о неустойчивости критического слоя.

Опишем кратко изменение СЗ с ростом волнового числа k на левом фрагменте рис. 4, где представлена первая траектория. При малых k значение c₂ отрицательно, и при k = 0.01 равно c₂ = -24.667. С ростом волнового числа модуль c₂ уменьшается, и при k в диапазоне 0.72 ≤ k ≤ 0.74 c₂(k) переходит через ноль. Далее c₂(k) растет с увеличением k вплоть до значения c₂ = 0.5699 (при k = 12.7), и затем с увеличением k значение c₂ уменьшается до нуля (при k ≈ 179), становясь отрицательным при k > 179. Таким образом, новое граничное условие подавляет неустойчивые возмущения на малых волновых числах. Однако и в данном случае существуют два типа неустойчивости течения: бароклинная неустойчивость и неустойчивость критического слоя. Это естественно следует и из других фрагментов рис. 4.

На рис. 5 представлены главные траектории для узкого фронта (Bu = 0.1) при а = 0.5 и прежних остальных параметрах. В данном случае расчетов результаты более наглядно, чем ранее, показывают, что в задаче II происходит стабилизация неустойчивых возмущений не только при малых, но и при больших волновых числах.

 

Рис. 5. Две первые траектории СЗ при Bu = 0.1, R = 10, Pr = 1, а = 0.5: а) задача I; б) задача II. Точки на траекториях показывают значения СЗ при k = 0.1 (4), k = 1 (3), k = 10 (2), k = 100 (1)

 

Важным параметром неустойчивых возмущений в рамках данной работы является инкремент роста kc₂ максимально растущего возмущения при различных граничных условиях. На рис. 6 даны инкременты роста для случая широкого фронта при различных значениях параметра a, построенные на основе первых траекторий задач I и II.

 

Рис. 6. Инкременты роста неустойчивых возмущений, построенные на основе первых траекторий, для различных граничных условий при Bu = 0.0001, R = 10, Pr = 1 и a = 1 (сплошные линии), a = 0.5 (штриховые линии): a) задача I; б) задача II

 

На рис. 6а мы видим два ярких максимума на кривых зависимости инкремента роста от волнового числа, как для а = 1, так и для а = 0.5. Данные максимумы свидетельствуют о максимально растущих возмущениях вследствие бароклинной неустойчивости (первый максимум) и неустойчивости критического слоя (второй максимум). (Описание данных двух видов неустойчивости для аналогичных задач дано, например, в [Кузьмина и др., 2019].) При a = 1 значения kc₂ в максимумах равны: kc₂ = 36.7 при k = 200 и kc₂ = 24.4 при k = 1450. При a = 0.5 значения kc₂ в максимумах следующие: kc₂ = 24.8 при k = 125 и kc₂ = 17.6 при k = 1325.

Отметим также, что инкременты роста на двух кривых рис. 6а спадают с ростом k после второго максимума, но являются положительными даже при k = 2500. Ограничение по волновому числу k = 2500 определялось с точки зрения целесообразности выполнения большого числа вычислений: хотя метод расчета практически применим при сколь угодно больших значениях волнового числа k, но для его реализации в таких случаях требуется высокая точность для представления десятичных чисел с очень большой мантиссой и, следовательно, значительное реальное время для необходимого объема вычислений.

Во втором случае (рис. 6б) на кривой для a = 1 наблюдается только один, первый максимум при k = 200. Максимальное значение kc₂ в этом случае составляет 38.3. При a = 0.5 в точке первого максимума kc₂ = 25.3 (k = 125), в точке второго максимума значение kc₂ = 25.3 (k = 405).

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

4. Обсуждение результатов

Здесь мы кратко представили результаты анализа влияния граничных условий на неустойчивость геострофического течения с параболическим вертикальным профилем общего вида. Результаты анализа свидетельствуют о том, что граничные условия прилипания (см. результаты расчетов задачи II) стабилизируют длинноволновые неустойчивые возмущения и сужают диапазон неустойчивых коротковолновых возмущений. Однако все типы неустойчивости течения, такие как бароклинная неустойчивость, неустойчивость критического слоя, а также новая неустойчивость, характеризующаяся фазовой скоростью, превышающей максимальную скорость течения (см. задачу I и подробное описание видов неустойчивости в работах [Кузьмина и др., 2019, 2023]), возникают и при граничных условиях прилипания, но в более узком диапазоне изменения параметров задачи. Этот факт является подтверждением того, что в рамках уравнения потенциального вихря с учетом вертикальной диффузии импульса и массы можно описывать при различных граничных условиях применительно к океану вихри, обусловленные неустойчивостью разного типа, и крупномасштабные интрузии океанского термоклина [Zhurbas, 2018; Кузьмина и др., 2018].

Отметим также, что вопрос о проблеме учета трения в задачах описания квазигеострофической динамики не является тривиальным. Действительно, если в уравнении (6) не учитывать диффузию массы, то это уравнение будет являться уравнением второго порядка. В таком случае на горизонтальных границах слоя следует принять, прежде всего, равенство нулю возмущения вертикальной скорости. Для дополнительного граничного условия с целью более эффективного учета влияния трения на возмущения необходимо в таком случае отказаться от классической параметризации трения. Так, в работе [McWilliams, 1989] влияние трения описывалось на основе четвертой, а не второй, производной возмущений скорости по вертикальной координате z. Такой подход усложняет и вид решения исходного уравнения, так как без учета диффузии массы исходное, классическое уравнение потенциального вихря второго порядка приводится к уравнению четвертого порядка. Насколько такой подход является обоснованным для описания реальных ситуаций в океане, требует отдельного исследования.

5. Выводы

Проведен анализ влияния граничных условий на неустойчивость геострофического зонального течения конечного поперечного масштаба с вертикальным параболическим профилем скорости общего вида в ограниченном по вертикали слое. Модель основана на уравнении потенциального вихря в квазигеострофическом приближении с учетом вертикальной диффузии массы и импульса. Уравнение и граничные условия сводились к спектральной задаче на собственные значения типа Орра–Зоммерфельда. Рассматривались два вида условий на горизонтальных границах слоя: 1) равенство нулю возмущений вертикальной скорости и потоков плавучести (задача I); 2) равенство нулю возмущений вертикальной скорости и возмущений горизонтальных скоростей (задача II). Для одних и тех же значений физических параметров уравнения рассчитывались траектории собственных значений и инкременты роста неустойчивых возмущений для граничных условий «1» и «2». На основе сравнения результатов получено, что использование граничных условий «2», которые включают условия прилипания, способствует стабилизации длинноволновых неустойчивых возмущений и сужает диапазон неустойчивых коротковолновых возмущений. Отмечается, однако, что все типы неустойчивости течения, полученные при условиях «1», такие как бароклинная неустойчивость, неустойчивость критического слоя, а также новая неустойчивость, характеризующаяся фазовой скоростью, превышающей максимальную скорость течения, возникают и при использовании граничных условий прилипания, но в более узком диапазоне изменения физических параметров задачи. Представлено также краткое обсуждение вопроса о проблеме учета трения в задачах описания квазигеострофической динамики в океане.

Благодарности

Авторы признательны анонимным рецензентам за полезные замечания.

Финансирование

Работа поддерживалась бюджетным финансированием Института океанологии им. П.П. Ширшова РАН (тема FMWE -2024-0015).

作者简介

N. Kuzmina

Shirshov Institute of Oceanology of the RAS

编辑信件的主要联系方式.
Email: kuzmina@ocean.ru
俄罗斯联邦, Nakhimovsky Prosp., 36, Moscow, 117997

S. Skorokhodov

Federal Research Center “Computer Science and Control” of the RAS

Email: sskorokhodov@gmail.com
俄罗斯联邦, str. Vavilov, 44, Moscow, 119333

N. Zhurbas

Shirshov Institute of Oceanology of the RAS

Email: kuzmina@ocean.ru
俄罗斯联邦, Nakhimovsky Prosp., 36, Moscow, 117997

D. Lyzhkov

Shirshov Institute of Oceanology of the RAS

Email: kuzmina@ocean.ru
俄罗斯联邦, Nakhimovsky Prosp., 36, Moscow, 117997

参考

补充文件