This tool allows you to compute the excitation functions
of the Earth rotation χ1, χ2, χ3
(according to the "Euler-Liouville" formalism) and to compare
them to the geophysical excitation functions, as far as these later
ones are available. Comparison is done through visual plot and computation
of the correlation coefficients.
- The observed excitation functions (χx, χy,
χz) are computed from the pole coordinates (x,y) and
length of day changes ΔLOD of the IERS C04 series (sampling
of 1 day, fluctuations > 6 days) according to the equations : χx
+ i χy=(x-i y) + i/σc d(x-iy)/dt where
σc is the Chandler angular frequency : (σc
= 2π/T ( 1 + i / 2 Q), T Chandler period, Q quality factor) and
χz=ΔLOD (actually the true adimensionel axial
excitation function is ΔLOD / LOD with LOD =86400 s TAI). The
equatorial excitation function χx + i χy
is computed according to an algorithm such as the one introduced by
Wilson and Vicente (1981,1985,2002) .
- Comparison can be done with :
- Atmospheric Angular Momentum Functions provided
by the Special Bureau for Atmosphere of the IERS (NCEP-NCAR reanalysis
time series) (http://files.aer.com/aerweb/AAM)
from 1962 to 2016/1/1 (see readme).
From 2016/1/1 up to the current week these data are
completed by those of the operational NCEP model
(
ftp://ftp.cpc.ncep.noaa.gov/long/aam/nmc).
In March 2019 we add the
operational ECWMF series
as computed by GFZ, daily updated and spanning the period 1976-one day back from now
- Oceanic Angular Momentum (OAM) functions of the ECCO model : file ECCO_50yr.chi
(10 day values from 1962 to 1993 daily interpolated at 0hUTC) + file ECCO_kf080
(daily values from 1993 to last current year, assimilation of altimetric measurements of the sea surface, NCEP/NCAR atmospheric forcing),
provided
by the Special Bureau for Oceans of IERS (responsible : Richard Gross).
In March 2019 we add the
operational MPIOP series,
as computed by GFZ from MPIOM (Max-Planck-Institute Ocean Model), daily updated and spanning the period 1976-one day back from now
- Hydrological Angular Momentum (HAM) functions (since March 2019).
HAM series of the LSDM model,
as estimated by GFZ, daily updated and spanning the period 1976-one day back from now
- Sea Level Angular Momentum (SLAM) functions (since March 2019).
SLAM series of the ECMWF-MPIOM-LSDM coupled model,
as estimated by GFZ, daily updated and spanning the period 1976-one day back from now
- For the observed axial excitation the effect of zonal gravitational
tides is removed according to the model recommanded by IERS Conventions 2010.
- The atmospheric and oceanic angular momentum functions are filtered
and sampled before comparison. By default the pressure term is associated
with oceans reacting as "Inverted Barometer" (IB) in front
of the pressure variations : this is a realist approximation for variations
larger than 10 days (for rough pressure term click the button "Non
IB").
- Geodetic Excitation functions are derived from C04 pole coordinates using the Wilson (1985) digital filter.
- The equatorial excitation function is based upon the knowledge
of the Chandler term period T and its quality factor Q. As these parameters
are affected by large uncertainties, we let you the possibility to
tune them within the allowed bands (426 <T<439 days; 50<Q<200).
The excitation function in the Chandler frequency band can be spoiled
by bad choice of these parameters. Outside this band the choice of
these parameter is not critic.
- Since March 2019 the user can fix by himself the coefficients of the effective angular momentum functions.
Thanks to Henryk Dobslaw and Robert Dill for fruitful discussion pertaining
to GFZ angular momentum series. Thanks to David Salstein (Atmospheric and Environmental Research, Boston)
and Olivier de Viron (Institut de Physique du Globe de Paris) for their
respective advices. C. Bizouard
|