A millennial-scale record of Holocene tsunamis on the Kronotskiy Bay coast, Kamchatka, Russia

Deposits from as many as 50 large tsunamis during the last 7000 years are preserved on the Pacific coast of the Kamchatka Peninsula near the mouth of the Zhupanova River, southern Kronotskiy Bay. These deposits are dated and correlated using Holocene marker tephra layers. The combined, preserved record of tsunami deposits and of numerous marker tephras on Kamchatka offers an unprecedented opportunity to study tsunami frequency. For example, from the stratigraphy along southern Kronotskiy Bay, we estimate frequency of large tsunamis (>5 m runup). In the last 3000 years, the minimum frequency is about one large tsunami per 100 years, and the maximum about one large tsunami per 30 years; the latter frequency occurred from about 0 to 1000 A.D. This time interval corresponds to a period of increased seismicity and volcanic activity that appears to be recorded in many places on the Kamchatka Peninsula.     

Evidence for fullerenes in solid bitumen from pillow lavas of Proterozoic age from Mítov (Bohemian Massif, Czech Republic)

Andesitic pillow lavas containing biogenic, solid bitumen (SB) are a constituent of a Neoproterozoic volcanosedimentary sequence (Teplá-Barrandian unit, Bohemian Massif) in the Mítov area of the Czech Republic. A black shale formation that is crosscut by these andesitic basalts is 565 Ma old. Carbon disulfide extracts of two powdered samples of SB contain 0.2 and 0.3 ppm of C₆₀, respectively, as determined by high-pressure liquid chromatography. The peak assignment based on retention time is fully supported by high-resolution electron ionization mass spectrometry (EI-MS). No C₇₀ was detected, nor was C₆₀ found in two other SB samples from this locality. Other investigated carbonaceous samples from Bohemia (coals and anthracites of Upper Paleozoic age and anthraxolite, graphitoids, and graphite of Upper Proterozoic age) did not contain fullerenes at concentrations above the detection limit of 0.01 ppm. The absence of C₆₀ in these samples was confirmed by EI-MS. The proposed mechanism of fullerene formation involves a primary algal phase, generation of a hydrocarbonaceous mixture in the course of thermal evolution of the sedimentary series, and their high-temperature transformation related to the extrusion of basalt. An important feature for fullerene conservation was the enclosure of fullerenes in SB with a structure similar to glasslike carbon, where the fullerene was protected against oxidation.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.     

Trapped Xe and I-Xe ages in aqueously altered CV3 meteorites

Twenty-two dark inclusions (DIs) from Allende (18), Leoville (2), Vigarano (1) and Efremovka (1) were studied by the I-Xe method. All except two of these DIs (Vigarano 2226 and Leoville LV2) produce well-defined isochrons, and precise I-Xe ages. The Allende DIs formed a tight group about 1.6 Ma older than Shallowater (4.566 ± 0.002 Ga), about 5 Ma older than four previously studied Allende CAIs. Most of the dark inclusions require trapped Xe with less ¹²⁹Xe (or more ¹²⁸Xe) than conventional planetary Xe (well restricted in composition by Q-Xe or OC-Xe). Studies of an irradiated/unirradiated DI pair from Allende demonstrate that the ¹²⁸Xe/¹³²Xe ratio in trapped is normal planetary, so that a ¹²⁹Xe/¹³²Xe ratio below planetary seems to be required. Yet, this is not possible given constraints on ¹²⁹Xe evolution in the early solar system. Trends among all of the Allende DIs suggest that an intimate mixture of partially decayed iodine and Xe formed a pseudo trapped Xe component enriched in both ¹²⁹Xe and ¹²⁷I, and subsequently in ¹²⁸Xe after n-capture during reactor irradiation. Enrichment in radiogenic ¹²⁹Xe, but with a ¹²⁹Xe/¹²⁷I ratio less than that observed in the iodine host phase, places closure of this trapped mixture ≥13 Ma after precipitation of the major iodine-bearing phase. Because the I-Xe isochron is a mixing line between iodine-derived and trapped Xe (pseudo or not), I-Xe ages, given by the slope of this mixing line, are not compromised by the presence of pseudo trapped Xe, and the precision of the I-Xe ages is given by the statistics of the line fit.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.