Generalized Pan-European Geological Database for Shallow Geothermal Installations

The relatively high installation costs for different types of shallow geothermal energy systems are obstacles that have lowered the impact of geothermal solutions in the renewable energy market. In order to reduce planning costs and obtain a lithological overview of geothermal potentials and drilling conditions, a pan-European geological overview map was created using freely accessible JRC (Joint Research Centre) data and ArcGIS software. JRC data were interpreted and merged together in order to collect information about the expenditure of installing geothermal systems in specific geological set-ups, and thereby select the most economic drilling technique. Within the four-year project of the European Union’s Horizon 2020 Research and Innovation Program, which is known as “Cheap-GSHPs” (the Cheap and efficient application of reliable Ground Source Heat exchangers and Pumps), the most diffused lithologies and corresponding drilling costs were analyzed to provide a 1 km × 1 km raster with the required underground information. The final outline map should be valid throughout Europe, and should respect the INSPIRE (INfrastructure for SPatial InfoRmation in Europe) guidelines.


Introduction
Nowadays, geothermal energy is one of the most seminal renewable energy sources, due to its high potential and multiple uses.With the aim of both reducing the overall costs of shallow geothermal systems and improving their installation safety, a European project is recently undertaken, under the Horizon 2020 EU framework program for Research and Innovation.The "Cheap-GSHPs" (the Cheap and efficient application of reliable Ground Source Heat exchangers and Pumps) project (http://cheap-gshp.eu/)involves 17 partners among nine European countries: Belgium, France, Germany, Greece, Ireland, Italy, Romania, Spain, and Switzerland.The project is financed for four years, for the period between June 2015 and May 2019.In order to achieve the planned targets, a holistic approach is adopted, where all of the involved elements in shallow geothermal activities are integrated.Since the technical feasibility, total performance, and installation costs are affected enormously by underground properties, it is indispensable to have detailed information about these parameters.Based on information provided by the partners, drilling costs differ from country to country as a function of the maturity of the market, the soils and near-surface geology, the hydrology, and the competitiveness of the drilling companies.The general aim of the project is to reduce the overall cost up to 20%.In order to reach that goal, installation techniques on the construction site are reviewed and rectified where applicable.Also, the probe and backfilling materials will be examined in terms of cost-saving adjustments, with at least the same or better performance than standard materials.In order to reduce the drilling costs, new helicoidal and coaxial GSHEs (ground source heat exchangers) are developed within the project.These types of GSHE are expected to reduce installation costs, since they can be installed at a much shallower depth than standard double-U probes.As for double-U probes, the drilling costs reach almost 40% of total installation costs of a new geoexchange system (Figure 1) [1].Belgium, France, Germany, Greece, Ireland, Italy, Romania, Spain, and Switzerland.The project is financed for four years, for the period between June 2015 and May 2019.In order to achieve the planned targets, a holistic approach is adopted, where all of the involved elements in shallow geothermal activities are integrated.Since the technical feasibility, total performance, and installation costs are affected enormously by underground properties, it is indispensable to have detailed information about these parameters.Based on information provided by the partners, drilling costs differ from country to country as a function of the maturity of the market, the soils and near-surface geology, the hydrology, and the competitiveness of the drilling companies.The general aim of the project is to reduce the overall cost up to 20%.In order to reach that goal, installation techniques on the construction site are reviewed and rectified where applicable.Also, the probe and backfilling materials will be examined in terms of cost-saving adjustments, with at least the same or better performance than standard materials.In order to reduce the drilling costs, new helicoidal and coaxial GSHEs (ground source heat exchangers) are developed within the project.These types of GSHE are expected to reduce installation costs, since they can be installed at a much shallower depth than standard double-U probes.As for double-U probes, the drilling costs reach almost 40% of total installation costs of a new geoexchange system (Figure 1) [1].Since it is very important to keep drilling costs low from an economic point of view, the drillers and planners need to have access to reliable data about the underground.This is an important challenge, as stratigraphy can be widely varied between different locations.Therefore, it is recommended to choose the drilling technique based on the physical properties of the underground to avoid extra financial burden.
To extrapolate these issues on a European scale, it is also important to provide a homogenous data set that is valid within all of the participating project partner countries to provide local predictions about expenditure.Additionally, the resolution of underground specific data sets has to be accurate.Soil and rock properties can change very rapidly within a small area.Within this project, newly developed heat basket-type GSHEs should reach around 15 m depth [2].Hence, the challenge is to map small-scale, eventually vertical, underground variability versus a European-wide demand of data, in combination with relevant data for the improvements of the economic factors of the installation of shallow geothermal heating and cooling systems.At a later stage of the "Cheap-GSHPs" project, the final map will be implemented into the DSS (Decision Support System) of the project's homepage.Therefore, interested users can select their property's location, and the DSS will pick up the required information for its deciding algorithm.The system will recommend the best combination of GSHE, heat pump, and drilling technique, through also taking into account building-specific properties (size, insulation, materials) and the best cost-benefit ratio for stakeholders.Since it is very important to keep drilling costs low from an economic point of view, the drillers and planners need to have access to reliable data about the underground.This is an important challenge, as stratigraphy can be widely varied between different locations.Therefore, it is recommended to choose the drilling technique based on the physical properties of the underground to avoid extra financial burden.
To extrapolate these issues on a European scale, it is also important to provide a homogenous data set that is valid within all of the participating project partner countries to provide local predictions about expenditure.Additionally, the resolution of underground specific data sets has to be accurate.Soil and rock properties can change very rapidly within a small area.Within this project, newly developed heat basket-type GSHEs should reach around 15 m depth [2].Hence, the challenge is to map small-scale, eventually vertical, underground variability versus a European-wide demand of data, in combination with relevant data for the improvements of the economic factors of the installation of shallow geothermal heating and cooling systems.At a later stage of the "Cheap-GSHPs" project, the final map will be implemented into the DSS (Decision Support System) of the project's homepage.Therefore, interested users can select their property's location, and the DSS will pick up the required information for its deciding algorithm.The system will recommend the best combination of GSHE, heat pump, and drilling technique, through also taking into account building-specific properties (size, insulation, materials) and the best cost-benefit ratio for stakeholders.

Standard Drilling Techniques and Developments
The most common drilling techniques are auger, rotary, and down-the-hole (DTH) hammer drilling.Auger is usually used for shallow drilling in clayey geological conditions with intercalations of sand and/or gravel layers.The spiral of the auger promotes the cuttings to the surface.The geological conditions define the drilling bit at the bottom of the auger.In primary clayey conditions, it is more feasible to use a drilling bit that executes sharping and cutting (e.g., drag bit, chevron bit).When using auger drilling in areas with more consolidated material, the drilling bit should be equipped with cutting teeth that have a more rupturing impact [3].
Another very common drilling technique is rotary drilling, which is based on a rotating drill stem and a drilling bit at the end of the drill rods.The drill rods are flushed with a drilling fluid, which can be water, water-based mud, air, or a mixture of air, water, and a foaming agent.The drilling fluid is ejected through nozzles at the drilling bit, and is recovered between the drill stem and borehole wall.The promoted cuttings have to settle out in a mud tank before the circulation system recirculates the fluid.The settled-out cuttings are declared hazardous waste, and have to be disposed of properly.However, rotary drilling also uses air as a drilling fluid if drilling in consolidated and/or rocky formations.
In general, there are several bits available for rotary drillings, whereas the most regular ones are tricone roller bits.Compared with other roller bits, they drill faster, with only small deviations along the vertical axis.Depending on a low, moderate, or high level of rock hardness, different designs are used for the cutting teeth, the angle of the cones compared with the vertical axis, the offset of the cones, and the dimensioning and robustness of the different bearings [3,4].
The DTH hammer system is used especially for hard, consolidated sediments and/or rocks, and provides a faster and more economical way to penetrate the underground compared with conventional drilling techniques.This technique combines hammering with the rotation of the bit.The pneumatic-driven hammer is located at the bottom of the drill stem.The efficiency of the DTH hammer depends on the air pressure that the above-ground compressor is able to supply.The compressed air is also used as drilling fluid, and promotes the drill cuttings to the surface.Small quantities of water and foaming agent can be added to the compressed air, using air water mist as drilling fluid.If dry compressed air is used as drilling fluid, a mud tank and waste disposal do not need to be provided on site [5].
Project participant HYDRA S.R.L. developed standard, easy drill technology that consists of high-tech drilling equipment with special drilling casings as drilling rods, coupled with a particular extractable drill bit that allows drilling boreholes with a diameter between 101-152 mm.The casing is designed to play a dual role: drilling tool, and casing to prevent the collapse of the borehole.This double function has been possible due to the special extractable drilling bit that can be removed at the end of the perforation, leaving the hollow passage for geothermal probes.To accomplish the bit extraction operation, a particular tool called a fishing tool is necessary.The fishing tool is connected to the winch of the machine, and is then dropped down into the hole.Then, it automatically connects to the drill bit, unlocking it from the casing.Rewinding the winch, the fishing tool will drag the drill bit out of the hole, leaving a hollow hole in the ground.Compared with other traditional drilling systems, a standard easy drill can lead to a cost reduction in terms of time, as fewer operations need to be completed.
Within the project, a new designed shallow heat exchanger (diameter 260-275 mm) was developed that allows lower drilling depths around 15 m, which reduces the time required for drilling and the amount of tools (casings, rods) used.At the same time, the helicoidal heat exchanger has a smaller diameter compared with the standard product (e.g., REHAU Helix diameter = 350 mm).This lowers the fuel consumption, as less material has to be removed.
To install the new designed GSHE, the cost-saving standard easy drill technology was modified to the enlarged easy drill technology (EEDT) within this project.The new design consists of 1500-mm long tubes (Figure 2a) with an external diameter of 356 mm.On the outer surface of the tube, a metal spiral with an external diameter of 450 mm has been welded.A drag/chevron-type drill-bit-to-loose has been designed with a low-cost manufacturing approach, and an unlocking system has been designed to unfasten the bit at the end of the drilling (Figure 2b) [6].
Geosciences 2018, 8, x FOR PEER REVIEW 4 of 15 tube, a metal spiral with an external diameter of 450 mm has been welded.A drag/chevron-type drill-bit-to-loose has been designed with a low-cost manufacturing approach, and an unlocking system has been designed to unfasten the bit at the end of the drilling (Figure 2b) [6].

Data Acquisition
The main idea was to identify a European-wide valid data set with the best possible resolution.The data acquisition should provide freely available and digital data with an opportunity for modifications.The working base was selected from several data sources that are on the market with different data quality levels.The first data set was from the European 1G-E project OneGeology-Europe, which created a harmonized data model of geological maps with a scale of 1:1,000,000.The project followed the INSPIRE (INfrastructure for SPatial InfoRmation in Europe) guidelines, but a small amount of countries did not participate [7], which limited the level of accuracy and comprehensiveness.Another potential data set dealing with thermal conductivities of the underground would have been the outcomes of the European ThermoMap project [8].Within this project, a pan-European outline map providing the very shallow geothermal potential (vSGP) was developed and expressed within ThermoMap-MapViewer [9].The pedological data set for creating the pan-European outline map was provided by the ESDAC's (European Soil Data Centre) TXSRFDO (Dominant surface textural class) for STU (soil typological unit) or by national datasets for different test areas on a small-scale level.[10] The last data sets in the collection were directly from the JRC (Joint Research Centre) of the ESDAC, which is organized under the umbrella of the European Commission.The ESDAC provides several thematic maps for soil-related data that cover Europe almost completely.There is the opportunity to access and download the European Soil Database (ESDB) v2.0 after registration.There are several groups within ESDB, which contain 73 soil-classifying attributes in total.Dominant value maps based on 1 km × 1 km raster data are selectable, as well as additionally corresponding purity and confidence level maps [10][11][12].The data covers the EU28 states (excluding Iceland and Cyprus) plus Switzerland, Norway, and the Balkans.As shallow geothermal systems could be installed within unconsolidated as well as consolidated material, the attribute group called 'parent material' was chosen.This attribute group provided information about the most common material in

Data Acquisition
The main idea was to identify a European-wide valid data set with the best possible resolution.The data acquisition should provide freely available and digital data with an opportunity for modifications.The working base was selected from several data sources that are on the market with different data quality levels.The first data set was from the European 1G-E project OneGeology-Europe, which created a harmonized data model of geological maps with a scale of 1:1,000,000.The project followed the INSPIRE (INfrastructure for SPatial InfoRmation in Europe) guidelines, but a small amount of countries did not participate [7], which limited the level of accuracy and comprehensiveness.Another potential data set dealing with thermal conductivities of the underground would have been the outcomes of the European ThermoMap project [8].Within this project, a pan-European outline map providing the very shallow geothermal potential (vSGP) was developed and expressed within ThermoMap-MapViewer [9].The pedological data set for creating the pan-European outline map was provided by the ESDAC's (European Soil Data Centre) TXSRFDO (Dominant surface textural class) for STU (soil typological unit) or by national datasets for different test areas on a small-scale level.[10] The last data sets in the collection were directly from the JRC (Joint Research Centre) of the ESDAC, which is organized under the umbrella of the European Commission.The ESDAC provides several thematic maps for soil-related data that cover Europe almost completely.There is the opportunity to access and download the European Soil Database (ESDB) v2.0 after registration.There are several groups within ESDB, which contain 73 soil-classifying attributes in total.Dominant value maps based on 1 km × 1 km raster data are selectable, as well as additionally corresponding purity and confidence level maps [10][11][12].The data covers the EU28 states (excluding Iceland and Cyprus) plus Switzerland, Norway, and the Balkans.As shallow geothermal systems could be installed within unconsolidated as well as consolidated material, the attribute group called 'parent material' was chosen.This attribute group provided information about the most common material in one spatial location.The attribute PAR-MAT-DOM1 is included in this group, which contains the SGDBE (Soil Geographical Database) codes and values for the "major group code for the dominant parent material of the STU (Soil Typological Units)".This major group code is called PARMADO1, and it is outlined in Table 1.The raster data were downloaded from the ESDB (https://esdac.jrc.ec.europa.eu/content/european-soil-database-v2-raster-library-1kmx1km) to set up a new geodatabase (*.gdb), and to build an ArcGIS project, focusing on geology, drillability, and shallow geothermal systems.Figure 3 shows an already existing working database for further creations of itemized maps covering almost completely the margins of the European Union.Some more peripheral areas, such as Iceland and Cyprus, are not covered within this database.
Geosciences 2018, 8, x FOR PEER REVIEW 5 of 15 one spatial location.The attribute PAR-MAT-DOM1 is included in this group, which contains the SGDBE (Soil Geographical Database) codes and values for the "major group code for the dominant parent material of the STU (Soil Typological Units)".This major group code is called PARMADO1, and it is outlined in Table 1.Anthropogenic deposits The raster data were downloaded from the ESDB (https://esdac.jrc.ec.europa.eu/content/european-soil-database-v2-raster-library-1kmx1km) to set up a new geodatabase (*.gdb), and to build an ArcGIS project, focusing on geology, drillability, and shallow geothermal systems.Figure 3 shows an already existing working database for further creations of itemized maps covering almost completely the margins of the European Union.Some more peripheral areas, such as Iceland and Cyprus, are not covered within this database.For the ArcGIS project, there is a single shapefile provided by the ESDAC that contains all 73 soil-classifying attributes with their unique codes and values.The cell size of the grid is maintained constant, with 1 km × 1 km.For the purpose of the Cheap-GSHPs project, only two codes were extracted from the shapefile for further processing: PARMADO and PARMADO1.

Data Modification
Within a first modification, defined classes were qualified after the level of rock hardness.These descriptions were done in close collaboration with a drilling machine producer and an applied drilling specialist that were participating within the Cheap-GSHPs project.This evaluation is essentially to choose the most practical and economical drilling technique, and to avoid exorbitant costs.In order to unify the major groups of PAR-MAT-DOM1, the sediment classes were structured after their assumed degree of hardness, as shown in Table 2. Within a certain major group of PAR-MAT-DOM1, the drilling time per meter (using auger, tricone, DTH hammer, etc.) is a function of the level of rock hardness.By applying this classification, predictions are easier to perform when making simplified statements about underground properties.Finally, the end user of the DSS should receive a three-color code that depends on the underground conditions.As a result of this advantage of knowledge, they are able to receive a more detailed offer from a drilling company.
To be more precise, the data set was extended in a second step.The attribute group PAR-MAT-DOM1 was originally built up out of 212 attributes and their corresponding codes from the attribute code PAR-MAT-DOM (Appendix A), the dominant parent material.There, the corresponding PARMADO code consisted of four digits, where the first defined the associated PARMADO1 code.As an example, the transition from PAR-MAT-DOM to PAR-MAT-DOM1 is explained in Table 3.
As the stability, and consequently the drillability, of a borehole depends on the underground's hardness and grain size, the first approach-to define all unconsolidated materials as slightly consolidated-is not detailed enough to declare drillability classes.In order to provide more exact information, the major groups #5-unconsolidated deposits (alluvium, weathering residuum, and slope deposits), #6-unconsolidated glacial deposits/glacial drift, #7-eolian deposits, and #9-anthropogenic deposits of PAR-MAT-DOM1 were discarded.Furthermore, these attributes were reinterpreted within the PARMADO codes #6000 and #9300 in order to discriminate between codes that contain mostly sandy, clayey, or gravely material.The classification is now focused on the material's grain size (distribution).Also, additional properties for unconsolidated materials such as water content, bulk density, thermal conductivity, and field capacity could be assigned more specifically to a certain material class.

Hardness Map
With the focus on the state of hardness, Table 2 was integrated into ArcGIS 10.3 software.The attribute tables of the corresponding shapefile were modified to receive a three-colored map, where red stands for intensively consolidated, orange stands for moderately consolidated and green stands for slightly consolidated lithology conditions.This provides a first hint about the drillability of the material, and which drilling technique will be the most suitable for a given area.A first stage drillability map is displayed in Figure 4.

Hardness Map
With the focus on the state of hardness, Table 2 was integrated into ArcGIS 10.3 software.The attribute tables of the corresponding shapefile were modified to receive a three-colored map, where red stands for intensively consolidated, orange stands for moderately consolidated and green stands for slightly consolidated lithology conditions.This provides a first hint about the drillability of the material, and which drilling technique will be the most suitable for a given area.A first stage drillability map is displayed in Figure 4.  Within the final map, which was used as a working base for the DSS of the Cheap-GSHPs project, unconsolidated sediments were distinguished between primarily gravely, sandy, or clayey material.The choice of the drilling equipment differs most from the lithology and its geotechnical properties, which are very disparate between the above-mentioned unconsolidated sediments.Further, according to VDI 4640 [13], the thermal conductivity values are different, and this crucially affects the heat extraction.Also, the type of GSHE is specified by the underground's parameter.For example, areal collector systems are usually not installed in an underground of igneous rock or metamorphic rock.However, probes or heat baskets can be installed through drillings.Therefore, the PAR-MAT-DOM attributes table was re-organized (Table 4) to provide a solution to discriminate these materials from each other, and further predict intending drilling activities at the best possible rate.The result is the attribute group FAU_PAR-MAT-CON, with its codes and values called PARMAFAU.The final attributes group FAU_PAR-MAT-CON contains eight classes: The PARMAFAU value "no information" represents mainly lakes and rivers, metropolises, and peripheral areas such as orogens and fjords.However, undefined anthropogenic material and waste is also added to this attribute.The values "igneous & metamorphic rocks" and "consolidated sedimentary rocks" remain unchanged, and keep their PARMADO codes."Unconsolidated material (undefined)" represents soft material, and does not differentiate between sand, clay, or gravel as the main component.All of the remaining codes and values of PARMADO were analyzed and redefined according to their grain size and expected occurrence, and were described best with a value of PARMAFAU.Finally the numerous code of PARMAFAU was transferred to the base raster grid with corresponding cell locations, and imported to ArcGIS to generate a new shapefile and produce the final version of the outline map.The final map for the dataset FAU_PAR-MAT-CON is displayed in Figure 5.

Discussion
As with every map, Figure 5 has certain limitations in accuracy when making predictions about the lithology in deeper depth regions.In fact, it cannot be ruled out that technical parameter of the lithology might change with the depth.If there are no data from outcrops, the declining accuracy has to be noted.In general, maps try to reflect a three-dimensional system on a two-dimensional image.The map provides only information about the topmost material, and does not provide information about its vertical extent.If lithology changes with depth, the drilling equipment has to be modified, corresponding to a change in hardness.Another parameter that would be very helpful for predicting heat extraction rates and defining drilling requirements is the hydrogeological information.Groundwater presence has a positive effect on the heat extraction rate of many soils and rocks, as thermal conductivity values increase with higher water content [14,15].Furthermore, drillers and planners have to be aware of the regional groundwater situation, as confined groundwater could lead to problems during the drilling process, which causes high costs for the client or drilling company.On the other side, several technical procedures have to be implemented when drilling in areas with groundwater present.During the approval process for geothermal applications, water protection areas have to be worked out to avoid grounds for rejection, as in many countries, drilling in such areas is forbidden by law.Another important geological aspect that is not covered by the FAU_PAR-MAT-CON data set is the appearance of swellable anhydrite.These sections can lead to extreme costly events of damage (e.g., Staufen im Breisgau, Germany), as anhydrite increases its volume by about 61% if it comes into contact with (ground) water [16].With this data set, which was generated within the Cheap-GSHPs project, an allocation of thermal conductivities for each FAU_PAR-MAT-CON class could be made according to VDI 4640, part 1, Table 1.However, this is only a rough overview, as the rock-specific thermal conductivity values cover a wide range, and the water saturation also significantly affects the thermal conductivity values.

Discussion
As with every map, Figure 5 has certain limitations in accuracy when making predictions about the lithology in deeper depth regions.In fact, it cannot be ruled out that technical parameter of the lithology might change with the depth.If there are no data from outcrops, the declining accuracy has to be noted.In general, maps try to reflect a three-dimensional system on a two-dimensional image.The map provides only information about the topmost material, and does not provide information about its vertical extent.If lithology changes with depth, the drilling equipment has to be modified, corresponding to a change in hardness.Another parameter that would be very helpful for predicting heat extraction rates and defining drilling requirements is the hydrogeological information.Groundwater presence has a positive effect on the heat extraction rate of many soils and rocks, as thermal conductivity values increase with higher water content [14,15].Furthermore, drillers and planners have to be aware of the regional groundwater situation, as confined groundwater could lead to problems during the drilling process, which causes high costs for the client or drilling company.On the other side, several technical procedures have to be implemented when drilling in areas with groundwater present.During the approval process for geothermal applications, water protection areas have to be worked out to avoid grounds for rejection, as in many countries, drilling in such areas is forbidden by law.Another important geological aspect that is not covered by the FAU_PAR-MAT-CON data set is the appearance of swellable anhydrite.These sections can lead to extreme costly events of damage (e.g., Staufen im Breisgau, Germany), as anhydrite increases its volume by about 61% if it comes into contact with (ground) water [16].With this data set, which was generated within the Cheap-GSHPs project, an allocation of thermal conductivities for each FAU_PAR-MAT-CON class could be made according to VDI 4640, part 1, Table 1.However, this is only a rough overview, as the rock-specific thermal conductivity values cover a wide range, and the water saturation also significantly affects the thermal conductivity values.

Conclusions
The freely available, INSPIRE conformal, data sets collected by the European Commission via the Joint Research Centre provide a powerful working base for different kind of underground-describing outline maps.The newly created dataset FAU_PAR-MAT-CON provides two major parameters.First, it provides the drillability and hardness of the ground, which are important for defining the major cost factor of the installation of shallow geothermal systems: the drilling costs.Second, in the next step, certain lithologies can be assigned to the literature values of thermal conductivities, which allows a classification in terms of the proposed heat extraction rate.The data sets of the JRC can also be very helpful in other areas of pan-European research as well; they are not only useful for shallow geothermal issues.All of the parameters are needed within the DSS tool of the Cheap-GSHPs project in order to provide the best results for the stakeholders' planning of a shallow geothermal system.Nevertheless, given that the data and maps described here are presented at a 1 km × 1 km scale, successfully planning and installing GSHE systems will require consultation with experts who are knowledgeable in the underground properties at the site being considered.Our data and maps are intended for use in the first stage in such planning.Notwithstanding the above, this work gives a first idea of what kind of natural circumstances could occur at a certain location, and whether these could cause any difficulties during the drilling process, which could affect the total cost of the installation.Additionally, local regulations and legislations should be considered before the planning phase.

Figure 1 .
Figure 1.European average costs of a standard geothermal installation outside the house (+in-house located heat pump [1]).

Figure 1 .
Figure 1.European average costs of a standard geothermal installation outside the house (+in-house located heat pump [1]).

Figure 5 .
Figure 5. Map for the dominant parent material for a classification due to hardness.The map allows assertions about the level of hardness and drillability.

Figure 5 .
Figure 5. Map for the dominant parent material for a classification due to hardness.The map allows assertions about the level of hardness and drillability.

Table 1 .
Soil Geographical Database (SGDBE) values PARMADO1: Major group code for the dominant parent material.

Table 1 .
Soil Geographical Database (SGDBE) values PARMADO1: Major group code for the dominant parent material.

Table 2 .
Assumed hardness levels of PARMADO1 code.

Table 3 .
Relationship between codes of PAR-MAT-DOM and PAR-MAT-DOM1.

Table 3 .
Relationship between codes of PAR-MAT-DOM and PAR-MAT-DOM1.