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% C.L. 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% C.L. in both PADE-I, and PADE-II. Although from the information criteria analysis, namely, AIC, BIC, DIC, the non-flat Λ-Cold Dark Matter model remains the best choice, however, in the light of DIC, PADE parametrizations are still appealing.