Search Results (193)

Comprehensive land improvement causes strong disturbances of land use patterns in the short term, resulting in changes in landscape structure and function. This study adopts the moving window method and semi-variation function to explore the spatial scale effect of landscape pattern metrics in the comprehensive land consolidation project area of Baimahu Farm, and the spatial variability and homologous ecological processes. The results showed that: (1) patch density, largest patch index, area-weighted average shape index, contagion, and division index all showed obvious scale effects, and the suitable first and second scale domains in the study area are 5–7 m and 35–40 m, respectively, and 5 m is the most suitable grain size for the study of landscape pattern change. (2) The block basis ratio of the semi-variogram of the six landscape level indices begins to stabilize at the window radius of 210 m. This scale can reflect the spatial variability of the landscape pattern in the study area and is the most suitable analysis range. (3) The fragmentation degree of paddy fields as landscape matrix decreased and the landscape dominance degree increased in the comprehensive land improvement; the degree of fragmentation of irrigated land and agricultural land for facilities increased, the aggregation of land for construction increased, the dominance degree of the pond surface decreased, and the overall landscape diversity of each mosaic decreased; the landscape heterogeneity of ditches, rural roads, forest and grassland corridors was weakened, and the ecosystem service function was weakened. (4) The trend of increased fragmentation, simplification of landscape types, and decreased diversity presented by the landscape pattern clearly indicates that the landscape pattern of the study area has been seriously damaged to some extent under the influence of human activities. This damage not only has a direct negative impact on the local ecological environment, but also poses a potential threat to the sustainable development of the region. Full article

From a global ecological management perspective, as a core tree species in the mountain ecosystem of Xinjiang, the study of the spatial distribution characteristics of Picea schrenkiana var. tianschanica is crucial for maintaining the ecological balance in the Tianshan region. This study focuses on the western section of the Tianshan mountains in Xinjiang and employs the variogram analysis technique to explore the spatial heterogeneity of Picea schrenkiana var. tianschanica biomass. Successively, the study implements ordinary kriging, multivariate linear regression, the random forest algorithm, and an innovative random forest residual kriging method to conduct a spatial interpolation analysis of Picea schrenkiana var. tianschanica biomass in the target area. The results indicate that the biomass of Picea schrenkiana var. tianschanica exhibits moderate spatial autocorrelation, with its distribution pattern being influenced by a combination of topography, climate, and soil conditions. After comparing multiple spatial interpolation methods, it is found that the hybrid model combining regression analysis and kriging, delivers the best performance (R² = 0.642, RMSE = 40.18, RMSPE = 44.6). This model not only significantly improves the prediction accuracy, but also provides an intuitive and accurate spatial distribution map of Picea schrenkiana var. tianschanica biomass in the western section of the Tianshan mountains which reveals the global ecological importance of Picea schrenkiana var. tianschanica in an intuitive and accurate way, providing valuable scientific evidence and practical guidance for the field of international ecological protection and resource management. Full article

Analyzing geological profiles is of great importance for various applications such as natural resource management, environmental assessment, and mining engineering projects. This study presents a novel geostatistical approach for subsurface geological profile interpolation using a fractional kriging method enhanced by random forest regression. Using bedrock elevation data from 49 boreholes in a study area in southeast China, we first use random forest regression to predict and optimize variogram parameters. We then use the fractional kriging method to interpolate the data and analyze the variability. We also compare the proposed model with traditional methods, including linear regression, K-nearest neighbors, ordinary kriging, and fractional kriging, using cross-validation metrics. The results indicate that the proposed model reduces prediction errors and enhances spatial prediction reliability compared to other models. The MSE of the proposed model is 25% lower than that of ordinary kriging and 10% lower than that of fractional kriging. In addition, the execution time of the proposed model is slightly higher than other models. The findings suggest that the proposed model effectively captures complex subsurface spatial relationships, offering a reliable and precise solution for performing spatial interpolation tasks. Full article

Geographic patterns determine geogenic radon factors that, changing over the territory, form spatial structures of different scales associated with regional and local variations. The study of these structures is important for assessing the possibility of using limited data to predict geogenic radon potential. Our research focuses on the study of the physical properties of soils (moisture, soil density, porosity and void ratio) in the Kuznetsk coal basin. Their variations are studied using statistical methods, a variogram cloud and spatial autocorrelation of data. Soil moisture and porosity have the greatest variability in space and with depth. We conclude that the assessment of geogenic radon predictors requires consideration of the variation coefficient and autocorrelation indices at different scales. Based on the variability of humidity and the fairly homogeneous nature of the studied soils (loams), to assess the radon hazard, it is necessary to study the influence of climatic conditions, since the permeability of the environment for radon will be determined by soil moisture. With the predominance of substantially clayey soils, it is necessary to study the content of ²²⁶Ra in the upper horizons, since it is assumed that radon is predominantly diffusely transferred, in which its role is dominant. Full article

Global sensitivity analysis (GSA) of the revised universal soil loss equation (RUSLE) factors is in its infancy but is crucial to rank the importance of each factor in terms of its non-linear impact on the soil erosion rate. Hence, the goal of this study was to perform a GSA of each factor of RUSLE for a soil erosion assessment in southern Bahia, Brazil. To meet this goal, three non-linear topographic factor (LS factor) equations alternately implemented in RUSLE, coupled with geographic information system (GIS) software and a variogram analysis of the response surfaces (VARSs), were used. The results showed that the average soil erosion rate in the Pardo River basin was 25.02 t/ha/yr. In addition, the GSA analysis showed that the slope angle which is associated with the LS factor was the most sensitive parameter, followed by the cover management factor (C factor) and the support practices factor (P factor) (CP factors), the specific catchment area (SCA), the sheet erosion (m), the erodibility factor (K factor), the rill (n), and the erosivity factor (R factor). The novelty of this work is that the values of parameters m and n of the LS factor can substantially affect this factor and, thus, the soil loss estimation. Full article

There is growing international interest in assessing population exposure to radiofrequency electromagnetic fields, especially those generated by mobile-phone base stations. The work presented here is an experimental study in which we assess exposure to radiofrequency electromagnetic fields in a university environment, where there is a site with mobile-phone antennas and where a large number of people live on a daily basis. The data were collected with a personal exposure meter in two samplings, one walking at ground level and the other using an aerial vehicle at a height higher than the buildings. The geo-referenced electric-field data were subjected to a process in which a theoretical model was adjusted to the experimental variograms, and heat maps were obtained using kriging interpolation. The research carried out is of great relevance, since it provides detailed measurements of the electromagnetic radiation levels both at ground level and at significant heights, using innovative methodologies such as the use of drones. Furthermore, the results obtained allow for contextualizing the exposures in relation to international safety limits, highlighting the importance of rigorous monitoring in everyday environments. Full article

The agro-pastoral ecotone of Inner Mongolia, one of China’s most ecologically vulnerable regions, requires careful evaluation and prediction of landscape ecological risks to improve its environment and support sustainable development. Our study built a model to assess the landscape ecological risks from 1990 to 2020 using land use data from Google Earth Engine. We examined the changes in landscape ecological risks and their driving factors through spatial autocorrelation analysis and geographic detectors. Future ecological risks from 2025 to 2040 were predicted using the multi-criteria evaluation-cellular automata-Markov model. Results revealed a declining trend in both disturbance and loss intensity across land use types, with the overall ecological risk index also decreasing. Higher risk areas were concentrated in the east and southwest, while lower risks were observed in the north and center. Temperature and precipitation are key natural factors, while the impact of Gross Domestic Product (GDP), a human factor, on ecological risk is increasing and surpassed natural influences in 2015 and 2020. In the future, the highest risk areas will remain in the southwest and northeast. This study provides detailed evidence and guidance for ecological safety and sustainable development in the agro-pastoral ecotone of Inner Mongolia. Full article

Uncertainty assessment is an essential part of modeling and mapping the spatial variability of key soil properties, such as texture. The study aimed to compare sequential Gaussian simulation (SGS) and turning bands simulation (TBS) for assessing the uncertainty in unknown values of the textural fractions accounting for their compositional nature. The study area was a forest catchment (1.39 km²) with soils classified as Typic Xerumbrepts and Ultic Haploxeralf. Samples were collected at 135 locations (0.20 m depth) according to a design developed using a spatial simulated annealing algorithm. Isometric log-ratio (ilr) was used to transform the three textural fractions into a two-dimensional real vector of coordinates ilr.1 and ilr.2, then 100 realizations were simulated using SGS and TBS. The realizations obtained by SGS and TBS showed a strong similarity in reproducing the distribution of ilr.1 and ilr.2 with minimal differences in average conditional variances of all grid nodes. The variograms of ilr.1 and ilr.2 coordinates were better reproduced by the realizations obtained by TBS. Similar results in reproducing the texture data statistics by both algorithms of simulation were obtained. The maps of expected values and standard deviations of the three soil textural fractions obtained by SGS and TBS showed no notable visual differences or visual artifacts. The realizations obtained by SGS and TBS showed a strong similarity in reproducing the distribution of isometric log-ratio coordinates (ilr.1 and ilr.2). Overall, their variograms and data were better reproduced by the realizations obtained by TBS. Full article

Essential soil nutrients are dynamic in nature and require timely management in farmers’ fields. Accurate prediction of the spatial distribution of soil nutrients using a suitable sampling density is a prerequisite for improving the practical utility of spatial soil fertility maps. However, practical research is required to address the challenge of selecting an optimal sampling density that is both cost-effective and accurate for preparing digital soil nutrient maps across regional extents. This study examines the impact of sampling density on spatial prediction accuracy for a range of soil fertility parameters over a regional extent of 8303 km² located in eastern India. Surface soil samples were collected from 1024 sample points. The performance of six levels of sampling densities for spatial prediction of 14 soil properties was compared using ordinary kriging. From the sample points, randomization was used to select 224 points for validation and the remaining 800 for calibration. Goodness-of-fit for the semi-variograms was evaluated by R² of model fit. Lin’s concordance correlation coefficient (CCC) and root mean square error (RMSE) were evaluated through independent validation as spatial prediction accuracy parameters. Results show that the impact of sampling density on prediction accuracy was unique for each soil property. As a common trend, R² of model fit and CCC scores improved, and RMSE values declined with the increasing sampling density for all soil properties. On the other hand, the rate of gain in the accuracy metrics with each increment in the sampling density gradually decreased and ultimately plateaued. This indicates that there exists a sampling density threshold beyond which the extra effort on additional sampling adds less to the spatial prediction accuracy. The findings of this study provide a valuable reference for optimizing soil nutrient mapping across regional extents. Full article

Groundwater, a vital resource for providing drinking water to populations, must be managed sustainably to ensure its availability and quality. This study aims to assess the groundwater quality in the Bourgogne-Franche-Comté region (~50,000 km²) of France and identify the processes responsible for its variability. Data were extracted from the Sise-Eaux database, resulting in an initial sparse matrix comprising 8723 samples and over 100 bacteriological and physicochemical parameters. From this, a refined full matrix of 3569 samples and 22 key parameters was selected. The data underwent logarithmic transformation before applying principal component analysis (PCA) to reduce the dimensionality of the dataset. The analysis of the spatial structure, using both raw and directional variograms, revealed a categorization of parameters, grouping major ions according to the regional lithology. Bacteriological criteria (Escherichia coli and Enterococcus) displayed strong spatial variability over short distances, whereas iron (Fe) and nitrates showed intermediate spatial characteristics between bacteriology and major ions. The PCA allowed the creation of synthetic maps, with the first seven capturing 80% of the information contained in the database, effectively replacing the individual parameter maps. These synthetic maps highlighted the different processes driving the spatial variations in each quality criterion. On a regional scale, the variations in fecal contamination were found to be multifactorial, with significant influences captured by the first four principal components. The 22 parameters can be grouped into six categories based on their spatial and temporal variations, allowing for the redefinition of a resource management and monitoring strategy that is adapted to the identified spatial patterns and processes at the regional scale, while also reducing analytical costs. Full article

Groundwater flow and transport models are essential tools for assessing and quantifying the migration of organic contaminants at polluted sites. Uncertainties in the hydrodynamic and transport parameters of the aquifer have a significant effect on model predictions. Uncertainties can be quantified with advanced sensitivity methods such as Sobol’s High Dimensional Model Reduction (HDMR) and Variogram Analysis of Response Surfaces (VARS). Here we present the application of VARS and HDMR to assess the global sensitivities of the outputs of a transient groundwater flow model of the Gállego alluvial aquifer which is located downstream of the Sardas landfill in Huesca (Spain). The aquifer is subject to the tidal effects caused by the daily oscillations of the water level in the Sabiñánigo reservoir. Global sensitivities are analyzed for hydraulic heads, aquifer/reservoir fluxes, groundwater Darcy velocity, and hydraulic head calibration metrics. Input parameters include aquifer hydraulic conductivities and specific storage, aquitard vertical hydraulic conductivities, and boundary inflows and conductances. VARS, HDMR, and graphical methods agree to identify the most influential parameters, which for most of the outputs are the hydraulic conductivities of the zones closest to the landfill, the vertical hydraulic conductivity of the most permeable zones of the aquitard, and the boundary inflow coming from the landfill. The sensitivity of heads and aquifer/reservoir fluxes with respect to specific storage change with time. The aquifer/reservoir flux when the reservoir level is high shows interactions between specific storage and aquitard conductivity. VARS and HDMR parameter rankings are similar for the most influential parameters. However, there are discrepancies for the less relevant parameters. The efficiency of VARS was demonstrated by achieving stable results with a relatively small number of simulations. Full article

Due to their complex morphology, karst terrains are particularly more fragile and vulnerable to environmental damage compared to most natural systems. Their hydraulic properties, such as their transmissivity (T) and spatial variability, can be relevant for understanding groundwater flow and, consequently, for the sustainable management of water resources. The application of geostatistical methods allows for spatial interpolation and mapping based on observations combined with uncertainty quantification. Direct measurements of T are typically scarce, while those of the specific capacity (S_c) are more frequent. We established a linear and spatial relationship between the logarithms of T and S_c measured in 174 wells in a semi-arid karst region in northeastern Brazil. These relationships were used to construct a cross-variogram, whose Linear Model of Coregionalization proved valid. The values and the cross-variogram of logT and logS_c were used to generate interpolations over 2554 values of logS_c, which did not spatially coincide with logT. We used ordinary co-kriging (CO-OK) and conditional sequential Gaussian co-simulation (CO-SGS) to generate the interpolations. The cross-variogram of logT and logS_c, when considering 174 wells, was isotropic with an exponential structure, a nugget effect of approximately 20% of the sill, and a range of 5 km. Cross-validation indicated an optimal number of 10 neighboring wells used in CO-OK, and we used 500 stochastic realizations in CO-SGS, which were then used to generate maps of logT estimates, deviations derived from the interpolations, and probabilistic scenarios. The resulting transmissivity maps are relevant for the design of groundwater management strategies, including stochastic approaches where the transmissivity realizations can be used to parameterize multiple executions of numerical flow models. Full article

Geostatistical prediction methods are increasingly used in earth sciences and engineering to improve upon our knowledge of attributes in space and time. During mining activities, it is very important to have an estimate of any contamination of the soil and groundwater in the area for environmental reasons and to guide the reclamation once mining operations are finished. In this paper, we present the geostatistical analysis of the water content in certain pollutants (Cd and Mn) in a group of mines in Northern Greece. The monitoring points that were studied are 62. The aim of this work is to create a contamination prediction map that better represents the values of Cd and Mn, which is challenging based on the small sample size. The correlation between Cd and Mn concentration in the groundwater is investigated during the preliminary analysis of the data. The logarithm of the data values was used, and after removing a linear trend, the variogram parameters were estimated. In order to create the necessary maps of contamination, we employed the method of ordinary Kriging (OK) and inversed the transformations using bias correction to adjust the results for the inverse transform. Cross-validation shows promising results ($\rho =65\%$ for Cd and $\rho =52\%$ for Mn, RMSE = 25.9 ppb for Cd and RMSE = 25.1 ppm for Mn). As part of this work, the Spartan Variogram model was compared with the other models and was found to perform better for the data of Mn. Full article

This article proposes a novel methodology for estimating metallurgical copper recovery, a critical feature in mining project evaluations. The complexity of modeling this nonadditive variable using geostatistical methods due to low sampling density, strong heterotopic relationships with other measurements, and nonlinearity is highlighted. As an alternative, a copula-based conditional quantile regression method is proposed, which does not rely on linearity or additivity assumptions and can fit any statistical distribution. The proposed methodology was evaluated using geochemical log data and metallurgical testing from a simulated block model of a porphyry copper deposit. A highly heterotopic sample was prepared for copper recovery, sampled at 10% with respect to other variables. A copula-based nonparametric dependence model was constructed from the sample data using a kernel smoothing method, followed by the application of a conditional quantile regression for the estimation of copper recovery with chalcocite content as secondary variable, which turned out to be the most related. The accuracy of the method was evaluated using the remaining 90% of the data not included in the model. The new methodology was compared to cokriging placed under the same conditions, using performance metrics RMSE, MAE, MAPE, and R². The results show that the proposed methodology reproduces the spatial variability of the secondary variable without the need for a variogram model and improves all evaluation metrics compared to the geostatistical method. Full article

This paper presents a Kriging interpolation method for constructing a database of atmospheric refractivity in Korea. To collect as much data as possible for the interpolation, meteorological data from 120 regions of Korea, including both land and sea areas, are examined. Then, the normalized atmospheric refractivity N_n is calculated for an altitude of 0 m, because the refractivity tends to vary depending on the altitude of the observation site. In addition, the optimal variogram model to obtain the spatial correlation required for the Kriging method is investigated. The estimation accuracy of the Kriging interpolation is compared with that of the inverse distance weighting (IDW) and the bi-linear methods. The average and maximum estimation errors when using the Kriging method are 0.24 and 1.32, respectively. The result demonstrates that the Kriging method is more suitable for the interpolation of atmospheric refractivity in Korea than the conventional methods. Full article