Research Publications: UM-SG-RS-2018-02
Evaluation of statistical methods for quantifying fractal scaling in water-quality time series with irregular sampling.
Zhang, Q; Harman, CJ; Kirchner, J
Source:Hydrology and Earth System Sciences 22(2):1175-1192
This article is freely available online. You can use the DOI number to find it through the journal's website or through a search engine.
River water-quality time series often exhibit fractal scaling, which here refers to autocorrelation that decays as a power law over some range of scales. Fractal scaling presents challenges to the identification of deterministic trends because (1) fractal scaling has the potential to lead to false inference about the statistical significance of trends and (2) the abundance of irregularly spaced data in water-quality monitoring networks complicates efforts to quantify fractal scaling. Traditional methods for estimating fractal scaling in the form of spectral slope (beta) or other equivalent scaling parameters (e.g., Hurst exponent) - are generally inapplicable to irregularly sampled data. Here we consider two types of estimation approaches for irregularly sampled data and evaluate their performance using synthetic time series. These time series were generated such that (1) they exhibit a wide range of prescribed fractal scaling behaviors, ranging from white noise (beta = 0) to Brown noise (beta = 2) and (2) their sampling gap intervals mimic the sampling irregularity (as quantified by both the skewness and mean of gap-interval lengths) in real water-quality data. The results suggest that none of the existing methods fully account for the effects of sampling irregularity on beta estimation. First, the results illustrate the danger of using interpolation for gap filling when examining autocorrelation, as the interpolation methods consistently underestimate or overestimate beta under a wide range of prescribed beta values and gap distributions. Second, the widely used Lomb-Scargle spectral method also consistently under-estimates beta. A previously published modified form, using only the lowest 5% of the frequencies for spectral slope estimation, has very poor precision, although the overall bias is small. Third, a recent wavelet-based method, coupled with an aliasing filter, generally has the smallest bias and root-mean-squared error among all methods for a wide range of prescribed beta values and gap distributions. The aliasing method, however, does not itself account for sampling irregularity, and this introduces some bias in the result. Nonetheless, the wavelet method is recommended for estimating beta in irregular time series until improved methods are developed. Finally, all methods' performances depend strongly on the sampling irregularity, highlighting that the accuracy and precision of each method are data specific. Accurately quantifying the strength of fractal scaling in irregular water-quality time series remains an unresolved challenge for the hydrologic community and for other disciplines that must grapple with irregular sampling.
Related Research Project(s) Funded by Maryland Sea Grant:
'Related Research Project(s)' link to details about research projects funded by Maryland Sea Grant that led to this publication. These details may include other impacts and accomplishments resulting from the research.
'Maryland Sea Grant Topic(s)' links to related pages on the Maryland Sea Grant website.