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.     

An autochthonous Avalonian basement source for the latest Ordovician Brenton Pluton in the Meguma terrane of Nova Scotia: U–Pb–Hf isotopic constraints and paleogeographic implications

International Journal of Earth Sciences - The Ediacaran–Ordovician Meguma Supergroup was thrust over Avalonia basement prior to the intrusion of post-Acadian, ca. 370 Ma, S-type...     

Atmospheric Correction Over Coastal Turbid Waters of Bay of Bengal Using OCEANSAT-I Ocean Colour Monitor (OCM) Data

Space-borne ocean-colour remote sensor-detected radiance is heavily contaminated by solar radiation backscattered by the atmospheric air molecules and aerosols. Hence, the first step in ocean-colour data processing is the removal of this atmospheric contribution from the sensor-detected radiance to enable detection of optically active oceanic constituents e.g. chlorophyll-a, suspended sediment etc. In standard atmospheric correction procedure for OCEANSAT-1 Ocean Colour Monitor (OCM) data, NIR bands centered at 765 and 865 nm wavelengths were used for aerosol characterization. Due to high absorption by water molecules, ocean surface in these two wavelengths acts as dark background, therefore, sensor detected radiance can be assumed to have major contribution from atmospheric scattering. For coastal turbid waters this assumption of dark surface fails due to the presence of highly scattering sediments which causes sufficient water-leaving radiance in NIR bands and lead to over-estimation of aerosol radiance resulting in negative water leaving radiance for λ < 700 nm. In the present study, for the turbid coastal waters in the northern Bay of Bengal, the concept of spatial homogeneity of aerosol and water leaving reflectance has been applied to perform atmospheric correction of OCAEANSAT-1 OCM data. The results of the turbid water atmospheric correction have also been validated using in-situ measured water-leaving radiance. Comparison of satellite derived water-leaving radiance for five coastal stations with in-situ measured radiance spectra, indicates an improvement over the standard atmospheric correction algorithm giving physically realistic and positive values. Root Mean Square Error (RMSE) between the in-situ measured and satellite derived water leaving radiance for wavelengths 412 nm, 443 nm, 490 nm, 512 nm and 555 nm was found to be 1.11, 0.718, 0.575, 0.611 and 0.651%, respectively, using standard atmospheric correction procedure. By the use of spatial homogeneity concept, this error was reduced to 0.125, 0.173, 0.176, 0.225, and 0.290 and the correlation coefficient arrived at 0.945, which is an improvement over the standard atmospheric correction procedure.     

A review of the NEMURO and NEMURO.FISH models and their application to marine ecosystem investigations

The evolution of the North Pacific Ecosystem Model for Understanding Regional Oceanography (NEMURO) family of models to study marine ecosystems is reviewed. Applications throughout the North Pacific have shown the models to be robust and to be able to reproduce 1D, 2D and 3D components of nutrient, carbon cycle and biogeochemical cycles as well as aspects of the lower trophic levels ecosystem (phyto- and zooplankton). NEMURO For Including Saury and Herring, an extension that includes higher trophic levels, can be run uncoupled or coupled to NEMURO. In the uncoupled mode, the growth and weight of an individual fish is computed using plankton densities simulated by NEMURO but with no feedback between fish consumption and plankton mortality. In the coupled mode, the feeding, growth and weight of a representative fish are computed, and prey removals due to feeding by fish appear as mortality terms on the prey. The NEMURO family of models continues to evolve, including effects of the microbial loop and iron limitation at lower trophic levels, and full life cycle, multi-species and multi-generational simulations at higher trophic levels. We outline perspectives for future end-to-end modeling efforts that can be used to study marine ecosystems in response to global environmental change.     

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.     

Propagation of plane waves in an unbounded elastic medium with Cauchy's initial stress

Summary Cauchy theory of initial stress has been applied to investigate the influence of uniform initial stress on the propagation of plane waves in an unbounded elastic medium. It is found that if the initial principal stresses are unequal, the velocity depends on the direction in which the wave propagates, whereas when the initial principal stresses are equal, the velocity is independent of the direction of wave propagation.     

Time and frequency domains response analyses of April 2015 Greece’s earthquake in the Nile Delta based on GNSS-PPP

Earthquake response characteristic and sensitivity of 1-Hz Global Navigation Satellite System - precise point positioning (GNSS-PPP) for a seismic response in time and frequency domains are the main objectives of this study. The Delta and Cairo areas, Egypt, experiencing the shock of the Greece earthquake occurring on April 16, 2015 (18:50 hours), is presented in this study. The results of this study reveal that the seismic wave effect is very small and can be neglected inside the Delta region and the strong motion occurred in an upward direction for the whole monitoring area. Furthermore, the time-frequency analysis illustrates the ability of the time and frequency domains to use and analyze the surface motion based on a 1-Hz sampling frequency and to detect the small- and high-strength motions of seismic waves.     

Comparing data-driven landslide susceptibility models based on participatory landslide inventory mapping in Purwosari area,Yogyakarta, Java

There are different approaches and techniques for landslide susceptibility mapping. However, no agreement has been reached in both the procedure and the use of specific controlling factors employed in the landslide susceptibility mapping. Each model has its own assumption, and the result may differ from place to place. Different landslide controlling factors and the completeness of landslide inventory may also affect the different result. Incomplete landslide inventory may produce significance error in the interpretation of the relationship between landslide and controlling factor. Comparing landslide susceptibility models using complete inventory is essential in order to identify the most realistic landslide susceptibility approach applied typically in the tropical region Indonesia. Purwosari area, Java, which has total 182 landslides occurred from 1979 to 2011, was selected as study area to evaluate three data-driven landslide susceptibility models, i.e., weight of evidence, logistic regression, and artificial neural network. Landslide in the study area is usually affected by rainfall and anthropogenic activities. The landslide typology consists of shallow translational and rotational slide. The elevation, slope, aspect, plan curvature, profile curvature, stream power index, topographic wetness index, distance to river, land use, and distance to road were selected as landslide controlling factors for the analysis. Considering the accuracy and the precision evaluations, the weight of evidence represents considerably the most realistic prediction capacities (79%) when comparing with the logistic regression (72%) and artificial neural network (71%). The linear model shows more powerful result than the nonlinear models because it fits to the area where complete landslide inventory is available, the landscape is not varied, and the occurence of landslide is evenly distributed to the class of controlling factor.     

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.