Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.     

Magnetostratigraphic correlation of the Oxfordian–Kimmeridgian boundary

A magnetic polarity pattern for Boreal and Sub-Boreal ammonite zones of the Upper Oxfordian to Lower Kimmeridgian was established and confirmed in four British sections, including the proposed Global Boundary Stratotype Section and Point (GSSP) on the Isle of Skye (Scotland) to define the base of the international Kimmeridgian Stage. A coeval pattern for Sub-Mediterranean ammonite zones was compiled from seven sections in Poland, one German section and multi-section composites from France and Spain. The mean paleopole for the European Craton (excluding Spain) at the Oxfordian–Kimmeridgian boundary is 74.2°N, 181.3°E (Α₉₅ = 3.8°). The common magnetic polarity scale enables inter-correlation of ammonite subzones among these three faunal provinces and to the marine magnetic-anomaly M-Sequence. The proposed GSSP at the base of the Pictonia baylei Zone is near the base of an extended interval dominated by reversed polarity, which is interpreted to be Chron M26r. This GSSP level projects to the lower to middle part of the Epipeltoceras bimammatum Subzone, which is the middle subzone of this E. bimammatum Zone in the Sub-Mediterranean standard zonation. In contrast, the traditional placement of the Oxfordian–Kimmeridgian boundary in that Sub-Mediterranean standard zonation (base of Sutneria platynota Zone) is at the base of Chron M25r, or nearly 1 million years younger.     

Meltwater discharge through the subglacial bed and its land‐forming consequences from numerical experiments in the Polish lowland during the last glaciation

Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady‐state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice‐sheet‐thickness‐modulated, down‐ice‐decreasing hydraulic heads. In relation to non‐glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice–bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo‐ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice–bed interface. Owing to the abundance of low‐permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd.     

An attempt to test Ambartsumian's idea of the origin of the galaxy . II. location of galaxies within clusters and galaxy position angles

The location of galaxies in 377 rich Abell galaxy clusters is discussed. We compared the distributions of galaxies in the sample containing all galaxies with the sample compiled from the 20 brightest objects. Counts in circular sectors with angle width equal to 30î show isotropy. Only in the case of BM I clusters and the coordinate system related to the cluster major axis did we find anisotropic distributions. We investigated also the distributions of galaxy position angles within clusters exhibiting isotropy. The structure position angles for both samples have been studied as well. They appeared to be different. The difference in location of galaxies in the case of cD clusters, as well as the more elliptical shape and different position angles of samples containing bright galaxies is observed.     

Is Newton's iteration faster than simple iteration for transformation between geocentric and geodetic coordinates?

Summary Two iterative algorithms for transformation from geocentric to geodetic coordinates are compared for numerical efficiency: the well known Bowring's algorithm of 1976, which employs the method of simple iteration, and the recent (1989) algorithm by Borkowski, which employs the Newton-Raphson method. The results of numerical tests suggest that the simple iteration method implemented in Bowring's algorithm executes approximately 30% faster than the Newton-Raphson method implemented in Borkowski's algorithm. Only two iterations of each algorithm are considered. Two iterations are sufficient to produce coordinates accurate to the comparable level of 1E-9 m, which exceeds the requirements of any practical application. Therefore, in the class of iterative methods, the classical Bowring's algorithm should be the method of choice.     

The New Algorithm for Fast Probabilistic Hypocenter Locations

The spatial location of sources of seismic waves is one of the first tasks when transient waves from natural (uncontrolled) sources are analysed in many branches of physics, including seismology, oceanology, to name a few. It is well recognised that there is no single universal location algorithm which performs equally well in all situations. Source activity and its spatial variability in time, the geometry of recording network, the complexity and heterogeneity of wave velocity distribution are all factors influencing the performance of location algorithms. In this paper we propose a new location algorithm which exploits the reciprocity and time-inverse invariance property of the wave equation. Basing on these symmetries and using a modern finite-difference-type eikonal solver, we have developed a new very fast algorithm performing the full probabilistic (Bayesian) source location. We illustrate an efficiency of the algorithm performing an advanced error analysis for 1647 seismic events from the Rudna copper mine operating in southwestern Poland.     

On the variability and spectral distortion of fluorescent iron lines from black hole accretion discs

We investigate the properties of fluorescent iron lines that arise as a result of the illumination of a black hole accretion disc by an X-ray source located above the disc's surface. We study in detail the light-bending model of the variability of the lines, extending previous work on the subject. We indicate that the bending of photon trajectories to the equatorial plane (a distinct property of the Kerr metric) is the most feasible effect underlying the reduced variability of the lines observed in several objects. A model involving an X-ray source with a varying radial distance, located within a few central gravitational radii around a rapidly rotating black hole, close to the disc's surface, may explain both the elongated red wing of the line profile and the complex variability pattern observed in MCG–6-30-15 by XMM–Newton . We also point out that illumination by radiation that returns to the disc (following the previous reflection) contributes significantly to the formation of the line profile in some cases. As a result of this effect, the line profile always has a pronounced blue peak (which is not observed in the deep minimum state in MCG–6-30-15), unless the reflecting material is absent within the innermost 2–3 gravitational radii.     

Short-Period Secular Variations (SPSV) of the geomagnetic field recorded in highly scattered palaeomagnetic records of Holocene lake sediments from North Poland

A simple method of calculating running means over different intervals within palaeomagnetic records has been used to identify short-period secular variations of the geomagnetic field. The records from Polish lake sediments reveal variations in declination with a period between 290 and 372 years. The periodicity of variations in the inclination has been estimated to be between 522 and 670 years.