Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Long-term seismic noise investigations on Shikotan Island: First results

In 2003–2004, long-term seismic noise observations were launched on Shikotan Island (Lesser Kuril Range) based on the “Shikotan” dormant regional seismic station. The geological and geophysical data on the registration area are reported. Information about the equipment and its technical specifications is given. The precursors to the strongest local earthquakesthat occurred in the Shikotan Island region in January 2005–March 2007 are identified.     

Permeability development in vesiculating magmas: implications for fragmentation

 Fragmentation, or the "coming apart" of magma during a plinian eruption, remains one of the least understood processes in volcanology, although assumptions about the timing and mechanisms of fragmentation are key parameters in all existing eruption models. Despite evidence to the contrary, most models assume that fragmentation occurs at a critical vesicularity (volume percent vesicles) of 75–83%. We propose instead that the degree to which magma is fragmented is determined by factors controlling bubble coalescence: magma viscosity, temperature, bubble size distribution, bubble shapes, and time. Bubble coalescence in vesiculating magmas creates permeability which serves to connect the dispersed gas phase. When sufficiently developed, permeability allows subsequent exsolved and expanded gas to escape, thus preserving a sufficiently interconnected region of vesicular magma as a pumice clast, rather than fully fragmenting it to ash. For this reason pumice is likely to preserve information about (a) how permeability develops and (b) the critical permeability needed to insure clast preservation. We present measurements and calculations that constrain the conditions (vesicularity, bubble size distribution, time, pressure difference, viscosity) necessary for adequate permeability to develop. We suggest that magma fragments explosively to ash when and where, in a heterogeneously vesiculating magma, these conditions are not met. Both the development of permeability by bubble wall thinning and rupture and the loss of gas through a permeable network of bubbles require time, consistent with the observation that degree of fragmentation (i.e., amount of ash) increases with increasing eruption rate. Received: 5 July 1995 / Accepted: 27 December 1995     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

LREE distribution in perovskite,apatite and titanite from South West Ugandan xenoliths and kamafugite lavas

Summary In-situ microprobe LREE analyses of perovskite and titanite (La, Ce, Nd), and apatite (La, Ce), from SW Ugandan clinopyroxenite xenoliths and kamafugite lavas indicate that LREE distribution in these minerals is determined by a number of factors related to their different parageneses: In particular LREE content is affected by whether the LREE-bearing minerals have crystallised from metasomatic carbonate or from silicate (i.e. metasomatic or magmatic) melts in the mantle. In this situation LREE partition favours carbonate over silicate melts. Distribution of LREE in perovskite and apatite crystallised from magmatic mantle melts or mantle-derived lavas is chiefly determined by preference of LREE for perovskite > apatite > titanite. LREE zoning in perovskite is influenced by changes in melt structure: increasing melt polymerisation enhancing mineral_LREE/melt_LREE partition into perovskite rims in magmatic xenoliths; decreasing melt polymerisation depleting LREE in lava perovskite rims. This zoning is reinforced by perovskite competition with apatite for LREE: perovskite (cores/rims) co-crystallising with apatite is reduced in LREE. There are 37 instances of perovskitewith Ce below detection while La and Nd levels are normal. These occur in both xenoliths and lavas; in grain zones or whole grains. Likewise Ce alone of the LREE is below detection in six out of ten titanite analyses. These observations are interpreted as evidence for increased f_{O 2}, Ce⁴ + being excluded from these mineral structures. Recognition of these various processes can elucidate the interpretation of bulk rock and bulk mineral LREE signatures in kamafugite volcanism.

---

LREE Verteilung in Perovskit, Apatit und Titanit aus Xenolithen und kamafugitischen Laven Südwest-Ugandas

Zusammenfassung In-situ LREE Analysen von Perovskit und Titanit (La, Ce, Nd) und Apatit (La, Ce) aus Klinopyroxenit-Xenolithen und kamafugitischen Laven Südwest-Ugandas zeigen, daß die LREE Verteilung in diesen Mineralen durch eine Vielzahl von Faktoren, die mit Unterschieden in den Paragenesen zusammenhängen, bestimmt wird: Der LREE-Gehalt wird im besonderen davon bestimmt, ob die LREE-führenden Minerale aus metasomatischen Karbonat- oder aus (metasomatischen oder magmatischen) Silikatschmelzen im Mantel auskristallisierten. Dabei erfolgt die LREE Fraktionierung zu Gunsten der Karbonatschmelzen. Die LREE-Verteilung von Perovskit und Apatit, die aus magmatischen Mantelschmelzen oder -laven kristallisierten, wird vorrangig durch den bevorzugten Einbau der LREE in Perovskit > Apatit > Titanit kontrolliert. Der LREE Zonarbau von Perovskit wird durch die Änderungen der Schmelzstruktur beinflußt: Verstärkte Schmelzpolymerisation führt zu verstärkter Mineral_LFEE/Schmelze_LREE Fraktionierung in den Perovskiträndern magmatischer Xenolithe, eine Abnahme der Schmelzpolymerisation hingegen resultiert in einer Abreicherung der LREE in den Perovskiträndern. Diese Art der Zonierung wird durch den Wettbewerb von Perovskit mit Apatit um die LREE verstärkt. Perovskit (Kerne/Ränder), der mit Apatit gemeinsam auskristallisierte, ist ärmer an LREE. 37 Fälle, in denenCe nicht nachweisbar war, La und Nd aber in normaler Konzentration auftreten, wurden sowohl in den Xenolithen als auch in den Laven gefunden; und zwar entweder in Kornbereichen oder in ganzen Körnern. Vergleichsweise liegt Ce nur in sechs von zehn Titanitproben unterhalb der Nachweisgrenze. Diese Beobachtungen werden als Hinweise auf erhöhte SauerstoffFugazitäten, bei denen Ce^4– aus der Mineralstruktur ausgeschlossen wird, angesehen.Ein Verständnis dieser verschiedenen Prozesse kann zur besseren Interpretation von LREE Gesamtgesteins- und Gesamtmineral-Signaturen in Kamafugiten beitragen.

---

---

With 3 Figures     

Detection of a new interstellar molecular ion, H2COH+ (protonated formaldehyde)

A new interstellar molecular ion, H2COH+ (protonated formaldehyde), has been detected toward Sgr B2, Orion KL, W51, and possibly in NGC 7538 and DR21(OH). Six transitions were detected in Sgr B2(M). The 1(1,0)-1(0,1) transition was detected in all sources listed above. Searches were also made toward the cold, dark clouds TMC-1 and L134N, Orion (3N, 1E), and a red giant, IRC + 10216, without success. The excitation temperatures of H2COH+ are calculated to be 60-110 K, and the column densities are on the order of 10(12)-10(14) cm-2 in Sgr B2, Orion KL, and W51. The fractional abundance of H2COH+ is on the order of 10(-11) to 10-(9), and the ratio of H2COH+ to H2CO is in the range 0.001-0.5 in these objects. The values in Orion KL seem to be consistent with the "early time" values of recent model calculations by Lee, Bettens, & Herbst, but they appear to be higher than the model values in Sgr B2 and W51 even if we take the large uncertainties of column densities of H2CO into account. We suggest production routes starting from CH3OH may play an important role in the formation of H2COH+.     

A window to the operation of microplate tectonics in the Tethys Ocean: the geochemistry of the Samothrace granite,Aegean Sea

Summary The island of Samothrace, northeastern Aegean Sea, consists of five main geological units: (i) A basement unit consisting of low grade metamorphic rocks (metapelites, marbles, metavolcanic rocks, and a metaconglomerate); (ii) an ophiolitic complex with K-Ar hornblende date of 154 ± 7 and 155 ± 7 Ma; (iii) A granite intrusion with biotite K-Ar dates of 14.5 ± 0.3 and 14.5 ± 0.5 Ma, and a contact metamorphic event dated at 40.9 + 2.2 Ma; (iv) a unit of Cenozoic volcanic rocks: orogenic volcanism apparently occurred in two cycles with Upper Eocene tholeiitic to calc-alkaline volcanic rocks and post-Eocene high-K andesites to trachytes. (v) Quaternary clastic sedimentary rocks which occur around the peripheral parts of the island. The granitic intrusion is predominantly a hornblende-biotite granite, granodiorite or quartz monzonite, with porphyritic variants and mafic enclaves. The pluton is cut by granophyre, aplite and rare granodioritic veins. All lithological units of the Samothrace intrusion show smooth and continuous major element trends and similar chondrite- and Ocean Ridge Granite-normalized incompatible element profiles. ORG-normalized incompatible element contents of Hf, Zr, Sm are explained with fractionation close to the normalizing values Y and Yb contents combined with high K/Yb ratios; Rb and Th are significantly enriched relative to Nb and Ta. In Y-Nb and Rb-SiO₂ space most samples of the Samothrace granite, plot in the volcanic arc and the syn-collisional granite fields. In Y + Nb-Rb space they are equally distributed within and transgress these two domains. The geochemical and regional data suggest a subduction or collision environment but biotite mineral data do not support a collisional setting for magma genesis. The Samothrace granite was probaby associated with a post-collisional domain after the closure of the Axios section of the Tethys Ocean.

---

Ein Einblick in das Wirken von Mikroplattentektonik in der Tethys—Die Geochemie des Samothrake Granites, Agäisches Meer

Zusammenfassung Die Insel Samothrake in der nordöstlichen Ägäis besteht aus fünf geologischen Haupteinheiten: (i) einem schwach metamorphen Basement (Metapelite, Marmore, Metavulkanite und Metakonglomerate); (ii) einem ophiolithischem Komplex, der mit K-Ar Datierungen an Hornblende ein Alter von 154 ± 7 und 155 ± 7 Ma ergab; (iii) ein granitischer Intrusionskörper mit K-Ar Altern an Biotit von 14.0 ± 0.3 und 14.5 ± 0.5 Ma und einem kontaktmetamorphen Ereignis, das mit 40.9 ± 2.2 Ma datiert ist; (iv) eine Abfolge känozoischer Vulkanite, wobei der orogene Vulkanismus offensichtlich in zwei Zyklen ablief mit tholeiitischen bis kalkalkalischen Vulkaniten im oberen Eozän und high-K Andesiten bis Trachyten im post-Eozän; (v) quartären klastischen Sedimentgesteinen, die im Randbereich der Insel auftreten. Die Granitintrusion setzt sich hauptsächlich aus Hornblende-Biotitgraniten, Granodioriten oder Quarzmonzoniten mit teilweise porphyrischen und mafischen Enklaven enthaltenden Varietäten zusammen. Der Pluton wird von Granophyren, Apliten und seltener von granodioritischen Gängen durchschlagen. Alle lithologischen Einheiten der Samothrake Intrusion zeigen kontinuierliche Hauptelementtrends und ähnliche Chondrit und ORG-normalisierte inkompatible Elementprofile. Die Gehalte an den inkompatiblen Elementen Hf, Zr, Sm sind sehr ähnlich denen von ozeanischen Graniten (ORG). Die niedrigen Y und Yb-Gehalte und die hohen K/Yb Verhältnisse werden durch Fraktionierung erklärt. Rb und Th sind signifikant angereichert im Vergleich zu Nb und Ta. In Y-Nb und Rb-SiO₂ Diagrammen plotten die meisten Proben des Samothrake Granites im Feld der vulkanischen Inselbogen- und Synkollisionsgranite. Im Y + Nb-Rb Diagramm zeigt sich eine gleichmäßige und überlappende Verteilung. Die geochemischen und regionalen Daten weisen auf einen Subduktions- oder Kollisionsbereich hin, obwohl die Biotitzusammensetzungen nicht für eine Bildung der Magmen in einem Kollisionsbereich sprechen. Die Bildung des Samothrake Granites steht möglicherweise mit post-Kollisionstektonik nach dem Schließen der Axioszone in der Tethys in Zusammenhang.

---

---

With 7 Figures     

Chernobyl radiocaesium in a karst system,Marble Cave,Crimea

 Surface contamination with radioactive caesium introduced into the environment after the accident at the Chernobyl nuclear plant was high enough in the Crimean Mountains to allow using radiocaesium as an indicator of penetration of radioactive contamination into a karst system. Caesium concentrations have been studied in soils above Marble Cave, Tchatyrdag Plateau, in percolation waters and in sediments transported by percolation waters within the cave. Contamination of the daylight surface with ¹³⁷Cs is about 30 kqB m^–2 which is approximately 13 times higher than the density of global fallouts. Besides ¹³⁷Cs, almost all samples showed presence of ¹³⁴Cs with the ¹³⁷Cs/¹³⁴Cs ratio close to that of Chernobyl contaminations. Concentrations of ¹³⁷Cs range from 9 to 15 mBq l^–1 in the present percolation waters in the cave. In sediments related to percolation waters ¹³⁴Cs is detected in some samples besides ¹³⁷Cs, although the effect of ²²⁸Ac is not ruled out. Surprisingly, the highest concentrations of radiocaesium were measured in "old" sediments in the cave's lower series. These sediments are not associated with modern percolation and are represented by a clay/moonmilk alternating sequence deposited in an old dried cave lake. Moonmilk layers have higher caesium content than clay. It is assumed that Chernobyl caesium was transported into the cave with aerosols which were then deposited mainly in areas where condensation occurs. The sampling site is located just in the boundary between two microclimatic zones with a temperature gradient of 0.5  °C. Active condensation processes occur in this area. Caesium was not detected in another similar sampling site (old lake deposits) located within homogeneous microclimatic conditions. If the above interpretation is correct, these results show the geochemical significance of the aerosol-condensation mechanism of mass transport and localisation. Received: 1 June 1995 · Accepted: 4 December 1995     

The population of Croatia

Croatia is a Central European and Mediterranean state, located in contact with and under the influence of various spheres of civilization, which for centuries have penetrated and conflicted on and around the territory of Croatia. Such a position has resulted in various influences which did not always have a positive effect on the development of Croatia. The gradual narrowing and decrease of the Croatian ethnic territory, as well as the presence of national minorities in it, was also a result of the aforementioned position and outside influences (ulji 1993/1994).The demographic structure of Croatia indicates a series of specificities which were primarily conditioned by the historical development of Croatia and which is particularly expressed in constant emigration since the end of the nineteenth century, the relatively large direct and indirect losses to the population during and immediately after the First and Second World Wars, emigration as a type of population movement in all inter-census periods after 1945, the appearance of a natural decline and the aging of the population on almost one half of the state territory.