Zircon geochronology and geochemistry to constrain the youngest eruption events and magma evolution of the Mid-Miocene ignimbrite flare-up in the Pannonian Basin,eastern central Europe

A silicic ignimbrite flare-up episode occurred in the Pannonian Basin during the Miocene, coeval with the syn-extensional period in the region. It produced important correlation horizons in the regional stratigraphy; however, they lacked precise and accurate geochronology. Here, we used U–Pb (LA-ICP-MS and ID-TIMS) and (U–Th)/He dating of zircons to determine the eruption ages of the youngest stage of this volcanic activity and constrain the longevity of the magma storage in crustal reservoirs. Reliability of the U–Pb data is supported by (U–Th)/He zircon dating and magnetostratigraphic constraints. We distinguish four eruptive phases from 15.9 ± 0.3 to 14.1 ± 0.3 Ma, each of which possibly includes multiple eruptive events. Among these, at least two large volume eruptions (>10 km³) occurred at 14.8 ± 0.3 Ma (Demjén ignimbrite) and 14.1 ± 0.3 Ma (Harsány ignimbrite). The in situ U–Pb zircon dating shows wide age ranges (up to 700 kyr) in most of the crystal-poor pyroclastic units, containing few to no xenocrysts, which implies efficient recycling of antecrysts. We propose that long-lived silicic magma reservoirs, mostly kept as high-crystallinity mushes, have existed in the Pannonian Basin during the 16–14 Ma period. Small but significant differences in zircon, bulk rock and glass shard composition among units suggest the presence of spatially separated reservoirs, sometimes existing contemporaneously. Our results also better constrain the time frame of the main tectonic events that occurred in the Northern Pannonian Basin: We refined the upper temporal boundary (15 Ma) of the youngest counterclockwise block rotation and the beginning of a new deformation phase, which structurally characterized the onset of the youngest volcanic and sedimentary phase.     

Partitioning of Ru, Rh, Pd, Re, Ir, and Au between Cr-bearing spinel, olivine, pyroxene and silicate melts

A series of high temperature experiments was undertaken to study partitioning of several highly siderophile elements (HSE; Ru, Rh, Pd, Re, Os, Ir, Pt and Au) between Cr-rich spinel, olivine, pyroxene and silicate melt. Runs were carried out on a Hawaiian ankaramite, a synthetic eucrite basalt, and a DiAn eutectic melt, at one bar, 19 kbar, and 20 kbar, respectively, in the temperature range of 1200 to 1300°C, at oxygen fugacities between the nickel-nickel oxide (NNO) and hematite-magnetite (HM) oxygen buffers. High oxygen fugacities were used to suppress the formation of HSE-rich “nuggets” in the silicate melts. The resulting oxide and silicate crystals (<100 μm) were analyzed using both SIMS and LA-ICP-MS, with a spatial resolution of 15 to 50 μm. Rhenium, Au and Pd were all found to be incompatible in Cr-rich spinel (D_Re^sp/melt = 0.0012-0.21, D_Au^sp/melt = 0.076, D_Pd^sp/melt = 0.14), whereas Rh, Ru and Ir were all found to be highly compatible (D_Rh^sp/melt = 41-530, D_Ru^sp/melt = 76-1143, D_Ir^sp/melt = 5-22000). Rhenium, Pd, Au and Ru were all found to be incompatible in olivine (D_Re^oliv/melt = 0.017-0.073, D_Pd^oliv/melt = 0.12, D_Au^oliv/melt = 0.12, D_Ru^oliv/melt = 0.23), Re is incompatible in orthopyroxene and clinopyroxene (D_Re^opx/melt = 0.013, D_Re^cpx/melt = 0.18-0.21), and Pt is compatible in clinopyroxene (D_Pt^cpx/melt = 1.5). The results are compared to and combined with previous work on HSE partitioning among spinel-structured oxides, and applied to some natural magmatic suites to demonstrate consistency.     

Landscape planning in the Transcaucasia

We have developed the planning documents for model territories of a different dimension: the landscape program of Adzharia (Georgia), the landscape structure plan for the Lake Sevan watershed basin (Armenia), and the landscape plan for the Shirvan National Park (Azerbaijan). Variants of the solution of regional natural-ecological problems are proposed, and the measures for implementation of the target concepts are substantiated.     

SHAPED CHARGE—A POWERFUL EXPLOSIVE SURFACE SEISMIC SOURCE *

The directional effect of shaped charge is a well-known feature used for a long time in military weapons, oil well guns, and steel industry. This charge was successfully applied as a seismic energy source by Petrobrás during the past three years under different surface and geological conditions. Preliminary amplitude measurements taken with fixed gain shallow refraction instruments showed a consistent difference between conventional and shaped charges. Lately, a similar difference has been noted in deep reflection energy recorded with digital binary gain instruments as well as in deep oil well velocity surveys. Direct comparisons along more than 50 km of multiple coverage field reflection shooting are in agreement with these results and have proved the practical advantage of this source as compared to conventional dynamite. This source has been used since 1971 in routine seismic operation in the Amazon jungle with 300 gram unit charges distributed in small and large shot arrays increasing substantially the coverage and halving the cost at a higher record quality. A large amount of production seismic field work has been carried out in several other areas attesting the successful application of the shaped charges.     

Utilization of data bank on grid square basis for regional planning

This paper deals with results of case studies using actual magnetic tape records of Statistics on Grid Square Basis and Digital National Land Information which have been compiled since 1965 by governmental agencies on the same grid square basis from census and map reading covering the whole country. The number of items on the records has been increasing, including land information and socio-economic statistics which seem to be useful for land development/conservation planning and regional planning. However, in reality a majority of local governments have not used to date these records enough. In this study the combination analysis of two or more items is made in order to make more effective use of the records for these plannings. The number of grid squares under fixed conditions are calculated and the distribution maps are drawn by means of line printer or X-Y plotter. The themes are as follows: 1. Population distribution and population growth and land use type by landform classification, 2. Proportion of built-up area in city planning area, 3. Selection of possible space for housing development, 4. Division of a given area based on natural conditions for daily life into three areas, namely suitable areas, possible but rather difficult areas, and unsuitable areas. 5. Selection of suitable areas for sightseeing and recreation. The items used for combination analysis are 1970 and 1975 Population Census, ground elevation, angle of inclination, value of relief, landform classification, surface geology classification, soil classification, land use classification, and city planning area. The data of land value and distance from nearest station are produced by map reading for the test area in this study. From these combinations an objective and detailed distribution can be determined. Finally, several examples of the studies are shown.     

Distribution of Mantle Up-welling Determined from Plate Motions: A Case for Large-scale Bénard Convection

—?The number and geometric distribution of putative mantle up-welling centers and the associated convection cell boundaries are determined from the lithospheric plate motions as given by the 14 Euler poles of the observational NUVEL model. For an assumed distribution of up-welling centers (called “cell-cores”) the corresponding cell boundaries are constructed by a Voronoi division of the spherical surface; the resulting polygons are called “Bénard cells.” By assuming the flow-kinematics within a cell, the viscous coupling between the flow and the plates is estimated, and the Euler poles for the plates are computed under the assumption of zero-net-torque. The positions of the cell-cores are optimized for the HS2-NUVEL1 Euler poles by a method of successive approximation (“subplex”); convergence to one of many local minima occurred typically after ～20,000 iterations. Cell-cores associated with the fourteen HS2-NUVEL1 Euler poles converge to a relatively small number of locations (8 to 10, depending on interpretation), irrespective of the number of convection cells submitted for optimized distribution (from 6 to 50). These locations are correlated with low seismic propagation velocities in tomography, uniformly occur within hotspot provinces, and may specifically be associated with the Hawaiian, Iceland, Reunion/Kerguelen (Indian Ocean), Easter Island, Melanesia/Society Islands (South Pacific), Azores/Cape Verde/Canary Islands, Tristan da Cunha (South Atlantic), Balleny Islands, and possibly Yellowstone hotspots. It is shown that arbitrary Euler poles cannot occur in association with mantle Bénard convection, irrespective of the number and the distribution of convection cells. Nevertheless, eight of the observational Euler poles – including the five that are accurately determined in HS2-NUVEL1 (Australia, Cocos, Juan de Fuca, Pacific, and Philippine) – are “Bénard-valid” (i.e., can be explained by our Bénard model). Five of the remaining six observational poles must be relocated within their error-ellipses to become Bénard-valid; the Eurasia pole alone appears to be in error by ～115°, and may actually lie near 40°N, 154°E. The collective results strongly suggest Bénard-like mantle convection cells, and that basal shear tractions are the primary factor in determining the directions of the plate motions as given by the Euler poles. The magnitudes of the computed Euler vectors show, however, that basal shear cannot be the exclusive driving force of plate tectonics, and suggest force contributions (of comparable magnitude for perhaps half of the plates) from the lithosphere itself, specifically subducting slab-pull and (continental) collision drag, which are provisionally evaluated. The relationship of the putative mantle Bénard polygons to dynamic chaos and turbulent flow is discussed.