Biomechanical modeling of the facial soft tissue behavior is needed in aesthetic or maxillo-facial surgeries where the simulation of the bone displacements cannot accurately predict the visible outcome on the patient’s face. Because these tissues have different nature and elastic properties across the face, depending on their thickness, and their content in fat or muscle, individualizing their mechanical parameters could increase the simulation accuracy. Using a specifically designed aspiration device, the facial soft tissues deformation is measured at four different locations (cheek, cheekbone, forehead, and lower lip) on 16 young subjects. The stiffness is estimated from the deformations generated by a set of negative pressures using an inverse analysis based on a Neo Hookean model. The initial Young’s modulus of the cheek, cheekbone, forehead, and lower lip are respectively estimated to be 31.0 kPa ± 4.6, 34.9 kPa ± 6.6, 17.3 kPa ± 4.1, and 33.7 kPa ± 7.3. Significant intra-subject differences in tissue stiffness are highlighted by these estimations. They also show important inter-subject variability for some locations even when mean stiffness values show no statistical difference. This study stresses the importance of using a measurement device capable of evaluating the patient specific tissue stiffness during an intervention.