Experiments on numerical modeling of intense convection

Results of numerical simulations using the WRF-ARW nonhydrostatic model are presented for eight episodes of intense convection over European Russia in the summer of 2007. The calculations were performed on four nested grids with horizontal grid meshes of 27, 9, 3, and 1 km. Convection was parametrized on the first two grids and explicitly resolved on the other two. It has been found that simulations on finer grids with explicit calculation of convective flows make it possible to reproduce heavy rainfalls and strong-wind zones in the areas of intense convection. A preliminary verification of the short-range predictions of convective systems shows that the maximum 12-h precipitation totals and the maximum winds at 10 m are close, in the order of magnitude, to the observed values. Prediction of convection centers is the weakest point. Difficulties in the model verification associated with the absence of data with high space-time resolution are discussed.     

Spectroscopic binaries near the North Galactic Pole paper 6: BD 33° 2206

Photoelectric radial-velocity measurements have confirmed O. C. Wilson’s finding that BD 33° 2206, the secondary star in the wide visual binary ADS 8470, is a spectroscopic binary. It has an eccentric orbit with a period of 100 days. Its γ-velocity is close to the constant radial velocity of the visual primary, confirming the physical association of the stars.     

Synchrotron spectra of short-period pulsars

A model with synchrotron radiation near the light cylinder is proposed to explain the observed spectra of short-period pulsars (P≤0.1 s). These spectra can be described if a power-law energy distribution of the emitting electrons with exponent γ=2–8 is assumed. For most pulsars, the peak frequency ν_m is below 10 MHz. The ν_m(γ) dependence is derived, and shows that the peak frequencies for pulsars with spectral indices α<1.5 may fall in the observable range. In particular, ν_m may be ν_m ～ 100 MHz for PSR J0751 + 1807 and PSR J1640 + 2224. The observed radio spectrum of Geminga (PSR J0633 + 1746) can be described by a synchrotron model with a monoenergetic or Maxwellian distribution of relativistic electrons and a small angle β between the spin axis and magnetic moment (β ～ 10°).     

Regularized derivatives of potential fields and their role in semi-automated interpretation methods

Evaluation of higher derivatives (gradients) of potential fields plays an important role in geophysical interpretation (qualitative and/or quantitative), as has been demonstrated in many approaches and methods. On the other hand, numerical evaluation of higher derivatives is an unstable process – it has the tendency to enlarge the noise content in the original data (to degrade the signal-to-noise ratio). One way to stabilize higher derivative evaluation is the utilization of the Tikhonov regularization. In the submitted contribution we present the derivation of the regularized derivative filter in the Fourier domain as a minimization task by means of using the classical calculus of variations. A very important part of the presented approach is the selection of the optimum regularization parameter – we are using the analysis of the C-norm function (constructed from the difference between two adjacent solutions, obtained for different values of regularization parameter). We show the influence of regularized derivatives on the properties of the classical 3D Euler deconvolution algorithm and apply it to high-sensitivity magnetometry data obtained from an unexploded ordnance detection survey. The solution obtained with regularized derivatives gives better focused depth-estimates, which are closer to the real position of sources (verified by excavation of unexploded projectiles).     

Analysis of the positive ionospheric response to a moderate geomagnetic storm using a global numerical model

Current theories of F-layer storms are discussed using numerical simulations with the Upper Atmosphere Model, a global self-consistent, time dependent numerical model of the thermosphere-ionosphere-plasmasphere-magnetosphere system including electrodynamical coupling effects. A case study of a moderate geomagnetic storm at low solar activity during the northern winter solstice exemplifies the complex storm phenomena. The study focuses on positive ionospheric storm effects in relation to thermospheric disturbances in general and thermospheric composition changes in particular. It investigates the dynamical effects of both neutral meridional winds and electric fields caused by the disturbance dynamo effect. The penetration of short-time electric fields of magnetospheric origin during storm intensification phases is shown for the first time in this model study. Comparisons of the calculated thermospheric composition changes with satellite observations of AE-C and ESRO-4 during storm time show a good agreement. The empirical MSISE90 model, however, is less consistent with the simulations. It does not show the equatorward propagation of the disturbances and predicts that they have a gentler latitudinal gradient. Both theoretical and experimental data reveal that although the ratio of [O]/[N₂] at high latitudes decreases significantly during the magnetic storm compared with the quiet time level, at mid to low latitudes it does not increase (at fixed altitudes) above the quiet reference level. Meanwhile, the ionospheric storm is positive there. We conclude that the positive phase of the ionospheric storm is mainly due to uplifting of ionospheric F2-region plasma at mid latitudes and its equatorward movement at low latitudes along geomagnetic field lines caused by large-scale neutral wind circulation and the passage of travelling atmospheric disturbances (TADs). The calculated zonal electric field disturbances also help to create the positive ionospheric disturbances both at middle and low latitudes. Minor contributions arise from the general density enhancement of all constituents during geomagnetic storms, which favours ion production processes above ion losses at fixed height under day-light conditions.     

Assessment of regional air pollution distribution by point cumulative semivariogram method at Erzurum urban center, TURKEY

In the assessment of air quality, regional distribution and dispersion with distance are important, together with the variations of pollutants in time. On this occasion, the point cumulative semi-variogram (PCSV) method is used in order to find simply regional distribution of pollutants of Erzurum urban centre. This method is based simply on the summation of square differences in air pollutant concentrations between different sites. Monthly regional variation maps of Erzurum are constructed by finding radius of influence (for SO₂, from 1000 m to 3500 m and, for TSP, 1000–2000 m) and PCSV scattering diagram data at different levels by using monthly average sulphur dioxide (SO₂) and total suspended particulate (TSP) matter concentrations in 2001–2002 winter season. Consequently, the air pollution distribution of Erzurum is assessed.     

Validation of an annual cycle simulation with a T42-L20 GCM

An annual cycle of an atmospheric general circulation model (AGCM) is presented. The winter and summer zonal averages of the atmospheric fields are compared with an observed climatology. The main features of the observed seasonal means are well reproduced by the model. One of the main discrepancies is that the simulated atmosphere is too cold, particularly in its upper part. Some other discrepancies might be explained by the interannual variability. The AGCM surface fluxes are directly compared to climatological estimates. On the other hand, the calculation of meridional heat transport by the ocean, inferred from the simulated energy budget, can be compared to transport induced from climatologies. The main result of this double comparison is that AGCM fluxes generally are within the range of climatological estimates. The main deficiency of the model is poor partitioning between solar and non-solar heat fluxes in the tropical belt. The meridional heat transport also reveals a significant energy-loss by the Northern Hemisphere ocean north of 45° N. The possible implications of model surface flux deficiencies on coupling with an oceanic model are discussed.This paper was presented at the International Conference on Modelling of Global Climate Change and Variability, held in Hamburg 11–15 September 1989 under the auspices of the Meteorological Institute of the University of Hamburg and the Max Planck Institute for Meteorology. Guest Editor for these papers is Dr. L. Dümenil     

Thermal evolution of an extensional detachment as constrained by organic metamorphic data and thermal modeling: Graz Paleozoic Nappe Complex (Eastern Alps)

Following Early Cretaceous nappe stacking, the Eastern Alps were affected by late-orogenic extension during the Late Cretaceous. In the eastern segment of this range, a Late Cretaceous detachment separates a very low- to low-grade metamorphic cover (Graz Paleozoic Nappe Complex, GPNC) above a low- to high-grade metamorphic basement. Synchronously, the Kainach Gosau Basin (KGB) collapsed and subsided on top of the section.Metamorphism of organic material within this section has been investigated using vitrinite reflectance data and Raman spectra of extracted carbonaceous material. In the southern part of the GPNC, vitrinite reflectance indicates a decrease in organic maturity towards the stratigraphic youngest unit. The remaining part of the GPNC is characterized by an aureole of elevated vitrinite reflectance values and Raman R2 ratios that parallels the margins of the GPNC. Vitrinite reflectance in the KGB shows a steep coalification gradient and increases significantly towards the western basin margin. The observed stratigraphic trend in the southern GPNC is a result of deep Paleozoic to Early Cretaceous burial. This maturity pattern was overprinted along the margins by advective heat and convective fluids during Late Cretaceous to Paleogene exhumation of basement rocks.During shearing, the fault zone was heated up to ca. 500 °C. This overprint is explained by a two-dimensional thermal model with a ramp-flat fault geometry and a slip rate of 1 to 1.5 cm/year during 5 Ma fault movement. The collapse basin above the detachment subsided in a thermal regime which was characterized by relaxing isotherms.