Calculation of the underground temperature resulting from heat injection/extraction into/from ground heat exchangers (GHEXs) with hourly variation is one of the most noteworthy challenges to address when simulating and designing a ground source heat pump (GSHP). In order to overcome this challenge, the authors introduce a method to calculate the underground temperature, by considering heat injection/extraction into/from GHEXs with hourly variation. The method applies the superposition of the infinite cylindrical source (ICS) solution and the infinite line source (ILS) solution to calculate the temperature change due to heat injection/extraction into/from the considered GHEX and other neighboring GHEXs, respectively. The calculation method also considers heat injection/extraction from GHEXs with different heat injection/extraction rates and is able to accommodate GHEXs with large diameters such as energy piles. The calculation method was evaluated by applying it to calculate the temperature variation of the heat carrier fluid in a GSHP system with energy piles.