Article of the Month - January 2020
     | 
  
		Detection of Service Pipes and the Risk of 
		Collapsing Sinkholes at the Lake of Constance in Switzerland Using 
		Ground Penetrating Radar (GPR)   
		Edi Meier, Inma Gutiérrez, Marco Baumann, Max 
		Bosshard, Rainer Heeb, Switzerland 
		
			
				
				   | 
				
				  | 
				
				   | 
			
			
				| Edi Meier  | 
				Inma Gutiérrez  | 
				Max Bosshard  | 
			
		
		
			
			This article in .pdf-format 
			(14 pages)
		This paper was presented at the FIG Working Week 
		2019 in Hanoi, Vietnam. The paper describes how Ground Penetrating Radar 
		is used to document the present underground condition of the water front 
		of the harbour promenade in Arbon (Lake of Constance) in Switzerland. 
		This nondestructive method gives as a result images of the subsurface 
		structures as well as of installations like service pipes.
			
		
						SUMMARY
		After a severe winter in Switzerland, two sinkholes opened suddenly 
		on the much used water front of the harbour promenade in Arbon (Lake of 
		Constance), following a time span with intensive wind that brought big 
		waves to the shore. The public building authority is responsible for the 
		safety of the boardwalk. They decided not only to search for further 
		sinkholes but also to document the present underground condition. This 
		will be used for the forthcoming restauration work. In order to get a 
		precise underground documentation of the area a Ground Penetrating Radar 
		(GPR) survey was executed. This nondestructive method gives as a result 
		images of the subsurface structures as well as of installations like 
		service pipes. The underground structure of the port facilities is 
		clearly shown in the radar images (reflectograms). For the city council 
		of Arbon this radar documentation is a helpful tool for the emergency 
		decisions as well as for the planning of the future lakeshore 
		constructions. 
		1.    EXTREME CONDITIONS IN THE DOCK OF ARBON
		1.1 Use of harbour during summertime
		The city of Arbon on the Swiss shore of the southern Lake of 
		Constance is famous for the “Summerdays”-festival. The festival 
		activities with openair concerts and temporary funfair installations 
		extend from the fairground plaza along the entire esplanade (Fig. 1).
		
		
		Fig. 1): Summerdays-festival in Arbon, Lake of Constance, Switzerland 
		[3]. 
		1.2 Situation in winter 2018 
		During winter 2018 extreme climatic conditions prevailed in the area 
		of the Lake Constance. At the beginning of February the lake water level 
		was higher than the since 1881 measured data (Fig. 2). Driven by a 
		strong north wind big waves lashed to the shore of the promenade. With 
		temperatures below zero the whole area was covered by a dangerous ice 
		shield (Fig. 3 left). 
		One month later, when the ice melted away, the ground under the 
		fairground plaza and under the pedestrian promenade suddenly collapsed 
		(Fig. 3 right). Only the massive retaining wall did not tumble down and 
		formed a bridge (Fig. 4 left). Apparently, the poor condition of the 
		shore ramp and the high water level - in combination with the strong 
		north wind - was the cause for the soil erosion. The lakewater could 
		penetrate through openings and cracks in the wall (Fig. 4 right) and 
		loosened the ground when freezing. 
		The city council was worried and assumed more hidden holes in the 
		underground (Bosshard M., 2018). This would be a major risk for the 
		coming "summerdays"-festival. 
		
		
		Fig. 2): Water level fluctuation of the Lake Constance since 1881. 
		Red curve = curve for the year 2018, green curve = mean water level, 
		black curve = min. water level, blue curve = max. water level.
		
		
		Fig. 3): Ice shield on the shore in February (left) [4]. Collapsed 
		ground after the melting of the ice in March 2018 (right). 
		2.    GPR INVESTIGATION 
		2.1 Task
		In March 2018 GPR investigations were carried out. The goal was to 
		detect hidden large holes in the ground that could cave in. At the same 
		time, the city council wanted to get a documentation of the underground 
		conditions for the future shore restoration. The measurements were 
		executed on the pedestrian promenade, the fairground plaza and on two 
		associated shore ramps. Parallel lines were recorded per zone. The 
		"radargrams" represent either a profile in depth or a profile 
		perpendicular to the ramp. Additionally 43 parallel profiles with a 
		spacing grid of 50 cm were executed, that allowed a 3-D analysis of the 
		underground. 
		
		
		Fig. 4): Collapsed site near the promenade (left). Holes between the 
		stones in the ramp (right). 
		2.2 Process of the measures
		Two areas were investigated: the pedestrian walkway and the 
		fairground plaza, both with associated ramp. The measured data were 
		collected in five parallel rows. During the first measurements, there 
		were still a few icy places on the shore. A GPR system with shielded 
		250- / 700-MHz double antennas was used (Fig. 5 above).
		During the entire series of measurements, the data was collected 
		simultaneously with the 700 MHz antenna as well as with the 250 MHz 
		antenna - in other words - twice. For a detailed analysis of the first 
		three meters depth, the data of the 700 MHz antennas were evaluated. In 
		order to locate deeper objects, the data of the 250 MHz antennas were 
		also analysed. In the upper zone, that is to a depth of three meters and 
		thus close to the surface, “radargrams” with a higher resolution are 
		possible using the 700 MHz antennas.
		A week later the two profiles on the shore ramp were executed. For 
		this purpose, two carriages were built, a guiding carriage for the quay 
		wall and a carriage for the ramp wall. They were fixed together with a 
		rope-system. 
		In the meantime the ice had melted away, but the sloping shore was 
		still very slippery. Because of the high waves and heavy rain, the field 
		work was done wearing a dry suit. 
		The guiding carriage was equipped with a precision trigger wheel. 
		This released a measurement every 5 cm interval. The carriage with the 
		radar antenna on the ramp wall was pulled with the same speed as the 
		guiding carriage. Thus, the metering of the ramp profiles is always 
		comparable to that of the profile F1 on the quay wall.
		The surface of the stones on the ramp wall was relatively uneven and 
		some obstacles had to be overcome. Because of the widely supported 
		suspension of the radar antenna, these obstacles gave less error signals 
		in the reflectograms than expected. 
		
		
		Fig. 5): GPR shots along the promenade: Three profiles on top of the 
		path, view towards the fairground plaza (top picture), two profiles on 
		the ramp wall, view in the opposite direction (pictures below). 
		2.3 Measuring method
		A radar pulse is emitted from the transmission antenna into the 
		ground, which is reflected at interfaces of rock packages or layers in 
		the underground (Fig. 6). Thereafter, the pulse returns to the receiving 
		antenna where it is recorded. The time that elapsed between transmission 
		and reception of the pulse provides information about the depth of the 
		object (reflector). The strength of a reflection serves as a clue to the 
		electrical properties of the material that caused the reflection.
		
		
		Fig. 6): Principle of GPR measurement (Meier E. et al., 2002)
		A measurement is always recorded as a profile. To do this, the echoes 
		of the transmitted radar pulses are recorded point by point along a line 
		(profile line). The reflection pattern is used to interpret the 
		condition of the subsoil. The resolution of the achieved GPR mapping 
		structures depends on the antenna frequency and the scan rate along the 
		profile line.
		2.4 Datapresentation and Dataanalysis 
		On top of the “radargrams” the length scale is given in meters, on 
		the left the signal propagation time in nanoseconds and on the right the 
		depth scale in meters. The depth specification is calculated from the 
		signal propagation time and the signal speed. The value of the depth 
		specification is only to be regarded as a guideline value. The zero 
		point of the vertical axes corresponds to the terrain surface. 
		The “radargrams” are checked for conspicuous structures, that is 
		diffraction patterns, so-called hyperbolas. They are caused by locally 
		limited strong reflectors such as pipes, manholes, cavities, etc. In 
		addition, linear structures are analysed, that can provide information 
		on geological stratifications in the underground. In order to specify 
		the depth of located objects or structures, the knowledge of the 
		propagation velocity of radar waves is necessary. For this purpose a 
		guideline value of v = 0.1 m/ns has been used in the “radargrams”. For 
		an exact depth specification, the depth scale would have to be 
		calibrated with the location of known pipes etc. or with data from 
		drilling holes. 
		2.5 Examples from the measurements on the promenade 
		The “radargrams” recorded on the pedestrian path show clear 
		structures up to about 4 m depth. Profile F3 was taken at a distance of 
		2 meters parallel to the quay wall. Figures 7 and 8 show the same 
		section as profile F3, but were recorded with two different frequencies. 
		While the 700 MHz antenna shows smaller structures (Fig. 7), the 250 MHz 
		antenna achieves greater depth penetration (Fig. 8). 
		
		
		Fig. 7): Part of profile F3, taken with the 700 MHz antenna.
		On both pictures is the same structure clearly visible: a structure 
		up to 1 m depth, that slopes to the left, is seen on the left half on 
		the “radargram” as well as the departure tracks of the construction 
		machines on the right half. 
		
		
		Fig. 8): Part of profile F3, taken with the 250 MHz antenna. 
		When interpreting these images, it must be ensured that they are not 
		to be regarded as a direct image of the underground. For example, a 
		highly reflective object can not only be seen vertically below the radar 
		antenna, but it is already recorded from an oblique angle from the 
		receiving antenna. In the “radargram” it is displayed perpendicular to 
		the profile. This leads to the above-mentioned “bow shape” (hyperbolas). 
		This can be seen on the “radargram” of profile F1 taken on the Quai 
		wall: there is a regular pattern of hyperbolas (Fig.9). Four transverse 
		reinforcing bars per meter create this pattern. The hyperbolas in the 
		“radargram” simulate a continuous arcuate support structure, that 
		certainly does not correspond to the reality. In fact, the “reinforcing 
		bars” are likely to be only a few cm in width.
		
		
		Fig. 9) Part of profile F1 on the quay wall, taken with 700 MHz 
		antenna. Reinforcing bars draw hyperbolas in the “radargram”. At the 
		apex of each hyperbola is a “reinforcing bar” in transverse direction to 
		the profile. 
		To correct these "artefacts" the geophysical processing method 
		"migration" is used. This method recalculates the hyperbola branches to 
		a circle, the actual origin of this object. However, migration only 
		makes sense if the hyperbolas cover underlying structures or if steep 
		structures are present. We recommend to keep the hyperbolas in the 
		“radargram” and to analyse them. The position and depth of the object 
		can be determined just as well without migration, because the objects 
		are always at the vertex of the hyperbola.
		In addition to cavities large stones or rubble can produce comparable 
		hyperbolas, as both air and stones have a different electrical 
		conductivity from the surrounding material. It is a fact, that the 
		electrical conductivity of the ground changes greatly at the apex of a 
		hyperbola. 
		The depth and size of an object detected in the “radargrams” can be 
		determined with a model calculation. For this purpose, a reflection 
		hyperbola is calculated and placed over the “radargram”. The shape of 
		the hyperbola is calculated from the “travelled” distance, the diameter 
		of the object and the mean velocity of the radar waves in the overlying 
		rock layer. In addition, the hyperbola shape is influenced by the angle 
		at which an object is cut. In our analysis we used in all “radargrams” 
		for the conversion of the run-time in depth the value of 0.1 m/ns. 
		Depending on the material that lies above the object, this depth 
		conversion must be adjusted.
		As an example, the analysis of the object in 2 m depth at profile 
		position 31 m in the radar profile F3 is shown in Fig.10. 
		
		
		Fig. 10): Model hyperbola adapted to the object in Fig. 8. The depth 
		scale on the right corresponds to the calculated speed of 0.07 m/ns 
		adapted in the model. 
		The model hyperbola gives a speed of 0.07 m/ns. The depth scale must 
		therefore be corrected at this point by a factor of 0.7. This slower 
		speed indicates that clayey material has been deposited over the object. 
		The model calculation of the value of the diameter of the object gave a 
		result of 2.5 m. However, it could also be an elongated object (pipe) 
		cut at a shallow angle. This would result in a comparable flat model 
		hyperbola. Since the object is not seen in the neighbouring profiles, 
		this interpretation is much more likely. A cavity washed out by the lake 
		can thus be excluded. 
		2.6 Analysis of the 3-D measurements on the fairground plaza 
		The 50cm line grid recorded on the fairground plaza allows both a 3-D 
		interpolation of the ground in the longitudinal and transverse 
		directions as well as a plane representation with depth (see Fig. 11). 
		These pictures are “put together” in the form of a film document. It can 
		be played on any video player program and stopped at any depth. In the 
		plane view, in 1.1 m depth, lines (=pipes) can be seen, which are shown 
		as a “Y-shape”. At a depth of 2 meters you can see a dense track in the 
		middle of the area and at the margin a deposit of landfill material is 
		assumed, because all the area is artificially filled up. On the basis of 
		these findings, dredging slits were executed to corroborate the 
		interpretation.
		
		
		
		
		Fig. 11): 3-D images: In the plane view, lines are visible as a 
		„Y-shape“ in 1.1 m depth. The lower picture shows the level at a depth 
		of 2 m. Here, especially on the left and right side of the area, there 
		are highly reflective zones. 
		The dredging slits to the left and right of the site did not show the 
		expected clayey landfill material. At first a large block of stone 
		appeared, which caused the hyperbola (Fig. 12, left and center). Below 
		the stone there followed potato-sized, round stones without finely 
		granulated material. Water immediately flowed through the large cavities 
		between the stones, (Fig. 12, right), which explains the change in the 
		reflection signals at a depth of 2 meters. This shows that there is a 
		rapid hydraulic connection to the lake water on the edge of the plaza 
		and that a lot of fine grained material from the original fill is washed 
		out. With great certainty it could be ruled out that the ground on the 
		plaza would suddenly collapse. At the most, slow subsidence can occur 
		when more fine grained material is washed out.
		
		Fig. 12): The Profile C7 has a distinct hyperbola (picture left). The 
		excavation brings out a large stone block (middle picture). The 
		remaining excavated material consists of round stones of 3 - 5 cm 
		diameter with little finely granulated material. The water flows 
		immediately (picture right).
		3.    CONCLUSION
		With the GPR method, the soil can be quickly “x-rayed”. However, the 
		interpretation is not always clear, as not the strength but the 
		electrical properties of the subsurface is shown. On the one hand, both, 
		air and large, dry stones will produce similar images. On the other 
		hand, water that penetrates into cavities can greatly modify the image. 
		Thus, even after 30 years of experience with GPR, we cannot rule out 
		misinterpretations.
		Probing with a dredge slit or with drill-holes has the advantage of 
		visually identifying the layer sequence at a particular location. But 
		using the GPR method will in short time and without destruction help you 
		to determine over the whole survey area, where you want to drill or 
		where to open up a dredge slit. 
		REFERENCES
		[1] Bosshard M., June 2018, Sanierung Ufermauer Arbon, unpublished 
		Protokoll Nr., Projekt Nr. 3100-0761, Arbon, Wälli Ingenieure
		[2] Meier E., Staubli P., Müller B. U., Stünzi J., Schubert E., 
		Dubois D., Juli 2002, „Georadar - der zerstörungsfreie Blick in den 
		Untergrund: Beispiel aus dem Naturschutzgebiet Zigermoss, Unterägeri/ZG 
		und der Deponie Riet Winterthur/ZH“, Bulletin für angewandte Geologie, 
		Volume 7, Nr. 1, S. 31-44, Winterthur.
		[3] Photo “Summerdays Arbon”
		
		https://www.summerdays.ch/galerien/208/
		[4] Photo “Eiswinter Arbon”
		
		http://www.arbon-online.ch/eiswinter2012/IMG_5623.html
		BIOGRAPHICAL NOTES
		Edi Meier is currently the managing director of the 
		engineering company Edi Meier + Partner AG, Winterthur, Switzerland. He 
		studied Geophysics at the Swiss Federal Institut of Technology, ETH 
		Zurich. Subsequently he worked as a manufacturer of seismic instruments 
		(Streckeisen Switzerland) for six years and founded his own engineering 
		company in 1987. His company is specialized in development and 
		construction of precision deformation measuring systems and in 
		geophysical services using Geoelectric and Georadar (GPR) tools.
		Inma Gutiérrez is a geologist and is currently 
		working at the engineering company Edi Meier + Partner AG, Winterthur, 
		Switzerland. She studied Geology at the University of Granada (Spain) 
		and completed a complementary Master in “Integrated water management” at 
		the University of Cádiz (Spain). After some experiences in the 
		geotechnical and water management branches, she is working since 2017 in 
		the geophysical branch, specializing in GPR.
		Marco Baumann is currently working as a Senior 
		Expert in the Section Water Resources, Office for the Environment OEN, 
		Canton TG, in Frauenfeld, Switzerland. After his Ph.D. Thesis in 1986 at 
		the Geological Institute, Swiss Federal Institut of Technology, ETH 
		Zürich, he worked in various engineering companies as a Geotechnical and 
		Scientific Consultant and Projectmanager and as a Geologist, 
		Hydrogeologist, Projectmanager and Senior Consultat in a geology 
		company. Since 1st August 1994 he is working in the Office for the 
		Environment OEN, Canton TG. First as head of section Water Resources 
		with the topics of integral water management, hydrological data, water 
		protection and water use and since 1st January 2015 with the added 
		topics of Natural hazards, flood protection, river engineering and as 
		expert in the competence center for use of geothermal energy.
		Max Bosshard studied civil engineering. After 
		graduating he worked as a research assistant at the Research Institute 
		for Hydraulic Engineering, Hydrology and Glaciology at the Swiss Federal 
		Institut of Technology, ETH Zurich. Subsequently, he joined as project 
		engineer some engineering companies in the Zurich area und worked in the 
		fields of tunneling, railway construction and hydraulic engineering. 
		Currently he is working at Wälli AG Ingenieure, Arbon, Switzerland. His 
		specialties at Wälli AG are constructions on the lake shore and in the 
		lake, such as port facilities, pipelines, footbridges and quay moorings 
		as well as river engineering and revitalization.
		 Rainer Heeb is a civil engineer FH. He 
		completed a postgraduate degree in business administration. Currently he 
		is the deputy head of the construction department of the city of Arbon, 
		Switzerland, and head of the civil engineering department.
		CONTACTS
		Edi Meier 
		Edi Meier + Partner AG
		SWITZERLAND
		Website: www.emp-winterthur.ch
		Inma Gutiérrez
		Edi Meier + Partner AG
		Website: www.emp-winterthur.ch
		Marco Baumann
		Kanton Thurgau
		Office fort he Environment
		Website: www.umwelt.tg.ch
		Max Bosshard
		Wälli AG Ingenieure
		Website: www.waelli.ch
		Rainer Heeb
		Stadt Arbon
		Website: www.arbon.ch