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.     

Meso-gamma scale forecasts using the nonhydrostatic model LM

Summary ?The nonhydrostatic model LM was developed for small scale operational predictions. Advances in computer development will give the possibility of operational models of a rather fine scale, which will cover the meso-gamma scale. The LM is currently applied at a scale of 7 km and an increase of the operational resolution to 2.5 km is planned for the next few years. Predictions of such high resolution require to abandon the hydrostatic assumption, which is used with most current operational weather prediction models. The LM was designed to cover all resolutions from 50 m to 50 km with an efficiency making it suitable for operational use. It is a fully elastic model, using second order centred finite differences. The time integration is done using the Klemp–Wilhelmson method, treating the slow modes by a larger time step than the fast modes. The vertical propagation of the fast waves is done implicitly. After describing the design of the LM, this paper gives examples of model predictions at the meso-γ scale. Some results of the current operational application at the resolution 7 km are presented. Deficiencies in the localisation of model generated precipitation are investigated using an idealised bell shaped mountain and applying different resolutions. In this way the convergence to the correct solution can be investigated. From these results it is concluded, that orographic filtering is necessary and the effect of such filtering on precipitation forecasts is investigated. Finally, the prediction of a squall line over northern Germany is shown in order to demonstrate the potential of the model in forecasting the meso-γ scale. Received May 15, 2001; revised September 21, 2001     

Accretion and exhumation at a Variscan active margin, recorded in the Saxothuringian flysch

The Saxothuringian flysch basin, on the north flank of the Central European Variscides, was fed and eventually overthrust by the northwestern, active margin of the Tepla-Barrandian terrane. Clast spectra, mineral composition and isotopic ages of detrital mica and zircon have been analyzed in order to constrain accretion and exhumation of rocks in the orogenic wedge. The earliest clastic sediments preserved are of early Famennian age (ca. 370?Ma). They are exposed immediately to the NW of the suture, and belong to the par-autochthon of the foreland. Besides ultramafic (?ophiolite) material, these rocks contain clasts derived from Early Paleozoic continental slope sediments, originally deposited at the NW margin of the Saxothuringian basin. These findings, together with the paleogeographic position of the Famennian clastics debris on the northwestern passive margin, indicate that the Saxothuringian narrow ocean had been closed by that time. Microprobe analyses of detrital hornblendes suggest derivation from the “Randamphibolit” unit, now present in the middle part of the Saxothuringian allochthon (Münchberg nappes). Detrital zircons of metamorphic rocks formed a little earlier (ca. 380?Ma) indicate rapid recycling at the tectonic front. The middle part of the flysch sequence (ca. early to middle Viséan), both in the par-autochthon and in the allochthon, contains abundant clasts of Paleozoic rocks derived from the northwestern slope and rise, together with debris of Cadomian basement, 500-Ma granitoids and 380?Ma (early Variscan) crystalline rocks. All of these source rocks were still available in the youngest part of the flysch (c. middle to late Viséan), but some clasts record, in addition, accretion of the northwestern shelf. Our findings permit deduction of minimum rates of tectonic shortening well in excess of 10–30?mm per year, and rates of exhumation of ca. 3?mm/a, and possibly more.     

Measurement and research method of metallogenic depths of Jiaojia orefield and Linglong orefield Shandong Province, China

Taking the contribution of the tectonic force to the total hydrostatic pressure into account, the author puts forward a new method on the calculation of the depth of petrogenesis and metallogenesis which is summarized as follows: first the tectonic added hydrostatic pressured P_s is subtracted from the total hydrostatic pressure P, then using their difference P_r, according to the general method the depth of petrogenesis and metallogenesis can be determined in consideration of lateral constraint. By the new method the following data on the depths of the metallogenesis are obtained: 2 243.6 m (No. I and No. II veins with metallogenic epoch of 105 Ma) and 1 632. 38 m (No. III vein with 105 Ma) for Jiaojia orefield, and 3 454.97 m (NE-trending zone with 213.2 Ma), 1 902.79 III (ENE-trending zone with 100.28 Ma), 1 090.97 m (NE-trending zone with 80.67 Ma) and 720.55 m (NNE-trending zone with 71.86 Ma) for Linglong orefield. Project supported by the Foundation of the State Planning Commission, China and the Foundation of the State Science and Technology Commission, China.