The Importance of Accounting for Atmospheric Effects in the Application of NDVI and Interpretation of Satellite Imagery Supporting Archaeological Research: The Case Studies of Palaepaphos and Nea Paphos Sites in Cyprus

This paper presents the findings of the impact of atmospheric effects when applied on satellite images intended for supporting archaeological research. The study used eleven multispectral Landsat TM/ETM+ images from 2009 until 2010, acquired over archaeological and agricultural areas. The modified Darkest Pixel (DP) atmospheric correction algorithm was applied, as it is considered one of the most simple and effective atmospheric corrections algorithm. The NDVI equation was applied and its values were evaluated before and after the application of atmospheric correction to satellite images, to estimate its possible effects. The results highlighted that atmospheric correction has a significant impact on the NDVI values. This was especially true in seasons where the vegetation has grown. Although the absolute impact on NDVI, after applying the DP, was small (0.06), it was considered important if multi-temporal time series images need to be evaluated and cross-compared. The NDVI differences, before and after atmospheric correction, were assessed using student’s t-test and the statistical differences were found to be significant. It was shown that relative NDVI difference can be as much as 50%, if atmosphere effects are ignored. Finally, the results had proven that atmospheric corrections can enhance the interpretation of satellite images (especially in cases where optical thickness of water vapour is minimized ≈ 0). This fact can assist in the detection and identification of archaeological crop marks. Therefore, removal of atmospheric effects, for archaeological purposes, was found to be of great importance in improving the image enhancement and NDVI values.


Introduction
Many applications of Remote Sensing Archaeology have used multi-temporal satellite imagery for the detection of sub-surface archaeological remains [1][2][3][4][5][6], monitoring of archaeological sites and their surroundings [7,8], or for risk assessment analysis [9].In cases of detection of buried archaeological remains, vegetation crop marks, which can be recognized from air/space, are often used.
Vegetation crop marks may be formed in areas where vegetation cover shallow sub-surface archaeological remains.In areas of archaeological interest, the soil tends to retain different percentages of moisture compared to the rest of the crops in the area without any archaeological remains.That means that if vegetation grows over a ditch, especially during periods when water stress develops, the crop growth is likely to be enhanced due to the fact that the topsoil holds more moisture than in the surrounding context.This phenomenon can be recorded from the air (satellite images, aerial photos) and is generally referred to as positive crop mark.However, in cases where there is not enough moisture in the retentive soil and there is a lack of available water for evapotranspiration (vegetation grown above building remains or compacted ground), the developed marks are characterized as negative crop marks and are generally less common than the positive ones [10].Furthermore, the evolution of crop marks can be a combination of different factors and parameters, such as, moisture availability, the availability of dissolved nutrients to the crops at crucial growing times and periods, for the plant, the crop, the soil type and finally the soil depth [11].
In such studies, vegetation indices have been widely used for the detection of archaeological remains through the enhancement of crop marks.The NDVI equation is the most known index applied in different archaeological areas of the world [12][13][14].Although satellite imagery and vegetation indices (e.g., NDVI) are used for the detection of crop marks, atmospheric corrections are often not applied in such studies.However, as Lillesand et al. [15] argued, satellite images need to be atmospherically corrected before being subjected to any post-processing techniques.The aim of this study is to validate atmospheric correction methodologies for the purposes of detecting archaeological remains and to assess their importance for crop marks detection.

Background
In the context of digital analysis of remotely sensed data, pre-processing is generally characterized by two types of data correction: (1) radiometric pre-processing, which addresses variations in the pixel intensities (digital numbers, DN); and (2) geometric correction, which addresses errors in the relative positions of pixels, mainly due to the sensor viewing geometry and terrain variations.Radiometric corrections are distinguished between calibrations, de-striping approaches, atmospheric corrections and removal of data errors or flaws.Radiometric correction is considered to be a more difficult technique than geometric correction since the distributions and intensities of these effects are often inadequately known.Despite the variety of techniques used to estimate the atmospheric effect, the atmospheric correction remains a difficult task in the pre-processing of image data.Furthermore, based on the target under investigation, its impact is an important issue in the retrieval of true and correct outcomes.
The effects of the atmosphere on spectral signatures and vegetation indices have become an important issue in relevant scientific literature since the 1980s (e.g., [16]).Atmosphere is a primary source of noise for accurate measurement of surface reflectance with optical remote sensing [17][18][19].Atmospheric effects are a result of molecular scattering and absorption, and influence the quality of the information extracted from remote sensing measurements, such as, vegetation indices.Such errors, caused by atmospheric effects, can increase the uncertainty up to 10%, depending on the spectral channel [20].Hadjimitsis et al. [21] highlighted the importance of considering atmospheric effects when several vegetation indices, such as NDVI were applied to Landsat TM/ETM+ images for agricultural applications.In their study a mean difference of 18% for the NDVI was recorded before and after the application of darkest pixel method.
Therefore, removal of the atmospheric effects is an important pre-processing step required in many remote sensing applications, since it is needed to convert the at-satellite spectral radiances of satellite imagery to their at-surface counterparts [22,23].Moreover, if an image is to be used for change detection and monitoring purposes, it is essential that adequate atmospheric and radiometric processing be applied, in order to bring all scenes to a common radiometric datum.
Several atmospheric correction algorithms are found in the literature, ranging from simple to sophisticated [24][25][26].Hadjimitsis et al. [27] classified these algorithms into the following two categories: • Category (A): Absolute corrections (corrections that lead to surface reflectance).This category can be subdivided into two sub-categories: image-based atmospheric corrections (for example, Darkest Pixel, Covariance Matrix method) and corrections using independent data for atmospheric optical conditions (including in situ measurements or historical records) using physical-based algorithms; • Category (B): Relative corrections (corrections that do not produce values of surface reflectance).
The advantage of the image-based algorithms is the fact that they do not require any ground data during satellite overpass [28].One of the simplest, fully-imaged based, is the Darkest Pixel (DP) algorithm, categorized in Category (A).The DP principle is based on the assumption that most of the signal reaching a satellite sensor from a dark object, such as, deep inland water bodies (as in the case of the Asprokremmos Dam followed in this study) is contributed to by the atmosphere at visible wavelengths and on the assumption that the near-infrared and middle infrared image data are free of atmospheric scattering effects [28].Therefore, the pixels from dark targets are indicators of the amount of upwelling path radiance in that band.The atmospheric path radiance adds to the surface radiance of the dark target, giving the target radiance at the sensor.The surface reflectance of the dark target, for the DP algorithm, is approximated to have zero surface reflectance.A modified adaptation of the DP method is to assume a known non-zero surface reflectance of the dark target based on ground truth data (e.g., spectroradiometric measurements).
Other than water bodies, atmospheric effects, account for the majority of the at-satellite measured radiance in the visible bands [29], and therefore such targets provide an opportunity to assess the effectiveness of the varying atmospheric correction methods [30,31].This has been shown by Hadjimitsis et al. [25], who provided a critical assessment of the effectiveness of most of the available algorithms using Landsat TM satellite imagery and in situ spectroradiometric measurements.It was found that the DP method provided a reasonable correction in Landsat TM bands 1, 2, and 3, at least for cloud-free scenes [25].It has been shown by other researchers [27,32,33] that the DP is the most effective atmospheric correction.The effectiveness of the DP in the NIR band must be taken into consideration in the cases where water vapor absorption is significant [27].Forster [34] argues that the water vapor thickness is very small for wavelengths less than 0.7 μm but should be taken into account for longer wavelengths such as NIR.
In many archaeological applications, there is a need to use, not only up-to date satellite imagery but, archived multi-temporal data as well.Therefore, the use of physical-based atmospheric correction (Category (A)) can be very difficult to apply, since auxiliary data may not exist for the case study.Moreover, image-based algorithms can be performed easier, even to archive data, regardless of satellite sensor and satellite characteristics.
In this study, an image-based atmospheric algorithm was applied using ground spectroradiometric data acquired over a water body (Asprokremmos Dam).In addition, meteorological data were acquired in order to retrieve optical thickness due to water vapor and to validate the atmospheric correction results.The results were compared to those calculated without removing atmospheric effects (apparent reflectance).For this reason NDVI diagram for a whole phenological cycle of barley crops was drawn and evaluated.Lastly, the results were assessed using visual interpretation of archaeological crop marks.

Spectral Reflectance (ρ)
The apparent reflectance model is used to correct the digital numbers of the satellite image from solar radiance and sun zenith angle effect.However, this model does not take into consideration any effects occurring from atmospheric scattering and absorption, or from topographic effects.For Landsat TM/ETM+ images calibration is performed based on Equation (1).All the necessary metadata were acquired for each image from the metadata header file: where: Lsat: apparent spectral radiance at satellite sensor (W•m 2 •ster 1 •μm 1 ), Lmaxλ and Lminλ: the calibration constants, (W•m 2 •ster 1 •μm 1 ), Qcalmax: the maximum quantized calibrated pixel value and Qcal is the digital number of the image.
Following the correction of the Digital Numbers of the images, from the solar radiance effect, the next step was to calibrate the images from the sun zenith angle effect as well.Therefore, to convert the at-satellite radiance values into at-satellite reflectance, using the solar irradiance at the top of the atmosphere, sun-earth distance correction and solar zenith angle Equation ( 2) was employed: where: ρ: spectral reflectance at the surface (%), d: earth-sun distance in Aus, ESUN λ : solar spectral irradiance on a surface perpendicular to the sun's ray outside of the atmosphere (W•m 2 •μm 1 ) and cosθ s : cosine of solar zenith angle.

Modified Darkest Pixel Atmospheric Correction Based on Spectroradiometric Measurements
The modified Darkest Pixel (DP) method [25] was applied in the current study.The surface radiance of the dark targets is assumed to have approximated zero surface radiance or reflectance [27].Instead of assuming L darkest to be zero value, the modified DP considers the 'true' ground radiance or reflectance value over dark targets as the L darkest .Therefore, the spectral reflectance is given by Equation (3): where: L darkest = Radiance at the darkest target.
The modified DP method applied in this study is based on the assumption that the dark target radiance is not equal to zero since dark targets such as inland water bodies (lakes, reservoirs, dams) are not fully black targets.Therefore, in situ measurements are needed in order to measure the "truth" reflectance of the darkest target.For the study area of Palaepaphos, a nearby large inland water body, the Asprokremmos Dam, was selected as the suitable selected dark target.Field spectroradiometric measurements at Asprokremmos Dam were acquired in order to retrieve the surface 'true' reflectance or radiance values (i.e., L darkest ).

Study Areas
For the purposes of this study, two cultivated areas from known archaeological areas (Arkalon-Marhcello, Sites 1-2), and one traditionally agricultural area (Mourokambos, Site 3), have been selected as indicated in Figure 1.Furthermore, the crop mark of the buried ancient amphitheatre of Nea Paphos (Site 4) was used for the evaluation of the importance of atmospheric corrections.Finally, the Asprokremmos Dam (Site 0) was used for the application of the modified DP algorithm.
Palaepaphos is situated on the SW coast of the island and was selected as a case study area.Palaepaphos is an extensive archaeological landscape and is listed in the World Heritage List of UNESCO since many important monuments, dating from the Chalcolithic period until today, are situated there.Its few visible secular and sepulchral monuments and its famous open-air sanctuary to an aniconic deity, known as Aphrodite, are scattered over an area of two square kilometers [35].
In this area, archaeological excavations have systematically been carried out, while in the past few years, systematic foot surveys and geophysical surveys have also been performed.Arkalon and Marhcello localities (Sites 1 and 2 respectively) are located in the Palaepaphos archaeological site and have attracted the interest of the archaeologists in recent years [14,[36][37][38].
The surroundings of Palaepaphos are also considered to be an important agricultural area for Cyprus.At the Mourokambos locality (Site 3), barley crops have been traditionally cultivated.It should be mentioned that during 2009-2010 (the period studied) in both sites (Palaepaphos and Mourokambos) barley crop was cultivated.
Moreover, Palaepaphos and Mourokambos areas are in close vicinity with the Asprokremmos Dam, the second largest dam in Cyprus (Site 0), which was used for the modified DP algorithm.
Nea Paphos (Site 4), according to written sources, was founded at the end of the 4th Century by Nicocles, the last king of Palaepaphos.In the beginning of the 3rd Century BC when Cyprus became part of the Ptolemaic Kingdom, which had its capital in Alexandria, Nea Paphos became the center of Ptolemaic administration in Cyprus.Until the end of the 2nd Century BC, the site became the political and economical centre of the region and therefore the Ptolemies made Nea Paphos the capital of Cyprus.
Nea Paphos regained some of its glamour during the Byzantine and Medieval periods but from the Venetian period and onwards, the coastal settlement of Nea Paphos was abandoned and the population moved further inland where the present town of Paphos was established [39].

Overall Methodology
For the objective of this study, Landsat TM/ETM+ multispectral images have been collected from USGS Glovis.The images were geometrically corrected and then, the apparent reflectance was calculated.In this step, the digital numbers (DN) of the images were converted into units of radiance and then reflectance.Following that, atmospherically corrected images were produced, using ground spectroradiometric measurements over the Asprokremmos Dam.These measurements were collected periodically from April to October (2010).
The next step included the calculation of the NDVI index over the Palaipaphos archaeological site.The NDVI values were compared both for the archaeological and non-archaeological site (agricultural areas).In both cases, the same kind of crop was cultivated (barley), while the areas have similar soil and climatic characteristics.Finally, interpretation of the results over the archaeological site of Nea Paphos and the Asprokremmos Dam was performed in order to evaluate atmospheric effects.For this purpose the global minimum and global maximum DN value for each band, for all images dataset, was calculated and then linear histogram enhancement was performed.

Resources
For the purposes of the study, multispectral Landsat TM/ETM+ satellite images have been used.Next, ground spectroradiometric data from the inland water of the Asprokremmos Dam have been acquired.Finally, meteorological data have been collected for the needs of the atmospheric corrections algorithms.A detail description of the data is provided in the followings chapters.

Satellite Imagery
With the opening of the Landsat archive by the USGS in 2009, a huge amount of images have become available to scientists [40].Eleven (11) Landsat TM/ETM+ satellite images acquired on 2009 and 2010 were used in this study, to cover the whole life cycle of the barleys crops as shown in Table 1.The ERDAS Imagine v.10 software was used for the pre-processing and post-processing of satellite imagery.

Field Spectroradiometric Data
In this study, a power-engine boat was used to support the in situ campaign in the Asprokremmos Dam in Paphos and a Global Position System Garmin GPS72 was used in order to store the selected sampling station [41].Sampling campaigns started in April 2010 and last until October 2010 including 16 sampling campaigns.A handheld GER-1500 field spectro-radiometer (spectral range covered by the instrument extends from 400 to 1,050 nm) equipped with a fiber optic probe was used, in order to retrieve the spectral signatures of the Asprokremmos Dam (see Figure 2).Reflectance was calculated as a ratio of the target radiance to the reference radiance.The target radiance value is the measured value taken on the surface water of the reservoir and the reference radiance value is the measured value taken on the standard Spectralon panel, representing the sun radiance, which reaches the earth surface-without atmospheric influence.Figure 3 show typical reflectance values of the dam.

Meteorological Data
Meteorological data were provided from the Meteorological Service of Cyprus, in order to search if any significant climatic differences have occurred during the satellite overpass.As it was found (Table 2), the relative humidity during the satellite overpass was similar (≈57-67%) for all days except 19/03/2010 (33%) and 15/06/2010 (77%).Relative humidity and temperature as shown by Forster [34] and by Hadjimitsis and Clayton [30] can be used to provide a measure of the equivalent mass of liquid water or water vapour thickness.Although DP algorithm can be applied without any auxiliary meteorological data, as a fully image-based technique, the authors used these data to investigate any possible water vapour absorption effect in atmospheric correction.According to Forster [34] and McClatchey et al. [42], optical thickness values of water vapor can be calculated based on relative humidity and temperature.In this case, the partial pressure of water vapour is calculated based at the assumption of ideal gas law.Absorption, due to water vapor for wavelengths lower than 0.7 μm, can be ignored [34].Therefore, the assumption that absorption of water vapor for visible bands of Landsat TM/ETM+ (Band 1: 0.45-0.52μm; Band 2: 0.52-0.60μm; Band 3: 0.63-0.69μm) is zero, was made for the objective of this study.However, optical thickness of water vapor, for the NIR band (Band 4: 0.76-0.90μm) for Landsat, needs to be calculated as indicated in Table 3.As it was found, the last three images (19/03/2010; 13/04/2010 and 15/06/2010) used in this study had the minimum optical thickness of the water vapor from all dataset.In the rest of the image dataset, the optical thickness ranged from 0.031 to 0.067.

Results
The inland water spectroradiometric measurements (Table 4) have shown that the reflectance values over the Asprokremmos Dam were very low and tend to be nearly zero, especially at the NIR band.In detail blue band has an average reflectance of 1.67% (±0.44%), green band 3.21% (±0.72%), red band 1.07% (±0.30) and the NIR band 0.09% (±0.03%).According to statistics, standard deviation of the results was less than 1% while for Red and NIR bands (which are often used for vegetations indices) the deviation was less than 0.30%.These results were used for the modified DP algorithm.
Modified DP algorithm was applied for removal of the atmospheric effects using the mean values of the spectroradiometric campaigns.The NDVI index was calculated for Sites 1-3 using atmospherically corrected satellite images and non-atmospherically corrected images.As is indicated in Figure 4, lower values of NDVI can be found in images where no atmospheric correction was applied.Modified DP algorithm tends to give higher NDVI values for all sites selected in this study.The average difference of NDVI values before and after the DP algorithm was calculated to be 0.10.However, the difference in relative values of NDVI before and after the atmospheric correction is estimated to be 25%.This is far less than using DP without any ground measurements (0.16 NDVI mean absolute difference or 32% relative difference).
The difference of NDVI values was assessed using the Student's t-test in order to evaluate the statistical significance of the differences between the two sample means: before and after the atmospheric correction.As it was found t observed for Sites 1-3 (t 1 = 5.205; t 2 = 5.064; t 3 = 14.587) was higher than the t statistical (t 1-3 = 2.228) for (n − 1 = 10) degrees of freedom and at a confidence level of 95%.This implies that NDVI values after the atmospheric correction have a significant statistical difference from the corresponding NDVI values, before the removal of the atmospheric effects.
In the case of Landsat TM or ETM+ satellite images acquired at 15/06/2010, where low water vapor optical thicknesses occurred, the difference of NDVI values before and after the modified DP algorithm is greater (≈0.15 absolute difference of NDVI values or 30% relative differences).Concerning DP algorithm, the difference was estimated at 0.30 absolute NDVI difference or 60% relative difference.These differences might be a result of the low impact of water absorption in the image reflectance.Such differences have also been reported in relevant literature.Indeed, Song and Woodcock [19] argue that, the NDVI calculated from apparent reflectance is the lowest while the NDVI values from different other algorithms (including DP algorithm) are much higher.Moreover, the observed pattern matches the theoretical analysis mentioned by Pinty and Verstraete [43], that atmospheric effects reduce NDVI.Special attention is given in the satellite images, acquired from September to May (No 6-10 in Figure 4), when barley crops were cultivated in the archaeological area of Palaepaphos.As is shown in Figure 4, significant differences might occur in the image interpretation, since NDVI values vary a lot in this period, if there is a need for time-series satellite interpretation.
However, it should be noticed that, although NDVI atmospheric corrected, provide the optimum results due to the removal of all the major radiometric effects, the NDVI values without taking into account the atmospheric effect, tend to give similar results with a variant offset of values in the Y-axes (NDVI values).This is a result of the NDVI equation, since it is a ratio of Red and NIR bands and therefore tends to self-calibrate atmospheric impact.Nevertheless, NDVI cannot minimize atmospheric effects like other indices, such as ARVI and SARVI.Indeed, such indices reduce the influence of the atmosphere by engaging the red-blue channel instead of the red one.ARVI and SARVI indices retrieve information regarding the atmosphere capacity, since compared to the red band, the blue band is much more easily scattered by the atmosphere particles [44].Evaluation of the importance of atmospheric effects for the detection of archaeological crop marks is a difficult task.In order to evaluate the impact of atmospheric correction at satellite interpretation, the authors selected a known archaeological crop mark, at the Nea Paphos archaeological site (Site 4) and a dark target (Site 0), for visual interpretation.The ancient amphitheatre of Nea Paphos (Figure 5) was selected due to its size, which is detectable from Landsat images (≈ 3 × 2 pixels in Vis and NIR bands).
The maximum and minimum histogram values of each band were estimated for all satellite images (before and after the atmospheric correction) in order to retrieve the max Global and min Global.Afterwards, a linear stretch was performed for each band based on these values.Using this simple technique, all scenes were corrected to a common radiometric range (see Figure 6).Table 6 indicates the results of the atmospheric correction for each band, at the surroundings of the Asprokremmos Dam.As is shown, the interpretation is improved after the atmospheric effects were removed.This is obvious for all bands and for the false-color images produced (NIR-R-B) by them.Atmospheric corrected images tend to enhance the image histogram and to improve both the final quality and the interpretation of the images.This is very important for remote sensing archaeological research since interpretation is considered to be a helpful tool for the detection of buried archaeological remains through crop marks.
Table 7 shows some of the results before and after the application of modified DP atmospheric correction, at the Nea Paphos archaeological site.Again, the interpretation qualities of the images are slightly improved after the removal of atmospheric effects.The Best photo interpretation quality was found at images with low water vapor optical thickness (19/03/2010; 13/04/2010 and 15/06/2010) for both sites (Nea Paphos and Asprokremmos Dam).On these days, the optical thickness of water vapor was minimum (less than 0.03).Nevertheless, for the rest of the images, where water vapor optical thickness was less than 0.05, the quality of the images after the atmospheric correction was improved.In the case of higher values (e.g., 14/07/2009, water vapor optical thickness = 0.063) the quality was not sufficiently improved.
Table 8 indicates some results before and after atmospheric correction, based on the DP algorithm.In this case, the Asprokremmos Dam was used as the dark target.This methodology assumes that the pixel of the dam with the lowest digital number (DN) or radiance in each band should actually have zero value and therefore its radiometric DN or radiance value will represent the atmospheric additive [27].As is shown, DP algorithm tends to enhance better archaeological crop marks in contrast with the non-corrected satellite image.Moreover, in the cases of low optical thickness of water vapor values (for example, as the case of the image acquired on 15/06/2010), this enhancement is maximized.However, during the linear stretch application (using the max Global and min Global ), it has been found from the image statistics that the image acquired at 15/06/2010 has a reflectance value of 33% for band 4, which corresponds to the lowest reflectance values from all multi-temporal imagery series.The global maximum value (from all images) was 57%.This caused an over-stretch enhancement to the image, i.e., presence of saturated pixels.However, the effect of the DP atmospheric correction, improves the image significantly, as shown in Figure 8 and discussed later.
Hence, it should be mentioned that DP can be used in cases where "true" dark target is visible in the image and water vapor optical thickness is minimized (≈ 0).Table 6.Interpretation of the results before and after the atmospheric correction at the Asprokremmos Dam.In order to evaluate the difference between NDVI values, before and after atmospheric correction, the relative NDVI values after atmospheric correction (percentage relative difference) was plotted against non-atmospheric corrected NDVI values, for all Sites, as shown in Figure 7.As was expected, for Site 0 (dam) a linear correlation (before and after atmospheric correction) was found, since a single value (average from red and NIR bands) was used.Nevertheless, for the rest of the sites, significant differences that are non-linear have been observed.Relative differences for low NDVI values (<0.40) can be as 50% more in cases of atmospheric corrected images in comparison with non-atmospheric images.For higher NDVI values (>0.40) relative differences tend to be smaller (≈20%) (see Table 9).Indeed, the effect of the atmosphere must be taken into consideration since the magnitude of atmospheric effects to the NDVI tends up to 50 % as shown in this study.In order to evaluate the effectiveness of the revised darkest pixel atmospheric correction, ground spectroradiometric measurements using the GER1500 field spectroradiometer from field campaigns were used.It is apparent that the difference of the NDVI results without atmospheric correction and ground measurements can be 15%.It is important to highlight that any omission of atmospheric effects may cause erroneous or non-qualitative outcomes regarding the presence of archaeological crop marks in the area of interest.Similarly, Figure 9 indicates the results obtained from an agricultural field near Paphos airport.Again, atmospheric corrected images match closely with spectroradiometer derived NDVI values.

Conclusions
Remote sensing archaeological investigations use not only up-to date satellite imagery, but also archived data.All these datasets should have, firstly, a radiometric correction, in order to retrieve valuable information.An important radiometric correction is the atmospheric correction which can be easily applied using image-based atmospheric corrections techniques, such as the simplest one, the DP algorithm.Most published atmospheric correction methods have been developed for specific targets, such as, coastal and inland waters or land, and their effectiveness for other targets and different applications needed first to be assessed [21,45,46].Indeed, an example of how important the atmospheric correction is, for archaeological research, especially when NDVI's are used for archaeological crops detection, is shown in this study.It is considered that omissions in considering the effects of the atmosphere when vegetation indices from satellite images are used, may lead to major discrepancies in the final estimation of evapotranspiration.The absence of such considerations as the effect of the atmosphere, when NDVI's are used in archaeological research, is demonstrated by this study.

Figure 1 .
Figure 1.Case study areas mentioned in the text.Site 0 is the Asprokremmos Dam, Sites 1 and 2 are the archaeological sites of Palaepaphos (Arkalon and Marchello localities), Site 3 is the Mourokambos agricultural field and Site 4 is the Nea Paphos archaeological site.

Figure 2 .
Figure 2. GER1500 field spectroradiometer (a) Reference spectroradiometric measurements on the white spectralon panel (b) Target spectroradiometric measurements using a fiber optic probe.

Figure 4 .
Figure 4. NDVI values for Site 1 and Site 2 using: (a) reflectance values and (b) after the modified DP algorithm reflectance values were applied.

Figure 5 .
Figure 5. Crop mark at the Nea Paphos archaeological site.

Figure 6 .
Figure 6.Global min and Global max diagram for histogram enhancement.

Figure 7 .
Figure 7. Relative NDVI values after atmospheric correction (modified DP) against NDVI values before atmospheric correction.

Figure 8 .
Figure 8.Comparison of NDVI values derived from satellite images (before and after atmospheric correction) and ground spectroradiometric measurements for Site 1.
Figure 8 shows the results from field campaign at Site 1 during March 2010 (05/03/2009; 12/03/2009 and 24/03/2009).As is shown, ground spectroradiometric results are closer to atmospheric corrected NDVI values (using the DP and modified DP) in relation to NDVI values before atmospheric correction.Modified DP atmospheric correction tends to be most accurate compared to DP algorithm.

Figure 9 .
Figure 9.Comparison of NDVI values derived from satellite images (before and after atmospheric correction) and ground spectroradiometric measurements for the agricultural site near Paphos airport.

Table 1 .
Satellite image acquisition dates.

Table 2 .
Meteorological data from the Paphos airport (Data provided from Meteorological Service of Cyprus).

Table 3 .
Equivalent mass of liquid water for the satellite images used in this study.

Table 7 .
Interpretation of the results before and after the modified DP atmospheric correction applied at the Nea Paphos crop mark (the archaeological crop mark in Nea Paphos has been circled, 1st row, 1st image from the left).

Table 8 .
Interpretation of the results before and after the DP atmospheric correction applied at the Nea Paphos crop mark.

Table 9 .
Relative NDVI differences for all sites, before and after atmospheric correction (min and max values for each site are underlined).