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     

An upgraded version of the Eta model

Upgrades implemented over a number of years in an open source version of the Eta model, posted at the CPTEC web site http://etamodel.cptec.inpe.br/, are summarized and examples of benefits are shown. The version originates from the NCEP’s Workstation Eta code posted on the NCEP web site http://www.emc.ncep.noaa.gov/mmb/wrkstn_eta, which differs from the NCEP’s latest operational Eta by having the WRF-NMM nonhydrostatic option included. Most of the upgrades made resulted from attention paid to less than satisfactory performance noted in several Eta results, and identification of the reasons for the problem. Others came from simple expectation that including a feature that is physically justified but is missing in the code should help. The most notable of the upgrades are the introduction of the so-called sloping steps, or discretized shaved cells topography; piecewise-linear finite-volume vertical advection of dynamic variables; vapor and hydrometeor loading in the hydrostatic equation, and changes aimed at refining the convection schemes available in the Eta. Several other modifications have to do with the calculation of exchange coefficients, conservation in the vertical diffusion, and diagnostic calculation of 10-m winds. Several examples showing improved performance resulting from the dynamics changes are given. One includes a case of unrealistically low temperatures in several mountain basins generated by a centered vertical advection difference scheme’s unphysical advection from below ground, removed by its replacement with a finite-volume scheme. Another is that of increased katabatic winds in the Terra Nova Bay Antarctica region. Successful forecast of the severe downslope zonda wind case in the lee of the highest peaks of the Andes is also shown, and some of the recent successful verification results of the use of the upgraded model are pointed out. The code is used at numerous places, and along with setup information it is available for outside users at the CPTEC Eta web site given above.     

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.     

Evaluation of a high-resolution regional climate simulation over Greenland

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.     

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.     

Tropospheric NO2 Columns over Northeastern North America： Comparison of CMAQ Model Simulations with GOME Satellite Measurements

We present comparisons of the NO2 regional Chemical Transport Model （CTM） simulations over North-eastern North America during the time period from May to September, 1998 with hourly surface NO2 observations and the NO2 columns retrieved from the GOME （Global Ozone Monitoring Experiment） satellite instrument. The model calculations were performed using the Mesoscale Meteorological Model 5 （MM5）, Sparse Matrix Operator Kernal Emissions （SMOKE）, and Community Multiscale Air Quality （CMAQ） modeling systems, using the emission data from the National Emissions Inventory （NEI） databases of 1996 （U.S.） and 1995 （Canada）. The major objectives were to assess the performance of the CMAQ model and the accuracy of the emissions inventories as they affected the simulations of this important short-lived atmospheric species. The modeled （NcMAQ） and measured （NGOME） NO2 column amounts, as well as their temporal variations, agreed reasonably well. The absolute differences （NcMAQ-NGOME） across the domain were between ±3.0×10^15 molecules cm^-2, but they were less than ±1.0×10^15 molecules cm^-2 over the majority （80%） of the domain studied. The overall correlation coefficient between the measurements and the simulations was 0.75. The differences were mainly ascribed to a combination of inaccurate emission data for the CTM and the uncertainties in the GOME retrievals. Of these, the former were the more easily identifiable.     

Ice-sheet growth and high-latitudes sea-surface temperature

The response of the LLN 2-D climate model to the insolation and CO₂ forcings during the Eemian interglacial is compared to reconstructions obtained from deep-sea cores drilled in the Norwegian Sea and in the North Atlantic. Both reconstructions and modeling results show a decrease of sea-surface temperature (SST) in the higher latitudes (70–75 °N zonal belt for the model and the Norwegian Sea for the proxy records), associated with a more moderate cooling at lower latitudes (50–55 °N and North Atlantic), at the middle of isotopic substage 5e, several millenia before the beginning of continental ice-sheet growth. Such a comparison between the simulated SST and ice volume of the Northern Hemisphere has been extended to the whole last glacial-interglacial cycle. The influence of the insolation forcing on SST and the shortcomings of the model due to its zonal character are discussed. Received: 6 July 1995/Accepted: 19 December 1995