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.     

河北平原水系密度与隐伏活动构造的关系

以1926年出版的比例尺为1：5万的顺直地形图作为基础图件，用地理信息系统技术对河北平原作水系密度计量分析，揭示隐伏在第四纪巨厚的沉积物之下活动构造的信息，研究结果证实水系线密度统计图与使用其它地球科学方法获得的成果，例如基底活动断块图，地貌类型图和第四纪厚度图都有很好的对应，特别与布格重力图对应更好，研究结果发现水系密度图与布格重力异常图在总体一致性较好的前提下，多处反向异常区均为近代强震发生的区域，提出这种地表水系密度统计研究的成果不仅可以揭示深部隐伏的活动构造，而且具有较明确的时间含义。     

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.     

Extracting flow structures from tracer data

A method of visualizing structures in closed chaotic flows out of homogenous particle distributions is presented in the example of models of a meandering jet. To this end, the system will be leaked or opened up by defining a region of the flow, so that a particle is considered to be escaped if it leaves this region. By applying this method to an ensemble of nonescaped tracers, we are able to characterize mixing processes by visualizing the converging and stretching filamentations (stable and unstable manifolds) in the flow without using additional mathematical tools. The possibility of applying the algorithm to analyze buoy data, and a comparison with the finite time manifolds are discussed.     

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.     

La vallée d'Ys : un paléoréseau hydrographique immergé en baie de Douarnenez (Finistère,France)

Interpretation of the recent high-resolution survey, CANADOU 2000, in the Bay of Douarnenez (Finistère, France) allowed us to restore the morphology of the substratum and the sedimentary filling of the bay. The Brioverian and Palaeozoic substratum reveals a well-defined network of incised valleys as results of successive emergence stages of the Bay during the Quaternary. Valleys join in a westward-widened mean valley, called Ys Valley. The present-day sedimentary fill of the bay of Douarnenez appears mainly controlled by the Holocene rise and the consecutive highstand. It comprises fluvial and estuarine deposits filling up incised valleys and marine sedimentation extending out of the incised valleys. To cite this article: G. Jouet et al., C. R. Geoscience 335 (2003).To cite this article: G. Jouet et al., C. R. Geoscience 335 (2003).     

Magnetic activity inside the auroral zones and field-aligned currents

Big perturbations of the magnetic field (amplitudes larger than 250 nT) are simply detected by subtracting the values of a model from the measurements of CHAMP satellite. Taking a full year of CHAMP data and organizing them in four subsets of three months length (spring, summer, autumn, winter), it is found that: (a) the two domains where such big perturbations mainly exist are limited, in both hemispheres, by a parallel of high latitude of the corrected geomagnetic coordinates system; (b) a conspicuous seasonal (annual) variation affects the density of the perturbations and is opposite in the two hemispheres. We hold that these perturbations are linked to the midday magnetic activity within the auroral zone, long ago described by one of us (Mayaud, 1956). The source of the perturbations observed at the satellite altitude would be field-aligned currents resulting from the penetration of the solar wind into the magnetospheric cusps. To cite this article: J.-L. Le Mouël et al., C. R. Geoscience 335 (2003).     

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