The obliquity of the ecliptic is essentially the angle of inclination of the Earth's equator with respect to the plane of its orbit, which to us, produces the apparent relative tilt of the Earth's polar axis which is the principal cause of the annual seasons.
This angle is approximately ±23.45 degrees and is the apparent maximum/minimum angle of the sun above or below the equator during the year as viewed from Earth. This value varies very slowly as the orbit of the Earth evolves over time in response to various internal and external influences, such as uneven mass distribution, tidal and gravitational perturbations of the sun, moon and other bodies, etc.
This is one of the fundamental values we need to know when computing the apparent position of a planet as viewed from the Earth at any given moment as it defines the instantaneous inclination of the Earth's orbit.
The ecliptic obliquity can also be used to determine the latitudes of the arctic and antarctic circles, however these values can vary slightly depending on the computational methods, the values of the applied constants, time scale, coordinates system, the applied geoid model, the accuracy of the applied nutational theory, among other factors, some unpredictable.
In the computation of the arctic/antarctic circles, the difference between applying the IAU 1980 nutation series or the 2000B series only amounts to a matter of inches difference on the ground. Due to the margin of error involved in such computations and the degree to which these margins overlap, this small difference is generally too insignificant to matter when computing the latitudes of the arctic/antarctic circles.
Laskar's formula for the mean obliquity theoretically spans 10,000 years, however, I have yet to find any reliable information as to the full validity span of either of the nutational series, so how accurately the computed nutations would apply to the mean obliquity 6000 years ago or in the future from now is still a bit hazy.
Over a period of several years, these values have varied as measurement-technologies improved their accuracy and consistency. They will continue to vary for the same reasons in the future, so the same predictions based on them will also vary to a small degree from generation to generation.
As for the theoretical accuracy, below are some excerpts taken from the original paper summary:
..... Excerpts from J. Laskar's original paper summary: Using the methods of the general theory given in (Laskar, 1985), we have analytically computed the differential system giving the secular variations of the orbital elements for the 8 major planets, at the order 2 with respect to the masses and up to degree 5 in the eccentricity-inclination variables with a relative precision of 10E-6. Relativistic and lunar perturbations are included. The entire system is integrated numerically over 10,000 years and then developed in Taylor expansion around J2000. We obtain new polynomial secular terms for the inner planets up to the power 10 of the time. Comparisons are made with Bretagon's theory VSOP82 and with the numerically integrated JPL ephemeris, DE102 (Newhall et al., 1983). The global accuracy is approximately 0.04202" / 1,000 yr for the inclination of the Earth. Using the theory of the rotation of the rigid Earth of Kinoshita (1977), we derive new formulas for the precessional quantities, up to the 10th power of t, and valid over 10,000 years. .....
PHP function to compute the mean obliquity in decimal degrees for any given JD argument
CPP function to compute the mean obliquity in decimal degrees for any given JD argument
VB.NET function to compute the mean obliquity in decimal degrees for any given JD argument
Jay Tanner - PHP Science Labs - 2018