Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

ASPI experiment: measurements of fields and waves on board the INTERBALL-1 spacecraft

The plasma-wave experiment ASPI (analysis of spectra of plasma waves and instabilities) on board the INTERBALL spacecraft is a combined wave diagnostics experiment. It performs measurements of the DC and AC magnetic field vector by flux-gate and searchcoil sensors, the DC and AC electric field vector by Langmuir double probes and the plasma current by Langmuir split probe. Preliminary data analysis shows the low noise levels of the sensors and the compatibility of new data with the results of previous missions. During several months of in-orbit operation a rich collection of data was acquired, examples of which at the magnetopause and plasma sheet are presented in second part of the paper.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of water level fluctuation on the spatial dynamic of the water chemistry at Lake Fertõ/Neusiedler See

Lake Fertõ/Neusiedler See is a steppe lake situated on the Hungarian-Austrian border with a surface area of 309 km². It is an extremely shallow lake with regulated outflow. 54% of the whole lake and 85% of the Hungarian part (75 km²) is covered by reed. There are numerous reedless areas (inner ponds) of variable size within the reed belt, which is enmeshed with canals connecting the inner ponds with the open water areas. Fundamental differences were recorded in the water chemistry of the three water type areas, open water, canals and inner ponds. The degree of these differences is basically determined by the distance from the open water areas and the water level of the actual years. An extremely high, an average and an extremely low water level year was examined, in order to present the spatial and temporal differences in the water chemistry of the investigated water bodies. Differences due to the water level fluctuation between years in the three different water types are discussed on the basis of electrical conductivity, cation, and anion concentration changes of the water. Our results support the concept of high instability of Lake Fertõ and describe habitats with remarkable temporal and spatial variability. The high influence of water level on habitat conditions at different sites of the lake recommend a carefull water level regulation strategy.     

A statistical analysis of methanol maser groups

An analysis of the flux densities of the 5₁-6₀ A ⁺ (6.7 GHz) and 2₀-3_?1 E (12.2 GHz) class II methanol maser lines in a large and homogeneous sample of maser sources has been carried out. For convenience, the maser lines were divided into three groups: group I contains spectral features for the lines most prominent in the 5₁-6₀ A ⁺ (6.7 GHz) transition, group II contains spectral features for the lines strongest in the 2₀-3_?1 E (12.2 GHz) transition, group III contains spectral features for which the velocities of the emission maxima of the two lines coincide. The same dependence was found for group II and group III: log S _6.7=(0.79±0.05)×log S _12.2+(0.79±0.05). The spectral features in group I do not obey this relation, and deviations from a linear dependence are considerably greater. It is suggested that methanol class II masers be divided into a subclass IIa, which has special conditions favoring 6.7 GHz masers, and a subclass IIb, which is comprised of the 12.2 GHz masers and those 6.7 GHz masers that necessarily accompany them under the same conditions.     

Relative variations of the physical parameters of variable extragalactic radio sources

A method for studying the physical conditions in compact components of extragalactic radio sources displaying variability on time scales of hundreds of days is proposed. The method can be used to estimate the relative variations of the magnetic-field strength and number density of relativistic electrons in superluminal jets from the cores of quasars and radio galaxies. Results are presented for the jets of the quasars 3C 120, 3C 273, 3C 279, and 3C 345. The energies of the magnetic field and relativistic particles in these objects are not in equipartition. As a rule, the magnetic-field strength decreases appreciably during the evolution of an expanding jet, while the number of relativistic electrons grows.     

Rainfall induced debris flows in pyroclastic deposits, Campania (southern Italy)

A combination of empirical and physically based hydrological models has been used to analyze historical data on rainfall and debris-flow occurrence in western Campania, to examine the correlation between rainfall and debris-flow events.

Rainfall data from major storms recorded in recent decades in western Campania were compiled, including daily series from several rain gauges located inside landslide areas, supplemented by hourly rainfall data from some of the principal storms.

A two-phase approach is proposed. During phase 1, soil moisture levels have been modelled as the hydrological balance between precipitation and evapotranspiration, on a daily scale, using the method of Thornthwaite [Geograph. Rev. 38 (1948) 55].

Phase 2 is related to the accumulation of surplus moisture from intense rainfall, leading to the development of positive pore pressures. These interactions take place on an hourly time scale by the “leaky barrel” (LB) model described by Wilson and Wiezoreck [Env. Eng. Geoscience, 1 (1995) 11]. In combination with hourly rainfall records, the LB model has been used to compare hydrological effects of different storms. The critical level of retained rain water has been fixed by the timing of debris-flow activity, related to recorded storm events.

New rainfall intensity–duration thresholds for debris-flow initiation in western Campania are proposed. These thresholds are related to individual rain gauge and assume a previously satisfied field capacity condition. The new thresholds are somewhat higher than those plotted by previous authors, but are thought to be more accurate and thus need less conservatism.     

Organic matter in a low productivity anoxic intraplatform basin in the Triassic Tethys

This paper presents isotopic, bulk geochemical and biomarker data measured on organic matter accumulated in a narrow extensional basin developed at the oceanward margin of the huge Triassic carbonate platform in the Alps–Appenines domain. The integration of isotope signatures, organic petrographical and biomarker evidence together with the composition of kerogen pyrolysates suggests immature organic matter predominantly of algal origin with a minor, but not negligible, higher plant derived and moderate bacterial contribution for the entire sequence. The mineral sources are dominated by platform-derived subtidal Dachstein limestone with a minor palaeosol input and a moderate contribution of autochtonous quartz. Nevertheless, parallel variations observed in the mineral content, as well as in the amount and the quality of the organic matter reflect variations in the palaeoenvironment. The increased humidity, existing in the period of the accumulation of the upper section of the sequence, led to the restriction of dolomitization. A slightly greater higher plant derived contribution, in this section, is evidenced by the composition of bitumen and the results on GC and GC/MS analyses on the non-aromatic hydrocarbon fraction of bitumen. Moreover, the climate-induced weathering enchanced the primary productivity and resulted in a pronounced increase in the TOC content. The average estimated value of the planktonic productivity is about three times higher for the calcite-rich sequence than the dolomite-rich one, being 44.2 and 15.3 tC_org/m²/Ma, respectively. The low to moderate planktonic productivity shows that anoxic conditions, observed for the entire succession, are a consequence of the stagnant water stratification rather than high planktonic productivity. Depth trends in the data measured on kerogens (HI, OI, δ¹³Corg values, composition of pyrolysate) together with the δ¹⁸O excursions and Δδ¹³C values appear to be controlled by sea-level fluctuations. Consistent with the high abundance of alkyltiophenes in kerogen pyrolysates, the high S_org/C ratios (ranging between 0.05 and 0.10) suggest the importance of natural sulfurization in the formation of the sulfur-rich type-II-S kerogen occurring in all of the samples.     

Dissipationless disk accretion

We consider disk accretion resulting purely from the loss of angular momentum due to the outflow of plasma from a magnetized disk. In this limiting case, the dissipation due to the viscosity and finite electrical conductivity of the plasma can be neglected. We have obtained self-consistent, self-similar solutions for dissipationless disk accretion. Such accretion may result in the formation of objects whose bolometric luminosities are lower than the flux of kinetic energy in the ejected material.