250 likes | 370 Views
NUMERICAL SIMULATION CLOUDS AND PRECIPITATION CAUSED CATASTROPHIC FLOODS ALONG THE ELBE RIVER IN AUGUST 2002. Krakovskaia S.V. , Palamarchuk L.V. and Shpyg V.M. Ukrainian Hydrometeorological Research Institute, Atmosphere Physics Department, Kiev, Ukraine KraSvet@antarc.icyb.kiev.ua.
E N D
NUMERICAL SIMULATION CLOUDS AND PRECIPITATION CAUSED CATASTROPHIC FLOODS ALONG THE ELBE RIVER IN AUGUST 2002 Krakovskaia S.V., Palamarchuk L.V. and Shpyg V.M. Ukrainian Hydrometeorological Research Institute, Atmosphere Physics Department, Kiev, Ukraine KraSvet@antarc.icyb.kiev.ua
Objective: • To carry out numerical experiments to investigate precipitation formation processes in the frontal cloud system caused catastrophic floods along the Elbe river in August 2002. Main tasks: • To fit combined model of cloudy troposphere (CMCT) to study summer clouds by varying intensity of microphysical precipitation formation processes. • To analyze mesostructure of thermodynamical characteristics with aim to find zones of instability and water vapor resource determined intense precipitation formation in mixed clouds. • To verify CMCT outputs by comparing with the precipitation measurements. • To estimate precipitation forecasting ability of CMCT.
COMBINED MODEL OF CLOUDY TROPOSPHERE (CMCT) CMCT is the combination of 3-D time-independent LAM and 1-D CRM with explicit microphysics, when initial thermodynamical characteristics in the vertical column where microphysics is calculated are continuously updated as it moves along horizontal axes over the initial point of 3-D domain at every time step (dt) in 1-D model on: Minus in equations means that the system moves opposite to the air mass displacement. At the same moment the total column movement in relation to initial point on 3-D domain (X0, Y0) will be: The speed uf,, vfof this movement could be constant and determined by synoptic charts or variable and calculated from the data of 3-D LAM as an average wind speed in layer up to height Zf, that corresponds to level n:
Ice supersaturation and wet-instability [pink] Vertical motions, cm/s
U-wind projection, m/s V-wind projection, m/s
Integral thermodynamic condensation rate as precipitation intensity (mm/h) and ice supersaturation (mm) as vapor reservoir 00 UTC 12.08.2002 06 UTC 12.08.2002
Integral thermodynamic condensation rate as precipitation intensity (mm/h) and ice supersaturation (mm) as vapor reservoir 12 UTC 12.08.2002 18 UTC 12.08.2002
Treks of 1D on 3D domainstarted 00 UTC 12.08.2002 Treks of 1D on 3D domainat 06 UTC 12.08.2002
1-D CLOUD RESOLVING MODEL 1. The kinetic equation of cloud droplets (k=1) size distribution function (f1): 2. The kinetic equation of cloud ice crystals (k=2) and raindrops (k=3) size distribution functions (fk): c2 = 1, c3 = 0 3. The equation of heat inflow:4. The equation of moisture inflow: 5. The state equation:
1-D CLOUD RESOLVING MODEL Cloud particles generation on CCN and IN is parameterized as follows: ifx > 0,ifx < 0 The rate of growth of individual particles due to condensation (sublimation) is: and due to gravitational collection by spherical CP of droplets: , k = 2, 3, whererk0aremin radii of CP; Е(r1,rk)arecoefficients of collection.
FREEZING PROCESSES IN 1-D CRM The processes of droplets (If1) and raindrops (If3) freezing and the amount of frozen water CP (If2) is parameterized as follows: , Cloud particles concentrations (Nk), their average sizes (r), ice and water contents (qk), intensity of precipitation (j): whererkminandrkmaxare minimum and maximum radii of CP. INTEGRAL CHARACTERISTICS OF MIXED CLOUD
5-min precipitation totals for Dresden 12 August, 00 UTC, to 13 August, 06 UTC 5-min precipitation totals for Doberlug 12 August, 00 UTC, to 13 August, 06 UTC
5-min precipitation totals for Chemnits 12 August, 00 UTC, to 13 August, 06 UTC 5-min precipitation totals for Cottbus 12 August, 00 UTC, to 13 August, 06 UTC
5-min precipitation totals for Cottbus 12 August, 06 UTC, to 13 August, 06 UTC
CONCLUSIONS • CMCT allowed to define zones with potentially intensive precipitation formation processes and to calculate probable precipitation intensity for desired points of 3-D domain. • Certain limitations in running CMCT on the concrete cloud system were defined: intensive updrafts with moisture advection resulted in constantly renewable water vapour resources in clouds which leaded to overestimated precipitation intensities and sums. • Varying of probabilities of activation of cloud nuclei by changing their parameters and even excluding of one of precipitation formation processes modified predictable precipitation sum on 10-20%. • It is necessary to define parameters of microphysical processes in dependence on thermodynamical conditions of cloud formation for getting meaningful precipitation forecast for concrete point and time. • Reliability of the model outputs directly depends on quantity and quality of inputs. Upper-air soundings minimum every 6 h are desirable to construct CMCT and to get reliable results.
3-D LAM results with 1-D CRM trek (digits near trek – tcloud, h) Precipitation rate (mm/h) Integral ice supersaturation (mm)
Initial vertical and time development of: (a) temperature (lines, oC), water (green) and ice supersaturations (shading, mg/kg); (b) vertical motions (cm/s) Vertical and time development of: (a) IC and (b) LWC (mg/kg) In the run with unmodified freezing