A constitutive model for mechanical and chemo‐mechanical compaction in sedimentary basins and finite element analysis

Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd.     

A probabilistic framework for improving reverse geocoding output

Reverse geocoding, which transforms machine‐readable GPS coordinates into human‐readable location information, is widely used in a variety of location‐based services and analysis. The output quality of reverse geocoding is critical because it can greatly impact these services provided to end‐users. We argue that the output of reverse geocoding should be spatially close to and topologically correct with respect to the input coordinates, contain multiple suggestions ranked by a uniform standard, and incorporate GPS uncertainties. However, existing reverse geocoding systems often fail to fulfill these aims. To further improve the reverse geocoding process, we propose a probabilistic framework that includes: (1) a new workflow that can adapt all existing address models and unitizes distance and topology relations among retrieved reference data for candidate selections; (2) an advanced scoring mechanism that quantifies characteristics of the entire workflow and orders candidates according to their likelihood of being the best candidate; and (3) a novel algorithm that derives statistical surfaces for input GPS uncertainties and propagates such uncertainties into final output lists. The efficiency of the proposed approaches is demonstrated through comparisons to the four commercial reverse geocoding systems and through human judgments. We envision that more advanced reverse geocoding output ranking algorithms specific to different application scenarios can be built upon this work.     

On the Leontovich boundary condition in geoelectromagnetism

The paper is dedicated to the 70th anniversary of the formulation of the Leontovich impedance boundary condition. This boundary condition is applied in geophysics for sounding of the lithosphere and diagnostics of the magnetosphere by electromagnetic waves. The use of the Leontovich condition is shown to open new possibilities for investigating the anharmonicity of MHD oscillations of the magnetosphere. The differential impedance equation, which is derived with the aid of the Leontovich boundary condition, can prove useful in the investigation of the electroconductivity of the lithosphere by the method of inductive sounding. It is shown within the framework of a simple model that the method of the Leontovich parabolic equation makes it possible to find a correction to the Leontovich boundary condition caused by the fact that, in general, the relation between the horizontal components of an electromagnetic field is not local. It is shown that the anharmonicity of MHD oscillations of the magnetosphere coupled with the nonlocality of the boundary condition on the Earth’s surface can lead to an apparent nonlinearity of the surface impedance calculated in accordance with the classical technique of magnetotelluric sounding.     

Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system

Simulations of flow for a discrete-fracture model in fractured porous rocks have gradually become more practical, as a consequence of increased computer power and improved simulation and characterization techniques. Discrete-fracture models can be formulated in a lower-dimensional framework, where the fractures are modeled in a lower dimension than the matrix, or in an equi-dimensional form, where the fractures and the matrix have the same dimension.     

Quantifying soil variability in GIS applications: II Spatial distribution of soil properties

Abstract

Application of Geographical Information Systems (GIS) to environmental problems requires a spatial delineation of soil properties and some a priori knowledge of their variability. In the U.S.A. estimated soil properties are available from soil surveys published by the National Cooperative Soil Survey. However, information on the extent and nature of their variability is generally lacking. To remedy the situation, variability of selected soil properties measured at farm and watershed scales was compared with published values using geostatistical methods. Results showed that the variability of bulk density and hydraulic conductivity can be estimated by comparing regularized variograms of measured and published values. An overlay method was proposed that combines measured and published delineations of soil properties into an overlay with the properties of both. The technique makes it possible to incorporate variability into a GIS analysis using data available from published soil surveys.     

Submarine Groundwater Discharge-Derived Nutrient Loads to San Francisco Bay: Implications to Future Ecosystem Changes

Submarine groundwater discharge (SGD) was quantified at select sites in San Francisco Bay (SFB) from radium (²²³Ra and ²²⁴Ra) and radon (²²²Rn) activities measured in groundwater and surface water using simple mass balance box models. Based on these models, discharge rates in South and Central Bays were 0.3?C7.4?m³?day^?1?m^?1. Although SGD fluxes at the two regions (Central and South Bays) of SFB were of the same order of magnitude, the dissolved inorganic nitrogen (DIN) species associated with SGD were different. In the South Bay, ammonium (NH ₄ ⁺ ) concentrations in groundwater were three-fold higher than in open bay waters, and NH ₄ ⁺ was the primary DIN form discharged by SGD. At the Central Bay site, the primary DIN form in groundwater and associated discharge was nitrate (NO ₃ ^? ). The stable isotope signatures (??¹⁵N_NO3 and ??¹⁸O_NO3) of NO ₃ ^? in the South Bay groundwater and surface waters were both consistent with NO ₃ ^? derived from NH ₄ ⁺ that was isotopically enriched in ¹⁵N by NH ₄ ⁺ volatilization. Based on the calculated SGD fluxes and groundwater nutrient concentrations, nutrient fluxes associated with SGD can account for up to 16?% of DIN and 22?% of DIP in South and Central Bays. The form of DIN contributed to surface waters from SGD may impact the ratio of NO ₃ ^? to NH ₄ ⁺ available to phytoplankton with implications to bay productivity, phytoplankton species distribution, and nutrient uptake rates. This assessment of nutrient delivery via groundwater discharge in SFB may provide vital information for future bay ecological wellbeing and sensitivity to future environmental stressors.     

Noble metals in the insoluble carbonaceous substance of black shales and ores: Direct atomic emission data

The way of extraction of insoluble carbonaceous substance (ICS) from black shales and technique of determination of noble metals (NM) in it by direct arc atomic emission analysis (AEA) were developed. The first results are reported on the determination of mineral-bound noble metals in ICS concentrates, which were extracted from the rocks and ores of the Sukhoi Log, Vernyi, Vasil’evskoe, Kuranakh, Daugyztau, Bakyrchik, Kamennoe, El’dorado, and other deposits.     

Liquefaction Identification Using IRS-1D Temporal Indices Data

One of the major after effect of Bhuj Earthquake which occurred on January 26, 2001 was wide spread appearance of liquefaction of soil in the Rann of Kachchh and the coastal areas of Kandla port covering an area of more than tens of thousands of kilometers. Remote sensing data products allow us to explore the land surface parameters at different spatial scales. In this work, an attempt has been made to identify the liquefied soil area using conventional indices from IRS-1D temporal images. The same has been investigated and compared with Class Based Sensor Independent (CBSI) spectral indices, while applying fuzzy based noise classification as soft computing approach using supervised classification. Seven spectral indices have been investigated to identify liquefied soil areas using temporal multi-spectral images. The result shows that the temporal variations can be accounted by using appropriate remote sensing based spectral indices. It is found that CBSI based TNDVI using temporal data yields the best results for identification of liquefied soil areas, while CBSI based SR gives best results for water body identification.