Mantle Convection Modeling with Viscoelastic/Brittle Lithosphere: Numerical Methodology and Plate Tectonic Modeling

--The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.     

Geo–ecology and management of sensitive montane landscapes

Montane (alpine) areas are generally of high value for nature conservation. Such environments and the habitats they support are dynamic and often fragile. They are vulnerable to disturbance from a range of human activities and are responsive to climate changes over short and long timescales. Biodiversity and conservation values are closely linked to geological history, geomorphological processes and soils, and it is crucial that management systems are based on understanding these links.
There are many similarities between the Cairngorm Mountains (Scotland), the Giant Mountains (Czech Republic) and Abisko Mountains (Sweden) in terms of geology, geomorphology, ecology, links with biodiversity and high conservation importance. Comparable pressures and management issues involve, to varying degrees, a history of human use and impacts from deforestation, pasturing, grazing, recreation and atmospheric pollution. Landscape change therefore involves a complex interplay between natural and anthropogenic factors. Managing such change requires better understanding of the geo–ecological processes involved and the factors that determine landscape sensitivity. This is illustrated through a simple framework and examples from the three areas. Comparison of landscape sensitivity between similar montane areas, but in different geographic locations and climatic environments, should allow more informed management planning and a precautionary approach in advance of further changes in human activity and from predicted global warming scenarios.     

The Role of Organic Matter During Metal Enrichment in Permian Kupferschiefer of the Rudna Mine,Southwest Poland

The bulk composition of organic matter and saturated and aromatic hydrocarbons extracted from 16 samples collected from two Kuperschiefer profiles in the Rudna mine,Southwest Poland has been analyzed to study the role of organic matter during base metal enrichment in the Kupferschiefer shale.The results indicated that the extract yields and saturated hydrocarbon yields decreased with increasing base metal contents.GC and GC/MS analyses indicated that n -alkanes and alkylated aromatic compounds were depleted and may have served as hydrogen donators for thermochemical sulfate reduction.The enrichment of base metal is closely connected with the destruction of hydrocarbons.     

Optimal Model for Geoid Determination from Airborne Gravity

Two different approaches for transformation of airborne gravity disturbances, derived from gravity observations at low-elevation flying platforms, into geoidal undulations are formulated, tested and discussed in this contribution. Their mathematical models are based on Green's integral equations. They are in these two approaches defined at two different levels and also applied in a mutually reversed order. While one of these approaches corresponds to the classical method commonly applied in processing of ground gravity data, the other approach represents a new method for processing of gravity data in geoid determination that is unique to airborne gravimetry. Although theoretically equivalent in the continuous sense, both approaches are tested numerically for possible numerical advantages, especially due to the inverse of discretized Fredholm's integral equation of the first kind applied on different data. High-frequency synthetic gravity data burdened by the 2-mGal random noise, that are expected from current airborne gravity systems, are used for numerical testing. The results show that both approaches can deliver for the given data a comparable cm-level accuracy of the geoidal undulations. The new approach has, however, significantly higher computational efficiency. It would be thus recommended for real life geoid computations. Additional errors related to regularization of gravity data and the geoid, and to accuracy of the reference field, that would further deteriorate the quality of estimated geoidal undulations, are not considered in this study.     

An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.     

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.     

Geoid determination using adapted reference field, seismic Moho depths and variable density contrast

 The traditional remove-restore technique for geoid computation suffers from two main drawbacks. The first is the assumption of an isostatic hypothesis to compute the compensation masses. The second is the double consideration of the effect of the topographic–isostatic masses within the data window through removing the reference field and the terrain reduction process. To overcome the first disadvantage, the seismic Moho depths, representing, more or less, the actual compensating masses, have been used with variable density anomalies computed by employing the topographic–isostatic mass balance principle. In order to avoid the double consideration of the effect of the topographic–isostatic masses within the data window, the effect of these masses for the used fixed data window, in terms of potential coefficients, has been subtracted from the reference field, yielding an adapted reference field. This adapted reference field has been used for the remove–restore technique. The necessary harmonic analysis of the topographic–isostatic potential using seismic Moho depths with variable density anomalies is given. A wide comparison among geoids computed by the adapted reference field with both the Airy–Heiskanen isostatic model and seismic Moho depths with variable density anomaly and a geoid computed by the traditional remove–restore technique is made. The results show that using seismic Moho depths with variable density anomaly along with the adapted reference field gives the best relative geoid accuracy compared to the GPS/levelling geoid. Received: 3 October 2001 / Accepted: 20 September 2002 Correspondence to: H.A. Abd-Elmotaal