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.     

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.     

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.     

Environmental impact assessment in the Moradabad industrial area (rivers Ramganga-Ganga interfluve), Ganga Plain,India

Sediment samples collected in the Moradabad area, lying in the interfluve of the Ganga and Ramganga Rivers, were analysed for heavy metals, after studying the geomorphology of the area. Geomorphologically, the area can be divided into three terraces - the T₀, T₁ and T₂ surfaces. The rivers on these three surfaces show varying amounts of pollution depending upon the input from industries and urban settlements. The Ramganga River on the T₀ surface shows the highest amount of pollution. However, the pollution levels in all these rivers show a downstream dilution effect. The characteristic feature of the vast interfluve area (T₂ surface) is the presence of several, independent basins which are closed and rarely interact with each other or with any river. The sediments are redistributed and redeposited within the basin itself, and thus these basins serve as sinks. The sediments of one such basin in the study area show significant concentrations of arsenic, chromium, copper, nickel, lead, zinc and organic carbon. The concentrations of heavy metals in such a basin will show exponential increases with time, because there is no activity to funnel out the sediments and dilute the effect of pollution. This increase will pose more threats, as ultimately it will make its way laterally and vertically through the sediments, thereby polluting groundwater.     

Strontium isotope evolution of Late Permian and Triassic seawater

The ⁸⁷Sr/⁸⁶Sr values based on brachiopods and conodonts define a nearly continuous record for the Late Permian and Triassic intervals. Minor gaps in measurements exist only for the uppermost Brahmanian, lower part of the Upper Olenekian, and Middle Norian, and only sparse data are available for the Late Permian. These 219 measurements include 67 brachiopods and 114 conodont samples from the Tethyan realm as well as 37 brachiopods and one conodont sample from the mid-European Middle Triassic Muschelkalk Sea. The Late Permian/Lower Triassic interval is characterized by a steep 1.3 × 10⁻³ rise, from 0.7070 at the base of the Dzhulfian to 0.7082 in the late Olenekian, a rate of change comparable to that in the Cenozoic. In the mid-Triassic (Anisian and Ladinian), the isotope values fall to 0.7075, followed again by a rise to 0.7081 in the Middle/Late Norian. The ⁸⁷Sr/⁸⁶Sr values decline again in the Late Norian (Sevatian) and Rhaetian to 0.7076.The sharp rise in the ⁸⁷Sr/⁸⁶Sr values during the Late Permian/Early Triassic was coincident with widespread clastic sedimentation. Because of the paucity of tectonic uplifts, the enhanced erosion may have been due to intermittent humid phases, during mainly an arid interval, coupled with the absence of a dense protective land plant cover following the mass extinction during the latest Permian. The apex of the ⁸⁷Sr/⁸⁶Sr curve at the Olenekian/Anisian boundary coincides with cessation of the large-scale clastic sedimentation and also marks the final recovery of land vegetation, as indicated by the renewed onset of coal formation in the Middle Triassic. The rising ⁸⁷Sr/⁸⁶Sr values from the Middle Carnian to the Late Norian coincide with the uplift and erosion of the Cimmeride-Indosinian orogens marking the closure of the Palaeotethys. The subsequent Rhaetian decline that continues into Jurassic (Pliensbachian/Toarcian boundary), on the other hand, coincides with the opening of the Vardar Ocean and its eastern continuation in the Izmir-Ankara Ophiolitic Belt.Samples from the Upper Muschelkalk are more radiogenic than the global trend. This may reflect separation of the basin from the open ocean. Due to strong meteoric influx from a large land mass in the north, the Germanic Basin became increasing brackish up section in the north and east, but because of the high evaporation rates, the salt content was not much reduced in the southern and central basin where a rich, but increasingly endemic, marine fauna survived.     

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     

Intercomparison of Stratospheric Chemistry Models under Polar Vortex Conditions

Several stratospheric chemistry modules from box, 2-D or 3-D models, have been intercompared. The intercomparison was focused on the ozone loss and associated reactive species under the conditions found in the cold, wintertime Arctic and Antarctic vortices. Comparisons of both gas phase and heterogeneous chemistry modules show excellent agreement between the models under constrained conditions for photolysis and the microphysics of polar stratospheric clouds. While the mean integral ozone loss ranges from 4–80% for different 30–50 days long air parcel trajectories, the mean scatter of model results around these values is only about ±1.5%. In a case study, where the models employed their standard photolysis and microphysical schemes, the variation around the mean percentage ozone loss increases to about ±7%. This increased scatter of model results is mainly due to the different treatment of the PSC microphysics and heterogeneous chemistry in the models, whereby the most unrealistic assumptions about PSC processes consequently lead to the least representative ozone chemistry. Furthermore, for this case study the model results for the ozone mixing ratios at different altitudes were compared with a measured ozone profile to investigate the extent to which models reproduce the stratospheric ozone losses. It was found that mainly in the height range of strong ozone depletion all models underestimate the ozone loss by about a factor of two. This finding corroborates earlier studies and implies a general deficiency in our understanding of the stratospheric ozone loss chemistry rather than a specific problem related to a particular model simulation.     

A Case Study Of An On-Ice Air Flow Over The Arctic Marginal Sea-Ice Zone

A case study of warm air advection over the Arctic marginalsea-ice zone is presented, based on aircraft observations with direct flux measurements carriedout in early spring, 1998. A shallow atmospheric boundary layer (ABL) was observed, which wasgradually cooling with distance downwind of the ice edge. This process was mainly connected with astrong stable stratification and downward turbulent heat fluxes of about 10–20 W m^-2, but wasalso due to radiative cooling. Two mesoscale models, one hydrostatic and the other non-hydrostatic,having different turbulence closures, were applied. Despite these fundamental differences betweenthe models, the results of both agreed well with the observed data. Various closure assumptions had amore crucial influence on the results than the differences between the models.Such an assumption was, for example,the parameterization of the surface roughness for momentum (z₀) and heat (z_T). This stronglyaffected the wind and temperature fields not only close to the surface but also within and abovethe temperature inversion layer. The best results were achieved using a formulation for z₀ that took intoaccount the form drag effect of sea-ice ridges together withz_T = 0.1z₀. The stability within theelevated inversion strongly depended on the minimum eddy diffusivity K_min. A simple ad hocparameterization seems applicable, where K_min is calculated as 0.005 timesthe neutral eddy diffusivity. Although the longwave radiative cooling was largest within the ABL, theapplication of a radiation scheme was less important there than above the ABL. This was related to theinteraction of the turbulent and radiative fluxes. To reproduce the strong inversion, it wasnecessary to use vertical and horizontal resolutions higher than those applied in most regional andlarge-scale atmospheric models.