Determining the Rate of Salinity of Persian Gulf Waters with the Aid of Satellite Images and Least Squares Method

As yet various methods have been used for determining the salinity rate of seas and oceans water. The current method of determining salinity rate of seas water has been field examination of various points of sea and determining its salinity rate. In the last decade, remote sensing satellite images have had high capability in determining sea waters salinity rate. Regarding that the present methods in remote sensing depend on the studied regions, therefore, the necessity of customization of these methods is felt. Fresh water springs due to impact on water salinity and temperature and also the environment physics and density like sound velocity are very significant and since coasts and islands of Persian Gulf are considered among arid and semi-arid regions and lack drinking water, access to fresh water springs has more significance. After studies performed, preparation of salinity rate observations and catching two series of proper images for felid data for complete coverage of the region, preprocessing and calibration was performed. For this purpose in turning the acquired radiance to reflection, ENVI software was used. The histogram of calibrated shades of gray rates in images was specified, so that reflection of each sample can be extracted from images. In this paper, the rate of least method efficiency in determining salinity rate of Persian Gulf waters was examined and finally identifying fresh water pits using remote sensing technique was done. The obtained results in the least squares methods after combining various bands of image with each other specified that combining 4 bands of 2, 3, 5 and 7 has the least standard deviation rate with training data and test, which is equal to 0.385 and 0.991978.


Introduction
In continental respect, Iran is locates on the arid and semiarid world belt and its rain average is 250 mm. according to the present statistics, 25% of Iran soil is located in extra dry regions, 40% in dry regions and 25% in semidry regions. Water sources are located in small part of the country area and should be transferred from these regions to two third of Iran area which is among shallow regions. Regarding the rate of water resources and consumption per capita, Iran is among countries which are in the group of countries faced with water shortage. This group includes countries which will face water physical shortage in 2025; this means that even with the highest possible efficiency in water consumption won't have enough water to supply their needs. About 25% of the world people including Iran are in this group.
Based on Falcon Mark indexes, National Organization index and International Institute of water management, Iran is in the threshold of standing in sever water crisis status. Based on mentioned indicators, for preserving its current status till 2025, Iran should be able to add 112% to its producible water resources that this rate regarding facilities and existing water resources seems impossible [1]. One of methods of increasing water resources in future is using present water resources in seas beds. Fresh water resources exist in some seas, so that the present underground aquifers present in lands could continue to seas and enter the sea as spring ( Figure 1). Persian Gulf is not exempt from this case, Persian Gulf bed is the continuation of Zagros Mountain so that a part of its anticlines is visible as islands. In geology and petrology respect, Zagros Mountains are calcareous that this matter helps appearance and existence of springs in Persian Gulf, even some of these underground aquifers have been transferred towards Arabia and have appeared as sweet water springs [2]. Exit of underground aquifers (SDG) in some cases is featured as fresh water resources in sea bed and creates some changes in water salinity, temperature and density. Change in water density could impact marine streams, creating sonic channels and waves features that with the aid of these factors, these springs could be locating. Fresh water springs existent in seas and oceans is among new and intact resources of fresh water. One of places which have been considered in this respect is Persian Gulf. Locating of fresh water springs or regions which have less salinity, not only makes possible access to fresh water, but also will considerably help preventing environmental pollutions. Regarding that about 60% of all the world humidification are in Middle East, one of problems of these machines is that the most modern humidification has 50% efficiency, it means they convert half of water they take and the rest enters sea water as bittern. The rate of salts and heavy metals existent in bittern is about 64 to 70 ppt. portion in thousand. This means that due to activity of humidification machines daily more than 1 million and 200 thousands m 2 bittern enters Persian Gulf. Persian Gulf due to its natural conditions doesn't have high speed of water full circulation, so for Persian Gulf water circulation to be able to have a full circulation and annihilate pollutions of Persian Gulf, 4 years' time is required. It is obvious that extracting fresh water springs or regions which have less salinity could help preventing from the sea environmental pollutions too [3] [4]. Also, desalination of waters which have less salinity will have less cost than saline waters, on this basis, one of the main objectives of this paper is to identify fresh water pits in saline environments with the aid of satellite images and remote sensing knowledge which besides having low cost comparing fresh water reactors is considerable in environment respect.

The Research History
Many studies have been performed bout identifying fresh water springs existent in seas abroad, but despite high significance of this issue, in arid and semiarid regions of Persian Gulf coasts and islands ,this matter has less been noticed. Here, some studies performed in the field of appearance of fresh water springs in seas and their physical specifications are pointed (Table 1).

Geographic Situation of Study Range
Persian Gulf is located in 24 to 30 degrees and 30 minutes northern width and 48 to 56 degrees and 25 minutes eastern width from Prime Meridian. This Gulf intercommunicates with Oman Sea through Strait of Hormuz and through it to high Seas. The studied region is located in east of Persian Gulf (Strait of Hormuz) and part of Oman Sea in 25 degrees and 40 minutes to 27 degrees and 10 minutes of northern width and 56 degrees to 57 degrees and 17 minutes of eastern length from Prime Meridian. Due to existence of salinity data from this region, this part has been selected as studied region [5].

Research Method
In this article, for achieving research objectives, the following stages were performed and its flowchart is specified in Table   This study's result shows correlation between sea water salinity rate and Landsat TM assessing bands. 1982 Khorram 1 They could use radar images ( microwaves) for determining seas water salinity rate 1982 Sonia et al.

2
In US, for the first time, he could register locating fresh water using electrical conduction in the form of an invention. In this method, a ship moved parallel to coastal line and using an electrical Conductometer which had been dropped out of the ship to the sea, identified fresh water springs.
1992 Laing 3 In Jordan country, he could identify the exit spot of fresh water springs using infrared imaging and field measuring of salinity and temperature factors in Dead Sea.

Akawwi 4
They showed that sea waters salinity rate is extractable using satellite optic images. 2008 Ahn et al.

5
They used radar images for determining waters salinity rate. In radar images, for determining sea waters salinity rate, measuring water surface temperature rate is required. Therefore, they could make a relation between sea waters salinity rate and surface temperature.  They used MODIS images and genetic algorithm for determining Taiwan Strait waters salinity rate. In this method, it was hypothesized that the previous research mathematics relation is not optimal and just supports a linear state. For this reason, they used genetic algorithm for determining the main relation between salinity and MODIS image bands.

The Research Findings
For identifying salinity rate of the region waters, Landsat 5 TM assessor of 1994 was used [6].
But regarding that the studied region didn't stand in one imaging scene, 2 series of images were used for complete coverage of the region. Two taken images with close distance have been selected to minimize errors due to time change. Meanwhile, it should be pointed out that due to two parts of studied region being located in various courses of satellite, this issue was unpredictable. False colored display of images has been shown in Figure 3 and Figure 4.   In the next stage, Landsat images were preprocessed and calibrated. After that, the rates of calibrated histogram shades of gray were obtained for images 3 and 4 ( Figure 5 and Figure 6). For Figure 3, major shades of gray rates have been spread in the range of 0 to 0.3. For Figure 4, major rates of shades of gray have been spread in the range of 0 to 0.25.

Salinity Data Taken from the Region
For determining and modeling salinity rate of the studied region pixels requires having salinity rate sample in the region, so 55 samples of sea water surface salinity which had been removed by ship in the region was used [7]. for finding the best parameters from Least Squares, training data was required, also for evaluating fitness rate performed, it is required that with other spots accuracy rate to be examined, for this purpose, form 55 existent regions, 40 spots were used for training and 15 spots for evaluating the results.
This part of work aims to create a relation between TM assessor spectral bands (Bi)  with salinity rate (S). For this purpose, 6 bands of Landsat images will be used. For this purpose, the general relation between bands and salinity rate will be the following relation: a a B a B a B a B a B a B = + * + * + * + * + * + * The aim of Least Squares will be acquiring unknowns from a 0 to a 6 that unknowns will be obtained according to what was stated in previous part. For evaluating the results accuracy on training spots and testing the estimated rate of salinity and real rate is obtained and from standard deviation criterion, error rate criterion will be calculated according the following equation: Regarding that there are 6 spectral bands and this question is propounded that which band will have the best result. So in this article, various permutations of 2, 3, 4, 5 and 6 will be used for modeling salinity so that the best combination to be obtained with the least error and finally using that combination and the region image, salinity rate of various pixels is calculated.

Binary Combinations
Various binary combinations were used among images bands and the rate of error and unknown parameters for each was calculated [8]. In Table 3, the error rate of training spots has been displayed for various combinations. In Table 4, the error rate of test spots has been displayed for various combinations. As it is observed, the rate of standard deviation of training data is less than test data that this issue is due to that their fitness is natural. Based on test data, the best results have been obtained for combination of two bands of TM assessors 5 and 7.

Triplex Combinations
Various triplex combinations were used among the image bands and the error rate and unknown parameters were calculated for each. In Table 5, the error rate of training spots has been displayed for various combinations. In Table 4, the error rate of training spot has been displayed for various combinations. Based on test data, the best results have been obtained for combination of 3 bands of TM assessors 3, 5 and 7.

Quadric Combinations
Various quadric combinations were used among image bands and the error rate and unknown parameters were calculated for each. In Table 6, the error rate of training spots has been displayed for various combinations. In Table 7, the error rate of test spots has been displayed for various combinations. Based on test data, the best results are for quadric combinations of TM assessors 2, 3, 5 and 7 bands.

Quintuple Combinations
Various quintuple combinations were used among image bands and the error rate and unknown parameters were calculated for each. In Table 8, the error rate of training spots has been displayed for various combinations. In Table 9, the error rate test spots  been obtained for quintuple combinations of TM assessors 1, 2, 3, 5 and 7. But, the obtained accuracy has been less than the best obtained accuracy for the best quadric combination, so the best bands which can determine salinity rate are: 2, 3, 5 and 7. Also, for

Salinity Map
For preparing salinity map, first the acquired equation between bands and salinity from Least Squares is according to the following relation: For identifying salinity rate of water regions, using ship file, mask region was created around these regions which is observed in Figure 7.
The obtained salinity map from Least Squares Methods has been displayed in Figure  8.

Conclusion
In this article, the main objective was to assess salinity rate of Persian Gulf water using satellite images and the Least Squares Method. So, regarding that there are 6 spectral bands in Landsat TM 5 images and that which band will have the best result, the best combination with least error has been calculated. The obtained results in the Least Squares method after combining the image various bands with each other, it was specified that the combination of 4 bands of 2, 3, 5 and 7 has the least standard deviation rate with test and training data, which is equal to 0.385 and 0.991978. These rates show the difference between the estimated salinity and real salinity rate of samples.