The present study deals with the stiffness and damping profiles of the leg joints during the ground-contact phase of hopping. A two-dimensional (sagittal plane) jumping model, consisting of four linked rigid segments and including the paired feet, shanks, thighs, and the head - arms - trunk segment, was developed. The segments were interconnected by damped torsional springs, representing the action of the muscles, tendons and ligaments across the joint and of the other joint tissues. A regressive function was used to express stiffness and damping, and included second-order dependence on angle and first-order dependence on angular velocity. By eliminating redundancies in the numerical solution using multicollinearity diagnostic algorithms, the model results revealed that the correct and sufficient nonlinearity for the joint stiffness is of the first order. Damping was found negligible. The stiffness profiles obtained were bell-shaped with a maximum near midstance and nonzero edge values. In predicting the joint moments, the obtained variable joint stiffnesses provided a closer agreement compared to a constant stiffness model. The maximal stiffness was found to be in linear correlation with the initial stiffness in each joint, providing support to the of muscles' preactivation strategy during the flight phase of hopping. All stiffnesses increased with increasing hopping frequency. The model presented provides an effective tool for future designing of artificial legs and robots and for the development of more accurate control strategies.