  We investigate the dark energy phenomenology in an extended parameter space
where we allow the curvature density of our universe as a free-to-vary
parameter. The inclusion of the curvature density parameter is motivated from
the recently released observational evidences indicating the closed universe
model at many standard deviations. Here we assume that the dark energy
equation-of-state follows the PADE approximation, a generalized parametrization
that may recover a variety of existing dark energy models. Considering three
distinct PADE parametrizations, labeled as PADE-I, SPADE-I and PADE-II, we
first constrain the cosmological scenarios driven by them using the joint
analyses of a series of recently available cosmological probes, namely,
Pantheon sample of Supernovae Type Ia, baryon acoustic oscillations, big bang
nucleosynthesis, Hubble parameter measurements from cosmic chronometers, cosmic
microwave background distance priors from Planck 2018 and then we include the
future Gravitational Waves standard sirens (GWSS) data from the Einstein
telescope with the combined analyses of these current cosmological probes. We
find that the current cosmological probes indicate a very strong evidence of a
dynamical dark energy at more than 99\% CL in both PADE-I, and PADE-II, but no
significant evidence for the non-flat universe is found in any of these
parametrizations. Interestingly, when the future GWSS data from the Einstein
telescope are included with the standard cosmological probes an evidence of a
non-flat universe is found in all three parametrizations together with a very
strong preference of a dynamical dark energy at more than 99\% CL in both
PADE-I, and PADE-II. Although from the information criteria analysis, namely,
AIC, BIC, DIC, the non-flat $\Lambda$-Cold Dark Matter model remains the best
choice, however, in the light of DIC, PADE parametrizations are still
appealing.
