Three-dimensional surface displacement field associated with the 25 April 2015 Gorkha, Nepal earthquake is derived from an integration of Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) measurements, with an extended SISTEM (Simultaneous and Integrated Strain Tensor Estimation From Geodetic and Satellite Deformation Measurements) approach (ESISTEM) proposed in this study. In ESISTEM approach, both surrounding InSAR and GPS measurements can be used as constraints in deriving surface displacements; while only surrounding GPS measurements are used in SISTEM approach. Besides the constraints from surrounding GPS measurements, the ESISTEM approach makes surrounding InSAR measurements available for constraining the derived deformations based on surface elastic theory for the first time. From the north to the south, derived surface displacement field shows prevailing southward horizontal deformations, and gradually varied vertical deformations ranging from −0.95 to 1.40 m within 120 km to the north of Kathmandu. This reveals that ruptures of Main Himalayan thrust (MHT) system were confined in subsurface and did not propagate to the Main Frontal Thrust (MFT) fault, in accordance with field investigation as well as geodetic and seismic studies. Relation between vertical deformations and earthquake-induced landslides is briefly discussed.