Water quality model of the Danube in its Bulgarian part

Models for predicting river water quality could assist protection and utilization of water resources under alternative management strategies. The model of the water quality characteristics of the Danube is based on an ecosystem approach. Considered are all waste water discharges and inflow of tributaries among the catchment area. The identification of model parameters is realized according to the method of integral transformation. For the performance of the model a computer program was elaborated. It is applied in predicting the water quality of the Danube considering almost 20 indicators for the abiotic and biotic parameters of river ecosystem.     

Hyperboloidal precession of a dynamically symmetric satellite. Construction of normal forms of the Hamiltonian

The Norma specialized program package, intended for normalization of autonomous Hamiltonian systems by means of computer algebra, is used in studies of small-amplitude periodic motions in the neighbourhood of regular precessions of a dynamically symmetric satellite on a circular orbit. The case of hyperboloidal precession is considered. Analytical expressions for normal forms and generating functions depending on frequencies of the system as on parameters are derived. Possible resonances are considered in particular. The 6th order of normalization is achieved. Though the intermediate analytical expressions occupy megabytes of computer's main memory, final ones are quite compact. Obtained analytical expressions are applied to the analysis of stability of small-amplitude periodic motions in the neighbourhood of hyperboloidal precession.     

On verification of the asymptotic model of the first kind

For the case of closed Friedmann models with the cosmological constant, an N - dependence is plotted, whereN is the conformal time of existence of the Universe, and is the ratio of the cosmological constant to its Einstein's value. The observational properties of the asymptotic model of the first kind (the so-called Al model) corresponding to the maximum of this dependence are analyzed. For the matter-dominated models, the maximum is achieved at the point = 1. Formulae for the calculation of the age of the Universe and for that of the photometric distance in the Al model are deduced. The observational properties of the A1 model are compared to the corresponding properties of the standard cosmological model ( = 0) which does not occupy any special distinguished position on the N - diagram. It is shown that from the standpoint of the modern observational cosmology these two models cannot be told one from another. However, the A1 model has better standing from the viewpoint of the strong wordings of the anthropic cosmological principle.     

Computation of additional components of the first-order ray approximation in isotropic media

This study shows that the use of the first-order additional components of the ray method in the seismic wave field modeling is easy and that it can bring a substantial improvement of the standard ray results obtained with the zero-order ray approximation only. For the calculation of a first-order additional component, spatial derivatives of the parameters of the medium and spatial derivatives of the zero-order ray amplitude term are necessary. The evaluation of the former derivatives is straightforward; the latter derivatives can be calculated approximately from neighboring rays by substituting the derivatives by finite differences. This allows an effective calculation of the first-order additional terms in arbitrary laterally varying layered media.The importance of the first-order additional terms is demonstrated by the study of individual higher-order terms of the ray series representing elementaryP andS elastodynamic Green functions for a homogeneous isotropic medium. The study shows clearly that the consideration of the first-order additional terms leads to a more substantial decrease of the difference between approximate and exact elementary Green functions than any other higher-order term. With this in mind, effects of the first-order additional terms on the ray synthetic seismograms for aVSP configuration are studied. It is shown that the use of the additional terms leads to such phenomena, unknown in the zero-order approximation of the ray method, like quasi-elliptical and transverse polarization of a singleP wave or longitudinal polarization of a singleS wave.     

Some properties of unstable slip on rough surfaces

In this paper we report results obtained from various friction experiments under direct and oblique shear loading conditions. We used four rock types of varying brittleness (quartzite, anhydrite, limestone, pyrophyllite) with different surface roughness. The observations concentrate on the time span several milliseconds before dynamic failure occurs. During this period a premonitory, unstable phase of slip (slip 2) occurs. This differs importantly from a premonitory, stable process (slip 1) with durations of hundreds of seconds. On smooth surfaces slip 2 is usually observed with ductile rocks and less reliably with brittle rocks. Slip 2 is mostly accompanied by acoustic emissions, which increase in rate of occurrence and in magnitude until the stick-slip event. Foreshocks are observed during approximately 50% of the slip 2 events on rough surfaces. Foreshocks far exceed the acoustic noise level, which is also prevalent before stick-slip events on rough surfaces. In the direct shear experiment, where two faults are being loaded simultaneously, in about 20% of the cases precursory slip 2 was observed on the opposite side on which the final stick-slip event occurred.     

A. International Symposium on Applications of Geodesy to Engineering May 13–17, Stuttgart, Germany

InformationsIAG Scientific Meetings in 1991

A. International Symposium on Applications of Geodesy to Engineering May 13–17, Stuttgart, Germany     

800-yr-long records of annual air temperature and precipitation over southern Siberia inferred from Teletskoye Lake sediments

A unique 800-yr-long record of annual temperatures and precipitation over the south of western Siberia has been reconstructed from the bottom sediments of Teletskoye Lake, Altai Mountains using an X-ray fluorescence scanner (XRF) providing 0.1-mm resolution timeseries of elemental composition and X-ray density (XRD). Br content appears to be broadly correlative with mean annual temperature variations because of changes in catchment vegetation productivity. Sr/Rb ratio reflects the proportion of the unweathered terrestrial fraction. XRD appears to reflect water yield regime and sediment flux. Sedimentation is rather continuous because annual clastic supply and deposited mass are the same. The artificial neural networks method was applied to convert annual sedimentary time-series of XRD, Br content, and Sr/Rb ratio to annual records of temperature and precipitation using a transfer function. Comparison of these reconstructed Siberian records with the annual record of air temperature for the Northern Hemisphere shows similar trends in climatic variability over the past 800 yr. Estimated harmonic oscillations of temperature and precipitation values for both historical and reconstructed periods reveal subdecadal cyclicity.     

Non-hyperbolic moveout inversion of wide-azimuth P-wave data for orthorhombic media

The azimuthally varying non‐hyperbolic moveout of P‐waves in orthorhombic media can provide valuable information for characterization of fractured reservoirs and seismic processing. Here, we present a technique to invert long‐spread, wide‐azimuth P‐wave data for the orientation of the vertical symmetry planes and five key moveout parameters: the symmetry‐plane NMO velocities, V⁽¹⁾_nmo and V⁽²⁾_nmo , and the anellipticity parameters, η⁽¹⁾, η⁽²⁾ and η⁽³⁾ . The inversion algorithm is based on a coherence operator that computes the semblance for the full range of offsets and azimuths using a generalized version of the Alkhalifah–Tsvankin non‐hyperbolic moveout equation. The moveout equation provides a close approximation to the reflection traveltimes in layered anisotropic media with a uniform orientation of the vertical symmetry planes. Numerical tests on noise‐contaminated data for a single orthorhombic layer show that the best‐constrained parameters are the azimuth ? of one of the symmetry planes and the velocities V⁽¹⁾_nmo and V⁽²⁾_nmo , while the resolution in η⁽¹⁾ and η⁽²⁾ is somewhat compromised by the trade‐off between the quadratic and quartic moveout terms. The largest uncertainty is observed in the parameter η⁽³⁾ , which influences only long‐spread moveout in off‐symmetry directions. For stratified orthorhombic models with depth‐dependent symmetry‐plane azimuths, the moveout equation has to be modified by allowing the orientation of the effective NMO ellipse to differ from the principal azimuthal direction of the effective quartic moveout term. The algorithm was successfully tested on wide‐azimuth P‐wave reflections recorded at the Weyburn Field in Canada. Taking azimuthal anisotropy into account increased the semblance values for most long‐offset reflection events in the overburden, which indicates that fracturing is not limited to the reservoir level. The inverted symmetry‐plane directions are close to the azimuths of the off‐trend fracture sets determined from borehole data and shear‐wave splitting analysis. The effective moveout parameters estimated by our algorithm provide input for P‐wave time imaging and geometrical‐spreading correction in layered orthorhombic media.     

Tin–polymetallic sulfide deposits in the eastern part of the Dachang tin field (South China) and the role of black shales in their origin

The Dafulou and Huile vein and stratabound cassiterite-sulfide deposits and sheeted ore veins at the Kangma cassiterite-sulfide deposit are located in the eastern part of the Dachang tin field. These deposits are hosted in a sedimentary sequence containing significant concentrations of organic matter in the form of Lower Devonian calcareous black shales and hornfels. These rocks together with the younger intrusion of Longxianggai granite (91DŽ Ma) actively participated in the formation of Sn-polymetallic deposits. The following three major stages have been distinguished in stratiform and vein-type orebodies at Dafulou, Huile and Kangma: stage I (cassiterite, pyrrhotite, arsenopyrite, tourmaline, carbonate), stage II - main sulfide stage (quartz, cassiterite, arsenopyrite, pyrrhotite, sphalerite, stannite, pyrite, carbonates) and stage III (native Bi, galena, electrum, sulfosalts). Stage IV (post-ore), recognized at Huile is represented by barren carbonates and zeolites. Whole rock geochemistry has revealed that at Dafulou, Bi and Cu correlate strongly with S, whereas V and Pb correlate well with C_org (organic carbon). The similar distribution patterns of selected elements in average slightly mineralized low-Ca black shales indicate a fluid composition similar for all deposits studied. Studies of graphitization of the organic matter in black shales adjacent to orebodies indicate that d₍₀₀₂₎ and FWHM (full width in half maximum)/peak height values gradually decrease in the following sequence: Dafulou deposit M Kangma deposit M Huile deposit. The pyrolysate of wall rocks at the Dafulou deposit is relatively enriched in asphaltenes and maltenes (55.6-72.0% of the pyrolysate) comparable with pyrolysate obtained from more distal black shales (19.2-28.5%). Typical GC-MS spectra of pyrolysate from distal black shales are dominated by alkanes in the n-C15 to n-C25 range, aromatic molecules being represented mostly by alkyl-naphthalenes. In contrast, only traces of aliphatic hydrocarbons in the n-C14 to n-C18 range and elemental sulfur were identified in pyrolysates from pyrrhotitized wall rocks. The earliest fluid inclusions of the studied system occur in the quartz-tourmaline-cassiterite assemblage of stage I at Dafulou. These inclusions are H₂O-CO₂-CH₄-rich, with 10 to 20 vol% of aqueous phase. P-T conditions of the trapping of inclusions are estimated to be up to 400 °C and 1.3 to 2.0 kbar (between 5.0 and 7.5 km under lithostatic pressure). In contrast, the presence of a low density gaseous CO₂-CH₄ phase indicates relatively low pressures during the formation of the breccia-type quartz-calcite-cassiterite-sulfide mineralization (stage II), when P-T conditions probably reached approx. 380 to 400 °C and 0.6 kbar (up to 6 km under hydrostatic pressure). Fluid inclusion data and oxygen isotope thermometry indicate that cassiterite-sulfide ores of the main sulfide stage (stage II) formed from aqueous-carbonic fluid (CO₂/CH₄ =ᄺ) at temperatures of up to 390 °C at Dafulou and in a temperature range of 250 to 360 °C at Huile and 260 to 370 °C at Kangma. The '³⁴S values of sulfides from Dafulou range mostly between -1 and -6‰, whereas sulfides from the Kangma and Huile deposits are characterized by more negative '³⁴S values (between -8 and -11‰, and between -9 and -12‰, respectively). These data suggest that bacteriogenic sulfides of black shales were a dominant source of reduced sulfur for epigenetic (vein and replacement) mineralization. Oxygen isotopic compositions of five quartz-cassiterite pairs from Dafulou and Huile show a relatively narrow range of calculated oxygen isotope temperatures (250-320 °C, using the equation of Alderton 1989) and high '¹⁸O_fluid values between +8 and +10‰ (SMOW), which are in agreement with fluid derivation from and/or high temperature equilibration with the Longxianggai granite. The carbon and oxygen isotope composition of carbonates reflects variable carbon sources. Stage I calcite is characterized by narrow ranges of '¹³C (-7.0 to -9.5‰ PDB) and '¹⁸O (+15.0 to +17.5‰ SMOW). This calcite shows ubiquitous deformation, evidenced by intense development of twins. Fluid compositions calculated at 330 °C for the Dafulou and Huile deposits and at 270-300 °C for the Kangma deposit ('¹⁸O_fluid between +10.0 and +11.5‰ SMOW, '¹³C_fluid between -5.5 and -7.5‰ PDB), agree with fluid derivation from and/or equilibration with the peraluminous, high-'¹⁸O Longxianggai granite and suggest a significant influence of contact metasedimentary sequences (carbon derived from decomposition and/or alteration of organic matter of calcareous black shales). The '¹³ C values of organic matter from the Lower to Upper Devonian host rocks at the Dafulou deposit (-24.0 and -28.0‰) fit with a marine origin from algae. However, organic matter adjacent to the host rock-ore contact displays a slight enrichment in ¹³C. The organic carbon from the Huile and Kangma deposits is even more ¹³C enriched (-24.6 to -23.5‰). The most heavy '¹³ C values (-16.5‰) were detected in hornfels sampled at the contact of the Upper Devonian sediments with the Longxianggai granite. The '¹³C data broadly correlate with the degree of structural ordering (degree of graphitization) of organic matter, which indicates that both variables are related to thermal overprint.