Identifying spatial heterogeneity of coal resource quality in a multilayer coal deposit by multivariate geostatistics

This study presents a new geostatistical approach to characterization of the geometry and quality of a multilayer coal deposit using the data of seam thickness as a geometric property and the contents of ash, sodium, total sulphur, and the heating value as quality properties. A coal deposit in East Kalimantan (Borneo), Indonesia, which has a synclinal geological structure, was chosen as the study site. Semivariogram analysis clarified the strong dependence of heating value on ash content in the top and bottom parts of each seam and the existence of a strong correlation with sodium content over the sub-seams in the same location. The correlations between the geometry and quality of the seams were generally weak. A linear coregionalization model was used to derive the spatial correlation coefficients of two variables at each scale component from the single- and cross-semivariogram matrices. Because the data were correlated spatially in the same seam or over different seams, multivariate techniques (ordinary cokriging and factorial cokriging) were mainly used and the resultant spatial estimates were compared to those derived using a univariate technique (ordinary kriging). A factorial cokriging was effective to decompose the spatial correlation structures with different scales. Another important characteristic was that the sodium content shows distinct segregation: the low zones are concentrated near the boundary of the sedimentary basin, while the high zones are concentrated in the central part. The main component of sodium originates from the abundance of saline water. Therefore, it can be inferred that seawater had stronger effects on the coal depositional process in the central basin than in the border part. The geostatistical modeling results suggest that the thicknesses of all the major seams were controlled by the syncline structure, while the coal qualities chiefly were originated from the coal depositional and diagenetic processes.     

Optimization of area–volume–elevation curve using GIS–SRTM method for rainwater harvesting in arid areas

The major limitation in planning water harvesting is the lack of knowledge in the estimation of surface area and storage volume at any depth of dam reservoir. The area–volume–elevation (AVE) curve of a reservoir plays a key role in estimating the most suitable depth, optimum surface area and highest capacity of reservoir storage. The existing methods to estimate the AVE curve are costly and time-consuming and require laborious work. This study attempts to develop a method to optimize the AVE curve for earth dams, using the digital elevation model generated by the Shuttle Radar Topography Mission (SRTM) data, and integrate it with the geographic information system (GIS), known as the GIS–SRTM. The proposed method was tested using field data in the Western Desert of Iraq, which is an arid environment. Three constructed small earth dams were selected for this study. The AVE curves were extracted for Horan 2 (H2), Al-gara 2 (G2) and Al-gara 4 (G4) earth dams. Comprehensive analyses have been carried out to evaluate the performance of the AVE curves using the proposed GIS–SRTM method and the field data. From the comparison, the proposed GIS–SRTM method was able to produce reliable AVE curves with a relative error less than 20%. Additionally, the proposed method was less time-consuming and the AVE curves can be visualized immediately. The proposed GIS–SRTM method is relatively supportive in analyzing spatial data to select the optimal site for rainwater harvesting and prevent excessive evaporation losses.     

The rupture process during the 1999 Düzce, Turkey, earthquake from joint inversion of teleseismic and strong-motion data

Teleseismic and strong-motion data are inverted to determine the rupture process during the November 1999 Düzce earthquake in NW Turkey. The fault geometry, rise time and rupture velocity are determined from the aftershock distribution and preliminary inversions of the teleseismic data. Joint inversion of the teleseismic and strong-motion data is then carried out for the slip distribution. We obtain the strike 264°, dip 64°, rake −172°, seismic moment 5.0×10¹⁹ N m (M_w 7.1), and average stress drop 7 MPa. This earthquake was characterized by bilateral fault rupture and asymmetric slip distribution. Two asperities (areas of large slip) are identified, the eastern one being 1.5 times larger than the western one. The derived slip distribution is consistent with the aftershock distribution, surface rupture and damage. The point of rupture initiation in this Düzce earthquake coincided with the eastern tip of the aftershock distribution of the August 1999 Izmit earthquake.     

Radiocarbon age distribution of groundwater in the Konya Closed Basin, central Anatolia, Turkey

Annual abstraction of 2.6?×?10⁹ m³ of groundwater in the 53,000 km² Konya Closed Basin of central Turkey has caused a head decline of 1 m/year over the last few decades. Therefore, understanding the hydrogeology of this large endorheic basin, in a semi-arid climate, is important to sustainable resource management. For this purpose, the groundwater’s radiocarbon age distribution has been investigated along a 150-km transect parallel to regional flow. Results show that the groundwater ranges in age from Recent at the main recharge area of the Taurus Mountains in the south, to about 40,000 years around the terminal Salt Lake located in the north. In this predominantly confined flow system, radiocarbon ages increase linearly by distance from the main recharge area and are in agreement with the hydraulic ages. The mean velocity of regional groundwater flow (3 m/year) is determined by the rate of regional groundwater discharge into the Salt Lake. Calcite dissolution, dedolomitization and geogenic carbon dioxide influx appear to be the dominant geochemical processes that determine the carbon isotope composition along the regional flow path. The groundwater’s oxygen-18 content indicates more humid and cooler paleorecharge. A maximum drop of 5°C is inferred for the past recharge temperature.     

Volcanic gaps due to oblique consumption of aseismic ridges

The present-day consumption of oceanic ridges and other buoyant rises and fragments at circum-Pacific subduction zones, and presumably elsewhere, are closely related to existing gaps in volcanism. Examples are the gaps associated with the Nazca, Juan Fernandez, Cocos, Marcus-Necker and Louisville ridges. The buoyancy of these ridges breaks the continuity of the subducted plate, which may lead to reduced water supply required for melting of magma, and therefore create temporary volcanic gaps. The oblique consumption of these ridges causes the gap to migrate with time. This mechanism may be useful in interpreting time-space patterns of past volcanic chains associated with subduction in terms of the consumption of the disruptive oceanic plateaus and ridges.     

Adaptive state estimation of groundwater contaminant boundary input flux in a 2-dimensional aquifer

In many circumstances involving heat and mass transfer issues,it is considered impractical to measure the input flux and the resulting state distribution in the domain.Therefore,the need to develop techniques to provide solutions for such problems and estimate the inverse mass flux becomes imperative.Adaptive state estimator(ASE)is increasingly becoming a popular inverse estimation technique which resolves inverse problems by incorporating the semi-Markovian concept into a Bayesian estimation technique,thereby developing an inverse input and state estimator consisting of a bank of parallel adaptively weighted Kalman filters.The ASE is particularly designed for a system that encompasses independent unknowns and/or random switching of input and measurement biases.The present study describes the scheme to estimate the groundwater input contaminant flux and its transient distribution in a conjectural two-dimensional aquifer by means of ASE,which in particular is because of its unique ability to efficiently handle the process noise giving an estimation of keeping the relative error range within 10%in 2-dimensional problems.Numerical simulation results show that the proposed estimator presents decent estimation performance for both smoothly and abruptly varying input flux scenarios.Results also show that ASE enjoys a better estimation performance than its competitor,Recursive Least Square Estimator(RLSE)due to its larger error tolerance in greater process noise regimes.ASE's inherent deficiency of being slower than the RLSE,resulting from the complexity of algorithm,was also noticed.The chosen input scenarios are tested to calculate the effect of input area and both estimators show improved results with an increase in input flux area especially as sensors are moved closer to the assumed input location.     

Assessing factors affecting drought,earthquake, and flood risk perception: empirical evidence from Bangladesh

Natural Hazards - Understanding household disaster risk perception is crucial to formulate and apply disaster risk reduction strategies. Using survey data from 300 households from three highly...