An insight into the surface velocity of Inylchek Glacier and its effect on Lake Merzbacher during 2006–2016 with Landsat time-series imagery

Mountain glacier is one of the extremely sensitive indicators for climate change, and its surface motion distribution and corresponding variation are valuable information for understanding ice mass exchange and glacier dynamics. This paper presents the long-term ice velocity distributions of Inylchek Glacier in the Tianshan region by pixel-tracking algorithm with time-series Landsat imagery acquired during 2006–2016. Then the monitored ice motion fields of Inylchek Glacier were carefully analyzed and revealed a generally similar spatial distribution characteristic. Most of the ice of the North Inylchek Glacier remains in a stagnant state except for the upstream part, but a relatively high velocity of 20–40 cm/day with an RMSE of 3 cm/day was observed on most part of the South Inylchek Glacier, except for the slow-moving glacier terminus. We also state the glacier dynamics around Lake Merzbacher and their possible effect on its glacier lake outburst flood (GLOF) risk. Besides, the surface velocity distribution on South Inylchek Glacier surface during the ablation period from 2014 to 2016 was also established and also compared with annual velocity. The corresponding difference yields that there is a positive relation between ice motion and temperature variation. Therefore, the time-series ice surface motion yielded by the Landsat imagery thus could provide us an efficient and low-cost way to analyze the current state and changes in glaciers, thanks to the continuous and regular spaceborne observations provided by the Landsat satellites.     

Hypsometric analysis of the micro-watersheds with different management practices located on Shivalik foothills

In this paper, the hypsometric curves and integrals of four neighboring micro-watersheds draining into Patiala-Ki-Rao stream which is situated in the Shivalik foothills of district SAS Nagar in the Punjab state (India) has been presented to access and compare the erosion regimes under different management practices. Area-elevation ratio method has been used to compute hypsometric curve and integral values for each micro-watershed through ArcGIS 10.3 and Microsoft Excel. The 9-year data of runoff and sediment yield for all these micro-watersheds under different management practices has been analyzed for their effect on land cover and soil quality. Thus, the results of present study are very useful for comparing, planning, implementing, and controlling soil erosion in similar watersheds.     

Evolution of the shape parameter in the Weibull distribution for brittle rocks under uniaxial compression

Applying the statistical damage theory based on the Weibull distribution to describe rock deformation and failure processes is an important development in rock mechanics. The shape parameter of the Weibull distribution, m, determines the basic shape of the distribution curve; additionally, it also represents a physical characteristic which can be applied when constructing rock constitutive models. To study the evolution of m during rock failure when applying the Weibull distribution to rock mechanics, uniaxial cyclic loading tests of shale specimens were conducted and previous rock mechanics experiments under different temperatures and loading rates were reviewed. The results indicate that m varied with the accumulation of damage but was almost constant between the volume expansion point and the peak strength point of each specimen. Combined with previous studies about the accelerated failure behavior of rocks, we conclude that between the volume expansion point and the peak strength point, the mechanical behavior of the rock fracture process did not change significantly. Based on the characteristics of m at different damage stages during the rock failure process, ranges of m values at different damage stages are proposed. The conclusions reached in this study may be used as an important reference for theoretical research on rock mechanics.     

ISPH numerical simulation of tsunami generation by submarine landslides

The present investigation focused on the numerical simulation of the gravity currents of non-Newtonian fluids by means of the incompressible smoothed particle hydrodynamics (ISPH) approach. Firstly, to solve the rheological properties of the non-Newtonian Bingham model by ISPH method, the multi-viscodensity approach has been introduced. Then, this methodology has been used to simulate the annular viscometer and landslide deformation test cases. Through simulating these test cases, the viscoplastic behavior of the non-Newtonian viscosity and propagation of tsunami waves due to underwater landslide movement have been observed. Numerical results were in good agreement with the theoretical and experimental studies and showed that this methodology can be used to investigate precisely the generation and propagation of tsunami surface waves.     

Geodetic observation crucial to sea level monitoring

The relative rate of rise of the sea levels measured by a tide gauge is made of a sea and a land component. The first is usually restricted to the global short-term effect of melting icecaps and expansion of water mass due to global temperature change. The second is often limited to the regional long-term effects of glacial isostatic adjustment (GIA). Sometimes, the regional subsidence, due to compaction and ground water withdrawal, is considered. Here we show as this assumption of regional subsidence fails to represent the relative sea level patterns of Sandy Hook, NJ, and The Battery, NY, as well as of Venezia Punta Della Salute, Venezia II, Trieste and Trieste II. The subsidence of the tide gauge instrument may only be addressed by the precise monitoring of the tide gauge vs. a Global Navigation Satellite System (GNSS) antenna, even if the GNSS tracking is only recent and not yet very accurate. The relative sea level records are much more complicated than what is thought.     

Assessment of spatial distribution of contaminants and their levels in soil and water resources of Calabar,Nigeria using geophysical and geological data

The contamination levels of soils and water resources in Calabar, Nigeria have been investigated using resistivity (vertical electrical sounding and electrical resistivity tomography), geochemical analyses of soil and water resources and textural analysis. Sixty randomly sited VES sites were investigated in two seasons while ERT investigations were performed along four profiles. The geochemical investigations were spread across seasons in order to track seasonal changes in physico-chemical parameters: hydrogen ion concentration (pH), electrical conductivity, total dissolved solids, chloride ion (Cl^?), nitrate ion (\( {\text{NO}}_{ 3}^{ - } \)), bicarbonate (\( {\text{HCO}}_{ 3}^{ - } \)), sulphate ion (\( {\text{SO}}_{ 4}^{2 - } \)), calcium ion (Ca²⁺), sodium ion (Na⁺), potassium ion (K⁺) and magnesium ion (Mg²⁺). Additionally, concentrations of ammonium, aluminium and nitrite ions in soils were determined. Results show that ionic concentrations in the sand-dominated soils and water are within permissible limits and baseline standards. The resistivities follow known trends in the area. However, at the central waste disposal site, a localised thin (< 5 m), low resistivity (< 15 Ωm) anomaly suspected to be due to contamination by leachates was observed. Comparatively, the contaminated area is also characterised by marginal increase in ionic concentrations. Strong attenuation capacities of overlying and adjoining clay/lateritic sediments and optimal design of the waste dump site probably reduced the spread of contaminants. The contaminated zone need to be closely monitored so that it does not extend to the aquifers. Hence, all strategies presently being used in managing wastes in Calabar should be sustained.     

A micro–macro investigation of the capillary strengthening effect in wet granular materials

Wet granular materials are three-dimensionally simulated by the discrete element method with water bridges incorporated between particles. The water bridges are simplified as toroidal shapes, and the matric suction is constantly maintained in the material. A comparison with experimental tests in the literature indicates that the toroidal shape approximation may be one of the best choices with high practicability and decent accuracy. Mechanical behaviours of wet granular materials are studied by triaxial tests. Effects of particle size distributions and void ratios are investigated systematically in this study. The hydraulic limit of the pendular state is also discussed. It gives the capillary cohesion function which is not only determined by the degree of saturation but also positively correlated to relative density and particle size polydispersity and inversely proportional to mean particle size. Furthermore, the capillary strengthening effect is also analysed microscopically in aid of the Stress–Force–Fabric relationship, mainly in fabric anisotropy, coordination number and stress transmission pattern, which revealed the micro-mechanisms of the additional effective stress induced by capillary effect.     

GaRiRO: Gradient and residual integrated rank ordering of stations in rainfall monitoring network

This paper presents a new method that integrates gradient and residual values for rank ordering of stations in a monitoring network (GaRiRO). The innovation is derived from the fact that the parameter (dependent variable) gauged through the monitoring network is modelled using independent variables that influence its measured quantity. And the dependent variable exhibit non-stationary spatial gradient with respect to the independent variables, particularly in complex terrain. GaRiRO technique was developed to prioritize the rain gauge stations for optimizing the existing network and selection of the best locations for relocation or installation of gauges. Although initially aimed to assist hydrologists with a ranking scheme for rain gauge stations, it can be applied to any environmental, meteorological or hydro-meteorological monitoring network. The new procedure is based on deriving gradient and residual value at each station by modeling the spatial relationship of dependent-independent variable using geographically weighted regression (GWR) technique. For the prospective stations with no record, the gradient value is estimated using GWR model and the residual value is derived from the residual map generated by applying kriging technique on the residual derived at all gauged locations. The method combines the decision factor with analytical strength of GIS for prioritizing the stations which results in limited number of trials for installation or relocation of gauges to yield optimized network configuration.     

Towards a Digital Earth: using archetypes to enable knowledge interoperability within geo-observational sensor systems design

Earth System Science (ESS) observational data are often inadequately semantically enriched by geo-observational information systems to capture the true meaning of the associated data sets. Data models underpinning these information systems are often too rigid in their data representation to allow for the ever-changing and evolving nature of ESS domain concepts. This impoverished approach to observational data representation reduces the ability of multi-disciplinary practitioners to share information in a computable way. Object oriented techniques that are typically employed to model data in a complex domain (with evolving domain concepts) can unnecessarily exclude domain specialists from the design process, invariably leading to a mismatch between the needs of the domain specialists, and how the concepts are modelled. In many cases, an over simplification of the domain concept is captured by the computer scientist. This paper proposes that two-level modelling methodologies developed by health informaticians to tackle problems of domain specific use-case knowledge modelling can be re-used within ESS informatics. A translational approach to enable a two-level modelling process within geo-observational sensor systems design is described. We show how the Open Geospatial Consortium’s (OGC) Observations & Measurements (O&M) standard can act as a pragmatic solution for a stable reference-model (necessary for two-level modelling), and upon which more volatile domain specific concepts can be defined and managed using archetypes. A rudimentary use-case is presented, followed by a worked example showing the implementation methodology and considerations leading to an O&M based, two-level modelling design approach, to realise semantically rich and interoperable Earth System Science based geo-observational sensor systems.     

Using landscape indicators and Analytic Hierarchy Process (AHP) to determine the optimum spatial scale of urban land use patterns in Wuhan,China

Quantifying land use patterns and functions is critical for modeling urban ecological processes, and an emerging challenge is to apply models at multiple spatial scales. Methods of determining the optimum scale of land use patterns are commonly considered using landscape metrics. Landscape metrics are quantitative indicators for analyzing landscape heterogeneity at the landscape level. In this study, due to their widespread use in urban landscape analyses and well-documented effectiveness in quantifying landscape patterns, landscape metrics that represent dominance, shape, fragmentation and connectivity were selected. Five metrics include Patch Density, Contagion, Landscape Shape Index, Aggregation Index and Connectivity. Despite a wide application of landscape metrics for land use studies, the majority mainly focuses on the qualitative analysis of the characteristics of landscape metrics. The previous models are limited in exploring the optimum scale of land use patterns for their lack of quantitation. Therefore, taking the City of Wuhan as an example, the land use unit was treated as a patch, and the landscape pattern metrics at different spatial scales were calculated and compared so as to find the optimum one. Furthermore, a mathematical model of landscape metrics was proposed to quantify the scale effect of urban land use patterns, generating a complementary tool to select the optimum scale. In addition, Analytic Hierarchy Process (AHP) was introduced to determine the respective weights of the chosen landscape metrics in this model. Fractal dimension was ultimately applied to verify the chosen optimum scale of our study region. The results indicated that 60 m is confirmed to be the optimum scale for capturing the spatial variability of land use patterns in this study area.