Modeling global net ecosystem exchange is essential to understanding and quantifying the complex interactions between the Earth’s terrestrial ecosystems and the atmosphere. Emphasizing the inter-relatedness between the global net ecosystem exchange, global sea surface temperature, and atmospheric $${\text {CO}}_{2}$$
CO
2
levels, intuitively suggests that all three systems may exhibit collective environmental memory. Motivated by this, we explicitly identified a collective memory function and showed a similar non-Markovian stochastic behavior for these systems exhibiting superdiffusive behavior in short time intervals. We obtained the values of the memory parameter, $$\mu$$
μ
, and the characteristic frequencies, $$\nu$$
ν
, for global net ecosystem exchange (GNEE) ($$\mu = 0.94\pm 0.03, \nu = 0.67\pm 0.08/mo.$$
μ
=
0.94
±
0.03
,
ν
=
0.67
±
0.08
/
m
o
.
), global sea surface temperature (GSST) ($$\mu = 0.68\pm 0.11, \nu = 0.30\pm 0.18/mo.$$
μ
=
0.68
±
0.11
,
ν
=
0.30
±
0.18
/
m
o
.
), and atmospheric $${\text {CO}}_{2}$$
CO
2
($$\mu = 0.78\pm 0.08, \nu =0.66\pm 0.13/wk.$$
μ
=
0.78
±
0.08
,
ν
=
0.66
±
0.13
/
w
k
.
). The values of the memory parameter are within the range, $$0< \mu < 1$$
0
<
μ
<
1
, and thus all three systems are in the superdiffusive regime. We emphasize, further, that these results were consistent with our previous analyses at the ecosystem level (i.e. Great Barrier Reef) suggesting scale invariance for these phenomena. Thus, the observed superdiffusive behavior operating at different scales suggests universality of the collective memory function for these systems.