Identification of heterogeneous aquifer transmissivity using an AE-based method

Determination of hydraulic head, H, as a function of spatial coordinates and time, in ground water flow is the basis for aquifer management and for prediction of contaminant transport. Several computer codes are available for this purpose. Spatial distribution of the transmissivity, T(x,y), is a required input to these codes. In most aquifers, T varies in an erratic manner, and it can be characterized statistically in terms of a few moments: the expected value, the variance, and the variogram. Knowledge of these moments, combined with a few measurements, permits one to estimate T at any point using geostatistical methods. In a review of transmissivity data from 19 unconsolidated aquifers, Hoeksema and Kitanidis (1985) identified two types of the logtransmissivity Y= ln(T) variations: correlated variations with variance sigma2Yc and correlation scale, I(Y), on the order of kilometers, and uncorrelated variations with variance sigma2Yn. Direct identification of the logtransmissivity variogram, Gamma(Y), from measurements is difficult because T data are generally scarce. However, many head measurements are commonly available. The aim of the paper is to introduce a methodology to identify the transmissivity variogram parameters (sigma2Yc, I(Y), and sigma2Yn) using head data in formations characterized by large logtransmissivity variance. The identification methodology uses a combination of precise numerical simulations (carried out using analytic element method) and a theoretical model. The main objective is to demonstrate the application of the methodology to a regional ground water flow in Eagle Valley basin in west-central Nevada for which abundant transmissivity and head measurements are available.     

Spontaneous Formation of Magnetic Flux Concentrations in Stratified Turbulence

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNSs) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in the formation of large-scale inhomogeneous magnetic structures from an initially uniform magnetic field. Earlier investigations of this instability in DNSs of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma ?? are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L ²/?? _t, where ?? _t is the turbulent diffusivity and L is the scale of the domain. DNSs are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.     

Structure of the Wollaston Lake Reflector (Trans-Hudson Orogen, Canada) from reflection AVO analysis: Fractured diabase intrusion, fluids, or silicified shear zone?

A nearly 160-km long Wollaston Lake Reflector (WLR) observed in seismic reflection profile S2b of the 1994 Lithoprobe Trans-Hudson Orogen transect (THOT) in northern Saskatchewan (Canada) is among the most spectacular and well-recorded features imaged within the crystalline crust. Based on modeling of its normal-incidence reflectivity, the observed bright spot reflector was originally interpreted as a series of tabular diabase intrusions. In order to further elucidate its structure, we reprocessed line S2b and analysed the WLR for the Amplitude Variations with Offset (AVO). By contrast to conventional (approximate) AVO analysis, we used the exact Zoeppritz equations and considered a thin-layer (mixed positive and negative polarities) reflectivity. The results suggest two possible interpretations of the WLR: 1) the reflector caused by a massive mafic intrusion as suggested earlier, in which case the intruded rocks should have anomalous Poisson's ratios of σ ≥ 0.33, and 2) the WLR represents a silicified shear zone, with only moderate (e.g.,  5–10%) alteration of the host rock and σ ≈ 0.2. Although both of these models may to some extent co-exist within the WLR, its brightness, sharpness, great lateral extent and smooth shape favour the second interpretation.In both models, a fractured fluid-filled zone within a major crustal detachment should have played a key role in the formation of the WLR. The association of the reflector with laterally- and depth-migrating fluids is also supported by magneto-telluric measurements of crustal conductivity beneath the WLR. Analogies from the studies of the Kola Superdeep Borehole (Russia), where free or metamorphic fluids were found at comparable depths, also suggest that fluids may contribute to WLR structure and formation.     

Zinc-rich and copper-bearing amphiboles from sulfide-free ore occurrences of the Pelagonian Massif,Republic of North Macedonia

Mineralogy and Petrology - This paper describes specific features of isomorphism of unusual amphiboles containing up to 23 wt% ZnO and up to 1.3 wt% CuO from sulfide-free...     

Comparative geochemical study of the cherty rocks of the Taukha terrane (Sikhote-Alin) and its paleogeodynamic significance

Acta Geochimica - The article presents the results of a comparative geochemical study of cherts and siliceous–clayey rocks composing the tectono-sedimentary complexes of various structural...     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part II: Pre-critical slowness match

Part II of this paper is a direct continuation of Part I, where we consider the same types of orthorhombic layered media and the same types of pure-mode and converted waves. Like in Part I, the approximations for the slowness-domain kinematical characteristics are obtained by combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. In Part I, the wide-angle ray was set to be the critical ray (‘critical slowness match’), whereas in Part II we consider a finite long offset associated with a given pre-critical ray (‘pre-critical slowness match’). Unlike the critical slowness match, the approximations in the pre-critical slowness match are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. Moreover, for the pre-critical slowness match, there is no need to distinguish between the high-velocity layer and the other, low-velocity layers. The form of the approximations in both critical and pre-critical slowness matches is the same, where only the wide-angle power series coefficients are different. Comparing the approximated kinematical characteristics with those obtained by exact numerical ray tracing, we demonstrate high accuracy. Furthermore, we show that for all wave types, the accuracy of the pre-critical slowness match is essentially higher than that of the critical slowness match, even for matching slowness values close to the critical slowness. Both approaches can be valuable for implementation, depending on the target offset range and the nature of the subsurface model. The pre-critical slowness match is more accurate for simulating reflection data with conventional offsets. The critical slowness match can be attractive for models with a dominant high-velocity layer, for simulating, for example, refraction events with very long offsets.     

Earth mantle conductivity beneath the Ukrainian territory

Studia Geophysica et Geodaetica - Magnetovariation methods, which are applicable to study the mantle conductivity, require long lasting registration of natural magnetic field variations. Such data...     

Pore to pore validation of pore network modelling against micromodel experiment results

Pore network modelling (PNM) has been widely used to study the multiphase flow and transport in porous media. Although a number of recent papers discussed the PNM validation on core-scale parameters such as permeability, relative permeability and capillary pressure; quantitative predictive potential of PNM on pore by pore basis has rarely been studied. The aim of this paper is to present a direct comparison between PNM simulations and corresponding micro-model experiments at the same scale and the same geometry. A number of well-defined and constrained two-phase flow in porous medium experimental scenarios were utilized to validate the physics solving part in PNM (filling rules, capillary and viscous pressure). This work validates that a dynamic pore network flow solver can predict two-phase flow displacements for these experiments for drainage situations at both pore and plug scales. A glass-etched micro-model is used to quantify the accuracy of a dynamic PNM solver on pore and core levels. Two-phase drainage micro fluidic experiments at different flow conditions are performed on micro-models. PNM simulations are performed on the same pattern and flow conditions as used in micro-model experiments. The two-phase distribution extracted from experiment images is registered onto rsults of PNM simulations for direct pore to pore comparison. Pore-scale matching level is found at around 75 % for all three test cases. The matching level of core-scale parameters such as S _{w c} and oil-phase permeability varies from case to case; the relative error to micro-model experiment measurements varies from 15 to 60 %. Possible reasons leading to discrepancies on core-scale parameters are discussed: missing considerations during validation of the combination of uncertainty in both simulator input parameters and experiments are seen as the principal factors.     

Northern Adriatic phytoplankton response to short Po River discharge pulses during summer stratified conditions

A 275‐km‐long transversel Northern Adriatic profile from the mouth of the Po River (Italian Adriatic coast) to the Kvarner region (Croatian coastal island area) was investigated in three successive case studies in August 2008, 2009 and 2010. The short Po River pulses in August result in the surface advection of riverine water, nutrients and phytoplankton from the western to the eastern side of the Adriatic. This surface spreading exhibits inter‐annual variability depending on the riverine discharge in the preceding period. The Po River discharge pulse in August 2010 in particular resulted in an extraordinary tongue‐like advection of riverine water, nutrients, and phytoplankton towards the Eastern Adriatic coast. The phenomenon was detected using both satellite imagery and classical oceanographic measurements. In the advective water, toxic dinoflagellates were most abundant in August 2010, when the influence of the Po was greatest.