Received 10 March 2016; accepted 24 May 2016; published 27 May 2016
The ongoing economic crisis, dependence on imported energy and the need to tackle climate change have forced the EU to promote more efficient energy production. The Energy Efficiency Directive (2012/27/EU) endorses the application of innovative technology to increase the use of renewable energy systems (RES) in Europe. The need to adopt more RES in all EU countries is vital, as the EU has a commitment to reduce greenhouse gas emissions from 85% to 90% below 1990 levels by 2050  . Finland, located in the north of Europe between the 60th and 70th degrees of latitude, is part of the EU and has set a target of 38% for the contribution of RES to gross overall energy consumption by 2020. In 2012, RES accounted for 35.1% of the overall energy consumption of Finland, 3.7% of which was produced by heat pumps  . One option to increase renewable energy production in Finland is to utilise groundwater reservoirs for energy production. Although groundwater has been widely used for decades as an RES, for instance in China  , Canada  and in Europe  , groundwater utilisation for energy purposes is still a new innovation in Finland, despite the locally significant groundwater energy potential  .
The main aquifers for water supply purposes are situated on shallow Quaternary eskers or ice-marginal end moraine complexes in Finland. The natural groundwater temperature is typically higher than the air temperature in cold regions  -  . This is due to the three reasons: firstly, in winter the snow cover functions as an insulator, preventing cold air conduction into the subsurface layers. Secondly, the change in the state of water releases latent heat into the soil when frost is formed and thirdly, frost acts as an insulator, reducing the flow of cold meltwater into deeper soil layers in early spring, when the melting of snow begins  -  . Finland’s mean air temperature was approximately 2.3˚C during the time period from 1981 to 2010  and average groundwater temperatures varied from 3.0˚C in northern to 6.6˚C in southern parts of the country  -  . The reason for groundwater temperature variations is mainly climatological. Urbanisation has elevated groundwater temperatures by several degrees    , resulting in an increased heating capacity of groundwater in urbanised compared to rural areas  -  . Because of the groundwater low temperature, i.e. low enthalpy energy, it is essential to exploit heating energy from groundwater with heat pumps and/or heat exchangers. Cooling energy can also be exploited using a free cooling system, but groundwater needs to be pumped from the aquifer to the heat transfer system. Hence, groundwater energy utilisation systems, both for heating and cooling, require electricity or some other form of exterior power to work properly. Groundwater exploitation is highly regulated by the Water and Environment act in Finland. However, Finnish legislation has no binding operational limits for groundwater temperature variations during energy utilisation as exist, for example, in Germany  or in Switzerland  .
Groundwater energy utilisation requires site-specific hydrogeological knowledge and an ability to recognise the collective risks caused by variations in the groundwater temperature and level   . Groundwater energy usage has to remain at a thermally sustainable level  to achieve both the environmental and economic operational requirements. Banks  presented a risk assessment approach to sustainable groundwater thermal use for cooling. Groundwater energy utilisation may also have positive environmental influences; in the Netherlands, an ATES system is used for the remediation of chemically contaminated groundwater   . Using groundwater for heating, i.e. injecting cooled groundwater to the aquifer, may provide a solution to reduce the urbanisation impact, which raises groundwater temperatures. Replacing convectional oil heating systems with groundwater energy system will reduce the risk of oil leaks to the aquifer.
The long-term groundwater energy potential of aquifers situated in northern regions has remained undetermined. Neither has the latest information on building energy demands been taken into consideration in groundwater energy system design in Finland. The present study investigated whether groundwater with a natural temperature of 4.9˚C could allow thermally sustainable energy utilisation from a Quaternary aquifer for three large building complexes during an operational period of 50 years. The groundwater energy potential was investigated for: 1) an area with 20 detached houses; 2) area with 3 apartment buildings and 3) a shopping centre. The dimensioning of the groundwater pumping rate was defined by hourly-based information on building energy demand. This study combined the benefits of building energy simulations and hydrogeological modelling.
2. Area Description, Material and Methods
The project area, Karhinkangas, was chosen this study due to available geological and hydrological data and suitable location on the area of naturally cold groundwater temperature. Karhinkangas is located in western Finland, near the Gulf of Bothnia (Figure 1), approximately on the latitude of 64˚ passing through Central Finland. According to the Finnish Meteorological Institute, during the time period from 1981 to 2010, the mean annual air temperature was 3˚C to 4˚C, the mean yearly precipitation 450 mm to 550 mm and snow cover normally persisted for 135 to 140 days in winter.
2.1. Geology and Hydrogeology
The bedrock is dominated by plutonic and metamorphic rocks, including granites, gneisses and granodiorites, and is associated with the Svecofennian (1900 Ma) orogeny  . The bedrock is mainly covered by glacial and postglacial sediments  deposited during the Weichselian glacial stage and the Holocene. The Karhinkangas esker consists of glaciofluvial sand and gravel with a thickness of approximately 20 m  . The average natural groundwater temperature at a depth of approximately 15 m is 4.9˚C according two years of continuous temperature data recorded by the Geological Survey of Finland. Groundwater temperature measurements have been taken by data loggers from 18 monitoring wells.
Figure 1. Location map of the Karhinkangas area. Finland’s capital, Helsinki, is also shown. Basemap database© Esri, DeLorme, Navteq. With permission from Golder Associates global ESRI licence.
2.2. Building Energy Simulations: Reference Year
The energy demand of three types of building was simulated. The net heating power for a detached house with an area of 134 m2 and an apartment building with an area of 814 m2 was simulated using the IDA Indoor Climate and Energy (IDA-ICE) 4.1 dynamic simulation tool, and the heating and cooling power demands of a shopping centre with area of 15,000 m2 were simulated with the RIUSKA application. The simulation results were presented as the hourly-based power distribution during a one-year period. This period is named as the reference year, and describes the current Finnish climatic conditions according to Kalamees et al.  . The thermal insulation of buildings fulfils the minimum requirements of the Finnish Building Code, part C3  . Power for household hot water heating, distribution losses of space heating and domestic hot water were not included in the net power simulations for houses or apartments. These additional energy demands were calculated based on the National Building code, parts D3 and D5   . These codes provide the total energy demand of additional systems, not the hourly-based power distribution. The additional energy demand was adjusted to the hourly power distribution with a four-phase calculation process. First, the hourly energy profiles of distribution losses from space heating and domestic hot water were produced using the same percentual energy distribution as the net power simulations. Secondly, the hourly profile of hot water heating power was produced, noting that hot water heating consumption is 14% less in June and July than in other months in Finland  . Thirdly, additional power distribution profiles were added to the net power profile, and finally, the heating power demand of one detached house was multiplied by 20 and that of one apartment building was multiplied by 3.
2.3. Groundwater Flow and Heat Transport Modelling
The groundwater flow model had previously been completed using the three-dimensional finite differences code MODFLOW  . Heat transport was simulated using the analogy between solute and heat transport. MT3DMS  was used in the simulations. The relationship between the solute transport equation and heat transport was derived by Thorne et al.  . The analogous heat transport equation is of the form:
where θ is the volumetric water content; vi is the mean pore water velocity vector; qs is a source or sink term; ρs is density of solid; ρw is density of water; cs and cw are specific heat capacity of the solid and water, respectively; T is the temperature; is the effective thermal conductivity of the porous media; Ts is the temperature of the source; and Dij is the dispersion tensor  . Thermal diffusion coefficient, Dm, is given as  :
Groundwater flow model and its calibration and validation were conducted in the previous study by Paalijärvi and Okkonen  . Some general information about the flow model is given here but the details of the flow model and its development can be found at that publication. The hydraulic conductivity varies between 1.2 ´ 10−6 - 2.3 ´ 10−3 m/s and groundwater recharge between 0 - 1.1 ´ 10−9 m/s (Figure 2) being 24,328 m3/d for the whole area (24.5 km2).
The heat transport focused on the southern part of the aquifer (Figure 2) and the area of interest was only a few hundreds of square meters, thus, so called regional to local model conversion was carried out  in order to ease the computational effort. From the regional model hydraulic conductivities, recharge rates, ground surface elevations and bedrock are interpolated to local model. The regional flow model was at first run with a different pumping scenarios (see last paragraph of the introduction) to get the groundwater level variations. These groundwater level variations were then applied in the local model as a time variant boundary conditions (Figure 3). Local model was then used in the heat transfer simulations. The maximum head change due to pumping (all scenarios included) at the east and west boundaries were only 0.005 meters. In the north and south no-flux boundary conditions was used. The local model was conceptualised in a rectangular area with a length of 800 m and a width of 400 m. The grid cell size in the local model was 5 m × 5 m (Figure 3). Initial condition was 4.9˚C at the whole domain. Heat content from the recharge and the east and south (time variant specified head) boundaries were 4.9˚C.
Figure 2. Conceptual model of the study area; the local model area is marked as a black box on the southern part of the esker. Red dots are groundwater observation wells. Recharge varies between 0 - 5.5 × 10−9 m/s and 5.5 × 10−9 - 1.1 × 10−8 m/s in light grey and dark grey areas, respectively. Hydraulic conductivity varies between 1.2 × 10−6 - 4.6 × 10−4 m/s in the red zone, 4.6 × 10−4 - 9.3 × 10−4 m/s in the yellow zone, 9.3 × 10−4 - 1.4 × 10−3 m/s in the green zone, 1.4 × 10−3 - 1.9 × 10−3 m/s in the light blue zone and 1.9 × 10−3 - 2.3 × 10−3 m/s in the blue zone.
Figure 3. Regional and local grids. The well distances, groundwater flow direction, heating and cooling sides of the system and injected groundwater temperatures are also shown.
The minimum pumping well distance in the energy transport models was estimated using an analytical equation:
A daily time step was used and the total simulation time was 50 years. The average recharge was used, but the pumping rate was based on the power demand during the reference year and adjusted on a monthly basis. Thus, intra-annual changes were taken into account, but these changes were the same for the whole 50 years. In the MT3DMS simulations, a standard finite difference method with an upstream weighting scheme was used.
The soil thermal conductivity was calculated as the weighted mean thermal conductivity of the main minerals and water (soil pores are filled with water in the saturated zone). The mineralogy was determined at the laboratory of the Geological Survey of Finland, and the thermal conductivity of minerals was taken from Clauser and Huenges  . Bulk and fluid densities, specific heat capacities, porosity, thermal diffusion and dispersion coefficients were estimated based on literature values (see Table 1). Transverse and longitudinal dispersivities were 0.025.
2.4. Groundwater Flux, Heating and Cooling Power Calculations
The groundwater pumping requirements needed to achieve simulated heating and cooling power (reference year) were calculated on an hourly basis (8760 hours in a year) using the following equations   :
Groundwater flux for heating:
Groundwater flux for cooling:
where FH = flux of water (kg∙s−1) for heating and FC for cooling respectively; H is heating power (W) and C cooling power (W); ΔTin-out = temperature difference between incoming and outgoing water in the heat pump / heat exchanger - temperature drop in heating mode and temperature rise in cooling mode (K); cw = specific heat capacity of water (J∙kg−1∙K−1); COPH = coefficient of performance of the heat pump for heating (dimensionless);
Table 1. Heat transport parameters. Square brackets indicate the reference.
and COPC = coefficient of performance of the heat exchanger system for cooling (dimensionless). The groundwater flux unit kg∙s−1 was changed to l∙s−1, as the change has no real effects on the results and l s−1 is universally used to describe groundwater flow. In heating mode, groundwater utilisation down to a temperature of +1.4˚C was assumed, and ΔTin-out was consequently 3.5˚C, because the initial groundwater temperature was 4.9˚C. Based on previous studies i.e. Andersson  and Cruickshanks and Adsett  and noting the Nordic location of Karhinkangas, a maximum groundwater return temperature of 12˚C was used for cooling calculations. Hence, ΔTin-out for cooling was 7.1˚C.
The specific heat capacities of water were taken from Yaws  . A heat pump was used for heating and a heat exchanger for cooling applications. Based on the information presented by Allen et al., Bayer et al., Saner et al.    and the European Heat Pump Association, EHPA  , a COPH of 3.5 and COPC of 25 was used in this analysis. Based on the buildings energy simulation results (see paragraph 2.2) the heat pumps were dimensioned to achieve 60% of the peak design power requirement of a detached house and an apartment building and 50% of that for a shopping centre. Cooling power was dimensioned to achieve 100% from groundwater flux.
3.1. Building Energy and Groundwater Flow Demands: The Reference Year
The largest requirement for power and groundwater flow was for the shopping centre (Table 2). Table 2 also states that heat pumps were designed with a nominal heating power of 120 kW for the detached houses, 75 kW for the apartment buildings and 350 kW for the shopping center. The groundwater flow for the second peak (Table 2) is calculated from the peak hourly value and hence is not based on modelled energy consumption.
Groundwater flow requirements vary significantly between days, especially in the ATES system (Figure 4). The peak pumping rate occurs on 1st and 3rd August, when maximum cooling power is needed. For example, the maximum groundwater flow per hour is 20.98 m3/h for 20 detached houses and 121.08 m3/h for shopping centre (Table 2). The average pumping rate per day is 7.38 m3/d and the median 6.88 m3/d for 20 detached houses and 23.76 m3/h and 18.96 m3/h for shopping centre respectively. The largest groundwater demand for a day is 1572 m3, which is 6.5% of the modelled recharge value of the aquifer.
Table 2. Summary of the energy, power and groundwater flow calculations.
-not analysed; aPeak heating and cooling power denotes the peak heating/cooling power demands of the building; bHeat pump energy and peak heat pump power denote the heating energy and power respectively producible in a building by using a heat pump.
Figure 4. The distribution of groundwater flux demands in (a) 20 detached houses and (b) a shopping centre. The level parts indicate hours when heat pump is working at its maximum power.
3.2. Groundwater Flow and Heat Transport Modelling
The groundwater temperature decreased as groundwater was only used for heating in well-doublet scheme simulations (see Figure 5 and Figure 6). The thermal plume stretched to 300 m in approximately 30 months after pumping had started and groundwater temperatures achieved a steady state condition after approximately 2 years of operation. For example, the temperatures remained constant between 2.8˚C and 2.9˚C in the detached houses scenario and at 3.9˚C in the apartment buildings scenario at an observation point 300 m from the injection well (Figure 5 and Figure 6). The yearly temperature variation in the injection well was 0.15˚C in the detached house model and 0.5˚C in the apartment building model.
In the shopping centre scenario, the thermal plume is more mixed than the other two schemes, as groundwater is utilised both for heating and cooling (Figure 7 and Figure 8). Temperature variations reached a constant annual cycle after five years of operation on the cooling side and approximately after the first year on the heating side (Figure 7). In the steady state conditions, the maximum operating temperature is 3.2˚C and minimum 1.4˚C on the cooling side and the minimum is 5.7˚C and maximum 11.6˚C on the heating side respectively. At the observation point 300 meters downstream from the injection well, the temperature begins to decrease after 27 months of operation. The temperature reaches its minimum level of 2.2˚C after 60 months of operation and then increases to the constant level of 2.3˚C after approximately 100 months of operation (Figure 8).
3.3. Results of Energy Calculations
The groundwater temperature remained constant during pumping and hence had no effects on the heating power reservoir in the detached house and apartment building scenarios. In the shopping centre scenario, the peak months for heating and cooling power are January and August, respectively. Groundwater would provide over 20% more heating power in January and over 25% more cooling power in August when compared to the reference year (Figure 9). The groundwater flux has been retained at the level of the reference year; only ΔTin-out has been changed according the modelled groundwater temperatures in Figure 9.
4.1. Power and Groundwater Flux Design Authors and Affiliations
Specific information on the heating and cooling power demands of buildings provides a possibility to accurately design groundwater pumping and energy transfer systems. Accurate planning adjusts the size of the energy dis-
tribution system, e.g. the nominal power of the heat and water pump, which will optimize the building and electricity costs of the project. Table 2 indicates that dimensioning a heat pump to cover 50% to 60% of the peak design power of a building provides 97.5% to 98.5% of the total heating energy with the given groundwater flux values. Holopainen et al.  reported that the lowest lifecycle cost will be achieved if a heat pump is dimen-
Figure 7. Groundwater temperatures on the heating (a) and cooling (b) side in the shopping centre scenario between 1 to 10 years of operation.
Figure 9. Monthly percentual change in the groundwater energy potential compared to the reference year in the a) heating and b) cooling model for the shopping centre in selected years. 100% indicates the energy need of the reference year.
sioned to cover 50% of the peak design power of apartment buildings in Finland. Rosen et al.  stated that in Sweden, the economically most suitable option is to dimension heat pumps to cover 50% to 60% of the peak design power of individual houses. Holopainen et al.  and Rosen et al.  made theirs calculation for a closed loop geoenergy system, i.e. a system where energy is exchanged from the ground to the fluid inside the heat exchanger pipes. Our calculation confirmed this assumption for an open-loop system in Nordic countries. As our results are comparable to previous studies in the Nordic climate   , power calculations and hence groundwater pumping requirements are on a reliable level. The Nordic climate is characterised by significant air temperature differences between summer and winter. These temperature fluctuations cause changes in building energy demands and hence observed variations in the groundwater pumping rate. Even the ATES system cannot be designed for a relatively stable pumping scenario in the Nordic environment.
Comparing our results to those of Banks  and Ferguson and Woodbury  , who investigated the cooling effects of buildings on groundwater, it appears that the thermal effect of groundwater energy utilisation is less harmful when more heating than cooling power is needed in buildings. Groundwater cooling can reduce microbial growth in groundwater. In general, warm groundwater provides more a suitable environment for harmful thermophilic microbes such as faecal bacteria than cool groundwater  . Arola and Korkka-Niemi  reported elevated groundwater temperatures due to urbanisation in Finland. Groundwater energy utilisation in urbanised areas restores the elevated temperature and hence changes the thermal environment of aquifers towards natural conditions. Cool groundwater is also a benefit if groundwater is distributed to the communal water system. However, cooled groundwater can change the natural vegetation on groundwater discharges areas and may consequently form threat to endangered species which are protected by Directive 2006/118/EU in the EU. Remediation of light non-aqueous phase liquids (LNAP) chlorinated solvents from groundwater  and seawater intrusion to coastal aquifers  can be accelerated with temperature variations. Heating or cooling power can be extracted from flowing groundwater similarly with remediation, and hence increase the sustainability of groundwater use.
4.2. Energy Utilisation
Groundwater can form a thermally long-term renewable energy source for large buildings in areas of cold groundwater such as Finland. This research demonstrates that careful, interdisciplinary planning involving thermo-geologists and HVAC engineers can improve the sustainability and economic viability of geothermal energy utilization. Accurate planning reduces the environmental risks and the overall economics of energy system could be improved. It can be estimated that groundwater could be effectively utilized until the temperature of approximately 4˚C. This research also indicates that the long-term maximum energy potential of groundwater can be estimated when the natural groundwater temperature, geological and hydrogeological conditions of the aquifer and energy requirements of buildings are precisely known. The ATES system generates different thermal regimes for an aquifer and needs more detailed system planning than groundwater utilization for heating or cooling only.
Groundwater energy utilization reduced the groundwater temperatures in the groundwater flow direction. In the Nordic environment, the groundwater temperature decreases due to the significantly larger heating than cooling energy requirements. In glaciofluvial esker formations, the size of groundwater thermal plume is dependent on site-specific thermo- and hydro-geological factors and can extend a few hundreds of meters from injection wells. Groundwater energy utilization may also have environmentally beneficial side effects; in urbanized areas, energy use could reduce the groundwater temperature to its natural level, the remediation of LNAP and/or seawater contaminated groundwater could be accelerated with the use of a heat pump or heat exchanger and replacing conventional oil heating systems with groundwater energy utilization will reduce the soil and groundwater contamination risk. The methods and results of this study represent well the Finnish geological environment and are applicable to similar environments globally.
The authors wish to thank Kirsti Korkka-Niemi, Veli-Pekka Salonen, Pirjo Tuomi and Jukka Hellen for their support and advice and Roy Siddall for language revision. Special thanks go to colleagues at the Turku office of Golder Associates Oy.
This research was funded by Golder Associates Oy and the K.H. Renlund Foundation.
 European Commission (2011) Communication from the Commission to the European Parliament, the Council, the European Economic and Social Committee and the Committee of the Regions. Energy Roadmap 2050, European Commission. Brussels.
 Banks, D. (2009) Thermogeological Assessment of Open-Loop Well-Doublet Schemes: A Review and Synthesis of Analytical Approaches. Hydrogeology Journal, 17, 1149-1155.
 Arola, T., Eskola, L., Hellen, J. and Korkka-Niemi, K. (2014) Mapping the Low Entalphy Geothermal Potential of Shallow Quaternary Aquifers in Finland. Geothermal Energy, 2, 9.
 Parnachev, V.P., Banks, D., Berezovsky, A.Y. and Garbe-Schonberg, D. (1999) Hydrochemical Evolution of Na-SO4-Cl Groundwaters in a Cold, Semi-Arid Region of Southern Siberia. Hydrogeology Journal, 7, 546-560.
 Rosen, B., Gabrielsson, A., Fallsvik, J., Hellstrom, G. and Nilsson, G. (2001) System for varme och kyla ur mark—En nulagesbeskrivning (Systems for Heating and Cooling from the Ground—A Status Report). Varia 511, Statens Geotekniska Institut, Lindkoping.
 McKenzie, J.M., Voss, C.I. and Siegel, D.I. (2007) Groundwater Flow with Energy Transport and Water Ice-Phase Change: Numerical Simulations, Benchmarks, and Application to Freezing in Peat Bogs. Advances in Water Resources, 30, 966-983.
 Tietavainen, H., Tuomenvirta, H. and Venalainen, A. (2010) Annual and Seasonal Mean Temperatures in Finland during the Last 160 Years Based on Gridded Temperature Data. International Journal of Climatology, 30, 2247-2256.
 Backman, B., Lahermo, P., Vaisanen, U., Paukola, T., Juntunen, R., Karhu, J., Pullinen, A., Rainio, H. and Tanskanen, H. (1999) Geologian ja ihmisen toiminnan vaikutus pohjaveteen. Seurantatutkimuksen tulokset vuosilta 1969-1996 (The Impact of Geology and Human Activities on Groundwater. The Results of the Follow-Up Investigation for the Years 1969-1996). The Geological Survey of Finland, 261.
 Oikari, H. (1981) Pohjaveden lampotila Etela-ja Keski-Suomessa vuosina 1975-1978 (Groundwater Temperature in Southern and Central Finland in 1975-1978). Master Thesis, University of Helsinki, Helsinki, 65 p.
 Kerl, M., Runge, N., Tauchmann, H. and Goldscheider, N. (2012) Hydrogeologisches Konzeptmodell von München: Grundlage für die thermische Grundwassernutzung (Conceptual Hydrogeological Model of the City of Munich, Germany, as a Basis for Geothermal Groundwater Utilisation). Grundwasser, 17, 127-135.
 Menberg, K., Blum, P., Shaffitel, A. and Bayer P. (2013) Long-Term Evolution of Anthropogenic Heat Fluxes into a Subsurface Urban Heat Island. Environmental Science & Technology, 47, 9747-9755.
 Allen, A., Milenic, D. and Sikora, P. (2003) Shallow Gravel Aquifers and the Urban “Heat Island” Effect: A Source of Low Enthalpy Geothermal Energy. Geothermics, 32, 569-578.
 Arola, T. and Korkka-Niemi, K. (2014) The Effect of Urban Heat Islands on Geothermal Potential: Examples from Quaternary Aquifers in Finland. Hydrogeology Journal, 22, 1953-1967.
 Zhu, K., Blum, P., Ferguson, G., Balke, K.-D. and Bayer, P. (2010) The Geothermal Potential of Urban Heat Islands. Environmental Research Letters, 5, Article ID: 044002.
 Haehnlein, S., Bayer, P. and Blum P. (2010) International Legal Status of the Use of Shallow Geothermal Energy. Renewable and Sustainable Energy Reviews, 14, 2611-2625.
 Bonte, M., Stuyfzand, P., Hulsmann, A. and van Beelen P. (2011) Underground Thermal Energy Storage: Environmental Risks and Policy Developments in the Netherlands and European Union. Ecology and Society, 16, 22.
 Brielmann, H., Griebler, C., Schmidt, S.I., Michel, R. and Lueders T. (2009) Effects of Thermal Energy Discharge on Shallow Groundwater Ecosystems. FEMS Microbial Ecology, 68, 273-286.
 De Keuleneer, F. and Renard, P. (2015) Can Shallow Open-Loop Hydrothermal Well-Doublets Help Remediate Seawater Intrusion? Hydrogeology Journal, 23, 619-629.
 Zuurbier, K.G., Hartog, N., Valstar, J., Post, V.E.A. and van Breukelen, BM. (2013) The Impact of Low-Temperature Seasonal Aquifer Thermal Energy Storage (SATES) Systems on Chlorinated Solvent Contaminated Groundwater: Modeling of Spreading and Degradation. Journal of Contaminant Hydrology, 147, 1-13.
 Paalijarvi, M. and Valjus, T. (2014) Karhinkankaan ja Sivakkokankaan pohjavesialueiden geologisen rakenteen selvitys 2009-2012 (Geological Structure Investigation of Karhinkangas and Sivakkokangas Aquifers in 2009-2012). The Geological Survey of Finland, Kokkola.
 Kalamees, T., Jylha, K., Tietavainen, H., Jokisalo, J., Ilomets, S., Hyvonen, R. and Saku, S. (2012) Development of Weighting Factors for Climate Variables for Selecting the Energy Reference Year According to the EN ISO 15927-4 Standard. Energy and Buildings, 47, 53-60.
 Koivuniemi, J. (2005) Lampiman kayttoveden mitoitusvirtaama ja lampotilakriteerit veden mikrobiologisen laadun kannalta kaukolammitetyissa asuinrakennuksissa [Dimensioning Standards of Domestic Hot Water Flow and Temperature Criteria according to Microbiological State of the Water in District Heated Houses]. Master Thesis. Helsinki University of Technology, Espoo.
 McDonald, M.G. and Harbaugh, A.W. (1988) A Modular Three-Dimensional Finite-Difference Ground-Water Flow Model. U.S. Geological Survey Techniques of Water-Resources Investigations, Book 6, Chap. A1, 586 p.
 Zheng, C. and Wang, P.P. (1999) MT3DMS, A Modular Three-Dimensional Multi-Species Transport Model for Simulation of Advection, Dispersion and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User’s Guide. U.S. Army Engineer Research and Development Center Contract Report SERDP-99-1, Vicksburg, 202 p.
 Thorne, D., Langevin, C.D. and Sukop, M.C. (2006) Addition of Simultaneous Heat and Solute Transport and Variable Fluid Viscosity to SEAWAT. Computers & Geosciences, 32, 1758-1768.
 Paalijarvi, M. and Okkonen, J. (2014) Karhinkankaan ja Sivakkokankaan vedenhankintatutkimukset ja pohjavedenvirtausmallinnus 2011-2014 (Groundwater Investigations and Groundwater Flow Model of Karhinkangas and Sivakkokangas Aquifers during the Period 2011-2014). Technical Report, Geological Survey of Finland, Kokkola, 64 p.
 Clauser, C. and Huenges, E. (2013) Thermal Conductivity of Rocks and Minerals. In: Ahrens, T.J., Ed., Rock Physics & Phase Relations: A Handbook of Physical Constants, American Geophysical Union, Washington DC, 105-126.
 Rau, G.C., Andersen, M.S. and Acworth, R.I. (2012) Experimental Investigation of the Thermal Dispersivity Term and Its Significance in the Heat Transport Equation for Flow in Sediments. Water Resources Research, 48, W03511.
 Andersson, O. (1994) Aquifer Thermal Energy Storages in Sweden—Experiences So Far and Market Potential. In: Kangas, M.T. and Lund, P.T. Eds., Thermal Energy Storage: Better Economy, Environment, Technology. Proceedings of the 6th International Conference on Thermal Energy Storage, Calorstock’94, Vol. 2, 22-25 August 1994, Espoo.
 Cruickshanks, F. and Adsett, E. (1994) Sussex Health Centre Aquifer Thermal Energy Storage. In: Kangas, M.T. and Lund, P.T., Eds., Thermal Energy Storage: Better Economy, Environment, Technology. Proceedings of the 6th International Conference on Thermal Energy Storage, Calorstock’94, Vol. 1, 22-25 August 1994, Espoo.
 Bayer, P., Saner, D., Bolay, S., Rybach, L. and Blum, P. (2011) Greenhouse Gas Emission Savings of Ground Source Heat Pump Systems in Europe: A Review. Renewable & Sustainable Energy Reviews, 16, 1256-1267.
 Saner, D., Juraske, R., Kübert, M., Blum, P., Helweg, S. and Bayer, P. (2010) Is It Only CO2 That Matters? A Life Cycle Perspective on Shallow Geothermal System. Renewable & Sustainable Energy Reviews, 14, 1798-1813.
 Holopainen, R., Vares, S., Ritola, J. and Pulakka, S. (2010) Maalammon ja viilennyksen hyodyntaminen asuinkerrostalon lammityksessa ja jaahdytyksessa (Exploitation of Ground Source Heat and Cooling in the Heating and Cooling of Apartment Buildings). VTT—Research Notes 2546, VTT Technical Research Centre of Finland, Espoo, 56 p.