Влияние орографических волн на общую циркуляцию и перенос озона в атмосфере тема диссертации и автореферата по ВАК РФ 25.00.29, кандидат физико-математических наук Коваль, Андрей Владиславович

• Коваль, Андрей Владиславович
• кандидат физико-математических наук
• 2011

Оглавление диссертации кандидат физико-математических наук Коваль, Андрей Владиславович

Введение

Глава 1 Текущее состояние исследований орографических волн

1.1 Теоретические исследования эффектов орографических 11 волн в средней и верхней атмосфере

1.2 Моделирование воздействия орографических волн

1.3 Наблюдения орографических гравитационных волн в 16 атмосфере

1.4 Влияние орографических волн на глобальный перенос озона

Глава 2 Динамическое и тепловое воздействие стационарных волн

2.1 Введение

2.2 Энергетические характеристики орографических волн во 20 вращающейся атмосфере

2.3 Параметризация воздействия орографических волн

2.4 Численная проверка параметризации.

2.5 Выводы

Глава 3. Влияние орографических волн на циркуляцию атмосферы

3.1 Введение

3.2 Модель средней и верхней атмосферы. 34 3.2.1 Численная реализация параметризации орографических волн

3.3 Среднезональные высотно-широтные распределения 37 скорости ветра и характеристик орографических волн

3.3.1 Распределения в январе месяце

3.3.2 Среднезональные распределения в июле 43 \ 3.3.3 Расчеты для апреля и октября

3.4 Широтно-долготные зависимости в разные сезоны

3.4.1 Глобальная структура зонального ветра и амплитуды орографических волн в январе - феврале

3.4.2 Структура ветра и волновых характеристик в июле - 63 августе

3.4.3 Широтно-долготные распределения в апреле 75 3.4.3 Структура ветра и волновых характеристик в октябре

3.5 Динамическое и тепловое воздействие орографических волн на среднюю атмосферу

3.5.1 Широтно-долготные распределения притоков тепла и 86 волновых ускорений в январе

3.5.2 Распределения в июле

3.5.3 Широтно-долготные распределения в апреле и октябре. 89 3.6 Выводы

Глава 4 Влияние орографических волн на глобальный перенос 93 озона.

4.1 Введение

4.2 Современные представления о глобальном переносе озона.

4.3 Влияние орографических волн на перенос озона общей 94 циркуляцией атмосферы

4.3.1 Потоки озона в январе

4.3.2 Потоки озона в июле 98 4.3.2 Перенос озона в апреле и октябре

4.3 Выводы

Введение диссертации (часть автореферата) на тему «Влияние орографических волн на общую циркуляцию и перенос озона в атмосфере»

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

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

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

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

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

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

В данной работе исследуется влияние мезосферных волн, генерируемых в тропосфере при обтекании неоднородного рельефа земной поверхности атмосферными потоками, на общую циркуляцию и тепловой режим вышележащих слоев атмосферы. Разработаны параметризации генерации орографических волн и их теплового и динамического воздействия на среднюю и верхнюю атмосферу. Парметризации включены в трехмерную нелинейную модель общей циркуляции средней и верхней атмосферы (МЦСВА). Эта модель позволяет рассчитывать циркуляцию атмосферы на всех высотах от тропосферы до термосферы. Модель учитывает все основные динамические и радиационные процессы, включая возможность генерации распространяющихся планетарных волн вблизи земной поверхности и расчеты с учетом переменного содержания углекислого газа в атмосфере. В качестве начального состояния используется безветренная атмосфера с заданным независящим от широты профилем температуры, который берется из модели США-86. Вертикальный поток массы в модели рассчитывается на каждом временном шаге сверху вниз из уравнения неразрывности с помощью конечных разностей на уровнях, расположенных между уровнями прогностических уравнений для горизонтальных составляющих скорости ветра и температуры. Естественным нижним граничным условием для приливных колебаний является требование обращения в нуль вертикальной скорости на поверхности Земли, записываемое через возмущение геопотенциала, обусловленное приливными колебаниями на нижней границе.

Актуальность темы диссертационного исследования

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

Цель диссертационной работы

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

Для достижения этой цели были решены следующие задачи:

1. Получение аналитических выражений для расчета динамических и энергетических характеристик орографических волн с учетом вращения атмосферы.

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

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

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

Научная новизна

Главной особенностью разработанного метода исследования орографических волн является использование новых математических выражений для расчета характеристик орографических волн с учетом вращения атмосферы. Уточнены поляризационные соотношения для случая стационарных волн. При использовании дисперсионного уравнения и поляризационных соотношений для стационарных волн во вращающейся атмосфере были получены уравнения волнового действия в неоднородной атмосфере, аналогичные полученным ранее для нестационарных волн. Эти уравнения решаются для получения параметров волн для различных атмосферных приложений. Указанные формулы получены в данной работе впервые и ранее при оценках эффектов орографических волн не использовались. \ ~ \

Практическая и научная значимость

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

Положения, выносимые на защиту

- полученные аналитические выражения для расчета динамических и энергетических характеристик орографических волн с учетом вращения атмосферы;

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

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

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

Личный вклад автора

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

Апробация и публикация работы

Работа выполнялась на физическом факультете СПбГУ. Тема диссертации включена в план работ кафедры физики атмосферы. Результаты работы докладывались на 9-ти российских\ и международных научных конференциях: на молодежной конференции «Физика и прогресс» (Санкт-Петербург, 2007 г.), международных симпозиумах стран СНГ «Атмосферная радиация и динамика» МСАРД-2009 и МСАРД-2011 (Санкт-Петербург, 2009 и 2011 гг.), международной конференции «25th International Laser Radar Conference» IRLC 25 (С-Петербург, 2010 г.), 12-м симпозиуме Международного Научного камитета по Солнечно-Земной Физике «SCOSTEP Symposium STP12» (Берлин, Германия, 2010 г.), школе молодых ученых РГГМУ (С-Петербург, 2010г.), 5-ой международной конференции "Атмосферная физика, климат и окружающая среда" (Санкт-Петербург, 2010 г.), 4-м семинаре по проблеме вертикального взаимодействия в системе атмосфера-ионосфера «4th IAGA/ICMA/CAWSES-II» (Прага, Чехия, 2011 г.), XI-ой международной научно-практической конференции

Фундаментальные и прикладные исследования, разработка и применение высоких технологий в промышленности" (Санкт-Петербург, 2011 г.). Также результаты исследований докладывались на семинарах кафедры физики атмосферы (2007 - 2011 гг.).

Работа выполнена при частичной поддержке Министерства образования и науки российской Федерации в рамках реализации мероприятия 1.1 Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (государственный контракт №02.740.11.0728 от 05 апреля 2010 г.) по направлению «Геофизика», мероприятия 1.3 федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (государственный контракт №П107 от 12 апреля 2010 г.) по направлению «Мониторинг и прогнозирование состояния атмосферы и гидросферы», а также Российского фонда фундаментальных исследований (Гранты РФФИ 07-05-00913-а (2007-2009 гг.) и 10-05-00719-а.(2012 - 2014 гг.)).

Основные результаты работы представлены в 9 публикациях, включая

1 в рецензируемом журнале. \

Структура и объем диссертации

Диссертация состоит из введения, четырех глав, заключения и списка литературы из 68 наимёнований. Работа содержит 110 страниц, 47 рисунков.

Заключение диссертации по теме «Физика атмосферы и гидросферы», Коваль, Андрей Владиславович

4.4 Выводы.

Для всех рассмотренных сезонов характерно изменение скорости потоков озона на 10 - 15%. Для каждого сезона есть свои особенности. В январе и июле наиболее сильные изменения потоков происходят в северном и южном (зимних) полушариях соответственно, что является следствием более активного распространения орографических волн именно в данных условиях (см. главу 3.3.1 и 3.3.2). В апреле выделяются две зоны генерации орографических волн: на высоких широтах южного полушария и на средних широтах северного. Именно на этих широтах наблюдаются отчетливые изменения вертикальных потоков в стратосфере. Для октября месяца характерно более сильное распространение орографических волн в северном полушарии, представленное на широтах наибольших горных систем, оно приводит к сильным изменениям потоков озона именно в этих областях.

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

Заключение

При выполнении диссертационной работы были получены следующие основные результаты:

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

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

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

4. Выполнены измерения изменений потоков озона в тропо - стратосфере при учете орографических волн.

По результатам научной работы можно сделать следующие выводы:

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

2. Орографические волны создают значительные притоки тепла (до 10-15 К/сут. на высотах 50 км), и ускорений среднего потока (до 20 м с" /сут), что может существенно влиять на общую циркуляцию, планетарные волны и тепловой рея^им средней атмосферы. \

3. В периоды солнцестояний наибольшее влияние орографических волн на общую / циркуляцию средней атмосферы происходит в зимних полушариях. В периоды равноденствия условия распространения орографических волн более равномерные в обоих полушариях. 4. Учет орографических волн приводит к изменениям вертикальных потоков озона, создаваемых общей циркуляцией атмосферы, на 10 - 15% на высотах от 20 до 60 км.

Список литературы диссертационного исследования кандидат физико-математических наук Коваль, Андрей Владиславович, 2011 год

1. Гаврилов Н.М. Основы численных моделей атмосферной динамики. // СПб., 2007, 96 стр.

2. Гаврилов H. М. Структура мезомасштабной изменчивости тропо-стратосферы по измерениям рефракции радиоволн на спутнике CHAMP. // Изв. РАН, Физика атмосферы и океана. 2007. Т. 43. No. 4. С. 492 -501.

3. Гаврилов H. М. Параметризация динамического и теплового воздействия установившихся внутренних гравитационных волн на среднюю атмосферу // Изв. РАН. Физика атмосферы и океана. 1989. Т. 25. №3.271-278.

4. Гаврилов H. М. Параметризация динамического и теплового воздействия установившихся внутренних гравитационных волн на среднюю атмосферу. // Изв. АН СССР, физика атмосферы и океана, т. 25, №3, с. 271 -278, 1989.

5. Гаврилов H. М. Распространение внутренних гравитационных волн в верхней атмосфере и их влияние на турбулентность и тепловой режим. // Автореферат канд. дисс. Л: ЛГУ, 1974, 10 с.

6. Груздев А.Н. Влияние озонного нагревания на динамику планетарных волн // Известия АН СССР. Физика атмосферы и океана. Т.21. №11. С. 1148-1159. 1985.

7. Госсард Э. Э., Хук У. X. Волны в атмосфере // М.: «Мир». 1978. 532 С.

8. Марчук Г.И. Численные методы в прогнозе погоды. Л.: Гидрометеоиздат, 356 с. 1967.

9. Погорельцев А.И. Генерация нормальных атмосферных мод стратосферными васциляциями // Изв. РАН. Физика атмосферы иокеана. Т. 43. №4. С.463-475. 2007.\

10. Погорельцев А.И., Суворова Е.В., Федулина И.Н., Ханна Э., Трехмернаяклиатическая модель распределении озона в средней атмосфере ///

11. Учейые записки. Научно-теоретический журнал. Вып. 10. С-Пб.: изд. РГГМУ. С. 463-475. 2007.

12. Belj aars A.C.M., A. Brown and N. Wood, A new parameterization of turbulent orographic form drag. // Q. J. R. Meteorol. Soc. 130, pp. 1327-1347. 2004.

13. Catry, B., Geleyn, J.-F., Bouyssel, F., Cedilnik, J., Broo, R., Derkova, M., Mladek, R. A new sub-grid scale lift formulation in a mountain drag parameterisation scheme. // Meteorologische Zeitschrift, V. 17, I. 2, Pp. 193208. 2008.

14. Chapman, S. and R.S.Lindzen, Atmospheric tides. // D.Reidel, Norwell, Mass. 201 pp. 1970.

15. Eckermann S. D., I.Hirota and W. K. Hocking. Gravity wave and equatorial wave morphology of the stratosphere derived from long-term rocket soundings. // Q. J. R. Meteorol. Soc. (1994), 121, pp. 149-186., 1994.

16. Eckermann, S. D., and P. Preusse (1999), Global measurements of stratospheric mountain waves from space. // Science. 1999. V. 286. 1534— 1537.

17. Egger, J. Time varying flow over mountains: Temperature perturbations at the surface. // Meteorologische Zeit., V. 18,1. 1, Pp. 101-106., 2009.

18. Fleming E.L., Sushil Chandra, M.R. Shoeberl, and J.J. Barnett, Monthly Mean Global Climatology of Temperature, Wind, Geopotential height, and Pressure for 0-120 km // NASA eypTinn^fl Memorandum 100697. pp. 87. 1988.

19. Forbes, J.M., Atmospheric tides, 1, Model description and results for the solar diurnal component // J. Geophys. Res. V. 87. P. 5222-5240. 1982.

20. Forget F., F. Hourdin, R. Fournier, Ch Hourdin and O. Talagrand, Improved General Circulation Models of the Martian atmosphere from the surface to above 80 km. // J. Geophys. Res. V. 4., No E10.,P. 24155-24175. 1999.

21. Fortuin J.P.F., Langematz U. An update on the global ozone clymatology and on concurrent ozone and temperature trends // Atmospheric Sensing and Modelling, Proc. SPIE 2311. P. 207-216. 1995.

22. Froehlich K., A. Pogoreltsev, and Ch. Jacobi. Numerical simulation of tides, Rossby and Kelvin waves with the COMMA-LIM model // Adv. Space Res. V. 32. No. 5. P. 863-868. 2003.

23. Fritts D. C., and Lu W. Spectral Estimates of Gravity Wave Energy and Momentum Fluxes. Part II: Parameterization of Wave Forcing and Variability // Journal of the Atmospheric Sciences. 1993. SO. 3695-3713.

24. Gavrilov, N. M. On the generation of internal gravity waves in the atmosphere by mesoscale turbulence // M.: Hydrometeoizdat Press. 1988. 7477.

25. Gabriel, A.D., D. Peters, I. Kichner, and H.-F. Graf. Effect of sonally asymmetric ozone on stratospheric temperature and planetary wave propagation // Geophys. Ras. Lett. V. 34. L06807. doi: 10.1029/2006GL028998, 2007.

26. Grisogono B. and L. Enger, Boundary-layer variations due to orographic-wave breaking in the presence of rotation. // Q. J. R. Meteorol. Soc., 130, pp. 2991-3014. 2004.

27. Hedin A. E.J. Geophys Res, V. 96, N. A2. P. 1159-1172. 1991. doi: 10.1029/90JA02125

28. Hines, C. O. A Modeling of Atmospheric Gravity Waves and Wave Drag Generated by Isotropic and Anisotropic Terrain // Journal of the Atmospheric Sciences. 1988. 45 309-322.

29. Jiang Q. and R. B. Smith, Gravity Wave Breaking in Two-Layer Hydrostatic Flow. // Journal if the Atm. Sci., V.60. pp. 1159-1172. 2003.

30. Jiang, J. H., D. L. Wu,\and S. D. Eckermann (2002), Upper Atmosphere Research Satellite (UARS) observation of mountain waves over the Andes. // J. Geophys. Res. 2002. y. 107. No. D20. 8273. doi:10.1029/2002JD002091.

31. Sensitivity, analysis method, and a case study. // J. Geophys. Res. 2002. V. 107. No. D23. 8178. doi:10.1029/2001JD000699.

32. Randel, W.J., and F. Wu. A stratospheric ozone profile data set for 19792005: Variability, trends, and comparisons with column ozone data. // J.Geophys. Res. V. 112. D06313, doi:10.1079/2006JD007339.

33. Sarin, V.B. , Forichon, M.Le Treut, H. Parameterization and influence of the orographic gravity-wave drag in the LMD-GCM // Mathematical and Computer Modelling, V. 24,1. 4., Pp. 71-84. 1996.

34. Strobel, D.F. Parametrization of the atmospheric heating rae from 15 to 120 rm ue to 02 and 03 absorption of solar radiation // J.Geophys. Res., V. 83. P. 6225-6230. 1978.

35. Scinocca J. F., and McFarlane N. A. The parametrization of drag induced by stratified flow over anisotrophic orography // Quarterly Journal of the Royal Meteorological Society. 2000. V. 126. No. 568. 2353-2393.

36. Scinocca J. F., and B.R. Sutherland, Self-Acceleration in the Parameterization of Orographic Gravity Wave Drag. // Journal if the Atm. Sci., V.67. No. 11. pp. 2537-2546. 2010.

37. Scorer, R. S. (1949), Theory of waves in the lee of mountains, // Quarterly Journal of the Royal Meteorological Society. 1949. V. 75. No. . 41- 56.

38. Shefov N. N., Pertsev N. N. Orographic disturbances of upper atmosphere emissions // Handbook for MAP. Urbana: SCOSTEP. 1984. V. 10. P. 171175.

39. Shimada, T., Kawamura, H., Wind-wave development under alternating wind jets and wakes induced by orographic effects. // Geophysical Research Letters, V. 33.1. 2. 2006.

40. Smith S., Baumgardner J., Mendillo M. Evidence of mesospheric gravity-waves generated by Orographic forcing in the troposphere // GeopHys. Res. Lett. 2009. V. 36. L08807. doi:10.1029/2008GL036936.

41. Steeneveld G.-J., C/ J. NAPPO, and A.A.M. Holstag, Estimation of Orographically Induced Wave Drag in the Stable Boundary Layer during the

42. CASES-99 Experimental Campaign. // Acta Geophysica, v. 57, no. 4, pp. 857-881 DOI: 10.2478/sl 1600-009-0028-3. 2009.

43. Swinbank R. and D.A. Orland, Compilation of wind data for the Upper Atmosphere Research (UARS) Reference Atmosphere Project // J.Geophys. Res. V. 108, NO. D19, 4615, doi:10.1029/2002JD003135, 2003.

44. Teixeira M. A. C. and P. M.A. Miranda, A linear model of gravity wave drag for hydrostatic sheared flow over elliptical mountains, // Q. J. R. Meteorol. Soc., 132, pp. 2439-2458. 2006.

45. Teixeira M. A. C., P. M.A. Miranda ,J. L.Argain and M. A. Valente, Resonant gravity-wave drag enhancement in linear stratified flow over mountains. // Q. J. R. Meteorol. Soc., 131, pp. 1795-1814. 2005.

46. Teixeira M. A. C. and P. M.A. Miranda, and R. M. Cardoso, Asymptotic gravity wave drag expressions for non-hydrostatic rotating flow over a ridge. // Q. J. R

47. Watanabe, S., Kawatani, Y., Tomikawa, Y., Miyazaki, K., Takahashi, M., Sato, K., General aspects of a T213L256 middle atmosphere general circulation model. // Journal of Geophysical Research D: Atmospheres, V. 113., I. 12, 2008.

48. Wood, N, Brown, A.R, Hewer, F.E, Parametrizing the effects of orography on the boundary layer: An alternative to effective roughness lengths. // Q. J. R, Meteorol.Soc. 127 (573), pp. 759-777. 2001.

49. Worthington R.M. Radar measurement of the effect of boundary-layer saturation on mountain-wave amplitude. // Meteorol Atmos Phys. V. 105,. pp. 29-35 DOI 10.1007/s00703-009-0032-9. 2009.