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.     

Dynamic nucleation process of shallow earthquake faulting in a fault zone

Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P ₁) are interrupted by the arrival of the second impulsive phase P₂ enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P₁ phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P₂ phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.     

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.     

Application of AVO Analysis to Seismic Data for Detection of Gas below Methane Hydrate Stability Zone in Nankai Trough Area

Abstract. For the purpose of development of methane hydrate, occurring in the deep marine subsurface, as a resource, the most important issue is to understand the methane hydrate system (generation, migration and accumulation) as well as to delineate the methane hydrate reservoir properties. We have applied the Amplitude Versus Offset (AVO) analysis to the seismic data acquired in the Nankai Trough, offshore Japan, in order to confirm the occurrence of gas just below the methane hydrate-bearing zone, assuming that gas will show a so-called Class-3 AVO response. Knowledge of the amount and occurrence of gas in the sediment below methane hydrate-bearing zone is one of the keys to understand the methane hydrate system.
We have utilized the qualitative analysis of AVO methodology to delineate how gas is located below the BSR, which is thought to be the reflection event from the interface between the methane hydrate-bearing zone and the underlying gas-bearing zone. In the region of MITI Nankai Trough Well PSW-3, we observe two BSRs separated by 25 ms. After AVO modeling using well data, we applied AVO attribute analysis and attribute crossplot analysis to the seismic data. Finally we applied an offset-amplitude analysis to CMP gather data at specific locations to confirm the results of AVO attribute analysis. The AVO analysis shows that there is very little gas located in the underlying sediment below methane hydrate-bearing zone. This result supports the fact that we could not obtain any clear evidence of gas occurrence just below the methane hydrate-bearing zone in the Nankai Trough well drilling.     

Detection of Methane Hydrate-Bearing Zones from Seismic Data

Abstract. The MITI Nankai Trough wells were drilled offshore Japan in the Tokai area in 1999 and 2000. The occurrence of methane hydrate was confirmed by various indicators in the borehole logs and from core data. These findings have a large impact on potential future Japanese energy resources and other related-scientific interests.
We first tried to find the methane hydrate-bearing zones using interval velocities derived from NMO velocity analysis. However, this analysis produced poor resolution. To achieve a more detailed delineation of the gas hydrate- and gas-bearing zones, we executed a seismic impedance inversion calibrated by the logs from two of the MITI Nankai Trough wells. Although these two wells are only about 90 m apart, we were able to produce an impedance section with fine detail by adopting a simple initial model and incorporating physical properties of the methane hydrate-bearing zones. The locations of the methane hydrate-bearing zones are readily apparent in the final section.