Indranil Sahoo,Joseph Guinness,Brian J. Reich
Abstract:Analysis of geostatistical data is often based on the assumption that the spatial random field is isotropic. This assumption, if erroneous, can adversely affect model predictions and statistical inference. Nowadays many applications consider data over the entire globe and hence it is necessary to check the assumption of isotropy on a sphere. In this paper, a test for spatial isotropy on a sphere is proposed. The data are first projected onto the set of spherical harmonic functions. Under isotropy, the spherical harmonic coefficients are uncorrelated whereas they are correlated if the underlying fields are not isotropic. This motivates a test based on the sample correlation matrix of the spherical harmonic coefficients. In particular, we use the largest eigenvalue of the sample correlation matrix as the test statistic. Extensive simulations are conducted to assess the Type I errors of the test under different scenarios. We show how temporal correlation affects the test and provide a method for handling temporal correlation. We also gauge the power of the test as we move away from isotropy. The method is applied to the near-surface air temperature data which is part of the HadCM3 model output. Although we do not expect global temperature fields to be isotropic, we propose several anisotropic models with increasing complexity, each of which has an isotropic process as model component and we apply the test to the isotropic component in a sequence of such models as a method of determining how well the models capture the anisotropy in the fields.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to conduct an isotropy test on spatial data on the sphere. Specifically, the author proposes a method based on spherical harmonic functions to test whether the spatial random field on the sphere satisfies the isotropy hypothesis. If this hypothesis does not hold, that is, there are differences in the correlations of the spatial random field in different directions, it will lead to errors in model prediction and statistical inference. Therefore, developing an isotropy test method suitable for spherical data is crucial for improving the accuracy and reliability of the model.
### Background of the Paper and Problem Definition
1. **Importance of Isotropy**:
- In geostatistical data analysis, it is usually assumed that the spatial random field is isotropic, that is, the correlation between observations only depends on the distance between them, not on the relative orientation.
- If this assumption does not hold, it will lead to deviations in model prediction and statistical inference.
2. **Special Characteristics of Spherical Data**:
- With the development of technology, many applications (such as atmospheric science) now need to process data covering the entire earth.
- The inherent correlation structure of spherical data is different from that of data in Euclidean space, so special methods are required to analyze spherical data.
3. **Limitations of Existing Methods**:
- Existing isotropy test methods are mainly applicable to data in Euclidean space, and usually assume stationarity as an alternative hypothesis.
- These methods cannot be directly applied to spherical data because there are no concepts of stationarity and geometric anisotropy on the sphere.
### Proposed Method
1. **Spherical Harmonic Function Representation**:
- The author uses spherical harmonic functions as orthogonal basis functions on the sphere and projects the data onto these basis functions.
- Under the isotropy hypothesis, the spherical harmonic coefficients are uncorrelated; while in the non - isotropy case, these coefficients are correlated.
2. **Test Statistic**:
- The author proposes using the maximum eigenvalue of the sample correlation matrix as the test statistic.
- The Type I error rate and power of this test in different types of scenarios are evaluated through simulation experiments.
### Application Examples
1. **HadCM3 Model Data**:
- The author applies this method to the near - surface air temperature data of the HadCM3 climate model, which is a daily temperature data set covering 2031 - 2035.
- The author constructs a series of anisotropic models with increasing complexity and tests the isotropic components of each model.
2. **Result Analysis**:
- The test results show that as the complexity of the model increases, the model's ability to explain the anisotropy of the temperature field gradually increases.
- Even for the most complex model, its isotropic component is still strongly rejected, indicating that there is indeed significant anisotropy in the temperature field.
### Conclusion
This paper successfully proposes an isotropy test method suitable for spherical data and verifies the effectiveness and practicality of this method through actual data. This method provides a new tool for processing spherical data and helps to improve the accuracy and reliability of geostatistical models.