A practical method for designing hybrid geothermal heat pump (GHP) systems that use closed-loop, earth heat exchangers installed in vertical boreholes is presented. The design difficulty with hybrid GHP systems is inherently an optimization problem that is best solved with a computer-based system simulation method. Many parameters can be optimized and there is no unique expression of the objective function. In this work the optimization problem is defined as balancing the annual thermal loads on the ground by minimizing the borehole heat exchanger length and supplemental equipment size. The supplemental equipment examined in this work has been limited to flat plate solar collectors.The design method for GHP systems was developed from results of 62 detailed computer simulations. Three dimensionless groups containing key GHP design parameters were identified using the Buckingham Pi Theorem, and correlated with a fitted surface equation. With typical design parameters available to a practitioner, the design method developed here can be used to estimate the total ground loop length for stand-alone GHP systems, along with the quantity of annual energy required to balance the annual ground loads. With additional input parameters also readily available to designers, the area of a solar collector array can be calculated, along with the corresponding reduced borehole heat exchanger length. Solar collector array area is calculated using the utilizability method.Units: SI