Crystal chemistry of selenates with mineral-like structures: V. Crystal structures of (H3O)2[(UO2)(SeO4)2(H2O)](H2O)2 and (H3O)2[(UO2)(SeO4)2(H2O)](H2O), new compounds with rhomboclase and goldichite topology

The crystal structures of two new compounds (H₃O)₂[(UO₂)(SeO₄)₂(H₂O)](H₂O)₂ (1, orthorhombic, Pnma, a = 14.0328(18), b = 11.6412(13), c = 8.2146(13) Å, V = 134.9(3) ų) and (H₃O)₂[(UO₂)(SeO₄)₂(H₂O)](H₂O) (2, monoclinic, P2₁/c, a = 7.8670(12), b = 7.5357(7), c = 21.386(3) Å, β = 101.484(12)°, V = 1242.5(3) ų) have been solved by direct methods and refined to R ₁ = 0.076 and 0.080, respectively. The structures of both compounds contain sheet complexes [(UO₂)(SeO₄)₂]^2? formed by cornershared [(UO₂)O₄(H₂O)] bipyramids and SeO₄ tetrahedrons. The sheets are parallel to the (100) plane in structure 1 and to (?102) in structure 2. The [(UO₂)(SeO₄)₂(H₂O)]^2? layers are linked by hydrogen bonds via interlayer groups H₂O and H₃O⁺. The sheet topologies in structures 1 and 2 are different and correspond to the topologies of octahedral and tetrahedral complexes in rhomboclase (H₂O₂)⁺[Fe(SO₄)₂(H₂O)₂] and goldichite K[Fe(SO₄)₂(H₂O)₂](H₂O)₂, respectively.     

Mineral assemblages of titanium-zirconium sands at the Central deposit,the east European platform

The study of mineral assemblages at the Central deposit allowed us to substantially refine the evolutionary model of the deposit and reveal the two main factors that control the variability of its mineralic space: (1) heterogeneity of the ore layer, consisting of a sublittoral bottom placer (the lower part) and a subaerial dune complex partly reworked in the course of a new transgression (the upper part), and (2) postore epigenetic alteration of the ore constituent of sands, which affected the quality of ore concentrates. The results obtained will be used in geological and technological mapping and development of the production program.     

Solar flares associated with singularities in the differential characteristics of magnetic fields

We analyze the flare events of October 28 and May 27, 2003 to examine the possible connection of solar flares with singularities in the differential characteristics of solar magnetic fields. The relation between flares and the behavior of F indicators (which describe structural singularities of the magnetic fields) introduced earlier is analyzed. It is shown that the generation of flares can be associated with self-intersections (or reconnection) of the F = 0 curves, where magnetic reconnection can occur. Consequently, the current sheets generated in such regions can be sources of energy release.     

Temporal variations in the position of the heliospheric equator

It is shown that the centroid of the heliospheric equator undergoes quasi-periodic oscillations. During the minimum of the 11-year cycle, the centroid shifts southwards (the so-called bashful-ballerina effect). The direction of the shift reverses during the solar maximum. The solar quadrupole is responsible for this effect. The shift is compared with the tilt of the heliospheric current sheet.     

Rotational and magnetic corrections to the g-mode frequencies in a one-dimensional MHD model

The central magnetic field and rotation of the solar radiative zone are responsible for corrections to the g-mode frequencies. Magnetogravitational spectra are calculated analytically in a simple one-dimensional MHD model that goes beyond the WKB approximation and avoid any cusp resonances that trap the wave within the radiative zone in the presence of a weak magnetic background. The calculations are compared with spacecraft observations of the 1% frequency shifts for candidate g-modes found in the SOHO GOLF experiment. The magnetic correction is the main contribution for a strong magnetic field satisfying the approximation used. It is shown that a constant magnetic field of 700 kG in the radiative zone provides the required frequency shift for the n = ?10 g-mode. The rotational correction, which is due to the Coriolis force in the one-dimensional model used, is much less than a percent (α_Ω ≤ 0.003).     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Geochemical Eccentricity of Ground Water Allied to Weathering of Basalts from the Deccan Volcanic Province,India: Insinuation on CO<Subscript>2</Subscript> Consumption

Analyses of 72 samples from Upper Panjhara basin in the northern part of Deccan Plateau, India, indicate that geochemical incongruity of groundwater is largely a function of mineral composition of the basaltic lithology. Higher proportion of alkaline earth elements to total cations and HCO₃>Cl + SO₄ reflect weathering of primary silicates as chief source of ions. Inputs of Cl, SO₄, and NO₃ are related to rainfall and localized anthropogenic factors. Groundwater from recharge area representing Ca + Mg–HCO₃ type progressively evolves to Ca + Na–HCO₃ and Na–Ca–HCO₃ class along flow direction replicates the role of cation exchange and precipitation processes. While the post-monsoon chemistry is controlled by silicate mineral dissolution + cation exchange reactions, pre-monsoon variability is attributable chiefly to precipitation reactions + anthropogenic factors. Positive correlations between Mg vs HCO₃ and Ca + Mg vs HCO₃ supports selective dissolution of olivine and pyroxene as dominant process in post-monsoon followed by dissolution of plagioclase feldspar and secondary carbonates. The pre-monsoon data however, points toward the dissolution of plagioclase and precipitation of CaCO₃ supported by improved correlation coefficients between Na + Ca vs HCO₃ and negative correlation of Ca vs HCO₃, respectively. It is proposed that the eccentricity in the composition of groundwater from the Panjhara basin is a function of selective dissolution of olivine > pyroxene followed by plagioclase feldspar. The data suggest siallitization (L < R and R _k) as dominant mechanism of chemical weathering of basalts, stimulating monosiallitic (kaolinite) and bisiallitic (montmorillonite) products. The chemical denudation rates for Panjhara basin worked out separately for the ground and surface water component range from 6.98 to 36.65 tons/km²/yr, respectively. The values of the CO₂ consumption rates range between 0.18 × 10⁶ mol//km²/yr (groundwater) and 0.9 × 10⁶ mol/km²/yr (surface water), which indicates that the groundwater forms a considerable fraction of CO₂ consumption, an inference, that is, not taken into contemplation in most of the studies.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.     

Coseismic responses and the mechanism behind MW 5.1 earthquake of March 14, 2005 in the Koyna–Warna region, India

An earthquake of M_w 5.1 occurred on March 14, 2005, in the seismically active Koyna–Warna region in western India, the site known for the largest reservoir triggered seismicity (RTS) in the world. For more than four decades, earthquakes with M  4.0 have occurred in this region at regular intervals. Impoundment of reservoirs and changes in lake levels can trigger earthquakes by two processes of stress modifications, namely direct loading effect of the reservoir and diffusion through various faults and fractures. In this paper we analysed the reservoir water level data at Koyna and Warna reservoirs prior to the occurrence of the March 14, 2005 earthquake, to explain the dominant mechanism behind its occurrence and its correlation with the observed coseismic changes. We conclude that the diffusion process, not the reservoir load effect, is the dominating mechanism triggering earthquakes in the region. The coseismic changes in deep well water levels sensitive to earth tides are found to be to the order of 1–12 cm.