Современные прогнозы погоды и климата основаны на численном решении уравнений гидродинамики атмосферы. В основе такого прогноза лежат гидродинамические модели атмосферы.
В настоящее время существуют, по крайней мере, три масштаба моделей, описывающих формирование погоды и климата с различным пространственным масштабом:
1) глобальные, которые используют в самых больших центрах погоды для прогноза крупномасштабных процессов и формирования динамических краевых условий к региональным прогнозам;
2) региональные, которые имею шаг около 3-10 километров, и служат для детального описания развития синоптических процессов над регионами, которые обслуживаются большинством центров погоды;
3) мелкомасштабные, с шагом от 100 до 500 метров, предназначенные для детального изучения процессов, формирующих погоду приземного и пограничного слоя.
Все модели описывают одну и ту же физику, но с разной детализацией и параметризациями явно не представляемых в модели, но важных физических процессов. В настоящее время наиболее полезны для профессиональной деятельности гидрометеорологов региональные модели. Они позволяют рассматривать все важнейшие опасные явления, такие, например, как наводнения всех видов, опасные конвективные явления, опасные температурные волны. Кроме того, региональные модели используются для статистического моделирования формирования и изменений регионального климата. Сложность этих моделей такова, что их разработкой заняты только большие коллективы высококвалифицированных работников. Но те модели, которые имеются в свободном доступе, дают возможность после изучения применять их для самых разнообразных региональных исследований, в частности, для восполнения недостающей климатологической информации.
Все гидродинамические модели атмосферы не лишены ошибок. Типов ошибок гидродинамического моделирования несколько. Это ошибки начальных данных, ошибки, связанные с постановкой граничных условий, ошибки в формулировке (упрощении) уравнений, ошибки, обусловленные выбором пространственного и временного разрешения модели, ошибки, связанные с выбором схем параметризации физических процессов подсеточного масштаба.
Все эти ошибки оказывают существенное влияние на качество моделирования. Это объясняет актуальность данного исследования, которое посвящено уменьшению влияния некоторых описанных выше ошибок на качество воспроизведения приземной температуры моделью WRF.
Целью данного исследования была оценка качества воспроизведения приземной температуры Санкт-Петербурга моделью WRF при различных конфигурациях модели и улучшение качества моделирования путём подбора оптимальной конфигурации.
Для достижения поставленной цели были определены следующие задачи:
1. Изучить мезомасштабную модель WRF.
2. Изучить интерпретационные языки программирования MATLAB и NCL для визуализации результатов моделирования.
3. Провести численные эксперименты с моделью WRF.
4. Произвести оценку качества моделирования и подобрать оптимальную конфигурацию модели для г. Санкт-Петербург.
Работа состоит из введения, пяти глав и заключения. В первой главе приводятся основные сведения о модели WRF. Во второй главе описывается процесс настройки модели и производится верификация результатов тестового моделирования. В третьей главе подбирается оптимальная комбинация схем параметризаций процессов подсеточного масштаба с целью улучшения качества моделирования. В четвёртой главе осуществляется подбор оптимальной конфигурации вычислительных сеток с целью увеличения качества моделирования. В пятой главе производится подбор эмпирических коэффициентов в схеме параметризации планетарного пограничного слоя также с целью улучшения качества моделирования. В конце работы приведён список использованных источников, состоящий из 44 наименований, а также два приложения.
В заключении сформулируем основные результаты и выводы по работе.
1. Во второй главе была подобрана вычислительная сетка, а также для того чтобы обеспечить устойчивость численного решения исходя из числа Куранта был подобран шаг по времени. Шаг по времени был равен 50 секундам.
2. Во второй главе был произведён тестовый запуск модели и произведена верификация результатов тестового моделирования.
3. В третьей главе по результатам верификации моделирования была подобрана оптимальная комбинация процессов подсеточного масштаба и тем самым улучшено качество моделирования приземной температуры воздуха в г. Санкт-Петербург.
4. В четвёртой главе было доказана неэффективность использования вложенных сеток в данной работе.
5. В пятой главе был подобран эмпирический коэффициент Ск, который используется для расчёта коэффициента турбулентного обмена импульсом и тем самым ещё немного улучшено качество моделирования приземной температуры в г. Санкт-Петербург.
Резюмируя всё вышеизложенное, в данной работе произведены оценки качества моделирования приземной температуры в г. Санкт-Петербург при различных конфигурациях модели и подобрана оптимальная конфигурация модели. Таким образом, модель оптимизирована для г. Санкт-Петербург и подготовлена для метеорологических исследований в данном регионе.
1. WRF-ARW V4: User’s Guide
2. https://www.ncdc.noaa.gov/data-access/model-data/model-datasets/global-forcast-system-gfs
3. http://meteo.ru/
4. Wicker, L. J., and W. C. Skamarock, 2002: Time-splitting methods for elastic models using forward time schemes. Mon. Wea. Rev., 130, 2088-2097
5. Наставление по службе прогнозов, 2017
6. Lin, Y.-L., R. D. Farley, and H. D. Orville, 1983: Bulk parameterization of the snow field in a cloud model. J. Climate Appl. Meteor., 22, 1065-1092.
7. Rutledge, S. A., and P. V. Hobbs, 1984: The mesoscale and microscale structure and organization of clouds and precipitation in midlatitude cyclones. XII: A diagnostic modeling study of precipitation development in narrow cloud-frontal rainbands. J. Atmos. Sci., 20, 2949-2972.
8. Tao, W.-K., J. Simpson, and M. McCumber 1989: An ice-water saturation adjustment, Mon. Wea. Rev., 117, 231-235.
9. Chen, S.-H., and W.-Y. Sun, 2002: A one-dimensional time dependent cloud model. J. Meteor. Soc. Japan, 80, 99-118.
10. Hong, S.-Y., J. Dudhia, and S.-H. Chen, 2004: A Revised Approach to Ice Microphysical Processes for the Bulk Parameterization of Clouds and Precipitation, Mon. Wea. Rev., 132, 103¬120.
11. Dudhia, J., S.-Y. Hong, and K.-S. Lim, 2008: A new method for representing mixed-phase particle fall speeds in bulk microphysics parameterizations. J. Met. Soc. Japan, 86, 33 -44.
12. Dudhia, J., 1989: Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model, J. Atmos. Sci., 46, 3077-3107.
13. Hong, S.-Y., H.-M. H. Juang, and Q. Zhao, 1998: Implementation of prognostic cloud scheme for a regional spectral model, Mon. Wea. Rev., 126, 2621-2639.
14. Thompson, G., P. R. Field, R. M. Rasmussen, and W. D. Hall, 2008: Explicit Forecasts of Winter Precipitation Using an Improved Bulk Microphysics Scheme. Part II: Implementation of a New Snow Parameterization. Mon. Wea. Rev., 136, 5095-5115.
15. Kain, J. S., and J. M. Fritsch, 1990: A one-dimensional entraining/ detraining plume model and its application in convective parameterization, J. Atmos. Sci., 47, 2784-2802.
16. Kain, J. S., and J. M. Fritsch, 1993: Convective parameterization for mesoscale models: The Kain-Fritcsh scheme, The representation of cumulus convection in numerical models, K. A. Emanuel and D.J. Raymond, Eds., Amer. Meteor. Soc., 246 pp.
17. Janjic, Z. I., 1994: The step-mountain eta coordinate model: further developments of the convection, viscous sublayer and turbulence closure schemes, Mon. Wea. Rev., 122, 927-945.
18. Betts, A. K., and M. J. Miller, 1986: A new convective adjustment scheme. Part II: Single column tests using GATE wave, BOMEX, and arctic air-mass data sets. Quart. J. Roy. Meteor. Soc., 112, 693-709.
19. Grell, G. A., and D. Devenyi, 2002: A generalized approach to parameterizing convection combining ensemble and data assimilation techniques. Geophys. Res. Lett., 29(14), Article 1693.
20. Mlawer, E. J., S. J. Taubman, P. D. Brown, M. J. Iacono, and S. A. Clough, 1997: Radiative transfer for inhomogeneous atmosphere: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res., 102 (D14), 16663-16682.
21. R.M. Goody, R. West, L. Chen, and D. Crisp. The correlated k-method for radiation calculations in nonhomogeneous atmosphere. J Quant Spectrosc Radiat Transfer, 42(6):539-550, 1989.
22. Cavallo, S. M., J. Dudhia and C. Snyder, 2011: A multi-layer upper boundary condition for longwave radiative flux to correct temperature biases in a mesoscale model. Mon. Wea. Rev., 139, 1952-1959.
23. Chou M.-D., and M. J. Suarez, 1999: A solar radiation parameterization for atmospheric studies. NASA Tech. Memo. 104606, 15, 40pp.
24. Chou, M. D., M. J. Suarez, X. Z. Liang, and M. M. H. Yan, 2001: A thermal infrared radiation parameterization for atmospheric studies. NASA Tech. Memo. 104606, 19, 68pp.
25. Matsui, T., S. Q. Zhang, W.-K. Tao, S. Lang, C. Ichoku, and C. Peters-Lidard, 2018: Impact of Radiation Frequency, Precipitation Radiative Forcing, and Radiation Column Aggregation on Convection-Permitting West African Monsoon Simulations. Clim. Dyn., 1-21, https://doi.org/10.1007/s00382-018-4187-2
26. Yang, P., L. Bi, B.A. Baum, K.-N. Liou, G.W. Kattawar, M.I. Mishchenko, and B. Cole, 2013: Spectrally consistent scattering, absorption, and polarization properties of atmospheric ice crystals at wavelengths from 0.2 to 100 pm. J. Atmos. Sci., 70, 330-347, doi:10.1175/JAS-D12-039.1
27. Радикевич В. М. Динамическая метеорология для океанологов, ЛГМИ, 1985, 157 c.
28. Hong, S.-Y., Y. Noh, and J. Dudhia, 2006: A new vertical diffusion package with an explicit treatment of entrainment processes. Mon. Wea. Rev., 134, 2318-2341.
29. Jimenez, P. A., and J. Dudhia, 2012: Improving the representation of resolved and unresolved topographic effects on surface wind in the WRF model. J. Appl. Meteor. Climatol., 51, 300- 316.
30. Lorente-Plazas, R., P. A. Jimenez, J. Dudhia, and J. P. Montavez, 2016: Evaluating and improving the impact of the atmospheric stability and orography on surface winds in the WRF model. Mon. Wea. Rev., 144, 2685-2693.
31. Wilson, T. H., and R. G. Fovell, 2018: Modeling the evolution and life cycle of radiative cold pools and fog. Weather and Forecasting. 33, 203-220.
32. Janjic, Z. I., 1996: The surface layer in the NCEP Eta Model, Eleventh Conference on Numerical Weather Prediction, Norfolk, VA, 19-23 August; Amer. Meteor. Soc., Boston, MA, 354-355. 32. Mellor and Yamada, 1982
33. Mellor, G. L., and T. Yamada, 1982: Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys., 20, 851-875.
34. Bougeault, P., and P. Lacarrere, 1989: Parameterization of orography-induced turbulence in a mesobeta-scale model. Mon. Wea. Rev., 117, 1872-1890.
35. Fairall, C., E. F. Bradley, J. Hare, A. Grachev, and J. Edson, 2003: Bulk parameterization of air-sea fluxes: Updates and verification for the COARE algorithm. J. Climate, 16, 571-591
36. Chen, F., and J. Dudhia, 2001: Coupling an advanced land-surface/ hydrology model with the Penn State/ NCAR MM5 modeling system. Part I: Model description and implementation. Mon. Wea. Rev., 129, 569-585.
37. Li, D., E. Bou-Zeid, M. Barlage, F. Chen, and J. A. Smith, 2013: Development and Evaluation of a Mosaic Approach in the WRF-Noah Framework. J. Geophys. Res., 118, 11918-11935.
38. Smirnova, T. G., J. M. Brown, and S. G. Benjamin, 1997: Performance of different soil model configurations in simulating ground surface temperature and surface fluxes. Mon. Wea. Rev., 125, 1870-1884.
39. Kusaka, H., H. Kondo, Y. Kikegawa, and F. Kimura, 2001: A simple single-layer urban canopy model for atmospheric models: Comparison with multi-layer and slab models. Bound.-Layer Meteor., 101, 329-358.
40. Kusaka, H. and F. Kimura, 2004: Coupling a single-layer urban canopy model with a simple atmospheric model: Impact on urban heat island simulation for an idealized case. J. Meteor. Soc. Japan, 82, 67-80.
41. Chen, F., M. Tewari, H. Kusaka, and T. T. Warner, 2006: Current status of urban modeling in the community Weather Research and Forecast (WRF) model. Joint with Sixth Symposium on the Urban Environment and AMS Forum: Managing our Physical and Natural Resources: Successes and Challenges, Atlanta, GA, USA, Amer. Meteor. Soc., CD-ROM. J1.4.41. Martilli et al., 2002
42. Martilli A, Clappier A, and Rotach M.W., 2002: An urban surface exchange parameterization for mesoscale models. Bound.-Layer Meteorol., 104, 261-304.
43. Salamanca, F., and A. Martilli, 2010: A new building energy model coupled with an urban canopy parameterization for urban climate simulations - part II. Validation with one dimension off-line simulations. Theor. Appl. Climatol., 99, 345-356.
44. Deardorff, J. W., 1972: Theoretical expression for the countergradient vertical heat flux. J. Geophys. Res., 77, 5900-5904.