Hydrochemical and Multivariate Statistical Evaluation of Heavy Metals in Shallow Alluvial Aquifers of North Brahmaputra Plain,India

The aim of this study was to display distribution and relationships of heavy elements in the unconfined, shallow alluvial aquifers of the lower Jia Bharali catchment and adjoining areas in central part of North Brahmaputra Plain (NBP), India using hydrochemical as well as multivariate statistical techniques such as principal component analysis (PCA) and cluster analysis. The original matrix was made up of 10 trace elements (As, Cd, Cu, Co, Cr, Fe, Mn, Pb, Ni and Zn) estimated from 50 shallow alluvial dug wells in both the wet and the dry season for a duration of 3 hydrological years (2008–2011). Except As, Cu and Zn all the other toxic metals in the shallow aquifers were found exceeding the WHO maximum permissible limits for drinking water. PCA extracted five varimax factors as geogenic, agricultural and anthropogenic explaining about 71.2% of the total variance in the wet season and 69.3% total variance in the dry season. Hierarchical cluster analysis classified the dug wells into two groups in the wet season and three groups in the dry season with respect to the heavy elements. The results emphasized the need for routine monitoring and management in order to avoid contamination of groundwater sources in the NBP with respect to the dissolved trace elements.     

Development and verification of the comprehensive model for physical properties of hydrate sediment

Natural gas hydrate is widely distributed all over the world and may be a potential resource in the near future, whereas hydrate dissociation during the development affects wellbore stability and drilling safety. However, the present modeling of hydrate reservoir parameters ignored the influence of effective stress and only considered the hydrate saturation. In this paper, a series of stress sensitivity experiments for the unconsolidated sandstone were carried out, and the influence of mean effective stress on physical parameters was obtained; a comprehensive model for the physical parameters of hydrate reservoir was developed subsequently. With the help of ABAQUS finite element software, the established comprehensive model was verified by the use of the wellbore stability numerical model of hydrate reservoir. The verification results show that ignoring the effect of mean effective stress on the parameters of hydrate formation aggravates the invasion of drilling fluid into the hydrate formation. Besides, ignoring the stress sensitivity of reservoir physical parameters will underestimate the wellbore instability during hydrate drilling, which will be a threat to the safety of gas hydrate drilling. At the end of the drilling operation, the maximum plastic strain of the model for considering and not considering stress sensitivity was 0.0145 and 0.0138, respectively. Therefore, the established comprehensive model will provide a theoretical support for accurately predicting the engineering geological disasters in hydrate development process.     

Assessing desertification by using soil indices

Desertification generally refers to land degradation in arid, semiarid, and dry semi-humid climatic zones. It involves five principal processes: vegetation degradation, water erosion, wind erosion, salinization and waterlogging, and soil crusting and compaction. The aim of this study is assessing desertification using soil criteria. For this purpose, nine indices including sodium absorption ratio (SAR), soil gypsum percentage, soil texture, the content of HCO₃ ^?, the percentage of the organic matter, electrical conductivity (EC), pH, the content of the soil sodium, and chloride were used. The soil samples were taken in the north of Zayandeh-Rood River in Isfahan province of Iran, using soil data randomly sampled in a depth of 0–20 cm. After assessing the normality of the samples using Kolmogorov-Smirnov test, indices were imported into GIS environment and interpolated with IDW and normal and discrete kriging methods for delineating soil characteristics maps based on MEDALUS model. In this model, the data were firstly changed from 100 to 200. Thus 100 and 200 are estimated as the best and worst quality, respectively. Then the final map of soil criteria has been created by geometric mean of its indicators. The results showed that the maximum area is related to the medium class of desertification and is equal to 44,746 ha. The areas of severe and very severe classes of desertification are equal to 30,949 and 351 ha, respectively. The results also revealed that the indices of the organic matter and soil gypsum percentage are the most influential indices which affect desertification phenomenon.     

A fine-granularity scheduling algorithm for parallel XDraw viewshed analysis

Viewshed analysis is widely used in many terrain applications such as siting problem, path planning problem, and etc. But viewshed computation is very time-consuming, in particular for applications with large-scale terrain data. Parallel computing as a mainstream technique with the tremendous potential has been introduced to enhance the computation performance of viewshed analysis. This paper presents a revised parallel viewshed computation approach based on the existing serial XDraw algorithm in a distributed parallel computing environment. A layered data-dependent model for processing data dependency in the XDraw algorithm is built to explore scheduling strategy so that a fine-granularity scheduling strategy on the process-level and thread-level parallel computing model can be accepted to improve the efficiency of the viewshed computation. And a parallel computing algorithm, XDraw-L, is designed and implemented taken into account this scheduling strategy. The experimental results demonstrate a distinct improvement of computation performance of the XDraw-L algorithm in this paper compared with the coarse-partition algorithm, like XDraw-E which is presented by Song et al. (Earth Sci Inf 10(5):511–523, 2016), and XDraw-B that is the basic algorithm of serial XDraw. Our fine-granularity scheduling algorithm can greatly improve the scheduling performance of the grid cells between the layers within a triangle region.     

Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Assessment of the Pacific decadal oscillation’s contribution to the occurrence of local torrential rainfall in north China

On 21–22 July 2012, torrential rains hit North China, with the daily precipitation record at Beijing station reaching 160.6 mm; this event is named the Beijing 7–21 case. This paper assesses the likelihood of the occurrence of local torrential rains, such as the Beijing 7–21 case, from the perspective of climate variability. In particular, the influence of the Pacific Decadal Oscillation (PDO) is assessed. There were five extreme events, with daily precipitation records equal to or larger than 160.6 mm, at Beijing station during the period 1951–2012; all of these events happened during negative phases of the PDO. The present analysis indicates that precipitation events more extreme than the Beijing 7–21 case should happen more than once per decade during negative phases of the PDO, but only about once every four decades during positive PDO phases. The negative phase of the PDO is found to be associated with a much greater probability of daily records of southerly winds in North China during summer. Strong southerly summer monsoons are deemed favorable for increasing the occurrence of local extreme rainfall over North China.     

Crossing Filaments

Solar filaments show the position of large-scale polarity-inversion lines and are used for the reconstruction of large-scale solar magnetic field structure on the basis of Hα synoptic charts for the periods that magnetographic measurements are not available. Sometimes crossing filaments are seen in Hα filtergrams. We analyze daily Hα filtergrams from the archive of Big Bear Solar Observatory for the period of 1999 – 2003 to find crossing and interacting filaments. A number of examples are presented and filament patterns are compared with photospheric magnetic field distributions. We have found that all crossing filaments reveal quadrupolar magnetic configurations of the photospheric field and presume the presence of null points in the corona.     

Mapping the thickness of sediments in the Ljubljana Moor basin (Slovenia) using microtremors

The Ljubljana Moor basin is characterized by moderate bedrock topography and thicknesses of Quaternary lacustrine and fluvial sediments ranging from 0 to 200 m. More than 65 boreholes which reached the bedrock were drilled in the area, but their distribution in the basin is very uneven and some data from the boreholes uncertain. There are also no data on S-velocity distribution within the basin, but seismic refraction measurements pointed out a rather uniform increase of P-velocity with depth, great impedance contrast with the bedrock and relatively small lateral velocity variations. The microtremor horizontal-to-vertical spectral ratio (HVSR) method was therefore applied as a complementary tool to seismic refraction survey to map the thickness of sediments. First, microtremors were measured at the locations of boreholes which reached the bedrock and the resonance frequencies determined. The inverse power relationship between the resonance frequency and the thickness of sediments was then determined from 53 data pairs. The quality of the correlation is moderate due to possible heterogeneities in sediments and possible 3D effects in some minor areas, but the obtained parameters correspond well to the values obtained in six other European basins. Secondly, a 16 km-long discontinuous seismic refraction profile was measured across the whole basin, leaving uncovered some larger segments where active seismic measurements were not possible. Microtremors were then measured at 64 locations along the same profile, using 250 m point spacing, without leaving any gaps. The frequency–thickness relationship was used to invert resonance frequencies to depths. These were first validated using the results of the seismic refraction survey, which showed good agreement, and finally used for interpolation in the segments of missing refraction data to obtain a continuous depth profile of the bedrock. The study has shown that the microtremor method can be used as a complementary tool for mapping the thickness of unconsolidated sediments also in areas characterized by moderate bedrock topography. As the input data are always to some extent uncertain, it is important to have a sufficiently large number of borehole data to establish a frequency–thickness relationship, as well as some additional independent geophysical information for its validation.     

Correlation between molecular absorption spectral slope ratios and fluorescence humification indices in characterizing CDOM

Ultraviolet–visible absorption spectral slope ratios S_R (slope in 275–295 nm divided by slope in 350–400 nm) and humification index (HIX, integrated fluorescence emission in 435–480 nm divided by that in 435–480 and 300–345 nm) were compared when characterizing chromophoric dissolved organic matter (CDOM) in three humic acids and 44 whole water samples. HIX increased with increasing pH for humic acids, while their S_R showed much more complicated dependencies on pH. There was a negative correlation between S_R and HIX. S_R increased in the order terrestrial coal/peat < terrestrial soil/river < seawater, while HIX increased in the order seawater < terrestrial soil/river < terrestrial coal/peat. The comparative study in this work indicates that terrestrially derived CDOM has higher HIX and lower S_R than marine CDOM. Investigators may potentially use these two indices to compare qualitatively the character of CDOM in different sources (e.g., terrestrial vs. marine).     

Position-domain integrity risk-based ambiguity validation for the integer bootstrap estimator

Integrity monitoring for ambiguity resolution is of significance for utilizing the high-precision carrier phase differential positioning for safety–critical navigational applications. The integer bootstrap estimator can provide an analytical probability density function, which enables the precise evaluation of the integrity risk for ambiguity validation. In order to monitor the effect of unknown ambiguity bias on the integer bootstrap estimator, the position-domain integrity risk of the integer bootstrapped baseline is evaluated under the complete failure modes by using the worst-case protection principle. Furthermore, a partial ambiguity resolution method is developed in order to satisfy the predefined integrity risk requirement. Static and kinematic experiments are carried out to test the proposed method by comparing with the traditional ratio test method and the protection level-based method. The static experimental result has shown that the proposed method can achieve a significant global availability improvement by 51% at most. The kinematic result reveals that the proposed method obtains the best balance between the positioning accuracy and the continuity performance.