Trend Analysis


Cutoff Frequency:
Output: Coefficients Graph
Data:       

Note:

  1. Paste daily mean data of a variable in the form to perform analysis. Each data line must contain year, month, day, and value of the variable. Data columns must be separated by one or more space characters. See the example of daily CO2 data from Hateruma station.
  2. The time sequence do not have to be continuous, but should be incremental. No missing value should be used.
  3. The default value 7.3 was used by Thoning et al. [J. Geophys. Res., 1989, 94(6):8549-8565] to remove short-term variations. You may use a small value (e.g., 0.5) to obtain inter-annual variation.
  4. The server may refuse a request if too many data are submitted.

Method:

The trend analysis uses Fast Fourier Transform with Thoning's exponential frequency filter (Thoning et al., 1989, J. Geophys. Res., 94(6):8549-8565) to remove short-term variations of time series data:

H(f) = exp(-c*(f/fc)p)

where fc is the cutoff frequency in cycles per year, c=ln(2) is the normalizing coefficient to make H(f)=0.5 when f=fc, and p=4.0 is used to modulate the decay of frequency.

Fast Fourier Transform requires a sequence of evenly spaced stationary data. Since a time series of CO2 data has trend and most likely has gaps of missing data, the equation

y = a + b*x + c1*cos(x) + d1*sin(x) + c2*cos(2x) + d2*sin(2x)

is used to remove trend and to fill gaps.

Internally, data of the first and the last years are patched to the beginning and the end respectively to reduce alias effect of end-points. The final results may be used to estimate variations of long-term and seasonal growth rates of green house gases.