Gaussian processes (GP) provide an elegant and model-independent method for extracting cosmological information from the observational data. In this work, we employ GP to perform a joint analysis by using the geometrical cosmological probes such as Supernova Type Ia (SN), Cosmic chronometers (CC), Baryon Acoustic Oscillations (BAO), and the H0LiCOW lenses sample to constrain the Hubble constant $$H_0$$
H
0
, and reconstruct some properties of dark energy (DE), viz., the equation of state parameter w, the sound speed of DE perturbations $$c^2_s$$
c
s
2
, and the ratio of DE density evolution $$X = \rho _\mathrm{de}/\rho _\mathrm{de,0}$$
X
=
ρ
de
/
ρ
de
,
0
. From the joint analysis SN+CC+BAO+H0LiCOW, we find that $$H_0$$
H
0
is constrained at 1.1% precision with $$H_0 = 73.78 \pm 0.84\ \hbox {km}\ \hbox {s}^{-1}\,\hbox {Mpc}^{-1}$$
H
0
=
73.78
±
0.84
km
s
-
1
Mpc
-
1
, which is in agreement with SH0ES and H0LiCOW estimates, but in $$\sim 6.2 \sigma $$
∼
6.2
σ
tension with the current CMB measurements of $$H_0$$
H
0
. With regard to the DE parameters, we find $$c^2_s < 0$$
c
s
2
<
0
at $$\sim 2 \sigma $$
∼
2
σ
at high z, and the possibility of X to become negative for $$z > 1.5$$
z
>
1.5
. We compare our results with the ones obtained in the literature, and discuss the consequences of our main results on the DE theoretical framework.