We present a new analytical model based on the finite line source that extends the steady state results for parallel horizontal pipes to the transient case and for any desired horizontal pipe layout. The analytical model is validated, when there is no freezing/thawing, by a 3D finite element numerical model. When the phase change is accounted for in the numerical model, the analytical model still provides good approximation to the ground temperature during the heating season and the heat extracted by the ground heat exchanger. However, summer ground temperature and thaw depth are overestimated by the analytical model. A case study for a typical building in Kuujjuaq (northern Canada) area is analyzed. The ground heat exchanger layout follows a spiral pattern characterized by three parameters: length L, depth D. and spacing S. The influence of each parameter on the amount of heat extracted from the ground and on the ground temperature at a control point is assessed. The results show that increasing depth D favors keeping the ground frozen at this depth and increases the amount of heat that may be extracted. Conversely, increasing S and/or L is beneficial for the amount of heat extracted, but it enhances the risk of thawing around the pipes. The model and case study provides useful ground heat exchanger design guidelines in cold regions for the double purpose of ground freezing and heat extraction. (C) 2011 Elsevier Ltd. All rights reserved.