Ground Penetrating Radar ( GPR ) Investigations for Architectural Heritage Preservation : The Case of Habib Sakakini Palace , Cairo , Egypt

A comprehensive Ground Penetration Radar (GPR) investigations and hazard assessment for the rehabilitation and strengthening of Habib Sakakini’s Palace in Cairo is presented herein, which is considered one of the most significant architectural heritage sites in Egypt. The palace located on an ancient water pond at the eastern side of Egyptian gulf besiding Sultan Bebris Al-Bondoqdary mosque is a place also called “Prince Qraja al-Turkumany pond”. That pond had been filled down by Habib Sakakini at 1892 to construct his famous palace in 1897. The integrated geophysical survey of the palace allowed the identification of several targets of potential archaeological and geotechnical engineering interest buried in fill and silty clay in the depth range between 100 700 cm. the methodological development focused on Multi-Fold (MF) Ground Penetrating Radar (GPR) imaging and subsurface characterization based on integrated velocity and attenuation analysis. Eight hundred sqm of Ground penetration Radar (GPR) profiling have been conducted to monitor the subsurface conditions. 600 meters are made in the surrounding area of the Palace and 200 sqm at the basement. The aim is to monitor the soil conditions beneath and around the Palace and to identify potential geological discontinuities, or the presence of faults and cavities. A suitable single and dual antenna are used (500 100 MHZ) is used to penetrate the desired depth of 7 meters (ASTM D6432). The GPR is used also detect the water table. At the building basement the GPR is used to identify the foundation thickness and soil-basement interface. As well as the inspection of cracks in some supporting columns, piers and masonry walls. The GPR also was used to investigate the floors and ceilings conditions and structural mapping. The results were validated by the geotechnical and structural surveys. All these results together with the seismic hazard analysis will be used for the complete analysis of the palace in the framework of the rehabilitation and strengthening works foreseen in a second stage.


Introduction
The modern architectural heritage of Egypt is rich, and extensively variable.It covers all kind of monumental structures like palaces, public buildings, residential and industrial buildings, bridges, springs, gardens and any other modern structure which falls within the definition of a monument and belongs to the Egyptian cultural heritage.
The majority of these buildings are private properties, which are either part of a traditional settlement, or isolated buildings incorporated into the wider area (such as towers and konaks) or into a completely changed and entirely modern area.
Preservation of the architectural heritage is considered a fundamental issue in the life of modern societies.In addition to their historical interest, cultural heritage build-ings are valuable because they contribute significantly to the economy by providing key attractions in a context where tourism and leisure are major industries in the 3rd millennium.The need of preserving historical constructions is thus not only a cultural requirement, but also an economical and developmental demand.
The study of historical buildings and other structures must be undertaken from an approach based on the use of modern technologies and science.The final aim must be to select and adequately manage the possible technical means needed to attain the required understanding of the morphology and the structural behavior of the construction and to characterize its repair needs.Modern requirements for an intervention include reversibility, unobtrusiveness, minimum repair, and respect of the original construction, as well as the obvious functional and structural requirements.Restoration operations complying with these principles require a scientific, multidisciplinary approach that comprehends historical understanding, modern non-destructive inspection techniques, and advanced experimental and computer methods of analysis.
The archaeological subsurface represents a potentially difficult problem for imaging targets with ground penetrating radar (GPR) systems.Structures within the ground may be very steeply dipping or exhibit large local variations in strike causing serious migration problems.Furthermore, the target may be located in a medium that is electrically conductive (e.g., clay-rich or waterlogged soil) such that the GPR wavelet is rapidly attenuated and the signal-to-noise ratio (SNR) in the resulting dataset is diminished.Recently, techniques used routinely by the seismic industry for imaging complex, low-SNR targets have been imported to the field of GPR acquisition.Specifically, use of three-dimensional (3-D) migration has been shown to significantly improve the resolution of the GPR target in regions of structural complexity, whilst the multi-offset (MO) method has been employed to boost SNR over otherwise poorly defined targets [1].When these methods are combined, the potential improvement to the image of an archaeological target may be considerable.An integrated 3-D, MO GPR acquisition was performed over an archaeological target, over a Romno-British villa at Ground well Ridge in Spring 2006.
Ground penetrating radar (GPR) is extensively used for a variety of applications.Among all the high-resolution geophysical methods, GPR has proven to be the most suitable for detection of karstic cavities and sinkholes, in a wide range of soil and rock conditions [2].However, one of the main limitations of GPR is the exact determination of the mean velocity of the electromagnetic waves, which is a key datum to estimate the depth of penetration into the ground.

Ground Penetrating Radar Survey
The GPR functions by sending high frequency electromagnetic waves into the from a transmitter antenna.Some of these waves are reflected back to the surface as they encounter changes in the dielectric permittivity of the matrix through which they are travelling, and are detected by a receiver antenna.The amplitude and twoway travel time of these reflections are recorded on a portable computer.This information is then used to construct a two-dimensional plot of horizontal distance versus travel time.Data collected in the field are stored on a portable computer for later analysis.A more complete and technical discussion of the method can be found elsewhere [3,4].
The effectiveness of GPR is controlled by the local soil conditions.GPR is most effective in locating buried objects in a homogenous soil matrix with a high electrical resistance.GPR is least effective in a heterogeneous environment with high electrical conductivity.A heterogeneous environment contributes to signal scattering and can result in insufficient depth of penetration and a "noisy" reflection (poor signal to noise ratio).A conductive environment can seriously inhibit depth of penetration due to conductive losses.Conductive loss is the result of the electromagnetic wave creating a conductive current in the soil medium, resulting in signal attenuation.
Although GPR survey can be performed in a number of ways, the method employed in this survey involves dragging the transmitter and receiver antennas together over the ground, called fixed offset reflection mode.The transmitter emits pulses at regular intervals along a transect which are picked up by the receiver.A laptop computer controls data collection and displays the data as a two dimensional profile.
The GPR is able to detect subsurface features whose electrical properties contrast with those of the surrounding soil.The GPR can detect human burials in several ways.It may detect the disturbed soil of the grave shaft, or a break in the natural stratigraphy or soil profile [5].It may also detect the coffin, bones, clothes and other articles in the burial.Reflections may be caused by air voids within the skull [6] or the coffin.It has also been suggested that the decomposition of bones may leach calcium salts into the surrounding soil for many years, which may change the electrical properties of the soil, making it visible to the radar.
The integrated geophysical survey of the palace allowed the identification of several targets of potential archaeological and geotechnical engineering interest buried in fill and silty clay in the depth range between 100 -700 cm. the methodological development focused on Multi-Fold (MF) Ground Penetrating Radar (GPR) imaging and subsurface characterization based on integrated velocity and attenuation analysis.
Information about geometry and location of targets is integrated by attenuation and velocity of radar waves, which provide direct evidence of varying electromagnetic properties (conductivity, dielectric constant) and allow a detailed characterization of surveyed volume based on physical characteristic of soil strata, buried foundations and materials.Applications of combined techniques can be successful in preliminary exploration of pre-historical sites, see [7] and of archaeological sites of later periods.The contrast in physical properties of the materials and the sensitivity of the different methods play a key role in applications of geophysical methods to archaeological and engineering surveys and determine the probability of success.The GPR is extensively applied to study archaeological sites, see [8][9][10][11][12].

Ground Penetration Radar (GPR) Data Acquisition and Analysis
ited to a finite bandwidth, an appropriate use of band limiting filtering may improve signal-to-noise without significantly altering the data [13].Taking into account the information obtained from the amplitude spectra of the raw data, a common band pass filter of 250 MHz to 750 MHz was applied to the whole set of 500 MHz profiles, and 150 MHz -300 MHz band-pass filter applied to the whole set of 250 MHz profiles to improve the signal quality.
Mala X3M Radar System coupled to shielded antenna of 500 MHz and 250 MHz central frequencies was used to accomplish the survey.Distances along the surveyed lines were accurately recorded using a measuring tape.The profile spacing was assumed according to the available space in the Palace, and the trace interval of 0.02 meters for the 500 MHz antenna profiles.About 100 profiles were conducted at the site (Figures 1 and 2).The radar antenna was moved along the lines, and the 2D profiles of a large number of periodic reflections were generated, thus producing a profile of the subsurface structure, with a fixed gain.Special care was taken to avoid artificial high-frequency noise that could deteriorate the radar signal.Time windows of 60 ns and 170 ns were used during data acquisition in order to receive information from sources located as deep as possible.The layout of the conducted profiles in basement and first floors are shown in (Figures 1(a) and (b)).
Since the radar velocities and the dielectric properties of the studied materials were unknown, two different techniques were used to estimate the mean velocities of the electromagnetic waves: 1) determination of the velocity that produces the best fit between the measured GPR reflections and the dips and depths of the geological structures identified with seismic surveys (e.g. a lithological boundary resulting in contrasting physical properties); 2) determination of the mean velocity of the material directly related to the geometry of the hyperbolic reflections originated by point sources.The second method used in the present study where a test was carried out on a known depth objects at the site.All the data were processed, modeled and interpreted using the software Reflex_W.4[14].In all the profiles, the position of antenna is represented in the horizontal axis, whereas depth is depicted in the vertical one.An example of data processing results is shown in (Figure 3).

GPR Data Processing
The analysis of GPR data is carried out by processing the data using different gains and filtering techniques.Gain is a value, by which raw data are multiplied, to enhance low amplitude reflections.Signal amplitude commonly decreases exponentially at increasing travel times (greater depths below surface).This is compensated by designing a time gain that increases the signal strength at greater travel times.Filtering is the use of mathematical processing algorithms to "clean" noises from the data and/or enhance certain characteristics of the data.This process was specifically customized for each profile.Data processing routine includes background noise removal, time-zero corrections and band-pass filtering applied to the acquired data.Given that, most of the energy is lim-

Interpretation of Geophysical Survey Results
The results of geophysical surveys of archaeological sites are generally presented graphically.This is done because anomalies of cultural origin are generally recognized by their pattern, rather than by their numeric values alone.
When rendered graphically, we can better recognize cultural and natural patterns and visualize the physical phenomena causing the detected anomalies.Interpretation of survey data must be a cooperative process involving both archaeological geophysicists and archaeologists that are familiar with the specific cultural context of the site being studied.An understanding of the geological context of the survey area is also very impor-  8).About 30 GPR profiles conducted at the area surround the palace.The target of these profiles is to depict the subsurface condition concerning the layering, water table, disturbed and collapsed soils and foundation locations.The soil su ce extended to about 7 meters.The fill consists of silt, clays and some stone fragment in some localities.The silt become more clayey below 4 meter depth as indicated from the amplitude of radar reflection.Shallow water table might be due to the seepage or the past history of the lagoon area at the site of the Palace.Some        Several targets of potential engineering interest were identified from the integrated GPR and geotechnical sury data and interpretations. The soil column of the site comprises heterogeneous materials as a fill for abou  The shallow ground water table, about 0.5 -1 meter, has a serious influence on the rigidity of t imposes humidity to the foundations and walls of the Palace. The composition of silt and clay soil for great depth might af  Many fractures detected in walls and concrete slabs that extended to about 60 cm in the w  The foundation type is shallow strip and spread stone foundations at depth of 2.5 to 3 meters on concret raft or mat. The wall masonry construction system is multiple leaf masonry wal  The added court in the eastern side is concrete structure.

Figure 1 .
Figure 1.(a) GPR layout at basement and area around the palace; (b) GPR layout at first floor of the palace.

Figure 2 .
Figure 2. Investigations of the Sakakini palace in Cairo.

Figure 3 .
Figure 3.Effect of the processing sequence on the raw GPR data.
at the Basement of the Ab profiles conducted at the basement floor.Palace out 50 GPR The target of these profiles is to depict the subsurface condition concerning the layering, disturbed and col-Copyright © 2012 SciRes.OJG lapsed soils and foundation locations.Figures 4-7 show the radar profile No. 7, 8, 22 and 29 and shows a disturbed soil and the expected locations of the foundation at depth of about 3 meters.The soil succession shows a fill below the floor exte 7. Results of GPR at Palace Surrounding ing nded to about 5 to 7 meters.The fill consists of silt, clays and some stone fragment in some localities.The silt become more clayey below 4 meter depth as indicated from the amplitude of radar reflection.Shallow water table might be due to the seepage or the past history of the lagoon area at the site of the Palace.The locations of the interpreted profiles in the basement floor are shown in (Figure

Figure 4 .
Figure 4. Interpreted GPR profile No. 7 at the basement area showing the disturbed soil at depth about 3 m.

Figure 5 .
Figure 5. Interpreted GPR profile No. 8 at the baseme nt area showing the expected foundation at depth about 2.5 to 3 m.

Figure 6 .
Figure 6.Interpreted GPR profile No. 22 at the basement area showing the expected disturbed soil due to foundation at about 3 m depth.

Figure 7 .
Figure 7. Interpreted GPR profile No. 29 at the basement area showing the expected foundation at depth about 3 m.

Figure 8 .
Figure 8. Interpreted locations of interrupted s (red

Figure 9 .
Figure 9. Interpreted GPR profile No. 56 at the basement area showing the expected disturbed soils near the tower.

Figure 10 .
Figure 10.Interpreted GPR profile No. 57 at the basemen

Figure 11 .
Figure 11.Interpreted GPR profile No. 61 a surround-t the ing area showing the expected foundations.

Figure 15 .
Figure 15.Interpreted locations of interrupted s (red oils circles) and expected foundation (green box), surrounding area.

Figure 16 .
Figure 16.Interpreted GPR profile No. 81 at the outside wall near elevator door.

Figure 17 .
Figure 17.Interpreted GPR profile No. 85 at the outside

Figure 19 .
Figure 19.Interpreted GPR profile No. 89 at the ground of the first floor.

Figure 20 .
Figure 20.Interpreted GPR profile No. 99 at the slab with fractures at Basement floor.ralheritage and underground, soil-foundation conditions and interactions, GWT and