A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

Dynamics of Saxothuringian subduction channel/wedge constrained by phase‐equilibria modelling and micro‐fabric analysis

Subduction and exhumation dynamics can be investigated through analysis of metamorphic and deformational evolution of associated high‐grade rocks. The Erzgebirge anticline, which forms at the boundary between the Saxothuringian and Teplá‐Barrandian domains of the Bohemian Massif, provides a useful study area for these processes owing to the occurrence of numerous meta‐basites preserving eclogite facies assemblages, and coesite and diamond bearing quartzofeldspathic lithologies indicating subduction to deep mantle depths. The prograde and retrograde evolution of meta‐basite from the Czech portion of the Erzgebirge anticline has been constrained through a combination of thermodynamic modelling and conventional thermobarometry. Garnet growth zoning indicates that the rocks underwent burial and heating to peak conditions of 2.6 GPa and at least 615 °C. Initial exhumation occurred with concurrent cooling and decompression resulting in the growth of amphibole and zoisite poikiloblasts overgrowing and including the eclogite facies assemblage. The development of clinopyroxene–plagioclase–amphibole symplectites after omphacite and Al‐rich rims on matrix amphibole indicate later heating at the base of the lower crust. Omphacite microstructures, in particular grain size analysis and lattice‐preferred orientations, indicate that the prograde evolution was characterized by a constrictional strain geometry transitioning into plane strain and oblate fabrics during exhumation. The initial constrictional strain pattern is interpreted as being controlled by competing slab pull and crustal buoyancy forces leading to necking of the subducting slab. The transition to plane strain and flattening geometries represents transfer of material from the subducting lithosphere into a subduction channel, break‐off of the dense slab and rebound of the buoyant crustal material.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.     

Biogenic volatile compound emissions from a temperate forest,China: model simulation

Emissions of biogenic volatile organic compounds (BVOC) were measured using a relaxed eddy accumulation (REA) technique on an above-canopy tower in a temperate forest (Changbai Mountain, Jilin province, China) during the 2010 and 2011 summer seasons. Solar global radiation and photosynthetically active radiation (PAR) were also measured. Based on PAR energy dynamic balance, an empirical BVOC emission and PAR transfer model was developed that includes the processes of BVOC emissions and PAR transfer above the canopy level, including PAR absorption and consumption, and scattering by gases, liquids, and particles (GLPs). Simulated emissions of isoprene and monoterpenes were in agreement with observations. The averages of the relative estimator biases for the flux were 39.3 % for isoprene, and 27.1 % for monoterpenes in the 2010 and 2011 growing seasons, with NMSE (normalized mean square error) values of 0.133 and 0.101, respectively. The observed and simulated mean diurnal variations of isoprene and monoterpenes in the 2010 and 2011 growing seasons were evaluated for the validation of the empirical model. Under observed atmospheric conditions, the sensitivity analysis showed that emissions of isoprene and monoterpenes were more sensitive to changes in PAR than to water vapor content or to the magnitude of the scattering factor. The emissions of isoprene and monoterpenes in the 2010 and 2011 growing seasons (from June to September) were estimated using this empirical model along with hourly observational data, with mean hourly emissions of 1.71 and 1.55 mg m^?2 h^?1 for isoprene, and 0.48 and 0.47 mg m^?2 h^?1 for monoterpenes in 2010 and 2011, respectively. As formaldehyde (HCHO) is considered as the main oxidation product of isoprene and monoterpenes, it is necessary to investigate the link between HCHO and BVOC emissions. GOME-2 HCHO vertical column densities (VCDs) can be used to estimate BVOC emission fluxes in the Changbai Mountain temperate forest.     

Initiation of large-volume silicic centers in the Yellowstone hotspot track: insights from H<Subscript>2</Subscript>O- and F-rich quartz-hosted rhyolitic melt inclusions in the Arbon Valley Tuff of the Snake River Plain

During the onset of caldera cluster volcanism at a new location in the Snake River Plain (SRP), there is an increase in basalt fluxing into the crust and diverse silicic volcanic products are generated. The SRP contains abundant and compositionally diverse hot, dry, and often low-δ¹⁸O silicic volcanic rocks produced through time during the formation of individual caldera clusters, but more H₂O-rich eruptive products are rare. We report analyses of quartz-hosted melt inclusions from pumice clasts from the upper and lower Arbon Valley Tuff (AVT) to gain insight into the initiation of caldera cluster volcanism. The AVT, a voluminous, caldera-forming rhyolite, represents the commencement of volcanism (10.44 Ma) at the Picabo volcanic field of the Yellowstone hotspot track. This is a normal δ¹⁸O rhyolite consisting of early and late erupted members (lower and upper AVT, respectively) with extremely radiogenic Sr isotopes and unradiogenic Nd isotopes, requiring that ~50 % of the mass of these elements is derived from melts of Archean upper crust. Our data reveal distinctive features of the early erupted lower AVT melt including: variable F concentrations up to 1.4 wt%, homogenous and low Cl concentrations (~0.08 wt%), H₂O contents ranging from 2.3 to 6.4 wt%, CO₂ contents ranging from 79 to 410 ppm, and enrichment of incompatible elements compared to the late erupted AVT, subsequent Picabo rhyolites, SRP rhyolites, and melt inclusions from other metaluminous rhyolites (e.g., Bishop Tuff, Mesa Falls Tuff). We couple melt inclusion data with Ti measurements and cathodoluminescence (CL) imaging of the host quartz phenocrysts to elucidate the petrogenetic evolution of the AVT rhyolitic magma. We observe complex and multistage CL zoning patterns, the most critical being multiple truncations indicative of several dissolution–reprecipitation episodes with bright CL cores (higher Ti) and occasional bright CL rims (higher Ti). We interpret the high H₂O, F, F/Cl, and incompatible trace element concentrations in the context of a model involving melting of Archean crust and mixing of the crustal melt with basaltic differentiates, followed by multiple stages of fractional crystallization, remelting, and melt extraction. This multistage process, which we refer to as distillation, is further supported by the complex CL zoning patterns in quartz. We interpret new Δ¹⁸O(Qz-Mt) isotope measurements, demonstrating a 0.4 ‰ or ~180 °C temperature difference, and strong Sr isotopic and chemical differences between the upper and lower AVT to represent two separate eruptions. Similarities between the AVT and the first caldera-forming eruptions of other caldera clusters in the SRP (Yellowstone, Heise and Bruneau Jarbidge) suggest that the more evolved, lower-temperature, more H₂O-rich rhyolites of the SRP are important in the initiation of a caldera cluster during the onset of plume impingement.     

Recursive algorithm for fast GNSS orbit fitting

Gaussian elimination is an efficient and numerically stable algorithm for estimating parameters and their precision. However, before estimating the parameters, it is often prudent to perform statistical tests to achieve the best fitting model. We use Gaussian elimination to select the best fitting model among candidate models. A succinct relationship between the weighted sum of squared residuals and the previous one is revealed by a volume formula. For quick parameter estimation and determination of weighted sum of squared residuals, a recursive elimination algorithm is proposed in the context of Gaussian elimination. In order to improve the model selection efficiency, the parameter estimation and the determination of the weighted sum of squared residuals are carried out in parallel using the proposed recursive elimination algorithm in which the improvement at each recursive stage is judged by the Bayesian information criterion. Ultimately, the computational complexity and numerical stability of the recursive elimination proposed are briefly discussed, and a GNSS orbit interpolation example is used to verify the results. It shows that the proposed recursive elimination algorithm inherits the numerical stability of the Gaussian elimination, and this algorithm can be used to examine the gain from the newly introduced parameter, dynamically assess the fitting model, and fix the optimal model efficiently. The optimal fitting model with the lowest information is very close to the real situation verified by checkpoints.     

A new method for specular and diffuse pseudorange multipath error extraction using wavelet analysis

Multipath remains one of the major challenges in Global Navigation Satellite System (GNSS) positioning because it is considered the dominant source of ranging errors, which can be classified into specular and diffuse types. We present a new method using wavelets to extract the pseudorange multipath in the time domain and breaking it down into the two components. The main idea is an analysis-reconstruction approach based on application of both continuous wavelet transform (CWT) and discrete wavelet transform (DWT). The proposed procedure involves the use of L1 code-minus-carrier (CMC) observable where higher-frequency terms are isolated as residuals. CMC residuals are analyzed by applying the CWT, and we propose the scalogram as a technique for discerning time–frequency variations of the multipath signal. Unlike Fourier transform, the potential of the CWT scalogram for examining the non-stationary and multifrequency nature of the multipath is confirmed as it simultaneously allows fine detection and time localization of the most representative frequencies of the signal. This interpretation of the CWT scalogram is relevant when choosing the levels of reconstruction with DWT, allowing accurate time domain extraction of both the specular and diffuse multipath. The performance and robustness of the method and its boundary applicability are assessed. The experiment was carried out using a receiver of Campania GNSS Network. The results are given in which specular multipath error is achieved using DWT level 7 approximation component and diffuse multipath error is achieved using DWT level 6 denoised detail component.     

Optimization and Parametrical Investigation to Assess the Reconstitution of Different Types of Indian Sand Using Portable Travelling Pluviator

By using the air pluviation technique, it is aimed to achieve the desired relative density with uniform void ratio throughout the specimen in order to maintain the homogeneity and to avoid the spatial variability. Further, in order to achieve the maximum deposition intensity, a systematic optimization study has been carried out rigorously in a test tank to determine the diameter of the orifice to be employed for the sieve plates of different porosity and the number of sieve plates to be installed in the diffuser sieve sets. The study has been conducted with four different patterns of sieves with different porosity to achieve a wide range of relative densities for four different uniformly graded Indian sands. The dynamic penetrometer which is considered to be one of the cost effective instruments has been efficiently used to determine the soil resistance at various locations of the test tank for every given height of fall in order to check the uniformity of placement density throughout the sand bed. The study reveals that the sand beds of different relative densities could be achieved using different patterns of diffuser sieves at optimum sand flow rate without compromising the uniformity. The effect of height of fall as well as porosity of diffuser sieves on the relative density of different sands has been studied in detail. The deposition intensity and the relative density obtained from the present study are compared with the values available in the literature.