Evaluating the suitability of different terrains for sustaining human settlements according to the local elevation range in China using the ASTER GDEM

The topographical suitability assessment of human settlements(SAHS) creates a solid foundation for regional population distribution and socio-economic development. Local elevation range(LER) is an important factor that can be used to assess the suitability of different terrains for sustaining human settlements. However, current digital elevation model(DEM)-based LER products suffer from some challenges typically because of their subjectively selected neighborhood scales and coarser spatial resolution. In this study, we initially determined the optimal statistical window and then calculated the appropriate LER with the finer resolution data of the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global DEM(ASTER GDEM) products for China. Then, the appropriate LER was used to evaluate the topographical SAHS and its correlations with the national gridded population distribution(1 km × 1 km) in 2010. The results show that the optimal statistical window for calculating a 1 arc-second(about 30 m) resolution GDEM LER for China can be determined using a 51 × 51 grid unit(width × height) within a rectangular neighborhood, corresponding to an area of about 2.34 km~2. Secondly, the LER values in the southern and western China were greater than those of the north and east, showing a trend which consistently reflects the general spatial features of landforms. Finally, the relationship between GDEM LER and population density was highly correlated with the R~2 value of 0.81. It showed that 85.22% of the Chinese population was located in areas where the LER is lower than 500 m. The topographically suitable area within China decreased from the southeastern coastal zone towards the northwestern inland areas due to transition from plains and basins to plateaus and mountains. The total area of moderate to high suitable level was 423.84 × 10~4 km~2, or 44.15% of the total land area, with 88.17% of the national population. Our study demonstrates the usefulness of appropriate LER in evaluating the topographical SAHS as well as its significant impact on population distribution.     

Analysis of sedimentary facies and depositional environments of Paleogene sequence of Cambay basin,India

A study on the sedimentary facies characterization and depositional environment interpretations for the K#Field (K-Oil Field) in Cambay petroleum basin of western onshore, India was conducted based on the sub-surface data from drilled wells, including well logs, borehole images, cores and the regional knowledge of the basin. In this work, an effort is made to integrate the current data from seismics and well logging, to study and analyze its depositional environments and establish the petroleum systems. The study regions for the present work are K45 and K48 blocks. The target strata includes 2 oil-bearing formations of Paleogene, which is about 3600 ft; they are M#Fm (M-Formation) of the Eocene and N#Fm (N-Formation) of Oligocene, subdivided into 11 zones. The sediment fill is mostly of Tertiary. The research attempts to decipher the oil - depositional source correlation problems of the basin. Sedimentary models were established referring to the core analysis, core photographs and well logs. Reservoir and heterogeneity study included reservoir lithology features, physical properties and pore structure features.Well facies analysis of oil well WELL-0297 and WELL-0129 was done and the results were analyzed for further drilling of new wells for oil and gas exploration. The study found that the Eocene, Oligocene, Miocene and Paleogene are fluviatile facies sand and mud interbed sediment with the thickness 2000-4000 ft, which are main oil-bearing formations in these areas. Studies concluded that the fluvial reservoirs of the K#Field are characterized by large variations from laterally extensive bodies with good interconnectedness and high net-to-gross ratios, multi-storey ribbon bodies with poor interconnections and low net-to-gross ratios.     

Felsic tuff from Rutland Island — A pyroclastic flow deposit in Miocene-sediments of Andaman-Java subduction complex

The bedded felsic tuff exposed in Rutland Island, Andaman, consists of two facies:

Resource allocation for regional earthquake risk mitigation: a case study of Tehran,Iran

This paper presents a new optimization model to help cities in seismically active developing countries decide (1) How much to spend on pre-earthquake mitigation versus waiting until after an event and paying for reconstruction or simply not rebuilding damaged buildings? (2) Which buildings to mitigate and how? and (3) Which buildings to reconstruct and how? It extends previously developed optimization models to consider the particular issues that arise in such countries. First, the model allows for the possibility that some damaged buildings will not be reconstructed immediately and keeps track of any lost building inventory. Second, buildings can be mitigated to, or when damaged, reconstructed to, any appropriate structural type and seismic design level. Finally, the model objectives include minimizing the chance of an extremely high death toll in any one earthquake and minimizing the average annual death toll across earthquakes. The model is illustrated through a case study analysis for Tehran, Iran.     

A semi-implicit method for incompressible three-phase flow in porous media

In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken.     

Genesis of the Paleoproterozoic Rare-Metal Granites of the Katugin Massif

We studied the petrography, mineralogy, and geochemistry of the Paleoproterozoic (2.06 Ga) granites of the Katugin massif (Stanovoy suture zone), which hosts the combined rare-metal Katugin deposit. Three groups of granites were distinguished: (1) biotite (Bt) and biotite–riebeckite (Bt–Rbk) granites of the western block of the massif; (2) biotite–arfvedsonite (Bt–Arf) granites of the eastern block; and (3) arfvedsonite (Arf), aegirine–arfvedsonite (Aeg–Arf), and aegirine (Aeg) granites of the eastern block. The Bt and Bt–Rbk granites of the first group are mainly metaluminous and peraluminous rocks with rather high CaO contents and the minimum F contents among the granites described here. It was suggested that the granites of this group could be derived from a source dominated by crustal rocks with a small addition of mantle materials. These granites probably crystallized from a metaluminous–peraluminous melt with elevated CaO and moderate F contents. Melts of such compositions are least favorable for the crystallization of ore minerals. The Bt–Arf granites of the second group are mainly peralkaline and show high contents of CaO and Y and low contents of Na₂O and F. A mixed mantle–crust source was proposed for the Bt–Arf granites. The initial melt of the Bt–Arf granites could have a peralkaline composition with elevated CaO content and moderate to high F content. The Arf, Aeg–Arf, and Aeg granites of the third group are enriched in ore mineral and were classified as peralkaline granites with very low CaO contents, elevated Na2O and F contents, and usually very high contents of Zr, Hf, Nb, and Ta. Based on the geochemical and isotopic data, it was supposed that the source of the granites of the third group could be derivatives of basaltic magmas produced in an OIB-type source with a minor addition of crustal material to the magma generation zone. It was suggested that the primary melt of this granite group could be a peralkaline CaO-poor and F-rich silicic melt, which is most favorable for the crystallization of ore minerals. Based on the analysis of the geochemical characteristics of the three granite groups and their relationships within the Katugin massif, a qualitative model of its formation was proposed. According to this model, the Bt and Bt–Rbk granites of the western block crystallized first, followed by the Bt–Arf granites of the eastern block and, eventually, the Arf, Aeg–Arf, and Aeg granites enriched in ore minerals.     

Deep seismic studies in the central Karelian Craton

On the basis of analysis and generalization of the extensive seismic information obtained over the last forty years, the modern concept of the deep crustal structure in the central Karelian Craton is developed. The tomographic models are reconstructed, the stratification of the Earth’s crust and geometry of seismic boundaries are refined, and the subsidence of the boundary between the crust and mantle down to a depth of 60 km is noted. The comparison and integrated interpretation of the sections studied by the common depth point (CDP), deep seismic sounding (DSS), and converted-wave methods together with the available geological data shows that the results derived by different seismic methods differ markedly. However, since various methods have revealed the different characteristics of the medium, the entire body of seismic evidence can be used for geological interpretation.     

Development of a GIS-based catastrophe theory model (modified DRASTIC model) for groundwater vulnerability assessment

This study developed a new paradigm for groundwater vulnerability assessment by modifying the standard DRASTIC index (DI) model based on catastrophe theory. The developed paradigm was called the catastrophe theory-based DI (CDI) model. The proposed model was applied to assess groundwater vulnerability to pollution index (GVPI) in Perak Province, Malaysia. The area vulnerability index was modeled by considering the DRASTIC multiple vulnerability causative factors (VCFs) obtained from different data sources. The weights and ranking of the VCFs were computed by using the inner fuzzy membership mechanism of the CDI model. The estimated vulnerability index values of the CDI model were processed in a geographic information system (GIS) environment to produce a catastrophe theory–DRASTIC groundwater vulnerability to pollution index (CDGVPI) map, which demarcated the area into five vulnerability zones. The produced CDGVPI map was validated by applying the water quality status–vulnerability zone relationship (WVR) approach and the relative operating characteristic (ROC) curve method. The performance of the developed CDI model was compared with that of the standard DI model. The validation results of the WVR approach exhibits 89.29% prediction accuracy for the CDI model compared with 75% for the DI model. Meanwhile, the ROC validation results for the CDI and DI models are 88.8% and 78%, respectively. The GIS-based CDI model demonstrated better performance than the DI model. The GVPI maps produced in this study can be used for precise decision making process in environmental planning and groundwater management.     

Partial derivatives of surface wave phase velocities for flat anisotropic models

Summary. Surface wave behaviour in flat anisotropic structures is first illustrated by performing an exact computation on a simple two-layer model. The variational procedure of Smith & Dahlen is then used to compute the partial derivatives of surface wave phase velocities with respect to the elastic parameters in more realistic earth models. Linear relationships between the partial derivatives for a general anisotropic structure and those for a transversely isotropic structure are derived. When considering waves propagating in a fixed direction, there are only four independent derivatives for Rayleigh waves, and two for Love waves. To avoid the lack of resolution in an inverse method, we propose to use physically constrained models. These results are illustrated by using a model with hexagonal symmetry and a symmetry axis oriented either vertically or horizontally. Quasi-Love- and quasi-Rayleigh-wave partial derivatives are computed for both axis orientations. Modes up to the second overtone and periods ranging between 45 and 130 s have been considered. Finally, anomalies of phase velocity are computed in an oceanic model made of 1/6 oriented olivine crystals with horizontal or vertical preferred orientations of the a -axis.     

Raman spectra of methane hydrate up to 86 GPa

High-pressure Raman studies of methane hydrate were performed using a diamond anvil cell in the pressure range of 0.1–86 GPa at room temperature. Raman spectra of the methane molecules revealed that new softened intramolecular vibration mode of ν ₁ appeared at 17 GPa and that the splitting of vibration mode of ν ₃ occurred at 15 GPa. The appearance of these two modes indicates that an intermolecular attractive interaction increases between the methane molecules and the host water molecules and between the neighboring methane molecules. These interactions might result in the exceptional stability of a high-pressure structure, a filled ice Ih structure (FIIhS) for methane hydrate, up to 40 GPa. At 40 GPa, a clear change in the slope of the Raman shift versus pressure occurred, and above 40 GPa the Raman shift of the vibration modes increased monotonously up to 86 GPa. A previous XRD study showed that the FIIhS transformed into another new high-pressure structure at 40 GPa. The change in the Raman spectra at 40 GPa may be induced by the transition of the structure.