A new design of the LAPS land surface scheme for use over and through heterogeneous and non-heterogeneous surfaces: Numerical simulations and tests

Numerical simulations and tests with the recently redesigned land–air parameterization scheme (LAPS) are presented. In all experiments, supported either by one-point micrometeorological, 1D or 3D simulations, the attention has been directed to: (1) comparison of simulation outputs, expressing the energy transfer over and through heterogeneous and non-heterogeneous surfaces, versus observations and (2) analysis of uncertainties occurring in the solution of the energy balance equation at the land–air interface. To check the proposed method for aggregation of albedo, “propagating hole” sensitivity tests with LAPS over a sandstone rock grid cell have been performed with the forcing meteorological data for July 17, 1999 in Baxter site, Philadelphia (USA). Micrometeorological and biophysical measurements from the surface experiments conducted over crops and apple orchard in Serbia, Poland, Austria and France were used to test the operation of LAPS in calculating surface fluxes and canopy environment temperatures within and above plant covers of different densities. In addition, sensitivity tests with single canopy covers over the Central Europe region and comparison against the observations taken from SYNOP data using 3D simulations were made. Validation of LAPS performances over a solid surface has been done by comparison of 2 m air temperature observations against 5-day simulations over the Sahara Desert rocky ground using 3D model. To examine how realistically the LAPS simulates surface processes over a heterogeneous surface, we compared the air temperature measured at 2 m and that predicted by the 1D model with the LAPS as the surface scheme. Finally, the scheme behaviour over urban surface was tested by runs over different parts of a hypothetical urban area. The corresponding 1D simulations were carried out with an imposed meteorological dataset collected during HAPEX-MOBILHY experiment at Caumont (France). The quantities predicted by the LAPS compare well with the observations and the various subcomponents of the scheme appear to operate realistically.     

The extreme value distribution of 5-min rainfall data at Belgrade

Summary For practical applications both the parent distribution of rainfall intensities and the distribution of their annual maxima are of interest. The relationship between these two distributions cannot be obtained from classical extreme value theory because of seasonal variation and serial correlation in the data. Mathematical results for the distribution of maxima in m-dependent sequences are presented to illustrate the effect of local dependence on the extreme value distribution. The average number of exceedances in a cluster is an important parameter in the relationship between the parent and the extreme value distribution. For 5-min rainfall data from Belgrade, quantiles of the annual maxima are overestimated by about 10 mm h^–1 if the effect of serial correlation is ignored. This bias can easily be removed by taking the local clustering of large rainfall intensities in a rainy spell into account.With 4 Figures     

Numerical case study of the Altai-Sayan lee cyclogenesis over East Asia

Summary Climatological studies show that the Altai-Sayan lee side is one of the major cyclogenesis areas in the Northern Hemisphere. In case of the Altai-Sayan lee cyclogenesis, the surface cyclone is generated when a primary cyclone is swept north of the mountains. In the mid-troposphere, a trough develops and finally turns into a cutoff low within 48 h. The main synoptic features are similar to those of Alpine cyclogenesis. Numerical simulations are performed to assess the effect of different representation of orography on the Altai-Sayan cyclogenesis. Two experiments are performed, a step-mountain (ETA) and an envelope orography (SGM) experiment. The ETA experiment produced the cyclogenesis in a way similar to that in the analysis both at the surface and at mid-troposphere. The SGM experiment failed in the simulation of the upper cutoff low. The difference in predicted pressure between the ETA and the SGM experiment shows a dipolar structure suggesting that the blocking effect of the mountains is essential in the development of the Altai-Sayan lee cyclogenesis.With 7 Figures     

Seeding agent dispersion within convective cloud as simulated by a 3-D numerical model

Summary The three-dimensional mesoscale cloud-resolving model ARPS (Advanced Regional Prediction System) was used to investigate the dispersal of an inert seeding agent within a cumulonimbus (Cb) cloud developing from two different initial states. In this paper, we stress the influence vortices in the cloud have on seeding agent dispersion. If a strong directional ambient wind shear is present in the lowest layer, a vortex pair formed at the flanks of the simulated cloud. Following the velocity field, a considerable amount of the injected seeding agent would be thrown out to the rear of the cloud, where both updrafts associated with vortices and downdrafts occurred. After a short time the agent was present only in the cloud periphery. If the Cb cloud developed under conditions where directional ambient wind shear did not exist, seeding agent dispersion would be quite different. In this case, almost all the seeding agent was transported into the main updraft, while the residence time of the agent within the cloud was longer due to the weaker cloud dynamics. Therefore, we must pay attention to whether or not the cloud contains vortices when we make the decision where to seed. This is necessary in order to minimize the loss of seeding material.     

Examination of physical processes of convective cell evolved from a MCS — Using a different model initialization

The present study is focused on examination of the physical processes of convective cell evolved from a MCS occurred on 4 November 2011 over Genoa, Italy. The Quantitative Precipitation Forecasts (QPF) have been performed using WRF v3.6 model under different configurations and cloud permitting simulations. The results indicate underestimation of the amount of precipitation and spatial displacement of the area with a peak 24-h accumulated rainfall in (mm). Our main objective in the research is to test the cloud model ability and performance in simulation of this particular case. For that purpose a set of sensitivity experiments under different model initializations and initial data have been conducted. The results also indicate that the merging process apparently alters the physical processes through low- and middle-level forcing, increasing cloud depth, and enhancing convection. The examination of the microphysical process simulated by the model indicates that dominant production terms are the accretion of rain by graupel and snow, probabilistic freezing of rain to form graupel and dry and wet growth of graupel. Experiment under WRF v3.6 model initialization has shown some advantage in simulation of the physical processes responsible for production and initiation of heavy rainfall compared to other model runs. Most of the precipitation came from ice-phase particles-via accretion processes and the graupel melting at temperature T₀ ≥ 0°C. The rainfall intensity and accumulated rainfall calculated by the model closely reflect the amount of rainfall recorded. Thus, the main benefit is to better resolve convective showers or storms which, in extreme cases, can give rise to major flooding events. In such a way, this model may become major contributor to improvements in weather analysis and small-scale atmospheric predictions and early warnings of such subscale processes.     

Spectral analysis of the “Koshava” wind

Summary “Koshava” is a gusty wind of changeable intensity, blowing from a south-easterly direction, over Serbia, Romania and Bulgaria. It is caused by the interaction between the synoptic circulation and the orography of the Carpathian and the Balkan mountains. This paper analyzes wind data measured at the Belgrade-Observatory during the longest period of consecutive days of “Koshava” which occurred from 14 January to 13 February 1972. Mean hourly wind speed data has been examined using spectral analysis. The power spectra are calculated using autocorrelation spectral analysis, the multi-taper method and wavelet transform. The maximum of which is about 122 h (5 days) corresponds to the time span of synoptic processes.     

Numerical Study of Winter Diurnal Convection Over the City of Krasnoyarsk: Effects of Non-freezing River,Undulating Fog and Steam Devils

We performed a numerical simulation of penetrative convection of an inversion-topped weakly stratified atmospheric boundary layer over urban terrain with a strong localized source of heat and moisture. With some simplifications, the case mimics the real environment of the Krasnoyarsk region in Russia where the non-freezing river Yenisei acts as a thermal and humidity source during winter, generating an undulating fog pattern along the river accompanied with scattered ‘steam devils’. An idealized full diurnal cycle was simulated using an unsteady Reynolds-averaged Navier–Stokes (RANS) three-equation algebraic flux model and the novel buoyancy-accounting functions for treating the ground boundary conditions. The results show a significant effect of the river on the net temperature and moisture distribution. The localized heat and moisture source leads to strong horizontal convection and marked non-uniformity of humidity concentration in the air. An interplay of several distinct large-scale vortex systems leads to a wavy pattern of moisture plumes over the river. The simulations deal with rare natural phenomena and show the capability of the RANS turbulence closure to capture the main features of flow and scalar fields on an affordable, relatively coarse, computational grid.     

A note on the performance of the multiply-upstream semi-Lagrangian advection schemes for one-dimensional nonlinear momentum conservation equation

Summary Two-time-level multiply-upstream semi-Lagrangian schemes were examined in the case of the self-advecting, one-dimensional nonlinear momentum conservation equation. The shock formation process was analyzed. It is pointed out that the shocks cannot be created in the truncated systems satisfying the Pudykiewicz, Benoit and Staniforth criterion.The numerical integrations were restricted to 12 h. It was shown that, at least in the sub-CFL range, increased complexity of the scheme can compensate reduced horizontal resolution. A considerable sensitivity of the schemes with respect to the time step was detected. In the super-CFL mode, several windows on various time scales were found within which the Pudykiewicz, Benoit and Staniforth criterion was satisfied. The time step of 1.44 times the maximum time step allowed by the CFL criterion was used in the semi-Lagrangian runs.The super-CFL, semi-Lagrangian solutions were diverging progressively from the sub-CFL ones as the forecasts advanced. This was also reflected in the energy spectra.Unacceptably large energy losses were encountered in the super-CFL, semi-Lagrangian runs. Most of these losses could be explained by the reduced mean wind speed, i.e., the amplitude of the zero wavenumber wave. At the same time, the energy content in the shorter waves increased. In a more complex model, such a situation would resemble a loss of zonal, and an increase of transient eddy kinetic energy.A trajectory error measure was defined as the maximum absolute value of the distance between the actual arriving point of the particle originating at the estimated departure point, and the grid point assumed to be the arrival point in the semi-Lagrangian procedure. In contrast to the sub-CFL regime, this measure could reach a considerable fraction of the grid distance in the computations with the super-CFL time steps.In the physical system considered, the trajectories are determined only by the velocities at the departure points. With the semi-Lagrangian schemes the distances traveled by the particles are estimated on the basis of the velocities at the points downstream with respect to the departure points. Thus, unless the solution is smooth (in space and time) on the scales of the extrapolation distances/times, the upstream extrapolation does not promise the convergence of the solution.With 16 Figures     

Changes in Components of the Water Balance in the Croatian Lowlands

Summary ?This paper deals with variations and trends in some components of the water balance: the soil water content; evaporation loss from the soil (from the surface and underlying layers); transpiration; ground water recharge and runoff. These components are calculated by means of the Palmer procedure. This analysis is based on data from Osijek, Croatia from this century (1900–1995). Besides the meteorological input parameters necessary for the water balance calculations, i.e. precipitation, temperature and relative humidity, the pedological characteristics of this area have also been taken into account. Fluctuations have been considered by means of the 11-year binomial filtered series and linear trends were tested by means of the Mann-Kendall rank test. For a closer look on the trends of water balance components, a progressive analysis of the time series was performed, too. The results show a significant increase in potential evapotranspiration and evapotranspiration and decrease in runoff and soil water content has occurred during the century. Received February 22, 1999/Revised August 3, 1999