The Scientific World Journal

Volume 2013 (2013), Article ID 259248, 9 pages

http://dx.doi.org/10.1155/2013/259248

## Understanding the Complexity of Temperature Dynamics in Xinjiang, China, from Multitemporal Scale and Spatial Perspectives

^{1}The Research Center for East-West Cooperation in China, The Key Laboratory of Geographic Information Science, Ministry of Education of PRC, East China Normal University, Shanghai 200241, China^{2}State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China

Received 14 April 2013; Accepted 28 May 2013

Academic Editors: H.-J. Wang and J. P. Wigneron

Copyright © 2013 Jianhua Xu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Based on the observed data from 51 meteorological stations during the period from 1958 to 2012 in Xinjiang, China, we investigated the complexity of temperature dynamics from the temporal and spatial perspectives by using a comprehensive approach including the correlation dimension (CD), classical statistics, and geostatistics. The main conclusions are as follows (1) The integer CD values indicate that the temperature dynamics are a complex and chaotic system, which is sensitive to the initial conditions. (2) The complexity of temperature dynamics decreases along with the increase of temporal scale. To describe the temperature dynamics, at least 3 independent variables are needed at daily scale, whereas at least 2 independent variables are needed at monthly, seasonal, and annual scales. (3) The spatial patterns of CD values at different temporal scales indicate that the complex temperature dynamics are derived from the complex landform.

#### 1. Introduction

The Earth’s climate system is a complex, interconnected system formed by the atmosphere, the oceans, and other bodies of water, land surface, snow, and ice cover together with all living organisms and linked by flows of energy and matter. To discover the complexity of climate change process, many concepts and methods, such as entropy, fractal, nonlinearity, chaos, wavelet, and artificial neural network, have been used by scholars [1–6]. The climate models used in the research on climate change have become more complex and are today believed to be able to provide fairly reliable predictions of future temperature ranges and climate developments. However, the models still do not cover a complete set of possible mechanisms and they include considerable uncertainties. Specifically, applying global scale simulation results to interpreting and predicting regional situations is challenging, and in fact its applicability is questionable [7].

Studies have suggested that the climatic process is a chaotic dynamic system, with nonlinearity as its basic characteristic; nevertheless, there are still many open questions on the complex system [3, 5, 8, 9] such as how to understand the complexity of spatial and temporal scales of the regional climatic system; the question has no satisfactory answer.

In the last 20 years, many studies have been conducted to evaluate climatic change in the arid and semiarid regions in northwestern China [10–14]. Some studies reached a conclusion that there was a visible climatic transition in the past half-century [15–18]. This transition was characterized by a temperature increase trend. However, the conclusion brought a question whether the increase trend is a regional response to global warming or merely a rising stage in the periodic dynamic process [19].

In order to understand the complexity of climatic dynamics in Xinjiang, China, based on observed data at 51 meteorological stations during the period from 1958 to 2012, this study investigated the temperature dynamics from multiple temporal scale and spatial perspectives by using a comprehensive approach including the correlation dimension (CD), classical statistics, and geostatistics.

#### 2. Study Area and Data

##### 2.1. Study Area

Located in the northwest of China, Xinjiang is a typical semiarid and arid area. It extends between 73°40′–96°23′E and 34°25′–48°10′N and covers an area of 166.04 × 10^{4} km^{2} (Figure 1). There are three mountain ranges in Xinjiang. From south to north, they are Kunlun, Tianshan, and Altay mountains. With their high elevations, these mountains block atmospheric circulations and create two vast desert basins in their rain-shadows between the mountains, that is, the Tarim Basin in the south and the Junggar Basin in the north. The Tianshan Mountain in the middle divides Xinjiang into the northern and southern parts. Northern Xinjiang has a continental arid and semiarid climate, with a mean temperature of −13°C in winter and 22.2°C in summer. Southern Xinjiang has a continental dry climate, with a winter mean temperature of −5.7°C and a summer mean temperature of 24.4°C. Annual precipitation is about 210 mm in northern Xinjiang while southern Xinjiang has less than 100 mm. Because of the dry climate, evaporation in Xinjiang is very strong with a mean annual pan evaporation between 1000 and 4500 mm, which is 500–1000 mm higher than other places at the same latitude in China [11, 20].

##### 2.2. Data

To ensure consistency and the longest continuous observation, data from 51 ground-based meteorological stations of the China Meteorological Administration (CMA) were used in this study. In order to understand the complexity of climate process from a multitemporal scale perspective, the daily, monthly, seasonal, and annual air temperature data from 1958 to 2012 were used for computation and analysis.

#### 3. Methods

In order to understand the complexity of temperature dynamics in Xinjiang, China, this paper conducted an integrative approach combining the correlation dimension (CD), classical statistics, and geostatistics method. Firstly, the CD value was computed to show the chaotic and fractal characteristics of temperature dynamics at different temporal scales. Secondly, correlations between the CD value with geographical location and elevation was showed by the correlation analysis and stepwise regression. Finally, the variogram and cokriging methods were used to reveal the spatial pattern of the CD values.

##### 3.1. Correlation Dimension

The correlation dimension (CD) is usually applied to analyze a time series and determine if it exhibits a chaotic dynamic characteristic [21, 22]. Consider , the time series of annual runoff, and suppose that it is generated by a nonlinear dynamic system with degrees of freedom. To restore the dynamic characteristic of the original system, it is necessary to construct an appropriate series of state vectors, , with delay coordinates in the -dimensional phase space according to the basic ideas initiated by Grassberger and Procaccia [23]: where is the embedding dimension and is an appropriate time delay.

The trajectory in the phase space is defined as a sequence of -dimensional vectors. If the dynamics of the system can be reduced to a set of deterministic laws, the trajectories of the system converge toward a subset of the phase space, which is called an “attractor.” Many natural systems do not conform with time to a cyclic trajectory. Some nonlinear dissipative dynamic systems tend to shift toward the attractors for which the motion is chaotic, that is, not periodic and unpredictable over long times. The attractors of such systems are called strange attractors. For the set of points on the attractor, using the G-P method [23], the correlation-integrals are defined to distinguish between stochastic and chaotic behaviors.

The correlation-integrals can be defined as follows: where is the surveyor's rod for distance, is the number of reference points taken from , and is the number of points, . The relationship between and is . is the Heaviside function, which is defined as

The expression counts the number of points in the dataset that are closer than the radius, , within a hypersphere of the radius, , and then divides this value by the square of the total number of points (because of normalization). As , the correlation exponent, , is defined as

It is apparent that the correlation exponent, , is given by the slope coefficient of versus . According to (, , can be obtained by the least squares method (LSM) using a log-log grid (as shown in Figure 2).

To detect the chaotic behavior of the system, the correlation exponent has to be plotted as a function of the embedding dimension (as shown in Figure 3).

If the system is purely random (e.g., white noise), the correlation exponent increases as the embedding dimension increases, without reaching the saturation value. If there are deterministic dynamics in the system, the correlation exponent reaches the saturation value, which means that it remains approximately constant as the embedding dimension increases. The saturated correlation exponent is called the correlation dimension (CD) of the attractor. The CD belongs to the invariants of the motion on the attractor. It is generally assumed that the CD equals the number of degrees of freedom of the system, and higher embedding dimensions are therefore redundant. For example, to describe the position of the point on the plane (two-dimensional system), the third dimension is not necessary because it is redundant. In addition, the CD value is often fractal and represented as a nonintegral dimension, which is typical for chaotic dynamical systems that are very sensitive to initial conditions.

The CD value provides the information regarding the dimension of the phase-space required for embedding the attractor. It is important for determining the number of dimensions necessary to embed the attractor and the number of variables present in the evolution of the process.

We used the previous correlation dimension method to analyze the chaotic and fractal characteristics for the temperature dynamics in this study.

##### 3.2. Correlation Analysis and Stepwise Regression

Correlation and regression analyses are the two commonly useful methods in various disciplines of geography [24], which were used to check the correlations between the CD value with geographical location and elevation in this study.

The correlation analysis is one of the most useful classical statistics, which is a statistical measurement of the correlationship between two variables. Possible correlations range from +1 to –1. A zero correlation indicates that there is no relationship between the variables. A negative correlation indicates that as one variable goes up, the other goes down. A positive correlation indicates that both variables move in the same direction together.

For the two variables, and , the correlation coefficient is calculated as where is the sample number; represents the value of for the sample ; represents the value of for the sample ; is the mean for all ; is the mean for all . Commonly, testing the significance of the correlation coefficient employs the distribution.

Stepwise regression can be achieved either by trying out one independent variable at a time and including it in the regression model if it is statistically significant, or by including all potential independent variables in the model and eliminating those that are not statistically significant, or by a combination of both methods. The multiple linear regression equation (MLRE) is as follows: where is dependent variable and is the coefficient of the independent variables (). In this study, the dependent variable is the CD value and the independent variables are elevation, latitude, and longitude.

##### 3.3. Geostatistics

Studies have shown that the parameters of temperature dynamics are typical regionalized variables, which are structural as well as stochastic [25, 26]. So its spatial variability can be analyzed by the geostatistics method [27, 28].

###### 3.3.1. The Variogram

The regionalized variable is regarded as the value of a variable at a location as a realization of a stochastic . This stochastic is assumed to be intrinsically stationary. The first is that the expected value of the stochastic, [], is constant for all . Secondly, the variance of the differences between the values of the variable at two different locations depends only on the lag vector separating the two locations and not on the absolute locations. In general, this variance may be a function of both the direction and length of the lag vector. If the regionalized variable is isotropic, the variogram is purely a function of the length of the vector which we denote by . Thus the relationship between values from different locations is described by the variogram as follows [27, 28]:

The variogram is estimated from variable values observed at sampled points, , . The method of estimator is the average of squared differences between observations separated by distance as follows: where indicates the magnitude of regionalized variable and is the total number of pairs of attributes that are separated by a distance .

###### 3.3.2. Kriging and Cokriging Methods

Based on the variogram, Kriging and cokriging can be used to estimate the values of regionalized variable at unsampled locations [29, 30].

Ordinary Kriging can mathematically be defined as given in the following: where is the estimated value and is the corresponding weight of each observation on the estimation. These weights are calculated to ensure that the estimator is unbiased and the estimation variance is a minimum. The nonbias condition requires that where is the variogram between sampled point and point , is the variogram between sampled point and estimated point, and is the Lagrange multiplier of minimum condition.

The general form of cokriging equations is where and are the primary and covariate (secondary) variables, respectively. In the cokriging method, the and are cross-correlated and the covariate contributes to the estimation of the primary variable. Generally, measuring the covariate is simpler than measuring the primary variable. For cokriging analysis, the cross variogram (or cross-variogram) should be determined in prior. Provided that there are points where both and have been measured, the cross-variogram is estimated by

#### 4. Results and Discussion

##### 4.1. The Chaotic Dynamic Process

Based on the meteorological data, we analyzed the chaotic dynamics with fractal characteristic for the temperature dynamics by using the G-P method [23].

Firstly, we randomly selected the time series of monthly data from 7 meteorological stations (i.e., Altay, Tacheng, Karamay, Urumqi, Turpan, Korla, and Hotan station) for a pilot study. The plots of correlation exponent versus embedding dimension were drawn as Figure 4.

Each plot in Figure 4 showed the gradual saturation process of the correlation exponent. It is evident that the correlation exponent increases with embedding dimension, , and a saturated correlation exponent, the correlation dimension of attractor, that is, CD, was obtained when .

Then, we repeated the previous process for the time series of daily, monthly, seasonal, and annual data in each meteorological station, and the results revealed that the correlation exponent reached the saturation value, which demonstrated that there is correlation dimension of the attractor in the temperature process at the each temporal scale (i.e., daily, monthly, seasonal, and annual scales).

Table 1 showed the correlation dimensions, that is, CD values, at daily, monthly, seasonal, and annual scales for 51 meteorological stations.

Because none of the CD values in Table 1 is integer, this indicated that the temperature process at each temporal scale is chaotic dynamic system with a fractal characteristic and is sensitive to the initial conditions.

##### 4.2. The Complexity of Temporal Scale

The last row in Table 1 showed the mean of CD values for every meteorological station at daily, monthly, seasonal, and annual scales.

The order of the MCD (2.5353 > 1.6397 > 1.4156 > 1.2995) reveals the complex order of the temperature dynamics at daily, monthly, seasonal, and annual scales; that is, the complexity of temperature dynamics decreases along with the increase of temporal scale. We think that the results accord with the facts, because the daily data series contains more details, and then followed by monthly data series, seasonal data series, and annual daily data series, respectively. Thus, we conclude that the temperature process at a smaller temporal scale is more complex than that at a larger temporal scale.

The MCD values also provided the information about the dimension of the phase-space required for embedding the attractor. Because all the CD values at daily scale are above 2, at least 3 independent variables are needed at to describe the dynamics of temperature process at daily scale. The same reason tells us that at least 2 independent variables are needed to describe the dynamics of temperature process at monthly, seasonal, and annual scales.

##### 4.3. The Effect of Geographical Location and Elevation

Table 1 showed that the CD values at different sites (the sites of meteorological station) are different at a same temporal scale (i.e., daily, monthly, seasonal, and annual scales). Maybe their spatial patterns are affected by the geographical location and elevation.

To reveal the correlation of the CD value with geographical location and elevation, we computed the correlation coefficients as in Table 2.

Table 2 showed that on the daily scale, the CD value positively correlates with elevation at the significant level of 0.05, whereas on the monthly scale, the CD value positively correlates both with elevation and latitude at the significant levels of 0.05 and 0.01, respectively. To verify the correlation, we used the stepwise regression analysis method to fit the multiple linear regression equations (MLREs) between the CD value with geographical location and elevation at daily and monthly scales, which are as in Table 3.

Table 3 told us that on the daily and monthly scales, the CD values are well explained by the geographical location and elevation at the significant levels of 0.006 and 0.000. The MLREs in Table 3 indicate that the site with higher elevation and latitude has a higher CD value. That is to say, the temperature dynamics at the site with higher elevation and latitude are of much higher complexity.

Though the MLREs in Table 3 well explained the relation between CD value with geographical location and elevation at daily and monthly scales, the CD value has no significant correlation with elevation, latitude, and longitude at the seasonal and annual scales. What is the reason for this?

Actually, beside the structural factor such as atmospheric circulation, the local temperature dynamics are also affected by the location, elevation, and other stochastic factors. Therefore, the CD value of temperature dynamic is a typical regionalized variable and its spatial pattern should be described by the variogram.

By using the aforementioned method for computing variogram, we fitted two variograms to describe the spatial variability of CD value at the seasonal and annual scales.

At seasonal scale, the spatial variability of CD value was well described by the variogram of Gaussian model as follows: where is the value of variogram, and is distance. The mean error and average standard error for model (13) are −0.0008275988 and 0.1726933, respectively.

At annual scale, the variogram of Gaussian model well described the spatial variability of CD value as follows: where and have the same meaning as in formula (13). The mean error and average standard error for model (14) are 0.0001671542 and 0.1709583, respectively.

Based on the previous models of variogram (13) and (14), choosing elevation and latitude as the two covariate variables, we used the aforementioned cokriging method to compute the interpolating of CD values at seasonal and annual scales.

Figure 5 presented the spatial pattern of CD values at seasonal scale, which showed that all the CD values are between 1.13 and 1.83. The higher values mainly distribute in the Tianshan, Kunlun, and Altun Mountains, which indicates that the temperature dynamics in these mountain areas are more complicated than other areas. The lower values mainly distribute in the Tarim Basin and the Hami Basin, which indicates that the complexity of the temperature dynamics in these basin areas is comparatively lower than other areas.

Figure 6 presented the spatial pattern of CD values at seasonal scale, which showed that all the CD values are between 1 and 1.51. Comparing it with Figure 5, the pattern of spatial distribution is a little different. The higher values mainly distribute in the Junggar Basin and part of the Altan, Kunlun, and Altun Mountains, whereas the lower values mainly distribute in the Tarim Basin, the Turpan Basin, and the Hami Basin.

Summarizing the results of Section 4.3, we came to the results at seasonal and annual scales as that the higher CD values mainly distribute on complex landform such as mountain areas, whereas the lower CD values mainly distribute on the comparative flat landform such as basin area. The results indicate that the complex temperature dynamics are derived from the complex landform.

#### 5. Conclusion

Summarizing the previous results, we elicited the conclusions as follows.(1)The integer CD values indicate that the temperature dynamics are a complex and chaotic system, which is sensitive to the initial conditions.(2)The order of the MCD (2.5353 > 1.6397 > 1.4156 > 1.2995) reveals the complex order of the temperature dynamics at daily, monthly, seasonal, and annual scales, that is, the complexity of temperature dynamics decreases along with the increase of temporal scale. To describe the temperature dynamics, at least 3 independent variables are needed at daily scale, whereas at least 2 independent variables are needed at monthly, seasonal, and annual scales.(3)The MLREs at the daily and monthly scales show that the site with higher elevation and latitude has a higher CD value, which indicates that the temperature dynamics at the site with higher elevation and latitude are of much higher complexity. (4)The results of the interpolating from cokriging method based on the variogram at seasonal and annual scales show that the higher CD values mainly distribute on complex landform such as mountain area, while the lower CD values mainly distribute on the comparative flat landform such as basin area. The results indicate that the complex temperature dynamics come from the complex landform.

#### Acknowledgments

This work was supported by National Basic Research Program of China (973 Program; no. 2010CB951003). The authors are grateful to the editor and referees whose comments helped them for improving the paper’s quality.

#### References

- J. W. Kantelhardt, D. Rybski, S. A. Zschiegner et al., “Multifractality of river runoff and precipitation: comparison of fluctuation analysis and wavelet methods,”
*Physica A*, vol. 330, no. 1-2, pp. 240–245, 2003. View at Publisher · View at Google Scholar · View at Scopus - S. Li, Q. Zhou, S. Wu, and E. Dai, “Measurement of climate complexity using sample entropy,”
*International Journal of Climatology*, vol. 26, no. 15, pp. 2131–2139, 2006. View at Publisher · View at Google Scholar · View at Scopus - J. A. Rial, R. A. Pielke Sr., M. Beniston et al., “Nonlinearities, feedbacks and critical thresholds within the earth's climate system,”
*Climatic Change*, vol. 65, no. 1-2, pp. 11–38, 2004. View at Publisher · View at Google Scholar · View at Scopus - J. H. Xu, Y. N. Chen, W. H. Li, M. H. Ji, S. Dong, and Y. L. Hong, “Wavelet analysis and nonparametric test for climate change in Tarim River Basin of Xinjiang during 1959–2006,”
*Chinese Geographical Science*, vol. 19, no. 4, pp. 306–313, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. H. Xu, Y. N. Chen, W. H. Li, Q. Nie, Y. L. Hong, and Y. Yang, “The nonlinear hydro-climatic process in the Yarkand River, northwestern China,”
*Stochastic Environmental Research and Risk Assessment*, vol. 27, no. 2, pp. 389–399, 2013. View at Google Scholar - J. H. Xu, Y. N. Chen, W. H. Li et al., “Combining BPANN and wavelet analysis to simulate hydro-climatic process—a case study of the Kaidu River, NW China,”
*Frontiers of Earth Science*, vol. 7, no. 2, pp. 227–237, 2013. View at Google Scholar - D. W. Pierce, T. P. Barnett, B. D. Santer, and P. J. Gleckler, “Selecting global climate models for regional climate change studies,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 106, no. 21, pp. 8441–8446, 2009. View at Publisher · View at Google Scholar · View at Scopus - T. N. Palmer, “A nonlinear dynamical perspective on climate prediction,”
*Journal of Climate*, vol. 12, no. 2, pp. 575–591, 1999. View at Google Scholar · View at Scopus - T. N. Palmer, “Predicting uncertainty in forecasts of weather and climate,”
*Reports on Progress in Physics*, vol. 63, no. 2, pp. 71–116, 2000. View at Publisher · View at Google Scholar · View at Scopus - H. Du, Z. Wu, M. Li, Y. Jin, S. Zong, and X. Meng, “Characteristics of extreme daily minimum and maximum temperature over Northeast China, 1961–2009,”
*Theoretical and Applied Climatology*, vol. 111, no. 1-2, pp. 161–171, 2013. View at Publisher · View at Google Scholar · View at Scopus - Q. Li, Y. Chen, Y. Shen, X. Li, and J. Xu, “Spatial and temporal trends of climate change in Xinjiang, China,”
*Journal of Geographical Sciences*, vol. 21, no. 6, pp. 1007–1018, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. J. Wang, Y. N. Chen, Z. S. Chen, and W. H. Li, “Changes in annual and seasonal temperature extremes in the arid region of China, 1960–2010,”
*Natural Hazards*, vol. 65, no. 3, pp. 1913–1930, 2013. View at Google Scholar - Q. Zhang, J. Li, Y. D. Chen, and X. Chen, “Observed changes of temperature extremes during 1960–2005 in China: natural or human-induced variations?”
*Theoretical and Applied Climatology*, vol. 106, no. 3-4, pp. 417–431, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. H. Xu, Y. W. Xu, and C. N. Song, “An integrative approach to understand the climatic-hydrological process: a case study of Yarkand River, Northwest China,”
*Advances in Meteorology*, vol. 2013, Article ID 272715, 9 pages, 2013. View at Publisher · View at Google Scholar - B. F. Li, Y. N. Chen, and X. Shi, “Why does the temperature rise faster in the arid region of northwest China?”
*Journal of Geophysical Research*, vol. 117, no. D16, Article ID D16115, 2012. View at Google Scholar - Y. F. Shi, Y. P. Shen, D. L. Li et al., “Discussion on the present climate change from warm-dry to warm-wet in Northwest China,”
*Quaternary Sciences*, vol. 23, no. 2, pp. 152–164, 2003 (Chinese). View at Google Scholar - Y. Shi, Y. Shen, E. Kang et al., “Recent and future climate change in northwest China,”
*Climatic Change*, vol. 80, no. 3-4, pp. 379–393, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. S. Wang, X. L. Fei, and F. Wei, “Further study of temperature change in northwest China in recent 50 years,”
*Journal of Desert Research*, vol. 28, no. 4, pp. 724–732, 2008 (Chinese). View at Google Scholar - J. Xu, Y. Chen, F. Lu, W. Li, L. Zhang, and Y. Hong, “The Nonlinear trend of runoff and its response to climate change in the Aksu River, western China,”
*International Journal of Climatology*, vol. 31, no. 5, pp. 687–695, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. Ling, H. Xu, J. Fu, Q. Zhang, and X. Xu, “Analysis of temporal-spatial variation characteristics of extreme air temperature in Xinjiang, China,”
*Quaternary International*, vol. 282, pp. 14–26, 2012. View at Publisher · View at Google Scholar · View at Scopus - B. Sivakumar, “Nonlinear determinism in river flow: prediction as a possible indicator,”
*Earth Surface Processes and Landforms*, vol. 32, no. 7, pp. 969–979, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. Xu, Y. Chen, W. Li, M. Ji, and S. Dong, “The complex nonlinear systems with fractal as well as chaotic dynamics of annual runoff processes in the three headwaters of the Tarim River,”
*Journal of Geographical Sciences*, vol. 19, no. 1, pp. 25–35, 2009. View at Publisher · View at Google Scholar · View at Scopus - P. Grassberger and I. Procaccia, “Characterization of strange attractors,”
*Physical Review Letters*, vol. 50, no. 5, pp. 346–349, 1983. View at Publisher · View at Google Scholar · View at Scopus - J. H. Xu,
*Mathematical Methods in Contemporary Geography*, Higher Education Press, Beijing, China, 2002. - G. Hudson and H. Wackernagel, “Mapping temperature using kriging with external drift: theory and an example from Scotland,”
*International Journal of Climatology*, vol. 14, no. 1, pp. 77–91, 1994. View at Google Scholar · View at Scopus - W. Z. Yue, J. H. Xu, H. J. Liao, and L. H. Xu, “Applications of spatial interpolation for climate variables based on geostatistics: a case study in Gansu Province, China,”
*Geographic Information Sciences*, vol. 9, no. 1-2, pp. 71–77, 2003. View at Google Scholar - P. A. Dowd, “The variogram and kriging: robust and resistant estimators,” in
*Geostatistics For Natural Resources Characterization*, G. Verly, Ed., pp. 91–106, D. Reidel Publishing Co., Amsterdam, The Netherlands, 1984. View at Google Scholar - A. G. Journel and C. Hiujbregts,
*Mining Geostatistics*, Academic Press, New York, NY, USA, 1978. - P. Goovaerts,
*Geostatistics For Natural Resources Evaluation*, Oxford University Press, New York, NY, USA, 1997. - J. A. Vargas-Guzmán and T.-C. J. Yeh, “Sequential kriging and cokriging: two powerful geostatistical approaches,”
*Stochastic Environmental Research and Risk Assessment*, vol. 13, no. 6, pp. 416–435, 1999. View at Google Scholar · View at Scopus