Numerical modeling of thunderstorm

Summary A brief review of numerical models of thunderstorms is presented. The attempts of modeling these clouds have to deal with an extremely broad range of microphysical and dynamic processes. Basic concepts of thunderstorm modeling, 1-D, 2-D and 3-D models are considered. Particular attention is given to the basic equations and to the capability of these models to predict characteristics of the tunderstorms. A brief treatment is also given of hail growth in thunderstorm models.

---

Zusammenfassung Im folgenden wird eine knappe Übersicht numerischer Gewittermodelle präsentiert. Bei der Modellierung von Gewitterwolken muß ein ausgesprochen weiter Bereich von mikrophysikalischen und dynamischen Prozessen beachtet werden. Grundkonzepte der Gewittermodellierung, 1-D-, 2-D- und 3-D-Modelle werden hierzu herangezogen. Besondere Aufmerksamkeit gilt dabei den Ausgangsgleichungen sowie den potentiellen Möglichkeiten dieser Modelle, Charakteristika der Gewitter vorherzusagen. Weiters wird auch das Hagelwachstum in Gewittermodellen kurz behandelt.

---

---

With 8 Figures     

Rujevac Sb-Pb-Zn-As polymetallic deposit,Boranja orefield,Western Serbia: native arsenic and arsenic mineralization

Rujevac is a low-temperature hydrothermal polymetallic Sb-Pb-Zn-As vein-type ore deposit, hosted within a volcanogenic-sedimentary zone situated in the Rujevac-Crvene Stene-Brezovica Diabase-Chert Formation (DCF) of the Podrinje Metallogenic District (PMD), Serbia. It is located several kilometers SE from the Boranja contact aureole, which is an integral part of the PMD in Western Serbia. Genetically related to the Tertiary granodioritic magma, the mineral assemblages are characterized by specific features. The mineral association of this deposit consists of sulfides, Pb-Sb(As) sulfosalts, native metals, oxides, hydroxides and gangue minerals. Chemical composition of the ore is very complex, where contents of valuable metals range as follows: Sb (0.17–24.31 wt.%), Zn (0.21–6.29 wt.%), Pb (0.15–6.33 wt.%), As (0.06–1.28 wt.%), Cd (25–747 ppm), Ag (7–408 ppm), Hg (13–473 ppm), and Tl (<1–29 ppm). Electron Probe Microanalyses (EPMA) of native arsenic from both the Rujevac and Stragari deposits showed contents of As up to 98.8 and 97.1 wt.%, with impurity contents of Sb up to 1.3 and 6.6 wt.%, and Tl up to 2 and 1.3 wt.%, respectively. Rhombohedral unit-cell parameters for native arsenic from Rujevac and Stragari deposits amount to: a?=?3.760(2), c?=?10.555(3) Å, V?=?129.23(7) ų and a?=?3.763(1), c?=?10.560(5) Å, V?=?129.48(8) ų, respectively. Mineral assemblages, deposition order and genesis of the Rujevac polymetallic deposit were also discussed in detail. Native arsenic mineralization here has been additionally compared with similar well-known global deposits.     

Magma hybridization in the Western Tatra Mts. granitoid intrusion (S-Poland, Western Carpathians)

In the Variscan Western Tatra granites hybridization phenomena such as mixing and mingling can be observed at the contact of mafic precursors of dioritic composition and more felsic granitic host rocks. The textural evidence of hybridization include: plagioclase?CK-feldspar?Csphene ocelli, hornblende- and biotite-rimmed quartz ocelli, plagioclase with Ca-rich spike zonation, inversely zoned K-feldspar crystals, mafic clots, poikilitic plagioclase and quartz crystals, mixed apatite morphologies, zoned K-feldspar phenocrysts. The apparent pressure range of the magma hybridization event was calculated at 6.1?kbar to 4.6?kbar, while the temperature, calculated by independent methods, is in the range of 810°C?770°C. U-Pb age data of the hybrid rocks were obtained by in-situ LA-MC-ICP-MS analysis of zircon. The oscillatory zoned zircon crystals yield a concordia age of 368?±?8?Ma (MSWD?=?1.1), interpreted as the age of magma hybridization and timing of formation of the magmatic precursors. It is the oldest Variscan magmatic event in that part of the Tatra Mountains.     

Spatial and temporal variability of winter and summer precipitation over Serbia and Montenegro

Summary The main characteristics of the spatial and temporal variability of winter and summer precipitation observed at 30 stations in Serbia and Montenegro were analysed for the period 1951–2000. The rainfall series were examined spatially by means of Empirical Orthogonal Functions (EOF) and temporally by means of the Mann-Kendall test and spectral analysis. The Alexandersson test was used to detect the inhomogeneity of the data set.The EOF analysis gave three winter and summer dominant modes of variations, which explained 89.7% and 70.4% of the variance, respectively. The time series associated with the first pattern showed a decreasing trend in winter precipitation. The spectral analysis showed a 16-year oscillation for the dominant winter pattern, around a 3-year oscillation for the dominant summer pattern, and a quasi-cycle of 2.5 years for the winter third pattern.     

Opening a new window to other worlds with spectropolarimetry

A high level of diversity has already been observed among the planets of our own Solar System. As such, one expects extrasolar planets to present a wide range of distinctive features, therefore the characterisation of Earth- and super Earth-like planets is becoming of key importance in scientific research. The Search (Spectropolarimetric Exoplanet AtmospheRe CHaracerisation) mission proposal of this paper represents one possible approach to realising these objectives. The mission goals of Search include the detailed characterisation of a wide variety of exoplanets, ranging from terrestrial planets to gas giants. More specifically, Search will determine atmospheric properties such as cloud coverage, surface pressure and atmospheric composition, and may also be capable of identifying basic surface features. To resolve a planet with a semi major axis of down to 1.4 AU and 30 pc distant Search will have a mirror system consisting of two segments, with elliptical rim, cut out of a parabolic mirror. This will yield an effective diameter of 9 m along one axis. A phase mask coronagraph along with an integral spectrograph will be used to overcome the contrast ratio of star to planet light. Such a mission would provide invaluable data on the diversity present in extrasolar planetary systems and much more could be learned from the similarities and differences compared to our own Solar System. This would allow our theories of planetary formation, atmospheric accretion and evolution to be tested, and our understanding of regions such as the outer limit of the Habitable Zone to be further improved.     

Origin of siderite mineralisation in Petrova and Trgovska Gora Mts., NW Dinarides

The Petrova and Trgovska Gora Mts. (Gora=Mountain) are Variscan basement units incorporated into the northwestern Dinarides during the Alpine orogeny. They host numerous siderite-quartz-polysulphide, siderite-chalcopyrite, siderite-galena and barite veins, as well as stratabound hydrothermal-replacement ankerite bodies within carbonates in non-metamorphosed, flysch-like Permo-Carboniferous sequences. The deposits have been mined for Cu, Pb, Ag and Fe ores since Medieval times. Fluid inclusion studies of quartz from siderite-polysulphide-quartz and barite veins of both regions have shown the presence of primary aqueous NaCl?CaCl₂±MgCl₂?H₂O±CO₂ inclusions. The quartz-sulphide stage of both regions show variable salinities; 2.7–26.2 wt% NaCl eq. for the Trgovska Gora region and 3.4–23.4 wt% NaCl eq. for the Petrova gora region, and similar homogenisation temperatures (100–230°C). Finally, barite is precipitated from low salinity-low temperature solutions (3.7–15.8 wt % NaCl equ. and 115–145°C). P-t conditions estimated via isochore construction yield formation temperatures between 180–250°C for the quartz-sulphide stage and 160–180°C for the barite stage, using a maximum lithostatic pressure of 1 kbar (cc. 3 km of overburden). The sulphur isotope composition of barite from both deposits indicates the involvement of Permian seawater in ore fluids. This is supported by the elevated bromium content of the fluid inclusion leachates (120–660 ppm in quartz, 420–960 ppm in barite) with respect to the seawater, indicating evaporated seawater as the major portion of the ore-forming fluids. Variable sulphur isotope compositions of galena, pyrite and chalcopyrite, between ?3.2 and +2.7‰, are interpreted as a product of incomplete thermal reduction of the Permian marine sulphate mixed with organically- and pyrite-bound sulphur from the host sedimentary rocks. Ore-forming fluids are interpreted as deep-circulating fluids derived primarily from evaporated Permian seawater and later modified by interaction with the Variscan basement rocks. ⁴⁰Ar/³⁹Ar data of the detrital mica from the host rocks yielded the Variscan age overprinted by an Early Permian tectonothermal event dated at 266–274 Ma. These ages are interpreted as those reflecting hydrothermal activity correlated with an incipient intracontinental rifting in the Tethyan domain. Nevertheless, 75 Ma recorded at a fine-grained sericite sample from the alteration zone is interpreted as a result of later resetting of white mica during Campanian opening/closure of the Sava back arc in the neighbouring Sava suture zone (Ustaszewski et al. 2008).     

Water losses from a reservoir built in karst: the example of the Boljunčica reservoir (Istria,Croatia)

The paper presents the case of the Boljunčica reservoir, which began operation in 1973. It is situated on the Istria peninsula (Croatia). This is a multipurpose reservoir which was built in order to protect the downstream area from flood, to store water for irrigation, and to control sediment transport. The reservoir is situated on the contact zone between water impermeable Eocene flysch and deep Eocene and Cretaceous limestone. The bottom of the reservoir is covered partly by both flysch and quaternary deposits. Water losses from the reservoir bottom are so large that the main service intended for the reservoir, the storage of water for irrigation, is impossible. After every intensive precipitation, which occurs often in this region, the reservoir fills very quickly. The problem is that its retention of water is very short, and lasts only a few days. The water volume of the reservoir at the spillway altitude of 93.00 m a. s. l. is about 6.5 × 10⁶ m³. Because of water losses from the reservoir bottom, the mean annual volume of water stored in it, during the period of 1977–2005 was only 0.5 × 10⁶ m³, which is less than 8% of the full reservoir volume. On the reservoir bottom, many new swallow holes opened through the sediment cover after each time it filled and emptied with water. Special attention is paid to the groundwater level analyses. Interdisciplinary analyses and investigations of hydrological and hydrogeological factors causing the formation of swallow holes and water losses from the Boljunčica reservoir are discussed. The example given in this paper explains one unsuccessful case of building a reservoir in Dinaric karst, caused mainly due to insufficient geological, hydrogeological and hydrological investigations. In order to prevent water losses from the Boljunčica reservoir, very complex and expensive work needs to be done, but its success regarding the reasonable reduction of water losses from the reservoir is in question.     

Slope instability mechanisms in dipping conglomerates over weathered marls: Bol landslide,Croatia

In February 2005, a landslide of significant dimensions occurred at Bol on the Island of Brač, in Croatia. The location of the landslide was in a geological structure of weakly bounded conglomerates lying above a flysch, in an area where numerous slides have been previously reported. In the paper the results of a geological engineering investigation of the Bol landslide are presented. The investigation results show that the landslide was triggered by a combination of improperly started earthworks, an increase in air temperature in days preceding landslide that caused snow melt, and abundant rainfall, but the primary cause of the landslide was weathering of the flysch. In other words, despite the fact that contact between the flysch and the conglomerates slopes is in approximate correspondence with the slope of the hillside, the landslide surface was not in contact between the conglomerates and marl, but deeper in the flysch, provoked by the penetration of the weathering front into the flysch. Numerical modelling results were able to reproduce this mechanism by incorporating strength degradation into the weathered zone, determined by field investigation. The results of the research can be used to predict potential slides not only in the area under study.     

Spatial and temporal analysis of dry spells in Croatia

Systematic statistical analysis of dry day sequences, which are defined according to 0.1, 1, 5 and 10 mm of precipitation-per-day thresholds, is performed on seasonal and yearly basis. The data analysed come from 25 Croatian meteorological stations and cover the period 1961–2000. Climatological features of the mean and maximum dry spell durations, as well as the frequency of long dry spells (>20 days) are discussed. The results affirm the three main climatological regions in Croatia, with the highlands exhibiting shorter dry spells than the mainland, and the coastal region exhibiting longer dry spells. The prevailing positive trend of both mean and maximal durations is detected during winter and spring seasons, while negative trend dominate in autumn for all thresholds. Positive field significant trends of mean dry spell duration with 5 and 10 mm thresholds are found during spring and the same is valid for annual maximum dry spell duration with a 10 mm threshold. It is found that the Discrete Autoregressive Moving Average (DARMA(1,1)) model can be used to estimate the probabilities of dry spells in Croatia that are up to 20–30 days long.     

Dependence of accumulated precipitation on cloud drop size distribution

Convective precipitation is the main cause of extreme rainfall events in small areas. Its primary characteristics are both large spatial and temporal variability. For this reason, the monitoring of accumulated precipitation fields (liquid and solid components) at the surface is difficult to carry out through the use of rain gauge networks or remote-sensing observations. Alternatively, numerical models seem to be the most powerful tool in simulating convective precipitation for various analyses and predictions. Due to a lack of comparisons between modelled and observed precipitation characteristics over a long period of time, we focus our research on comparisons between observations and three model samples of accumulated convective precipitation over a particular study area. We use a numerical cloud model with two model schemes involving the unified Khrgian–Mazin size distribution of cloud drops and a model scheme involving a monodisperse cloud droplet spectrum and the Marshall–Palmer size distribution for raindrops, respectively. For comparison, we have selected a study area with a sounding site. Our analysis shows that the model version with the Khrgian–Mazin size distribution exhibits a better agreement with the observed mean, median and range of extreme values of accumulated convective precipitation. Model simulations with the Khrgian–Mazin size distribution most closely match observations, with a correlation coefficient of 0.91. Use of the Marshall–Palmer size distribution, on the other hand, systemically underestimates the observed precipitation and has the lowest correlation coefficient among the methods, 0.83. Such an investigation is crucial to improve predictions of accumulated convective precipitation for various climatological and hydrological analyses and predictions.